source: LMDZ5/branches/IPSLCM5A2.1_ISO/libf/phyiso/sisvat/YOMCST_SISVAT.F90

Last change on this file was 3331, checked in by acozic, 7 years ago

Add modification for isotopes

  • Property svn:executable set to *
File size: 2.3 KB
Line 
1MODULE YOMCST_SISVAT
2
3IMPLICIT NONE
4                !
5! $Header$
6!
7!  ATTENTION!!!!: ce fichier include est compatible format fixe/format libre
8!                 veillez à n'utiliser que des ! pour les commentaires
9!                 et à bien positionner les & des lignes de continuation
10!                 (les placer en colonne 6 et en colonne 73)
11!
12!
13! A1.0 Fundamental constants
14      REAL,PARAMETER :: RPI=3.141592653589793238462643e0,               &
15     &                  RCLUM=299792458.,   RHPLA=6.6260755E-34,        &
16     &                  RKBOL=1.380658E-23, RNAVO=6.0221367E+23
17
18! A1.1 Astronomical constants
19      REAL,PARAMETER :: RDAY=86400.,REA=149597870000.,REPSM=0.409093,   &
20     &        RSIYEA=365.25*RDAY*2.*RPI/6.283076,                       &
21     &        RSIDAY=RDAY/(1.+RDAY/RSIYEA), ROMEGA=2.*RPI/RSIDAY
22
23! A1.1.bis Constantes concernant l'orbite de la Terre:
24      REAL,PARAMETER :: R_ecc=0.016715,R_peri=102.7,R_incl=23.441
25
26! A1.2 Geoide
27      REAL,PARAMETER :: RG=9.80665, RA=6371229.,R1SA=1.D0/RA
28
29! A1.3 Radiation
30      REAL,PARAMETER :: RSIGMA = 2.*RPI**5 * (RKBOL/RHPLA)**3           &
31     &                   * RKBOL/RCLUM/RCLUM/15.
32!     REAL RSIGMA,RI0
33
34! A1.4 Thermodynamic gas phase
35      REAL,PARAMETER :: R=RNAVO*RKBOL, RMD=28.9644,    RMO3=47.9942,    &
36     &                  RMV=18.0153,   RD=1000.*R/RMD, RV=1000.*R/RMV,  &
37     &                  RCPD=3.5*RD,   RCVD=RCPD-RD,   RCPV=4. *RV,     &
38     &                  RCVV=RCPV-RV,  RKAPPA=RD/RCPD, RETV=RV/RD-1.
39! A1.5,6 Thermodynamic liquid,solid phases
40      REAL,PARAMETER :: RCW=RCPV,      RCS=RCPV
41
42! A1.7 Thermodynamic transition of phase
43      REAL,PARAMETER :: RTT=273.16,    RLVTT=2.5008E+6,RLSTT=2.8345E+6, &
44     &                  RLMLT=RLSTT-RLVTT,             RATM=100000.
45
46! A1.8 Curve of saturation
47      REAL,PARAMETER :: RESTT=611.14,  RGAMW=(RCW-RCPV)/RV,             &
48     &                  RBETW=RLVTT/RV+RGAMW*RTT,                       &
49     &                  RGAMS=(RCS-RCPV)/RV,                            &   
50     &                  RBETS=RLSTT/RV+RGAMS*RTT,                       & 
51     &                  RGAMD=RGAMS-RGAMW, RBETD=RBETS-RBETW
52      REAL           :: RALPW,RALPS,RALPD
53
54!      RALPW = LOG(RESTT)+RBETW/RTT+RGAMW*LOG(RTT)
55!      RALPS = LOG(RESTT)+RBETS/RTT+RGAMS*LOG(RTT)
56!      RALPD = RALPS-RALPW
57   
58
59END MODULE YOMCST_SISVAT
Note: See TracBrowser for help on using the repository browser.