ALHevp = ALHC
IF ( aim_energPrecip .OR. useThSIce ) ALHevp = ALHC + ALHF
dEvp(J) = MAX( dEvp(J), -dShf(J)/ALHevp )
& - ( Shf0(J) + ALHevp*Evp0(J) )
& - ( SHF(J) + ALHevp*EVAP(J) )
& - ( dShf(J) + ALHevp*dEvp(J) )
_RL ALHevp, Fstb0, dTstb, dFstb
IF ( iceornot.GE.1 ) ALHevp = ALHC + ALHF