_RL latentHeatFluxMax (1:sNx,1:sNy)
latentHeatFluxMax(I,J) = 0.0 _d 0
latentHeatFluxMax(I,J) = lhSublim * recip_deltaTtherm *
latentHeatFluxMax(I,J) = ZERO
latentHeatFluxMaxMult(I,J,IT) = latentHeatFluxMax(I,J)*pFac
CALL DIAGNOSTICS_FILL(latentHeatFluxMax,