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

Diff of /trunk/Sources/phylmd/suphec.f

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

revision 169 by guez, Mon Sep 14 17:13:16 2015 UTC revision 198 by guez, Tue May 31 16:17:35 2016 UTC
# Line 11  module suphec_m Line 11  module suphec_m
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
# Line 38  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, save:: RCPV, RCVD, RCVV    real, parameter:: RCPV = 4. * RV ! Cpv, gas phase
46    
47      real, save:: RCVD, RCVV
48    
49    real, parameter:: RCPD = 7. / 2 * RD    real, parameter:: RCPD = 7. / 2 * RD
50    ! specific heat capacity for dry air, in J K-1 kg-1    ! specific heat capacity for dry air, in J K-1 kg-1
# Line 48  module suphec_m Line 54  module suphec_m
54    real, save:: RETV    real, save:: RETV
55    
56    ! A1.5, 6 Thermodynamic liquid, solid phases    ! A1.5, 6 Thermodynamic liquid, solid phases
57    REAL, save:: RCW, RCS  
58      REAL, parameter:: RCW = RCPV ! LIQUID PHASE Cw
59      real, save:: RCS
60    
61    ! A1.7 Thermodynamic transition of phase    ! A1.7 Thermodynamic transition of phase
62    REAL, save:: RLMLT    REAL, save:: RLMLT
# Line 77  contains Line 85  contains
85    
86      ! 2. DEFINE ASTRONOMICAL CONSTANTS      ! 2. DEFINE ASTRONOMICAL CONSTANTS
87    
     RSIYEA = 365.25*RDAY*2.*PI/6.283076  
     RSIDAY = RDAY/(1.+RDAY/RSIYEA)  
88      ROMEGA = 2.*PI/RSIDAY      ROMEGA = 2.*PI/RSIDAY
89    
90      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  
91      print '('' omega = '', E13.7, '' s-1'')', ROMEGA      print '('' omega = '', E13.7, '' s-1'')', ROMEGA
92    
93      ! 3. DEFINE GEOIDE.      ! 3. DEFINE GEOIDE.
# Line 103  contains Line 104  contains
104      ! 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.      ! 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.
105    
106      RCVD = RCPD-RD      RCVD = RCPD-RD
     RCPV = 4. * RV  
107      RCVV = RCPV-RV      RCVV = RCPV-RV
108      RETV = RV / RD - 1.      RETV = RV / RD - 1.
109      print *, 'Thermodynamics, gas'      print *, 'Thermodynamics, gas'
# Line 112  contains Line 112  contains
112      print *, "rv = ", RV, "J K-1 kg-1"      print *, "rv = ", RV, "J K-1 kg-1"
113      print '('' Cpd = '', e13.7)', RCPD      print '('' Cpd = '', e13.7)', RCPD
114      print '('' Cvd = '', e13.7)', RCVD      print '('' Cvd = '', e13.7)', RCVD
     print '('' Cpv = '', e13.7)', RCPV  
115      print '('' Cvv = '', e13.7)', RCVV      print '('' Cvv = '', e13.7)', RCVV
116      print '('' Rd/Cpd = '', e13.7)', RKAPPA      print '('' Rd/Cpd = '', e13.7)', RKAPPA
117      print '('' Rv / Rd - 1 = '', e13.7)', RETV      print '('' Rv / Rd - 1 = '', e13.7)', RETV
118    
     ! 6. DEFINE THERMODYNAMIC CONSTANTS, LIQUID PHASE.  
   
     RCW = RCPV  
     print *, 'Thermodynamic, liquid '  
     print '('' Cw = '', E13.7)', RCW  
   
119      ! 7. DEFINE THERMODYNAMIC CONSTANTS, SOLID PHASE.      ! 7. DEFINE THERMODYNAMIC CONSTANTS, SOLID PHASE.
120    
121      RCS = RCPV      RCS = RCPV
Legend:
Removed from v.169  
changed lines
  Added in v.198
  ViewVC Help
Powered by ViewVC 1.1.21