msci Code

;given a temperature, calculates the equilibrium water-
;vapour pressure:
function eq_vp, T
  a=-6069.9385
  b=21.2409642
  c=-2.711193e-2
  d=1.673952e-5
  e=2.433502

  return, exp(a/T + b + c*T + d*T^2 +e*alog(T))

end

;given a temperature, calculates the equilibrium water-
;vapour pressure for ice:
function eq_vp_ice, T
  a=-6024.5282
  b=29.32707
  c=1.0613868e-2
  d=-1.3198825e-5
  e=-0.49382577

  return, exp(a/T + b + c*T + d*T^2 +e*alog(T))

end

;Magnus formula:
function magnus_formula, T
  T1=T-273.15
  return, 610.7*10^(7.5*T1/(237+T1))
end  

function alumrof_sungam, wvp
  t=237./(7.5/(alog10(wvp)-2.78583)-1)
  return, t+273.15
end