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