/[lmdze]/trunk/phylmd/suphec.f
ViewVC logotype

Diff of /trunk/phylmd/suphec.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/libf/phylmd/suphec.f90 revision 70 by guez, Mon Jun 24 15:39:52 2013 UTC trunk/phylmd/suphec.f revision 254 by guez, Mon Feb 5 10:39:38 2018 UTC
# Line 1  Line 1 
1  module suphec_m  module suphec_m
2    
3      use nr_util, only: pi
4    
5    implicit none    implicit none
6    
7    ! A1.0 Fundamental constants    ! A1.0 Fundamental constants
8    REAL RPI    real, parameter:: RCLUM = 299792458. ! speed of light, m s-1
9    real, parameter:: RCLUM = 299792458.    real, parameter:: RHPLA = 6.6260755E-34 ! Planck constant, J s
   real, parameter:: RHPLA = 6.6260755E-34  
10    real, parameter:: KBOL = 1.380658E-23 ! Boltzmann constant, in J K-1    real, parameter:: KBOL = 1.380658E-23 ! Boltzmann constant, in J K-1
11    real, parameter:: NAVO = 6.0221367E23 ! Avogadro number, in mol-1    real, parameter:: NAVO = 6.0221367E23 ! Avogadro number, in mol-1
12    
13    ! A1.1 Astronomical constants    ! A1.1 Astronomical constants
14    REAL RSIYEA, RSIDAY, ROMEGA  
15      REAL ROMEGA
16    real, parameter:: RDAY = 86400.    real, parameter:: RDAY = 86400.
17    real, parameter:: REA = 149597870000.  
18    real, parameter:: REPSM = 0.409093    REAL, parameter:: RSIYEA = 365.25 * RDAY * 2. * PI / 6.283076
19      ! sideral year, in s
20    
21      REAL, parameter:: RSIDAY = RDAY / (1. + RDAY / RSIYEA) ! sideral day, in s
22    
23    ! A1.2 Geoide    ! A1.2 Geoide
24    real, parameter:: RG = 9.80665 ! acceleration of gravity, in m s-2    real, parameter:: RG = 9.80665 ! acceleration of gravity, in m s-2
25    real, parameter:: RA = 6371229.    real, parameter:: RA = 6371229.
26    
27    ! A1.3 Radiation    ! A1.3 Radiation
28    REAL RSIGMA    REAL, parameter:: rsigma = 2. * pi**5 * (kbol / rhpla)**3 * kbol / rclum**2 &
29           / 15.
30    
31    ! A1.4 Thermodynamic gas phase    ! A1.4 Thermodynamic gas phase
32    REAL, parameter:: R = NAVO * KBOL ! ideal gas constant, in J K-1 mol-1    REAL, parameter:: R = NAVO * KBOL ! ideal gas constant, in J K-1 mol-1
# Line 36  module suphec_m Line 42  module suphec_m
42    ! specific ideal gas constant for dry air, in J K-1 kg-1    ! specific ideal gas constant for dry air, in J K-1 kg-1
43    ! (factor 1e3: conversion from g to kg)    ! (factor 1e3: conversion from g to kg)
44    
45    real RCPV, RCVD, RCVV    real, parameter:: RCPV = 4. * RV
46      ! specific heat capacity at constant pressure of water vapor, in J K-1 kg-1
47    
48      real, save:: RCVD
49      real, save:: RCVV
50    
51    real, parameter:: RCPD = 7. / 2 * RD    real, parameter:: RCPD = 7. / 2 * RD
52    ! specific heat capacity for dry air, in J K-1 kg-1    ! specific heat capacity at constant pressure of dry air, in J K-1 kg-1
53    
54    real, parameter:: RMO3 = 47.9942    real, parameter:: RMO3 = 47.9942
55    REAL, parameter:: RKAPPA = RD/RCPD    REAL, parameter:: RKAPPA = RD/RCPD
56    real RETV    real, save:: RETV
57    
58    ! A1.5, 6 Thermodynamic liquid, solid phases    ! A1.5, 6 Thermodynamic liquid, solid phases
59    REAL RCW, RCS  
60      REAL, parameter:: RCW = RCPV ! LIQUID PHASE Cw
61      real, save:: RCS
62    
63    ! A1.7 Thermodynamic transition of phase    ! A1.7 Thermodynamic transition of phase
64    REAL RLMLT    REAL, save:: RLMLT
65    real, parameter:: RTT = 273.16    real, parameter:: RTT = 273.16
66    real, parameter:: RLVTT = 2.5008E+6  
67      real, parameter:: RLVTT = 2.5008E+6
68      ! specific latent heat of vaporization of water at triple point, in J kg-1
69    
70    real, parameter:: RLSTT = 2.8345E+6    real, parameter:: RLSTT = 2.8345E+6
71    real, parameter:: RATM = 100000.    real, parameter:: RATM = 1e5
72    
73    ! A1.8 Curve of saturation    ! A1.8 Curve of saturation
74    REAL RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS    REAL, save:: RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS
75    real, parameter:: RESTT = 611.14    real, parameter:: RESTT = 611.14
76    REAL RALPD, RBETD, RGAMD    REAL, save:: RALPD, RBETD, RGAMD
77    
78    save    private pi
79    
80  contains  contains
81    
# Line 73  contains Line 88  contains
88    
89      PRINT *, 'Call sequence information: suphec'      PRINT *, 'Call sequence information: suphec'
90    
     ! 1. DEFINE FUNDAMENTAL CONSTANTS  
   
     print *, 'Constants of the ICM'  
     RPI = 2.*ASIN(1.)  
     print *, 'Fundamental constants '  
     print '('' PI = '', E13.7, '' -'')', RPI  
     print '('' c = '', E13.7, ''m s-1'')', RCLUM  
     print '('' h = '', E13.7, ''J s'')', RHPLA  
   
91      ! 2. DEFINE ASTRONOMICAL CONSTANTS      ! 2. DEFINE ASTRONOMICAL CONSTANTS
92    
93      RSIYEA = 365.25*RDAY*2.*RPI/6.283076      ROMEGA = 2.*PI/RSIDAY
     RSIDAY = RDAY/(1.+RDAY/RSIYEA)  
     ROMEGA = 2.*RPI/RSIDAY  
94    
95      print *, 'Astronomical constants '      print *, 'Astronomical constants '
     print '('' day = '', E13.7, '' s'')', RDAY  
     print '('' half g. axis = '', E13.7, '' m'')', REA  
     print '('' mean anomaly = '', E13.7, '' -'')', REPSM  
     print '('' sideral year = '', E13.7, '' s'')', RSIYEA  
     print '('' sideral day = '', E13.7, '' s'')', RSIDAY  
96      print '('' omega = '', E13.7, '' s-1'')', ROMEGA      print '('' omega = '', E13.7, '' s-1'')', ROMEGA
97    
98      ! 3. DEFINE GEOIDE.      ! 3. DEFINE GEOIDE.
# Line 104  contains Line 103  contains
103    
104      ! 4. DEFINE RADIATION CONSTANTS.      ! 4. DEFINE RADIATION CONSTANTS.
105    
     rsigma = 2.*rpi**5 * (kbol/rhpla)**3 * kbol/rclum/rclum/15.  
106      print *, ' Radiation '      print *, ' Radiation '
107      print '('' Stefan-Bol. = '', E13.7, '' W m-2 K-4'')', RSIGMA      print '('' Stefan-Bol. = '', E13.7, '' W m-2 K-4'')', RSIGMA
108    
109      ! 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.      ! 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.
110    
111      RCVD = RCPD-RD      RCVD = RCPD-RD
     RCPV = 4. * RV  
112      RCVV = RCPV-RV      RCVV = RCPV-RV
113      RETV = RV / RD - 1.      RETV = RV / RD - 1.
114      print *, 'Thermodynamics, gas'      print *, 'Thermodynamics, gas'
     print '('' Ozone mass = '', e13.7)', RMO3  
115      print *, "rd = ", RD, "J K-1 kg-1"      print *, "rd = ", RD, "J K-1 kg-1"
116      print *, "rv = ", RV, "J K-1 kg-1"      print *, "rv = ", RV, "J K-1 kg-1"
117      print '('' Cpd = '', e13.7)', RCPD      print '('' Cpd = '', e13.7)', RCPD
118      print '('' Cvd = '', e13.7)', RCVD      print '('' Cvd = '', e13.7)', RCVD
     print '('' Cpv = '', e13.7)', RCPV  
119      print '('' Cvv = '', e13.7)', RCVV      print '('' Cvv = '', e13.7)', RCVV
120      print '('' Rd/Cpd = '', e13.7)', RKAPPA      print '('' Rd/Cpd = '', e13.7)', RKAPPA
121      print '('' Rv / Rd - 1 = '', e13.7)', RETV      print '('' Rv / Rd - 1 = '', e13.7)', RETV
122    
     ! 6. DEFINE THERMODYNAMIC CONSTANTS, LIQUID PHASE.  
   
     RCW = RCPV  
     print *, 'Thermodynamic, liquid '  
     print '('' Cw = '', E13.7)', RCW  
   
123      ! 7. DEFINE THERMODYNAMIC CONSTANTS, SOLID PHASE.      ! 7. DEFINE THERMODYNAMIC CONSTANTS, SOLID PHASE.
124    
125      RCS = RCPV      RCS = RCPV
# Line 140  contains Line 129  contains
129      ! 8. DEFINE THERMODYNAMIC CONSTANTS, TRANSITION OF PHASE.      ! 8. DEFINE THERMODYNAMIC CONSTANTS, TRANSITION OF PHASE.
130    
131      RLMLT = RLSTT-RLVTT      RLMLT = RLSTT-RLVTT
132      print *, 'Thermodynamic, trans. '      print *, 'Thermodynamic, transition of phase:'
133      print '('' Fusion point = '', E13.7)', RTT      print '('' Fusion point = '', E13.7)', RTT
     print '('' RLvTt = '', E13.7)', RLVTT  
134      print '('' RLsTt = '', E13.7)', RLSTT      print '('' RLsTt = '', E13.7)', RLSTT
135      print '('' RLMlt = '', E13.7)', RLMLT      print '('' RLMlt = '', E13.7)', RLMLT
136      print '('' Normal press. = '', E13.7)', RATM      print '('' Normal pressure = '', E13.7)', RATM
137    
138      ! 9. SATURATED VAPOUR PRESSURE.      ! 9. SATURATED VAPOUR PRESSURE.
139    
Legend:
Removed from v.70  
changed lines
  Added in v.254
  ViewVC Help
Powered by ViewVC 1.1.21