source: LMDZ.3.3/trunk/libf/dyn3d/q_sat.F @ 2

Last change on this file since 2 was 2, checked in by lmdz, 25 years ago

Initial revision
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.5 KB
Line 
1      FUNCTION q_sat(kelvin, millibar)
2c
3      IMPLICIT none
4c======================================================================
5c Autheur(s): Z.X. Li (LMD/CNRS)
6c Objet: calculer la vapeur d'eau saturante (formule Centre Euro.)
7c======================================================================
8c Arguments:
9c kelvin---input-R: temperature en Kelvin
10c millibar--input-R: pression en mb
11c
12c q_sat----output-R: vapeur d'eau saturante en kg/kg
13c======================================================================
14c
15      REAL q_sat, kelvin, millibar
16c
17      REAL r2es
18      PARAMETER (r2es=611.14 *18.0153/28.9644)
19c
20      REAL r3les, r3ies, r3es
21      PARAMETER (R3LES=17.269)
22      PARAMETER (R3IES=21.875)
23c
24      REAL r4les, r4ies, r4es
25      PARAMETER (R4LES=35.86)
26      PARAMETER (R4IES=7.66)
27c
28      REAL rtt
29      PARAMETER (rtt=273.16)
30c
31      REAL retv
32      PARAMETER (retv=28.9644/18.0153 - 1.0)
33c
34      REAL zqsat
35      REAL temp, pres
36C     ------------------------------------------------------------------
37c
38c
39      temp = kelvin
40      pres = millibar * 100.0
41c      write(*,*)'kelvin,millibar=',kelvin,millibar
42c      write(*,*)'temp,pres=',temp,pres
43c
44      IF (temp .LE. rtt) THEN
45         r3es = r3ies
46         r4es = r4ies
47      ELSE
48         r3es = r3les
49         r4es = r4les
50      ENDIF
51c
52      zqsat=r2es/pres * EXP ( r3es*(temp-rtt) / (temp-r4es) )
53      zqsat=MIN(0.5,ZQSAT)
54      zqsat=zqsat/(1.-retv  *zqsat)
55c
56      q_sat = zqsat
57c
58      RETURN
59      END
Note: See TracBrowser for help on using the repository browser.