/[lmdze]/trunk/Sources/phylmd/suphec.f
ViewVC logotype

Annotation of /trunk/Sources/phylmd/suphec.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 198 - (hide annotations)
Tue May 31 16:17:35 2016 UTC (8 years ago) by guez
File size: 4194 byte(s) 
Removed variables nk1 and nk in cv_driver and below. These arrays were
just equal to the constant minorig. (This is also the case in LMDZ.)

In cv_thermo, removed some variables which were copies of variables of
suphec_m. Changed some variables to named constants.
1 guez 3 module suphec_m
2    
3 guez 169 use nr_util, only: pi
4    
5 guez 3 implicit none
6    
7 guez 38 ! A1.0 Fundamental constants
8 guez 169 real, parameter:: RCLUM = 299792458. ! speed of light, m s-1
9     real, parameter:: RHPLA = 6.6260755E-34 ! Planck constant, J s
10 guez 70 real, parameter:: KBOL = 1.380658E-23 ! Boltzmann constant, in J K-1
11     real, parameter:: NAVO = 6.0221367E23 ! Avogadro number, in mol-1
12 guez 38
13     ! A1.1 Astronomical constants
14 guez 171
15     REAL ROMEGA
16 guez 40 real, parameter:: RDAY = 86400.
17 guez 38
18 guez 171 REAL, parameter:: RSIYEA = 365.25 * RDAY * 2. * PI / 6.283076
19     ! sideral year, in s
20    
21     REAL, parameter:: RSIDAY = RDAY / (1. + RDAY / RSIYEA) ! sideral day, in s
22    
23 guez 38 ! A1.2 Geoide
24 guez 47 real, parameter:: RG = 9.80665 ! acceleration of gravity, in m s-2
25 guez 40 real, parameter:: RA = 6371229.
26 guez 38
27     ! A1.3 Radiation
28 guez 169 REAL, parameter:: rsigma = 2. * pi**5 * (kbol / rhpla)**3 * kbol / rclum**2 &
29     / 15.
30 guez 38
31     ! A1.4 Thermodynamic gas phase
32 guez 70 REAL, parameter:: R = NAVO * KBOL ! ideal gas constant, in J K-1 mol-1
33     real, parameter:: MV = 18.0153 ! molar mass of water, in g mol-1
34 guez 47
35 guez 70 real, parameter:: RV = 1e3 * R / MV
36     ! specific ideal gas constant for water vapor, in J K-1 kg-1
37     ! (factor 1e3: conversion from g to kg)
38    
39     real, parameter:: MD = 28.9644 ! molar mass of dry air, in g mol-1
40    
41     real, parameter:: RD = 1e3 * R / MD
42 guez 47 ! specific ideal gas constant for dry air, in J K-1 kg-1
43 guez 70 ! (factor 1e3: conversion from g to kg)
44 guez 47
45 guez 198 real, parameter:: RCPV = 4. * RV ! Cpv, gas phase
46 guez 70
47 guez 198 real, save:: RCVD, RCVV
48    
49 guez 53 real, parameter:: RCPD = 7. / 2 * RD
50     ! specific heat capacity for dry air, in J K-1 kg-1
51    
52 guez 40 real, parameter:: RMO3 = 47.9942
53 guez 53 REAL, parameter:: RKAPPA = RD/RCPD
54 guez 169 real, save:: RETV
55 guez 38
56 guez 40 ! A1.5, 6 Thermodynamic liquid, solid phases
57 guez 38
58 guez 198 REAL, parameter:: RCW = RCPV ! LIQUID PHASE Cw
59     real, save:: RCS
60    
61 guez 38 ! A1.7 Thermodynamic transition of phase
62 guez 169 REAL, save:: RLMLT
63 guez 40 real, parameter:: RTT = 273.16
64     real, parameter:: RLVTT = 2.5008E+6
65     real, parameter:: RLSTT = 2.8345E+6
66 guez 88 real, parameter:: RATM = 1e5
67 guez 38
68     ! A1.8 Curve of saturation
69 guez 169 REAL, save:: RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS
70 guez 40 real, parameter:: RESTT = 611.14
71 guez 169 REAL, save:: RALPD, RBETD, RGAMD
72 guez 38
73 guez 169 private pi
74 guez 38
75 guez 3 contains
76    
77     SUBROUTINE suphec
78    
79 guez 53 ! From phylmd/suphec.F, version 1.2 2005/06/06 13:16:33
80     ! Initialise certaines constantes et certains paramètres physiques.
81 guez 3
82     !------------------------------------------
83    
84 guez 37 PRINT *, 'Call sequence information: suphec'
85 guez 17
86 guez 37 ! 2. DEFINE ASTRONOMICAL CONSTANTS
87 guez 17
88 guez 169 ROMEGA = 2.*PI/RSIDAY
89 guez 17
90 guez 37 print *, 'Astronomical constants '
91 guez 40 print '('' omega = '', E13.7, '' s-1'')', ROMEGA
92 guez 17
93 guez 40 ! 3. DEFINE GEOIDE.
94 guez 17
95 guez 40 print *, ' Geoide '
96     print '('' Gravity = '', E13.7, '' m s-2'')', RG
97     print '('' Earth radius = '', E13.7, '' m'')', RA
98 guez 17
99 guez 40 ! 4. DEFINE RADIATION CONSTANTS.
100 guez 17
101 guez 40 print *, ' Radiation '
102     print '('' Stefan-Bol. = '', E13.7, '' W m-2 K-4'')', RSIGMA
103 guez 3
104 guez 40 ! 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.
105 guez 17
106 guez 40 RCVD = RCPD-RD
107     RCVV = RCPV-RV
108 guez 70 RETV = RV / RD - 1.
109     print *, 'Thermodynamics, gas'
110 guez 40 print '('' Ozone mass = '', e13.7)', RMO3
111 guez 70 print *, "rd = ", RD, "J K-1 kg-1"
112     print *, "rv = ", RV, "J K-1 kg-1"
113 guez 40 print '('' Cpd = '', e13.7)', RCPD
114     print '('' Cvd = '', e13.7)', RCVD
115     print '('' Cvv = '', e13.7)', RCVV
116     print '('' Rd/Cpd = '', e13.7)', RKAPPA
117 guez 70 print '('' Rv / Rd - 1 = '', e13.7)', RETV
118 guez 17
119 guez 40 ! 7. DEFINE THERMODYNAMIC CONSTANTS, SOLID PHASE.
120 guez 17
121 guez 40 RCS = RCPV
122 guez 37 print *, 'thermodynamic, solid'
123 guez 40 print '('' Cs = '', E13.7)', RCS
124 guez 17
125 guez 40 ! 8. DEFINE THERMODYNAMIC CONSTANTS, TRANSITION OF PHASE.
126 guez 17
127 guez 40 RLMLT = RLSTT-RLVTT
128 guez 88 print *, 'Thermodynamic, transition of phase:'
129 guez 40 print '('' Fusion point = '', E13.7)', RTT
130     print '('' RLvTt = '', E13.7)', RLVTT
131     print '('' RLsTt = '', E13.7)', RLSTT
132     print '('' RLMlt = '', E13.7)', RLMLT
133 guez 88 print '('' Normal pressure = '', E13.7)', RATM
134 guez 17
135 guez 40 ! 9. SATURATED VAPOUR PRESSURE.
136 guez 17
137 guez 40 RGAMW = (RCW-RCPV)/RV
138     RBETW = RLVTT/RV+RGAMW*RTT
139     RALPW = LOG(RESTT)+RBETW/RTT+RGAMW*LOG(RTT)
140     RGAMS = (RCS-RCPV)/RV
141     RBETS = RLSTT/RV+RGAMS*RTT
142     RALPS = LOG(RESTT)+RBETS/RTT+RGAMS*LOG(RTT)
143     RGAMD = RGAMS-RGAMW
144     RBETD = RBETS-RBETW
145     RALPD = RALPS-RALPW
146 guez 17
147 guez 3 END SUBROUTINE suphec
148    
149     end module suphec_m
  ViewVC Help
Powered by ViewVC 1.1.21