_RL sigP, recT, tmpQ
tmpQ = E0*EXP(C1*(TA(J)-T0)/(TA(J)-T1))
QSAT(J)= QS1*tmpQ/(sigP-QS2*tmpQ)
RH(J) = QSAT(J)*C1*(T0-T1)*recT*recT*sigP/(sigP-QS2*tmpQ)
tmpQ = E0*EXP(C2*(TA(J)-T0)/(TA(J)-T2))
RH(J) = QSAT(J)*C2*(T0-T2)*recT*recT*sigP/(sigP-QS2*tmpQ)