source: LMDZ.3.3/trunk/libf/phylmd/suphec.F @ 346

Last change on this file since 346 was 2, checked in by lmdz, 25 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 6.4 KB
Line 
1      SUBROUTINE suphec
2C
3#include "YOMCST.h"
4#include "YOETHF.h"
5C      -----------------------------------------------------------------
6C
7C*       1.    DEFINE FUNDAMENTAL CONSTANTS.
8C              -----------------------------
9C
10      WRITE(UNIT=6,FMT='(''0*** Constants of the ICM   ***'')')
11      RPI=2.*ASIN(1.)
12      RCLUM=299792458.
13      RHPLA=6.6260755E-34
14      RKBOL=1.380658E-23
15      RNAVO=6.0221367E+23
16      WRITE(UNIT=6,FMT='('' *** Fundamental constants ***'')')
17      WRITE(UNIT=6,FMT='(''           PI = '',E13.7,'' -'')')RPI
18      WRITE(UNIT=6,FMT='(''            c = '',E13.7,''m s-1'')')
19     S RCLUM
20      WRITE(UNIT=6,FMT='(''            h = '',E13.7,''J s'')')
21     S RHPLA
22      WRITE(UNIT=6,FMT='(''            K = '',E13.7,''J K-1'')')
23     S RKBOL
24      WRITE(UNIT=6,FMT='(''            N = '',E13.7,''mol-1'')')
25     S RNAVO
26C
27C     ----------------------------------------------------------------
28C
29C*       2.    DEFINE ASTRONOMICAL CONSTANTS.
30C              ------------------------------
31C
32      RDAY=86400.
33      REA=149597870000.
34      REPSM=0.409093
35C
36      RSIYEA=365.25*RDAY*2.*RPI/6.283076
37      RSIDAY=RDAY/(1.+RDAY/RSIYEA)
38      ROMEGA=2.*RPI/RSIDAY
39c
40c exp1      R_ecc = 0.05
41c exp1      R_peri = 102.04
42c exp1      R_incl = 22.5
43c exp1      print*, 'Parametres orbitaux modifies'
44c ref      R_ecc = 0.016724
45c ref      R_peri = 102.04
46c ref      R_incl = 23.5
47      R_ecc = 0.016724
48      R_peri = 102.04
49      R_incl = 23.5
50c
51      WRITE(UNIT=6,FMT='('' *** Astronomical constants ***'')')
52      WRITE(UNIT=6,FMT='(''          day = '',E13.7,'' s'')')RDAY
53      WRITE(UNIT=6,FMT='('' half g. axis = '',E13.7,'' m'')')REA
54      WRITE(UNIT=6,FMT='('' mean anomaly = '',E13.7,'' -'')')REPSM
55      WRITE(UNIT=6,FMT='('' sideral year = '',E13.7,'' s'')')RSIYEA
56      WRITE(UNIT=6,FMT='(''  sideral day = '',E13.7,'' s'')')RSIDAY
57      WRITE(UNIT=6,FMT='(''        omega = '',E13.7,'' s-1'')')
58     S                  ROMEGA
59      write(unit=6,fmt='('' excentricite = '',e13.7,''-'')')R_ecc
60      write(unit=6,fmt='(''     equinoxe = '',e13.7,''-'')')R_peri
61      write(unit=6,fmt='(''  inclinaison = '',e13.7,''-'')')R_incl
62C
63C     ------------------------------------------------------------------
64C
65C*       3.    DEFINE GEOIDE.
66C              --------------
67C
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     S      RG
74      WRITE(UNIT=6,FMT='('' Earth radius = '',E13.7,'' m'')')RA
75      WRITE(UNIT=6,FMT='('' Inverse E.R. = '',E13.7,'' m'')')R1SA
76C
77C     -----------------------------------------------------------------
78C
79C*       4.    DEFINE RADIATION CONSTANTS.
80C              ---------------------------
81C
82c z.x.li      RSIGMA=2. * RPI**5 * RKBOL**4 /(15.* RCLUM**2 * RHPLA**3)
83      rsigma = 2.*rpi**5 * (rkbol/rhpla)**3 * rkbol/rclum/rclum/15.
84      RI0=1370.
85      WRITE(UNIT=6,FMT='('' ***        Radiation       ***'')')
86      WRITE(UNIT=6,FMT='('' Stefan-Bol.  = '',E13.7,'' W m-2 K-4''
87     S )')  RSIGMA
88      WRITE(UNIT=6,FMT='('' Solar const. = '',E13.7,'' W m-2'')')
89     S      RI0
90C
91C     -----------------------------------------------------------------
92C
93C*       5.    DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.
94C              ------------------------------------------
95C
96      R=RNAVO*RKBOL
97      RMD=28.9644
98      RMV=18.0153
99      RD=1000.*R/RMD
100      RV=1000.*R/RMV
101      RCPD=3.5*RD
102      RCVD=RCPD-RD
103      RCPV=4. *RV
104      RCVV=RCPV-RV
105      RKAPPA=RD/RCPD
106      RETV=RV/RD-1.
107      WRITE(UNIT=6,FMT='('' *** Thermodynamic, gas     ***'')')
108      WRITE(UNIT=6,FMT='('' Perfect gas  = '',e13.7)') R
109      WRITE(UNIT=6,FMT='('' Dry air mass = '',e13.7)') RMD
110      WRITE(UNIT=6,FMT='('' Vapour  mass = '',e13.7)') RMV
111      WRITE(UNIT=6,FMT='('' Dry air cst. = '',e13.7)') RD
112      WRITE(UNIT=6,FMT='('' Vapour  cst. = '',e13.7)') RV
113      WRITE(UNIT=6,FMT='(''         Cpd  = '',e13.7)') RCPD
114      WRITE(UNIT=6,FMT='(''         Cvd  = '',e13.7)') RCVD
115      WRITE(UNIT=6,FMT='(''         Cpv  = '',e13.7)') RCPV
116      WRITE(UNIT=6,FMT='(''         Cvv  = '',e13.7)') RCVV
117      WRITE(UNIT=6,FMT='(''      Rd/Cpd  = '',e13.7)') RKAPPA
118      WRITE(UNIT=6,FMT='(''     Rv/Rd-1  = '',e13.7)') RETV
119C
120C     ----------------------------------------------------------------
121C
122C*       6.    DEFINE THERMODYNAMIC CONSTANTS, LIQUID PHASE.
123C              ---------------------------------------------
124C
125      RCW=4218.
126      WRITE(UNIT=6,FMT='('' *** Thermodynamic, liquid  ***'')')
127      WRITE(UNIT=6,FMT='(''         Cw   = '',E13.7)') RCW
128C
129C     ----------------------------------------------------------------
130C
131C*       7.    DEFINE THERMODYNAMIC CONSTANTS, SOLID PHASE.
132C              --------------------------------------------
133C
134      RCS=2106.
135      WRITE(UNIT=6,FMT='('' *** thermodynamic, solid   ***'')')
136      WRITE(UNIT=6,FMT='(''         Cs   = '',E13.7)') RCS
137C
138C     ----------------------------------------------------------------
139C
140C*       8.    DEFINE THERMODYNAMIC CONSTANTS, TRANSITION OF PHASE.
141C              ----------------------------------------------------
142C
143      RTT=273.16
144      RLVTT=2.5008E+6
145      RLSTT=2.8345E+6
146      RLMLT=RLSTT-RLVTT
147      RATM=100000.
148      WRITE(UNIT=6,FMT='('' *** Thermodynamic, trans.  ***'')')
149      WRITE(UNIT=6,FMT='('' Fusion point  = '',E13.7)') RTT
150      WRITE(UNIT=6,FMT='(''        RLvTt  = '',E13.7)') RLVTT
151      WRITE(UNIT=6,FMT='(''        RLsTt  = '',E13.7)') RLSTT
152      WRITE(UNIT=6,FMT='(''        RLMlt  = '',E13.7)') RLMLT
153      WRITE(UNIT=6,FMT='('' Normal press. = '',E13.7)') RATM
154      WRITE(UNIT=6,FMT='('' Latent heat :  '')')
155C
156C     ----------------------------------------------------------------
157C
158C*       9.    SATURATED VAPOUR PRESSURE.
159C              --------------------------
160C
161      RESTT=611.14
162      RGAMW=(RCW-RCPV)/RV
163      RBETW=RLVTT/RV+RGAMW*RTT
164      RALPW=LOG(RESTT)+RBETW/RTT+RGAMW*LOG(RTT)
165      RGAMS=(RCS-RCPV)/RV
166      RBETS=RLSTT/RV+RGAMS*RTT
167      RALPS=LOG(RESTT)+RBETS/RTT+RGAMS*LOG(RTT)
168      RGAMD=RGAMS-RGAMW
169      RBETD=RBETS-RBETW
170      RALPD=RALPS-RALPW
171C
172C     ------------------------------------------------------------------
173c
174c calculer les constantes pour les fonctions thermodynamiques
175c
176      RVTMP2=RCPV/RCPD-1.
177      RHOH2O=RATM/100.
178      R2ES=RESTT*RD/RV
179      R3LES=17.269
180      R3IES=21.875
181      R4LES=35.86
182      R4IES=7.66
183      R5LES=R3LES*(RTT-R4LES)
184      R5IES=R3IES*(RTT-R4IES)
185C
186      RETURN
187      END
Note: See TracBrowser for help on using the repository browser.