source: LMDZ.3.3/branches/rel-LF/libf/phylmd/suphec.F @ 400

Last change on this file since 400 was 393, checked in by lmdzadmin, 22 years ago

Modifications de JLD sur la conservation de l'energie
On supprime les modifs de Pascale sur le cdrag, elles refroidissaient trop
l'atmosphere
LF
  • 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
RevLine 
[2]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
[393]125      RCW=RCPV
[2]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
[393]134      RCS=RCPV
[2]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.