source: trunk/src/dew_temp_airt_rel_humidity_jv.pro @ 81

Last change on this file since 81 was 50, checked in by pinsard, 13 years ago

replace TROPFLUX by PROJECT
  • Property svn:executable set to *
File size: 1.1 KB
Line 
1;+
2;
3; =============================
4; dew_temp_airt_rel_humidity_jv
5; =============================
6;
7; .. function:: dew_temp_airt_rel_humidity_jv(t, td)
8;
9;     :param t: air temperature in degC
10;     :param td: dew point temperature in degC
11;
12;
13; SEE ALSO
14; ========
15;
16; EXAMPLES
17; ========
18;
19; ::
20;
21;    IDL> t=+todo+
22;    IDL> td=+todo+
23;    IDL> result=dew_temp_airt_rel_humidity_jv(t, td)
24;
25;
26; TODO
27; ====
28;
29; hard coded directory - usage of ${PROJECT_ID}
30;
31; coding rules
32;
33; EVOLUTIONS
34; ==========
35;
36; - fplod 20101214T093615Z aedon.locean-ipsl.upmc.fr (Darwin)
37;
38;   * minimal header
39;
40; - pbk 2008
41;
42;   * creation
43;
44;-
45function dew_temp_airt_rel_humidity_jv, t,td
46
47if (not keyword_set(t)) then t=28            ; air temperature in degC
48if (not keyword_set(td)) then td=23           ; dew point temperature in degC
49es0=6.11        ; reference saturation vapor pressure (hPa)
50
51;; calculating fields.  First saturation vapor pressure (Es) and
52;; actual vapor pressure are calculated.  The ratio (E/Es)*100
53;; gives the relative humidity
54
55e=es0*10^(7.5*td/(237.7+td))
56es=es0*10^(7.5*t/(237.7+T))
57
58rh=(e/es)*100
59;;;print, 'RH = ',rh
60
61return, rh
62
63end
Note: See TracBrowser for help on using the repository browser.