*=*=*=*= qsat.html =*=*=*=*
subroutine qsat (dq, q, e, p, t, r) tf=273.16 if (t.gt.tf) go to 1 a = 21.874 b = 7.66 go to 2 1 a = 17.269 b = 35.86 2 e = 6.11*exp(a*(t-tf)/(t-b)) q = 0.622*r*e/(p-r*e) dq=a*p*q*(tf-b) dq=dq/((p-r*e)*(t-b)*(t-b)) return end