[FLASH-USERS] setting up a user BC

[turcotte]

Dear all,

I am trying to set-up an inflow entering the computational domain from the lower 
y-axis boundary in 2-D cylindrical coordinates.
Initially I used the same set-up as the WindTunnel example, modifying  
Grid_applyBCEdge.F90 to have a uniform inflow across the lower boundary.
For example, the density of the inflowing material is assigned like:
        case (DENS_VAR)
           do i = 1,guard
            dataRow(i)= max(rho_initial,sim_rho_ambient)
           end do

I need to implement a more complex BC however. To do so I define the inflow in  
Grid_bcApplyToRegionSpecialized.F90. 
My first step was to try to reproduce with Grid_bcApplyToRegionSpecialized the 
simple uniform inflow.
I expected that something like this:
         case (DENS_VAR)
            do i = 1,guard
               regionData(i,:,:,DENS_VAR)= max(rho_initial,sim_rho_ambient)
            end do
would yield the same inflow as before but it's not the case, the inflow is not 
uniform. There are variations that are aligned with block boundaries. 

All else is identical in the two formulations, such as EOS, physical parameters 
of the inflow, etc...
 
What am I missing?
 
Thank you,
Sylvain

---
Sylvain Turcotte
Physics Department    
Bishop's University