A boundary condition question: In cylindrical coordinates, how is locus corresponding to r=0 handled in relativistic hydrodynamical calculations? If someone could point me to the relevant stretch of code within FLASH, for example, that would be great, --many thanks! --John