[FLASH-USERS] Problem with Grid_bcApplyToRegion.F90 file

[Klaus Weide]

On Thu, 22 Mar 2018, María Fernanda Clever Uribe wrote:

> I tried Klaus' recommendation about the j index but I still don't get what
> I wish for.
> 
> 
> 
> At the moment, following the Grid_bcApplyToRegionSpecialized.F90 boundary
> condition from de ShuOsher problem, I am able to inject a wind but from all
> the boundary, and I am not able to only inject the wind from a small region.
>[.....] 
> Do I have to replace the regionSize defined in the Grid file and use
> blkLimitsGC instead, or maybe use endPoints?

I believe you don't want the boundary values depend on j (or any index 
variable) directly, but on the corresponding real coordinate values.

You should find several examples in the FLASH source where something 
like this is done.  Look for files

  Grid_bcApplyToRegionSpecialized.F90 ,

or perhaps just

 Grid_bcApplyToRegion.F90 ,

with the following properties:

 * There is a call (or several) to gr_extendedGetCellCoords.
 * That call (or one of them) applies to the "second direction".
   That is, look for something like

      call gr_extendedGetCellCoords(secondDir, 

   not

      call gr_extendedGetCellCoords(axis, 

 * See how the real coordinates that are retrieved with that call
   (which may be called secondCoord, yCoord, xCoords, or something else)
   are used in the code.

You should find several files - you'll probably want to narrow down the 
choice to those that match your intended use more closely and where you 
can follow the code.

Klaus