_RL sigP, recT, tmpQ
sigP = PS(1)
IF (SIG.GT.0.0) sigP=SIG*PS(J)
QSAT(J)= QS1*tmpQ/(sigP-QS2*tmpQ)
RH(J) = QSAT(J)*C1*(T0-T1)*recT*recT*sigP/(sigP-QS2*tmpQ)
RH(J) = QSAT(J)*C2*(T0-T2)*recT*recT*sigP/(sigP-QS2*tmpQ)