beta = 1. _d 0/(1. _d 0+alpha*bLoc)
tg(1) = ( cg(1)*tg(1) + fLoc*delT*beta
& / ( cg(1) + aLoc + bLoc*delT*beta
tSurf = ( tg(1) + alpha*fLoc ) * beta
eLoc = eg(1) + fLoc*delT*beta
cLoc = cg(1) + aLoc + bLoc*delT*beta
eg(1) = eLoc - (aLoc + bLoc*delT*beta)*tg(1)
& tmp2,(tmp1+alpha*fLoc)*beta
tSurf = alpha*fLoc * beta
eg(1) = eg(1) + aLoc*tg(2) + fLoc*delT*beta
& tmp1,(tmp1+alpha*fLoc)*beta
_RL aLoc, bLoc, cLoc, eLoc, fLoc, alpha, beta
beta = 1.D-11