_RL hght(Nrphys), slope,pr1,pr2,hpr1,hpr2
slope = ( hpr1-hpr2 )/( pr1-pr2 )
hght(L2) = hpr2 + ( pres(L2)-pr2 )*slope
_RL qpr1(irun), qpr2(irun), slope(irun)
slope(i) = 0
slope(i) = ( xlat(i) -watlats(jlat-1) )
slope(i) = 1
h2o_time_lat(i,k) = qpr1(i) + slope(i)*(qpr2(i)-qpr1(i))
slope(i) =(QPR1(i)-QPR2(i)) / (PR1(i)-PR2(i))
slope(i) = zero
slope(i) = ( XLAT(I) -ozlats(jlat-1) )
slope(i) = one
o3int1(i,k) = qpr1(i) + slope(i)*( qpr2(i)-qpr1(i) )