[FLASH-USERS] problems getting coordinates for face variables

[Ken Dere]

Hi,

I am trying to use the unsplit mesh mhd solver in a 2D geometry with
constant gravity in FLASH 3.2.  I am using

call  gr_extendedGetCellCoords(IAXIS, blockHandle, gr_myPE, FACES,
.true., xCoord,  sizesGC(IAXIS))

from  Grid_bcApplyToRegionSpecialized and get the following error
message when trying to get FACEX coordinates:

 pe,block,accblk,axis,j           0           1           1          
2          13 coords differ!!! T
  18000000000.000000000000000     .NE.  
17000000000.000000000000000       2000000000.0000000000000000     
3.81469726562500000000000000E-06

 Driver_abort called. See log file for details.
 Error message is [gr_extendedGetCellCoords] Two ways to calculate cell
coordinates give significantly different answers!

any help would be appreciated.

Ken Dere

---------------------------------------
Kenneth P. Dere
Research Professor of Solar Physics
Dept. of Computational and Data Sciences
MS 6A2
George Mason University
4400 University Drive
Fairfax VA 22030
703-993-4555
703-993-1992 FAX
kdere at gmu.edu