DTC =
STRDG1(ChNo)=DTC
& DHSDTC(ChNo)*DTC
ESATNW = ESATTC(CHNO)+DESDTC(CHNO)*DTC
DTC = ( Q0 + ALHX(CHNO)*ETURB(ChNo) - DHSDEA(CHNO)*DEA ) /
IF(TC(CHNO)+DTC .GT. TF .AND. SNOW(CHNO).GT.0.) THEN
DTC = ( Q0X + ALHX(CHNO)*ETURB(ChNo) - DHSDEA(CHNO)*DEA ) /
IF(TC(CHNO)+DTC .LT. TF) THEN
DTC=TF-TC(CHNO)
Q0SNOW=A11*DTC+A12*DEA
DTC = ( Q0 - A12*DEA ) / A11
EVAP(ChNo) = ETURB(ChNo) + DEDEA(ChNo)*DEA + DEDTC(ChNo)*DTC
DTC = ( Q0X - A12*DEA ) / A11
TC(ChNo) = TC(ChNo) + DTC
HLWUP(CHNO) = HLWTC + DHLWTC*DTC
& DEA, DTC, SNLEFT, Q0X, Q0SNOW,
IF(TC(CHNO)+DTC .GT. TF .AND. SNLEFT.GT.0.) THEN
DEA=(F0-A21*DTC)/A22
& DEDTC(ChNo)*DTC