[1990] | 1 | MODULE YOMCST_SISVAT |
---|
| 2 | |
---|
| 3 | IMPLICIT 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 | |
---|
| 59 | END MODULE YOMCST_SISVAT |
---|