24 |
IF (firstcall) THEN |
IF (firstcall) THEN |
25 |
PRINT *, 'suphec initialise les constantes du GCM' |
PRINT *, 'suphec initialise les constantes du GCM' |
26 |
firstcall = .FALSE. |
firstcall = .FALSE. |
27 |
! |
|
28 |
!* 1. DEFINE FUNDAMENTAL CONSTANTS. |
!* 1. DEFINE FUNDAMENTAL CONSTANTS. |
29 |
! |
|
30 |
WRITE(UNIT=6,FMT='(''0*** Constants of the ICM ***'')') |
WRITE(UNIT=6,FMT='(''0*** Constants of the ICM ***'')') |
31 |
RPI=2.*ASIN(1.) |
RPI=2.*ASIN(1.) |
32 |
RCLUM=299792458. |
RCLUM=299792458. |
43 |
RKBOL |
RKBOL |
44 |
WRITE(UNIT=6,FMT='('' N = '',E13.7,''mol-1'')') & |
WRITE(UNIT=6,FMT='('' N = '',E13.7,''mol-1'')') & |
45 |
RNAVO |
RNAVO |
46 |
! |
|
|
! |
|
47 |
!* 2. DEFINE ASTRONOMICAL CONSTANTS. |
!* 2. DEFINE ASTRONOMICAL CONSTANTS. |
48 |
! |
|
49 |
RDAY=86400. |
RDAY=86400. |
50 |
REA=149597870000. |
REA=149597870000. |
51 |
REPSM=0.409093 |
REPSM=0.409093 |
52 |
! |
|
53 |
RSIYEA=365.25*RDAY*2.*RPI/6.283076 |
RSIYEA=365.25*RDAY*2.*RPI/6.283076 |
54 |
RSIDAY=RDAY/(1.+RDAY/RSIYEA) |
RSIDAY=RDAY/(1.+RDAY/RSIYEA) |
55 |
ROMEGA=2.*RPI/RSIDAY |
ROMEGA=2.*RPI/RSIDAY |
56 |
! |
|
57 |
WRITE(UNIT=6,FMT='('' *** Astronomical constants ***'')') |
WRITE(UNIT=6,FMT='('' *** Astronomical constants ***'')') |
58 |
WRITE(UNIT=6,FMT='('' day = '',E13.7,'' s'')')RDAY |
WRITE(UNIT=6,FMT='('' day = '',E13.7,'' s'')')RDAY |
59 |
WRITE(UNIT=6,FMT='('' half g. axis = '',E13.7,'' m'')')REA |
WRITE(UNIT=6,FMT='('' half g. axis = '',E13.7,'' m'')')REA |
62 |
WRITE(UNIT=6,FMT='('' sideral day = '',E13.7,'' s'')')RSIDAY |
WRITE(UNIT=6,FMT='('' sideral day = '',E13.7,'' s'')')RSIDAY |
63 |
WRITE(UNIT=6,FMT='('' omega = '',E13.7,'' s-1'')') & |
WRITE(UNIT=6,FMT='('' omega = '',E13.7,'' s-1'')') & |
64 |
ROMEGA |
ROMEGA |
65 |
! |
|
66 |
!* 3. DEFINE GEOIDE. |
!* 3. DEFINE GEOIDE. |
67 |
! |
|
68 |
RG=9.80665 |
RG=9.80665 |
69 |
RA=6371229. |
RA=6371229. |
70 |
R1SA=SNGL(1.D0/DBLE(RA)) |
R1SA=SNGL(1.D0/DBLE(RA)) |
73 |
RG |
RG |
74 |
WRITE(UNIT=6,FMT='('' Earth radius = '',E13.7,'' m'')')RA |
WRITE(UNIT=6,FMT='('' Earth radius = '',E13.7,'' m'')')RA |
75 |
WRITE(UNIT=6,FMT='('' Inverse E.R. = '',E13.7,'' m'')')R1SA |
WRITE(UNIT=6,FMT='('' Inverse E.R. = '',E13.7,'' m'')')R1SA |
76 |
! |
|
|
! |
|
77 |
!* 4. DEFINE RADIATION CONSTANTS. |
!* 4. DEFINE RADIATION CONSTANTS. |
78 |
! |
|
|
! z.x.li RSIGMA=2. * RPI**5 * RKBOL**4 /(15.* RCLUM**2 * RHPLA**3) |
|
79 |
rsigma = 2.*rpi**5 * (rkbol/rhpla)**3 * rkbol/rclum/rclum/15. |
rsigma = 2.*rpi**5 * (rkbol/rhpla)**3 * rkbol/rclum/rclum/15. |
80 |
!IM init. dans conf_phys.F90 RI0=1365. |
!IM init. dans conf_phys.F90 RI0=1365. |
81 |
WRITE(UNIT=6,FMT='('' *** Radiation ***'')') |
WRITE(UNIT=6,FMT='('' *** Radiation ***'')') |
82 |
WRITE(UNIT=6,FMT='('' Stefan-Bol. = '',E13.7,'' W m-2 K-4'')') RSIGMA |
WRITE(UNIT=6,FMT='('' Stefan-Bol. = '',E13.7,'' W m-2 K-4'')') RSIGMA |
83 |
|
|
|
! |
|
84 |
!* 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE. |
!* 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE. |
85 |
! |
|
86 |
R=RNAVO*RKBOL |
R=RNAVO*RKBOL |
87 |
RMD=28.9644 |
RMD=28.9644 |
88 |
RMO3=47.9942 |
RMO3=47.9942 |
108 |
WRITE(UNIT=6,FMT='('' Cvv = '',e13.7)') RCVV |
WRITE(UNIT=6,FMT='('' Cvv = '',e13.7)') RCVV |
109 |
WRITE(UNIT=6,FMT='('' Rd/Cpd = '',e13.7)') RKAPPA |
WRITE(UNIT=6,FMT='('' Rd/Cpd = '',e13.7)') RKAPPA |
110 |
WRITE(UNIT=6,FMT='('' Rv/Rd-1 = '',e13.7)') RETV |
WRITE(UNIT=6,FMT='('' Rv/Rd-1 = '',e13.7)') RETV |
111 |
! |
|
|
! |
|
112 |
!* 6. DEFINE THERMODYNAMIC CONSTANTS, LIQUID PHASE. |
!* 6. DEFINE THERMODYNAMIC CONSTANTS, LIQUID PHASE. |
113 |
! |
|
114 |
RCW=RCPV |
RCW=RCPV |
115 |
WRITE(UNIT=6,FMT='('' *** Thermodynamic, liquid ***'')') |
WRITE(UNIT=6,FMT='('' *** Thermodynamic, liquid ***'')') |
116 |
WRITE(UNIT=6,FMT='('' Cw = '',E13.7)') RCW |
WRITE(UNIT=6,FMT='('' Cw = '',E13.7)') RCW |
117 |
! |
|
|
! |
|
118 |
!* 7. DEFINE THERMODYNAMIC CONSTANTS, SOLID PHASE. |
!* 7. DEFINE THERMODYNAMIC CONSTANTS, SOLID PHASE. |
119 |
! |
|
120 |
RCS=RCPV |
RCS=RCPV |
121 |
WRITE(UNIT=6,FMT='('' *** thermodynamic, solid ***'')') |
WRITE(UNIT=6,FMT='('' *** thermodynamic, solid ***'')') |
122 |
WRITE(UNIT=6,FMT='('' Cs = '',E13.7)') RCS |
WRITE(UNIT=6,FMT='('' Cs = '',E13.7)') RCS |
123 |
! |
|
|
! |
|
124 |
!* 8. DEFINE THERMODYNAMIC CONSTANTS, TRANSITION OF PHASE. |
!* 8. DEFINE THERMODYNAMIC CONSTANTS, TRANSITION OF PHASE. |
125 |
! |
|
126 |
RTT=273.16 |
RTT=273.16 |
127 |
RLVTT=2.5008E+6 |
RLVTT=2.5008E+6 |
128 |
RLSTT=2.8345E+6 |
RLSTT=2.8345E+6 |
135 |
WRITE(UNIT=6,FMT='('' RLMlt = '',E13.7)') RLMLT |
WRITE(UNIT=6,FMT='('' RLMlt = '',E13.7)') RLMLT |
136 |
WRITE(UNIT=6,FMT='('' Normal press. = '',E13.7)') RATM |
WRITE(UNIT=6,FMT='('' Normal press. = '',E13.7)') RATM |
137 |
WRITE(UNIT=6,FMT='('' Latent heat : '')') |
WRITE(UNIT=6,FMT='('' Latent heat : '')') |
138 |
! |
|
|
! |
|
139 |
!* 9. SATURATED VAPOUR PRESSURE. |
!* 9. SATURATED VAPOUR PRESSURE. |
140 |
! |
|
141 |
RESTT=611.14 |
RESTT=611.14 |
142 |
RGAMW=(RCW-RCPV)/RV |
RGAMW=(RCW-RCPV)/RV |
143 |
RBETW=RLVTT/RV+RGAMW*RTT |
RBETW=RLVTT/RV+RGAMW*RTT |
148 |
RGAMD=RGAMS-RGAMW |
RGAMD=RGAMS-RGAMW |
149 |
RBETD=RBETS-RBETW |
RBETD=RBETS-RBETW |
150 |
RALPD=RALPS-RALPW |
RALPD=RALPS-RALPW |
151 |
! |
|
|
! |
|
152 |
! calculer les constantes pour les fonctions thermodynamiques |
! calculer les constantes pour les fonctions thermodynamiques |
153 |
! |
|
154 |
RVTMP2=RCPV/RCPD-1. |
RVTMP2=RCPV/RCPD-1. |
155 |
RHOH2O=RATM/100. |
RHOH2O=RATM/100. |
156 |
R2ES=RESTT*RD/RV |
R2ES=RESTT*RD/RV |
160 |
R4IES=7.66 |
R4IES=7.66 |
161 |
R5LES=R3LES*(RTT-R4LES) |
R5LES=R3LES*(RTT-R4LES) |
162 |
R5IES=R3IES*(RTT-R4IES) |
R5IES=R3IES*(RTT-R4IES) |
|
|
|
163 |
ELSE |
ELSE |
164 |
PRINT *, 'suphec DEJA APPELE ' |
PRINT *, 'suphec DEJA APPELE ' |
165 |
ENDIF |
ENDIF |