*=*=*=*= qsat.html =*=*=*=*
subroutine qsat

subroutine qsat


       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