[FLASH-BUGS] spherical symmetry II

[JB Gallagher]

Dear developers,

as I wrote in my email from 5. Februar, I think there is a bug in FLASH2.3
concerning the use of 1d spherical coordinates (spherical symmetry).


>> ./setup sedov_sph -1d -auto
>> make
>  
>

produces code which crashes in
source/mesh/amr/paramesh2.0/quadratic_spherical/amr_prolong_gen_unk_fun.F90.
There is a loop around line number 440, starting with


>> do j = 1, 2*nguard+nyb
>  
>

the variables ypl,ypr,yp,dyp are changed here with index j = 1 to 9.
However, they are defined as scalars for ndim=1. As consequence, myGeom is
overwritten (set to CARTESIAN) which leads to an abort.

Can you reproduce this error? The same error occurs in my application of
flash, so it is an urgent problem for me.

Please help,

Peter

-- 
Brad Gallagher
ASCI Flash Center
jbgallag at flash.uchicago.edu