_RL HZ(irun),ARG0(irun),ARG1(irun),DX(irun)
ARG1(I) = 1. - X1(I)
IF ( Z(I) .LT. 0.013 ) ARG1(I) =
ARG1(I) = ARG1(I) * ( 1.+X0(I) )
ARG1(I) = ARG1(I) / ARG0(I)
ARG1(I) = LOG( ARG1(I) )
PSI2(I) = 2. * ARG0(I) + ARG1(I)
ARG1(I) = 1. - Y1(I)
IF( Z(I) .LT. 0.0065 ) ARG1(I) =
ARG1(I) = ARG1(I) * ( 1. + Y0(I) )
PSI2(I) = LOG( ARG1(I) )
ARG1(indx) = VZH(I)
HZ(I) = 1. / ARG1(I)
ARG1(I) = HZ(I)*HZ(I)*X0DEN(I)*X1DEN(I)
ARG1(I) = 0.5 * ARG1(I)
TEMP(I) = ARG1(I) + CM5 * ARG0(I)
ARG1(I) = LOG( Z2ZWM(I) )
PSI2(I) = TEMP(I) + CM6 * ARG1(I)
ARG1(I) = CM7 * Z2ZWM(I) / ( CM2 + Z2(I) )
ARG1(I) = Z2(I) * RZCM
PSI2(I) = TEMP(I) + ARG0(I) * ARG1(I)