RCINV = ( DUM1*LOG(DUM2) + ZK*ZLAI(ChNo) ) / DUM3
VPSIM(I) = LOG( VHZ(I) )
CM8 = 6. * LOG(CM8ARG)
ARG1(I) = LOG( ARG1(I) )
PSI2(I) = LOG( ARG1(I) )
ARG1(I) = LOG( Z2ZWM(I) )
& + 2.*BETA*GAMMA*(1.-GAMMA*LOG((1.+GAMMA)/GAMMA))
& ( rF(3) * LOG ( rF(3) / rF(2) ) /
& ( LOG ( zRef(i,j) / rF(2) ) +
zwln = LOG(hu/zref)
ztln = LOG(ht/zref)
& *( LOG( (exf_one + exf_two*x + xsq)
& *exf_two*LOG( exf_half*(exf_one+xsq) )
term2=10. _d 0*LOG(termP/atm_po)
t = SQRT(-2.0*LOG(xs)/xs)
WVSurf(I2,myThid) = (LOG(PSA(I2))-SIGL(K))*WVI(K-1,1)
SIGL(K) = LOG(SIG(K))
WVI(K,2)=(LOG(SIGH(K))-SIGL(K))*WVI(K,1)
zwln = LOG(zwd/zref)
ztln = LOG(zth/zref)
& ( LOG( (one + two*x + xsq)*(one+xsq)*.125 )
& + (one-stable)*( two*LOG(half*(one+xsq)) )
& ( LOG( (exf_one + exf_two*x + xsq)*
& ( exf_two*LOG(exf_half*(exf_one+xsq)) )