& rSphere, recip_rSphere, radius_fromHorizGrid,
_RL radius_fromHorizGrid
CALL WRITE_0D_RL( radius_fromHorizGrid, INDEX_NONE,
IF ( rSphere.NE.radius_fromHorizGrid ) THEN
tmpFac = rSphere / radius_fromHorizGrid
& radius_fromHorizGrid.NE.UNSET_RL ) THEN
rSphere = radius_fromHorizGrid
IF ( radius_fromHorizGrid.EQ.UNSET_RL ) THEN
radius_fromHorizGrid = rSphere
& rSphere, radius_fromHorizGrid, deepAtmosphere,
radius_fromHorizGrid= UNSET_RL