source: trunk/src/qsat.pro @ 152

Last change on this file since 152 was 80, checked in by pinsard, 13 years ago

progress on humidity processing
  • Property svn:keywords set to URL
File size: 792 bytes
Line 
1;+
2;
3; ========
4; qsat.pro
5; ========
6;
7; .. function:: qsat(T,P)
8;
9;     :param T: SST (degC)
10;     :param P: surface pressure (mb)
11;
12;     :returns: Air humidity at saturation (g/kg)
13;     
14; EXAMPLES
15; ========
16;
17; ::
18;
19;  IDL> T=+todo+
20;  IDL> P=+todo+
21;  IDL> result=qsat(T,P)
22;
23; TODO
24; ====
25;
26; coding rules
27;
28; complete example
29;
30; get rid of uppercase
31;
32; SEE ALSO
33; ========
34;
35; :func:`rh_to_spechum`
36;
37; EVOLUTIONS
38; ==========
39;
40; $Id$
41;
42; $URL$
43;
44; - fplod 20110811T131113Z aedon.locean-ipsl.upmc.fr (Darwin)
45;
46;   * add minimal header to pk release given today (!)
47;
48; - pbk 2008
49;
50;   * creation
51;
52;-
53function qsat, T, P
54;
55x=T
56es=6.112*exp(17.502*x/(x+241.))*(1.0007+3.46e-6*p)
57y=es*622/(p-.378*es)
58
59return, y
60
61end
Note: See TracBrowser for help on using the repository browser.