_RL HZ(irun),ARG0(irun),ARG1(irun),DX(irun)
ARG0(I) = 1. - X0(I)
IF ( Z0(I) .LT. 0.013 ) ARG0(I) =
ARG0(I) = ARG0(I) * ( 1.+X1(I) )
ARG1(I) = ARG1(I) / ARG0(I)
ARG0(I) = -DX(I) / ( 1. + X1(I)*X0(I) )
ARG0(I) = ATAN( ARG0(I) )
PSI2(I) = 2. * ARG0(I) + ARG1(I)
ARG0(I) = 1. - Y0(I)
IF( Z0(I) .LT. 0.0065 ) ARG0(I) =
ARG0(I) = ARG0(I) * ( 1. + Y1(I) )
ARG0(I) = (Z1(I) + 0.1) * (Z0(I) + 0.1)
ARG0(I) = CM3 + ARG0(I) * CM4
ARG0(I) = ( Z1(I) - Z0(I) ) / ARG0(I)
TEMP(I) = ARG1(I) + CM5 * ARG0(I)
ARG0(I) = 6.
ARG0(I) = YCM
PSI2(I) = TEMP(I) + ARG0(I) * ARG1(I)