/[lmdze]/trunk/libf/phylmd/suphec.f90
ViewVC logotype

Diff of /trunk/libf/phylmd/suphec.f90

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

revision 40 by guez, Tue Feb 22 13:49:36 2011 UTC revision 47 by guez, Fri Jul 1 15:00:48 2011 UTC
# Line 6  module suphec_m Line 6  module suphec_m
6    REAL RPI    REAL RPI
7    real, parameter:: RCLUM = 299792458.    real, parameter:: RCLUM = 299792458.
8    real, parameter:: RHPLA = 6.6260755E-34    real, parameter:: RHPLA = 6.6260755E-34
9    real, parameter:: RKBOL = 1.380658E-23    real, parameter:: RKBOL = 1.380658E-23 ! Boltzmann constant, in J K-1
10    real, parameter:: RNAVO = 6.0221367E+23    real, parameter:: RNAVO = 6.0221367E+23 ! Avogadro number, in mol-1
11    
12    ! A1.1 Astronomical constants    ! A1.1 Astronomical constants
13    REAL RSIYEA, RSIDAY, ROMEGA    REAL RSIYEA, RSIDAY, ROMEGA
# Line 16  module suphec_m Line 16  module suphec_m
16    real, parameter:: REPSM = 0.409093    real, parameter:: REPSM = 0.409093
17    
18    ! A1.2 Geoide    ! A1.2 Geoide
19    real, parameter:: RG = 9.80665    real, parameter:: RG = 9.80665 ! acceleration of gravity, in m s-2
20    real, parameter:: RA = 6371229.    real, parameter:: RA = 6371229.
21    
22    ! A1.3 Radiation    ! A1.3 Radiation
23    REAL RSIGMA    REAL RSIGMA
24    
25    ! A1.4 Thermodynamic gas phase    ! A1.4 Thermodynamic gas phase
26    REAL R, RD, RV, RCPD, RCPV, RCVD, RCVV    REAL, parameter:: R = RNAVO * RKBOL ! ideal gas constant, in J K-1 mol-1
27    real, parameter:: RMD = 28.9644    real RV, RCPD, RCPV, RCVD, RCVV
28      real, parameter:: RMD = 28.9644 ! molar mass of dry air, in g mol-1
29    
30      real, parameter:: RD = 1000. * R / RMD
31      ! specific ideal gas constant for dry air, in J K-1 kg-1
32    
33    real, parameter:: RMO3 = 47.9942    real, parameter:: RMO3 = 47.9942
34    real, parameter:: RMV = 18.0153    real, parameter:: RMV = 18.0153
35    REAL RKAPPA, RETV    REAL RKAPPA, RETV
# Line 96  contains Line 101  contains
101    
102      ! 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.      ! 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.
103    
     R = RNAVO * RKBOL  
     RD = 1000. * R / RMD  
104      RV = 1000.*R/RMV      RV = 1000.*R/RMV
105      RCPD = 3.5*RD      RCPD = 3.5*RD
106      RCVD = RCPD-RD      RCVD = RCPD-RD
# Line 107  contains Line 110  contains
110      RETV = RV/RD-1.      RETV = RV/RD-1.
111      print *, 'Thermodynamic, gas '      print *, 'Thermodynamic, gas '
112      print '('' Perfect gas = '', e13.7)', R      print '('' Perfect gas = '', e13.7)', R
     print '('' Dry air mass = '', e13.7)', RMD  
113      print '('' Ozone mass = '', e13.7)', RMO3      print '('' Ozone mass = '', e13.7)', RMO3
114      print '('' Vapour mass = '', e13.7)', RMV      print '('' Vapour mass = '', e13.7)', RMV
115      print '('' Dry air constant = '', e13.7)', RD      print '('' Dry air constant = '', e13.7)', RD
Legend:
Removed from v.40  
changed lines
  Added in v.47
  ViewVC Help
Powered by ViewVC 1.1.21