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

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

Parent Directory Parent Directory | Revision Log Revision Log

Revision 17 - (show annotations)
Tue Aug 5 13:31:32 2008 UTC (15 years, 9 months ago) by guez
File size: 5753 byte(s) 
Created rule for "compare_sampl_*" files in
"Documentation/Manuel_LMDZE.texfol/Graphiques/GNUmakefile".

Extracted "qcheck", "radiornpb", "minmaxqfi" into separate files.

Read pressure coordinate of ozone coefficients once per run instead of
every day.

Added some "intent" attributes.

Added argument "nq" to "ini_histday". Replaced calls to "gr_fi_ecrit"
by calls to "gr_phy_write_2d". "Sigma_O3_Royer" is written to
"histday.nc" only if "nq >= 4". Moved "ini_histrac" to module
"ini_hist".

Compute "zmasse" in "physiq", pass it to "phytrac".

Removed computations of "pftsol*" and "ppsrf*" in "phytrac".

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