call Simulation_locatePhotosphere(real(INOUT) :: rad_photosphere)
Loops through blocks and zones starting from the outer boundary of the simulation domain and integrates optical depth up to tau=2/3 to determine the location of the photosphere.
rad_photosphere : the photosphere radius is returned here
There is no useful default implementation. A simulation has to provide an implementation if needed.