[FLASH-USERS] Grid_bcApplyToRegionSpecialized

[turcotte]

Hello Dongwook and Klaus,

Many thanks for your suggestions.
I will do that right away.

Sylvain

>Hi Sylvain,
>
>I assume you're applying user BC only on the x=xmin boundary.
>If this is not true, and if you're also applying user BC on the
>y direction(s) with the current setups of x,y coordinates, this
>will not work for you, because your BC_DIR also varies.
>
>My first suggestion is to replace Grid_getCellCoords with
>gr_extendedGetCellCoords. This routine has an extended functionality
>that should be used internally within the Grid
>unit (i.e., should be called from Grid_*.F90) and can give you
>coordinate information even for remote blocks.
>
>I suggest your do loop look like (assuming userBC on x=xmin,
>i.e., axis=IAXIS, secondDir=JAXIS, thirdDir=KAXIS):
>(For more general cases, please have a look how secondDir
>and thirdDir are set in gr_bcApplyToOneFace.F90.)
>
>------------------------------------------------------
>sizeX = regionSize(axis)
>sizeY = regionSize(secondDir)
>sizeZ = regionSize(thirdDir)
>
>allocate(xCoord(sizeX))
>allocate(yCoord(sizeY))
>
>call gr_extendedGetCellCoords&
>      (axis, blockHandle, gr_myPE, CENTER, .true., xCoord, sizeX)
>
>call gr_extendedGetCellCoords&
>      (secondDir, blockHandle, gr_myPE, CENTER, .true., yCoord, sizeY)
>
>do i=1,guard
>    do j=1,sizeY
>       regionData(i,j,1:ke,VELY_VAR)=value(xCoord(i),yCoord(j))
>    enddo
>enddo
>
>---------------------------------------------------------
>
>Note again that if you're applying userBC to x=xmax, as well
>as y=ymin or y=ymax, then you have to be more careful in
>setting up your xCoord and yCoord in using axis, secondDir
>as well as proper offsets for face=HIGH.
>
>If you need a BC for that case, I can give you an example
>of Grid_bcApplyToRegionSpecialized.F90 I've implemented.
>
>Hope this helps,
>Dongwook
>
>
>----------------------------------------------------
>Dongwook Lee, Ph.D., Research Scientist
>ASC Center for Astrophysical Thermonuclear Flashes
>University of Chicago
>5640 S. Ellis, RI 484
>(773) 834-6830
>

Sylvain Turcotte
Physics Department    turcotte at apollo.ubishops.ca
Bishop's University   tel: (819) 822-9600 x 2287
Lennoxville (Quebec)  fax: (819) 822-9661
Canada  J1M 1Z7       office: J 11B
URL:  http://apollo.ubishops.ca/~turcotte