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