Last change
on this file since 16 was
12,
checked in by pinsard, 13 years ago
|
fix for graph in PDF
|
-
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 | ; |
---|
| 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 ${TROPFLUX_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 | ;- |
---|
[5] | 45 | function dew_temp_airt_rel_humidity_jv, t,td |
---|
| 46 | |
---|
| 47 | if (not keyword_set(t)) then t=28 ; air temperature in degC |
---|
| 48 | if (not keyword_set(td)) then td=23 ; dew point temperature in degC |
---|
| 49 | es0=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 | |
---|
| 55 | e=es0*10^(7.5*td/(237.7+td)) |
---|
| 56 | es=es0*10^(7.5*t/(237.7+T)) |
---|
| 57 | |
---|
| 58 | rh=(e/es)*100 |
---|
| 59 | ;;;print, 'RH = ',rh |
---|
| 60 | |
---|
| 61 | return, rh |
---|
| 62 | |
---|
| 63 | end |
---|
Note: See
TracBrowser
for help on using the repository browser.