/[lmdze]/trunk/libf/phylmd/suphec.f90
ViewVC logotype

Diff of /trunk/libf/phylmd/suphec.f90

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 17 by guez, Tue Aug 5 13:31:32 2008 UTC revision 47 by guez, Fri Jul 1 15:00:48 2011 UTC
# Line 2  module suphec_m Line 2  module suphec_m
2    
3    implicit none    implicit none
4    
5      ! A1.0 Fundamental constants
6      REAL RPI
7      real, parameter:: RCLUM = 299792458.
8      real, parameter:: RHPLA = 6.6260755E-34
9      real, parameter:: RKBOL = 1.380658E-23 ! Boltzmann constant, in J K-1
10      real, parameter:: RNAVO = 6.0221367E+23 ! Avogadro number, in mol-1
11    
12      ! A1.1 Astronomical constants
13      REAL RSIYEA, RSIDAY, ROMEGA
14      real, parameter:: RDAY = 86400.
15      real, parameter:: REA = 149597870000.
16      real, parameter:: REPSM = 0.409093
17    
18      ! A1.2 Geoide
19      real, parameter:: RG = 9.80665 ! acceleration of gravity, in m s-2
20      real, parameter:: RA = 6371229.
21    
22      ! A1.3 Radiation
23      REAL RSIGMA
24    
25      ! A1.4 Thermodynamic gas phase
26      REAL, parameter:: R = RNAVO * RKBOL ! ideal gas constant, in J K-1 mol-1
27      real RV, RCPD, RCPV, RCVD, RCVV
28      real, parameter:: RMD = 28.9644 ! molar mass of dry air, in g mol-1
29    
30      real, parameter:: RD = 1000. * R / RMD
31      ! specific ideal gas constant for dry air, in J K-1 kg-1
32    
33      real, parameter:: RMO3 = 47.9942
34      real, parameter:: RMV = 18.0153
35      REAL RKAPPA, RETV
36    
37      ! A1.5, 6 Thermodynamic liquid, solid phases
38      REAL RCW, RCS
39    
40      ! A1.7 Thermodynamic transition of phase
41      REAL RLMLT
42      real, parameter:: RTT = 273.16
43      real, parameter:: RLVTT = 2.5008E+6
44      real, parameter:: RLSTT = 2.8345E+6
45      real, parameter:: RATM = 100000.
46    
47      ! A1.8 Curve of saturation
48      REAL RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS
49      real, parameter:: RESTT = 611.14
50      REAL RALPD, RBETD, RGAMD
51    
52      save
53    
54  contains  contains
55    
56    SUBROUTINE suphec    SUBROUTINE suphec
57    
58      ! From phylmd/suphec.F,v 1.2 2005/06/06 13:16:33      ! From phylmd/suphec.F, v 1.2 2005/06/06 13:16:33
   
59      ! Initialise certaines constantes et parametres physiques.      ! Initialise certaines constantes et parametres physiques.
60    
61      use YOMCST, only: rpi, rclum, rhpla, rkbol, rnavo, rday, rea, repsm, &      !------------------------------------------
          rsiyea, rsiday,romega, rg, ra, r1sa, rsigma, r, rmd, rmo3, rmv, rd, &  
          rv, rcpd, rcvd, rcpv, rcvv, rkappa, retv, rcw, rcs, rtt, rlvtt, &  
          rlstt, rlmlt, ratm, restt, rgamw, rbetw, ralpw, rgams, rbets, ralps, &  
          rgamd, rbetd, ralpd  
     use yoethf  
62    
63      LOGICAL:: firstcall = .TRUE.      PRINT *, 'Call sequence information: suphec'
64    
65      !------------------------------------------      ! 1. DEFINE FUNDAMENTAL CONSTANTS
66    
67      IF (firstcall) THEN      print *, 'Constants of the ICM'
68         PRINT *, 'suphec initialise les constantes du GCM'      RPI = 2.*ASIN(1.)
69         firstcall = .FALSE.      print *, 'Fundamental constants '
70        print '('' PI = '', E13.7, '' -'')', RPI
71         !*       1.    DEFINE FUNDAMENTAL CONSTANTS.      print '('' c = '', E13.7, ''m s-1'')', RCLUM
72        print '('' h = '', E13.7, ''J s'')', RHPLA
73         WRITE(UNIT=6,FMT='(''0*** Constants of the ICM   ***'')')      print '('' K = '', E13.7, ''J K-1'')', RKBOL
74         RPI=2.*ASIN(1.)      print '('' N = '', E13.7, ''mol-1'')', RNAVO
75         RCLUM=299792458.  
76         RHPLA=6.6260755E-34      ! 2. DEFINE ASTRONOMICAL CONSTANTS
77         RKBOL=1.380658E-23  
78         RNAVO=6.0221367E+23      RSIYEA = 365.25*RDAY*2.*RPI/6.283076
79         WRITE(UNIT=6,FMT='('' *** Fundamental constants ***'')')      RSIDAY = RDAY/(1.+RDAY/RSIYEA)
80         WRITE(UNIT=6,FMT='(''           PI = '',E13.7,'' -'')')RPI      ROMEGA = 2.*RPI/RSIDAY
81         WRITE(UNIT=6,FMT='(''            c = '',E13.7,''m s-1'')') &  
82              RCLUM      print *, 'Astronomical constants '
83         WRITE(UNIT=6,FMT='(''            h = '',E13.7,''J s'')') &      print '('' day = '', E13.7, '' s'')', RDAY
84              RHPLA      print '('' half g. axis = '', E13.7, '' m'')', REA
85         WRITE(UNIT=6,FMT='(''            K = '',E13.7,''J K-1'')') &      print '('' mean anomaly = '', E13.7, '' -'')', REPSM
86              RKBOL      print '('' sideral year = '', E13.7, '' s'')', RSIYEA
87         WRITE(UNIT=6,FMT='(''            N = '',E13.7,''mol-1'')') &      print '('' sideral day = '', E13.7, '' s'')', RSIDAY
88              RNAVO      print '('' omega = '', E13.7, '' s-1'')', ROMEGA
89    
90         !*       2.    DEFINE ASTRONOMICAL CONSTANTS.      ! 3. DEFINE GEOIDE.
91    
92         RDAY=86400.      print *, ' Geoide '
93         REA=149597870000.      print '('' Gravity = '', E13.7, '' m s-2'')', RG
94         REPSM=0.409093      print '('' Earth radius = '', E13.7, '' m'')', RA
95    
96         RSIYEA=365.25*RDAY*2.*RPI/6.283076      ! 4. DEFINE RADIATION CONSTANTS.
97         RSIDAY=RDAY/(1.+RDAY/RSIYEA)  
98         ROMEGA=2.*RPI/RSIDAY      rsigma = 2.*rpi**5 * (rkbol/rhpla)**3 * rkbol/rclum/rclum/15.
99        print *, ' Radiation '
100         WRITE(UNIT=6,FMT='('' *** Astronomical constants ***'')')      print '('' Stefan-Bol. = '', E13.7, '' W m-2 K-4'')', RSIGMA
101         WRITE(UNIT=6,FMT='(''          day = '',E13.7,'' s'')')RDAY  
102         WRITE(UNIT=6,FMT='('' half g. axis = '',E13.7,'' m'')')REA      ! 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.
103         WRITE(UNIT=6,FMT='('' mean anomaly = '',E13.7,'' -'')')REPSM  
104         WRITE(UNIT=6,FMT='('' sideral year = '',E13.7,'' s'')')RSIYEA      RV = 1000.*R/RMV
105         WRITE(UNIT=6,FMT='(''  sideral day = '',E13.7,'' s'')')RSIDAY      RCPD = 3.5*RD
106         WRITE(UNIT=6,FMT='(''        omega = '',E13.7,'' s-1'')') &      RCVD = RCPD-RD
107              ROMEGA      RCPV = 4. *RV
108        RCVV = RCPV-RV
109         !*       3.    DEFINE GEOIDE.      RKAPPA = RD/RCPD
110        RETV = RV/RD-1.
111         RG=9.80665      print *, 'Thermodynamic, gas '
112         RA=6371229.      print '('' Perfect gas = '', e13.7)', R
113         R1SA=SNGL(1.D0/DBLE(RA))      print '('' Ozone mass = '', e13.7)', RMO3
114         WRITE(UNIT=6,FMT='('' ***         Geoide         ***'')')      print '('' Vapour mass = '', e13.7)', RMV
115         WRITE(UNIT=6,FMT='(''      Gravity = '',E13.7,'' m s-2'')') &      print '('' Dry air constant = '', e13.7)', RD
116              RG      print '('' Vapour constant = '', e13.7)', RV
117         WRITE(UNIT=6,FMT='('' Earth radius = '',E13.7,'' m'')')RA      print '('' Cpd = '', e13.7)', RCPD
118         WRITE(UNIT=6,FMT='('' Inverse E.R. = '',E13.7,'' m'')')R1SA      print '('' Cvd = '', e13.7)', RCVD
119        print '('' Cpv = '', e13.7)', RCPV
120         !*       4.    DEFINE RADIATION CONSTANTS.      print '('' Cvv = '', e13.7)', RCVV
121        print '('' Rd/Cpd = '', e13.7)', RKAPPA
122         rsigma = 2.*rpi**5 * (rkbol/rhpla)**3 * rkbol/rclum/rclum/15.      print '('' Rv/Rd-1 = '', e13.7)', RETV
123         !IM init. dans conf_phys.F90   RI0=1365.  
124         WRITE(UNIT=6,FMT='('' ***        Radiation       ***'')')      ! 6. DEFINE THERMODYNAMIC CONSTANTS, LIQUID PHASE.
125         WRITE(UNIT=6,FMT='('' Stefan-Bol.  = '',E13.7,'' W m-2 K-4'')')  RSIGMA  
126        RCW = RCPV
127         !*       5.    DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.      print *, 'Thermodynamic, liquid '
128        print '('' Cw = '', E13.7)', RCW
129         R=RNAVO*RKBOL  
130         RMD=28.9644      ! 7. DEFINE THERMODYNAMIC CONSTANTS, SOLID PHASE.
131         RMO3=47.9942  
132         RMV=18.0153      RCS = RCPV
133         RD=1000.*R/RMD      print *, 'thermodynamic, solid'
134         RV=1000.*R/RMV      print '('' Cs = '', E13.7)', RCS
135         RCPD=3.5*RD  
136         RCVD=RCPD-RD      ! 8. DEFINE THERMODYNAMIC CONSTANTS, TRANSITION OF PHASE.
137         RCPV=4. *RV  
138         RCVV=RCPV-RV      RLMLT = RLSTT-RLVTT
139         RKAPPA=RD/RCPD      print *, 'Thermodynamic, trans. '
140         RETV=RV/RD-1.      print '('' Fusion point = '', E13.7)', RTT
141         WRITE(UNIT=6,FMT='('' *** Thermodynamic, gas     ***'')')      print '('' RLvTt = '', E13.7)', RLVTT
142         WRITE(UNIT=6,FMT='('' Perfect gas  = '',e13.7)') R      print '('' RLsTt = '', E13.7)', RLSTT
143         WRITE(UNIT=6,FMT='('' Dry air mass = '',e13.7)') RMD      print '('' RLMlt = '', E13.7)', RLMLT
144         WRITE(UNIT=6,FMT='('' Ozone   mass = '',e13.7)') RMO3      print '('' Normal press. = '', E13.7)', RATM
145         WRITE(UNIT=6,FMT='('' Vapour  mass = '',e13.7)') RMV  
146         WRITE(UNIT=6,FMT='('' Dry air cst. = '',e13.7)') RD      ! 9. SATURATED VAPOUR PRESSURE.
147         WRITE(UNIT=6,FMT='('' Vapour  cst. = '',e13.7)') RV  
148         WRITE(UNIT=6,FMT='(''         Cpd  = '',e13.7)') RCPD      RGAMW = (RCW-RCPV)/RV
149         WRITE(UNIT=6,FMT='(''         Cvd  = '',e13.7)') RCVD      RBETW = RLVTT/RV+RGAMW*RTT
150         WRITE(UNIT=6,FMT='(''         Cpv  = '',e13.7)') RCPV      RALPW = LOG(RESTT)+RBETW/RTT+RGAMW*LOG(RTT)
151         WRITE(UNIT=6,FMT='(''         Cvv  = '',e13.7)') RCVV      RGAMS = (RCS-RCPV)/RV
152         WRITE(UNIT=6,FMT='(''      Rd/Cpd  = '',e13.7)') RKAPPA      RBETS = RLSTT/RV+RGAMS*RTT
153         WRITE(UNIT=6,FMT='(''     Rv/Rd-1  = '',e13.7)') RETV      RALPS = LOG(RESTT)+RBETS/RTT+RGAMS*LOG(RTT)
154        RGAMD = RGAMS-RGAMW
155         !*       6.    DEFINE THERMODYNAMIC CONSTANTS, LIQUID PHASE.      RBETD = RBETS-RBETW
156        RALPD = RALPS-RALPW
        RCW=RCPV  
        WRITE(UNIT=6,FMT='('' *** Thermodynamic, liquid  ***'')')  
        WRITE(UNIT=6,FMT='(''         Cw   = '',E13.7)') RCW  
   
        !*       7.    DEFINE THERMODYNAMIC CONSTANTS, SOLID PHASE.  
   
        RCS=RCPV  
        WRITE(UNIT=6,FMT='('' *** thermodynamic, solid   ***'')')  
        WRITE(UNIT=6,FMT='(''         Cs   = '',E13.7)') RCS  
   
        !*       8.    DEFINE THERMODYNAMIC CONSTANTS, TRANSITION OF PHASE.  
   
        RTT=273.16  
        RLVTT=2.5008E+6  
        RLSTT=2.8345E+6  
        RLMLT=RLSTT-RLVTT  
        RATM=100000.  
        WRITE(UNIT=6,FMT='('' *** Thermodynamic, trans.  ***'')')  
        WRITE(UNIT=6,FMT='('' Fusion point  = '',E13.7)') RTT  
        WRITE(UNIT=6,FMT='(''        RLvTt  = '',E13.7)') RLVTT  
        WRITE(UNIT=6,FMT='(''        RLsTt  = '',E13.7)') RLSTT  
        WRITE(UNIT=6,FMT='(''        RLMlt  = '',E13.7)') RLMLT  
        WRITE(UNIT=6,FMT='('' Normal press. = '',E13.7)') RATM  
        WRITE(UNIT=6,FMT='('' Latent heat :  '')')  
   
        !*       9.    SATURATED VAPOUR PRESSURE.  
   
        RESTT=611.14  
        RGAMW=(RCW-RCPV)/RV  
        RBETW=RLVTT/RV+RGAMW*RTT  
        RALPW=LOG(RESTT)+RBETW/RTT+RGAMW*LOG(RTT)  
        RGAMS=(RCS-RCPV)/RV  
        RBETS=RLSTT/RV+RGAMS*RTT  
        RALPS=LOG(RESTT)+RBETS/RTT+RGAMS*LOG(RTT)  
        RGAMD=RGAMS-RGAMW  
        RBETD=RBETS-RBETW  
        RALPD=RALPS-RALPW  
   
        ! calculer les constantes pour les fonctions thermodynamiques  
   
        RVTMP2=RCPV/RCPD-1.  
        RHOH2O=RATM/100.  
        R2ES=RESTT*RD/RV  
        R3LES=17.269  
        R3IES=21.875  
        R4LES=35.86  
        R4IES=7.66  
        R5LES=R3LES*(RTT-R4LES)  
        R5IES=R3IES*(RTT-R4IES)  
     ELSE  
        PRINT *, 'suphec DEJA APPELE '  
     ENDIF  
157    
158    END SUBROUTINE suphec    END SUBROUTINE suphec
159    
Legend:
Removed from v.17  
changed lines
  Added in v.47
  ViewVC Help
Powered by ViewVC 1.1.21