Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

suphel.F @ 4400

Last change on this file since 4400 was 879, checked in by Laurent Fairhead, 17 years ago
Suite de la bascule vers une physique avec thermiques, nouvelle convection, poche froide ... LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 7.1 KB

Line
1	!
2	! $Header$
3	!
4	SUBROUTINE suphel
5	C
6	#include "YOMCST.h"
7	#include "YOETHF.h"
8	cIM cf. JLD
9	LOGICAL firstcall
10	SAVE firstcall
11	c$OMP THREADPRIVATE(firstcall)
12	DATA firstcall /.TRUE./
13
14	IF (firstcall) THEN
15	PRINT*, 'suphel initialise les constantes du GCM'
16	firstcall = .FALSE.
17	ELSE
18	PRINT*, 'suphel DEJA APPELE '
19	RETURN
20	ENDIF
21	C -----------------------------------------------------------------
22	C
23	C* 1. DEFINE FUNDAMENTAL CONSTANTS.
24	C -----------------------------
25	C
26	WRITE(UNIT=6,FMT='(''0* Constants of the ICM *'')')
27	RPI=2.*ASIN(1.)
28	RCLUM=299792458.
29	RHPLA=6.6260755E-34
30	RKBOL=1.380658E-23
31	RNAVO=6.0221367E+23
32	WRITE(UNIT=6,FMT='('' * Fundamental constants *'')')
33	WRITE(UNIT=6,FMT='('' PI = '',E13.7,'' -'')')RPI
34	WRITE(UNIT=6,FMT='('' c = '',E13.7,''m s-1'')')
35	S RCLUM
36	WRITE(UNIT=6,FMT='('' h = '',E13.7,''J s'')')
37	S RHPLA
38	WRITE(UNIT=6,FMT='('' K = '',E13.7,''J K-1'')')
39	S RKBOL
40	WRITE(UNIT=6,FMT='('' N = '',E13.7,''mol-1'')')
41	S RNAVO
42	C
43	C ----------------------------------------------------------------
44	C
45	C* 2. DEFINE ASTRONOMICAL CONSTANTS.
46	C ------------------------------
47	C
48	RDAY=86400.
49	REA=149597870000.
50	REPSM=0.409093
51	C
52	RSIYEA=365.25RDAY2.*RPI/6.283076
53	RSIDAY=RDAY/(1.+RDAY/RSIYEA)
54	ROMEGA=2.*RPI/RSIDAY
55	c
56	c exp1 R_ecc = 0.05
57	c exp1 R_peri = 102.04
58	c exp1 R_incl = 22.5
59	c exp1 print*, 'Parametres orbitaux modifies'
60	c ref R_ecc = 0.016724
61	c ref R_peri = 102.04
62	c ref R_incl = 23.5
63	c
64	cIM 161002 : pour avoir les ctes AMIP II
65	cIM 161002 R_ecc = 0.016724
66	cIM 161002 R_peri = 102.04
67	cIM 161002 R_incl = 23.5
68	cIM on mets R_ecc, R_peri, R_incl dans conf_phys.F90
69	c R_ecc = 0.016715
70	c R_peri = 102.7
71	c R_incl = 23.441
72	c
73	WRITE(UNIT=6,FMT='('' * Astronomical constants *'')')
74	WRITE(UNIT=6,FMT='('' day = '',E13.7,'' s'')')RDAY
75	WRITE(UNIT=6,FMT='('' half g. axis = '',E13.7,'' m'')')REA
76	WRITE(UNIT=6,FMT='('' mean anomaly = '',E13.7,'' -'')')REPSM
77	WRITE(UNIT=6,FMT='('' sideral year = '',E13.7,'' s'')')RSIYEA
78	WRITE(UNIT=6,FMT='('' sideral day = '',E13.7,'' s'')')RSIDAY
79	WRITE(UNIT=6,FMT='('' omega = '',E13.7,'' s-1'')')
80	S ROMEGA
81	c write(unit=6,fmt='('' excentricite = '',e13.7,''-'')')R_ecc
82	c write(unit=6,fmt='('' equinoxe = '',e13.7,''-'')')R_peri
83	c write(unit=6,fmt='('' inclinaison = '',e13.7,''-'')')R_incl
84	C
85	C ------------------------------------------------------------------
86	C
87	C* 3. DEFINE GEOIDE.
88	C --------------
89	C
90	RG=9.80665
91	RA=6371229.
92	R1SA=SNGL(1.D0/DBLE(RA))
93	WRITE(UNIT=6,FMT='('' * Geoide *'')')
94	WRITE(UNIT=6,FMT='('' Gravity = '',E13.7,'' m s-2'')')
95	S RG
96	WRITE(UNIT=6,FMT='('' Earth radius = '',E13.7,'' m'')')RA
97	WRITE(UNIT=6,FMT='('' Inverse E.R. = '',E13.7,'' m'')')R1SA
98	C
99	C -----------------------------------------------------------------
100	C
101	C* 4. DEFINE RADIATION CONSTANTS.
102	C ---------------------------
103	C
104	c z.x.li RSIGMA=2. * RPI*5 RKBOL*4 /(15. RCLUM*2 RHPLA**3)
105	rsigma = 2.rpi5 (rkbol/rhpla)*3 rkbol/rclum/rclum/15.
106	cIM init. dans conf_phys.F90 RI0=1365.
107	WRITE(UNIT=6,FMT='('' * Radiation *'')')
108	WRITE(UNIT=6,FMT='('' Stefan-Bol. = '',E13.7,'' W m-2 K-4''
109	S )') RSIGMA
110	cIM init. dans conf_phys.F90 WRITE(UNIT=6,FMT='('' Solar const. = '',E13.7,'' W m-2'')')
111	cIM init. dans conf_phys.F90 S RI0
112	C
113	C -----------------------------------------------------------------
114	C
115	C* 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.
116	C ------------------------------------------
117	C
118	R=RNAVO*RKBOL
119	RMD=28.9644
120	RMO3=47.9942
121	RMV=18.0153
122	RD=1000.*R/RMD
123	RV=1000.*R/RMV
124	RCPD=3.5*RD
125	RCVD=RCPD-RD
126	RCPV=4. *RV
127	RCVV=RCPV-RV
128	RKAPPA=RD/RCPD
129	RETV=RV/RD-1.
130	WRITE(UNIT=6,FMT='('' * Thermodynamic, gas *'')')
131	WRITE(UNIT=6,FMT='('' Perfect gas = '',e13.7)') R
132	WRITE(UNIT=6,FMT='('' Dry air mass = '',e13.7)') RMD
133	WRITE(UNIT=6,FMT='('' Ozone mass = '',e13.7)') RMO3
134	WRITE(UNIT=6,FMT='('' Vapour mass = '',e13.7)') RMV
135	WRITE(UNIT=6,FMT='('' Dry air cst. = '',e13.7)') RD
136	WRITE(UNIT=6,FMT='('' Vapour cst. = '',e13.7)') RV
137	WRITE(UNIT=6,FMT='('' Cpd = '',e13.7)') RCPD
138	WRITE(UNIT=6,FMT='('' Cvd = '',e13.7)') RCVD
139	WRITE(UNIT=6,FMT='('' Cpv = '',e13.7)') RCPV
140	WRITE(UNIT=6,FMT='('' Cvv = '',e13.7)') RCVV
141	WRITE(UNIT=6,FMT='('' Rd/Cpd = '',e13.7)') RKAPPA
142	WRITE(UNIT=6,FMT='('' Rv/Rd-1 = '',e13.7)') RETV
143	C
144	C ----------------------------------------------------------------
145	C
146	C* 6. DEFINE THERMODYNAMIC CONSTANTS, LIQUID PHASE.
147	C ---------------------------------------------
148	C
149	RCW=RCPV
150	WRITE(UNIT=6,FMT='('' * Thermodynamic, liquid *'')')
151	WRITE(UNIT=6,FMT='('' Cw = '',E13.7)') RCW
152	C
153	C ----------------------------------------------------------------
154	C
155	C* 7. DEFINE THERMODYNAMIC CONSTANTS, SOLID PHASE.
156	C --------------------------------------------
157	C
158	RCS=RCPV
159	WRITE(UNIT=6,FMT='('' * thermodynamic, solid *'')')
160	WRITE(UNIT=6,FMT='('' Cs = '',E13.7)') RCS
161	C
162	C ----------------------------------------------------------------
163	C
164	C* 8. DEFINE THERMODYNAMIC CONSTANTS, TRANSITION OF PHASE.
165	C ----------------------------------------------------
166	C
167	RTT=273.16
168	RLVTT=2.5008E+6
169	RLSTT=2.8345E+6
170	RLMLT=RLSTT-RLVTT
171	RATM=100000.
172	WRITE(UNIT=6,FMT='('' * Thermodynamic, trans. *'')')
173	WRITE(UNIT=6,FMT='('' Fusion point = '',E13.7)') RTT
174	WRITE(UNIT=6,FMT='('' RLvTt = '',E13.7)') RLVTT
175	WRITE(UNIT=6,FMT='('' RLsTt = '',E13.7)') RLSTT
176	WRITE(UNIT=6,FMT='('' RLMlt = '',E13.7)') RLMLT
177	WRITE(UNIT=6,FMT='('' Normal press. = '',E13.7)') RATM
178	WRITE(UNIT=6,FMT='('' Latent heat : '')')
179	C
180	C ----------------------------------------------------------------
181	C
182	C* 9. SATURATED VAPOUR PRESSURE.
183	C --------------------------
184	C
185	RESTT=611.14
186	RGAMW=(RCW-RCPV)/RV
187	RBETW=RLVTT/RV+RGAMW*RTT
188	RALPW=LOG(RESTT)+RBETW/RTT+RGAMW*LOG(RTT)
189	RGAMS=(RCS-RCPV)/RV
190	RBETS=RLSTT/RV+RGAMS*RTT
191	RALPS=LOG(RESTT)+RBETS/RTT+RGAMS*LOG(RTT)
192	RGAMD=RGAMS-RGAMW
193	RBETD=RBETS-RBETW
194	RALPD=RALPS-RALPW
195	C
196	C ------------------------------------------------------------------
197	c
198	c calculer les constantes pour les fonctions thermodynamiques
199	c
200	RVTMP2=RCPV/RCPD-1.
201	RHOH2O=RATM/100.
202	R2ES=RESTT*RD/RV
203	R3LES=17.269
204	R3IES=21.875
205	R4LES=35.86
206	R4IES=7.66
207	R5LES=R3LES*(RTT-R4LES)
208	R5IES=R3IES*(RTT-R4IES)
209	C
210	RETURN
211	END

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format