delT = land_deltaT
tg(1) = ( cg(1)*tg(1) + fLoc*delT*beta
& / ( cg(1) + aLoc + bLoc*delT*beta
eLoc = eg(1) + fLoc*delT*beta
cLoc = cg(1) + aLoc + bLoc*delT*beta
eg(1) = eLoc - (aLoc + bLoc*delT*beta)*tg(1)
eg(1) = eg(1) + aLoc*tg(2) + fLoc*delT*beta
land_Pr_m_Ev(i,j,bi,bj) = dMsn/delT
& + delT*fLoc - land_Lfreez*mSnow + aLoc*tg(2)
& / ( cg(1)+aLoc + delT/alpha - aLoc*aLoc/(cg(2)+aLoc) )
tg(1) = cg(1)*tg(1) / ( cg(1)+aLoc + delT/alpha )
dMsn = ( fLoc+tg(1)/alpha )*delT*recip_Lfreez
& + delT*fLoc - land_Lfreez*dMsn
& + delT*(fLoc-bLoc*Tsurf)
_RL mSnow, dMsn, delT