/[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/Sources/phylmd/suphec.f revision 198 by guez, Tue May 31 16:17:35 2016 UTC trunk/phylmd/suphec.f revision 266 by guez, Thu Apr 19 17:54:55 2018 UTC
# Line 42  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, parameter:: RCPV = 4. * RV ! Cpv, gas phase    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, RCVV    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, save:: RETV    real, parameter:: RETV = RV / RD - 1.
57    
58    ! A1.5, 6 Thermodynamic liquid, solid phases    ! A1.5, 6 Thermodynamic liquid, solid phases
59    
# Line 61  module suphec_m Line 63  module suphec_m
63    ! A1.7 Thermodynamic transition of phase    ! A1.7 Thermodynamic transition of phase
64    REAL, save:: 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 = 1e5    real, parameter:: RATM = 1e5
72    
# Line 105  contains Line 110  contains
110    
111      RCVD = RCPD-RD      RCVD = RCPD-RD
112      RCVV = RCPV-RV      RCVV = RCPV-RV
     RETV = RV / RD - 1.  
113      print *, 'Thermodynamics, gas'      print *, 'Thermodynamics, gas'
     print '('' Ozone mass = '', e13.7)', RMO3  
114      print *, "rd = ", RD, "J K-1 kg-1"      print *, "rd = ", RD, "J K-1 kg-1"
115      print *, "rv = ", RV, "J K-1 kg-1"      print *, "rv = ", RV, "J K-1 kg-1"
116      print '('' Cpd = '', e13.7)', RCPD      print '('' Cpd = '', e13.7)', RCPD
# Line 127  contains Line 130  contains
130      RLMLT = RLSTT-RLVTT      RLMLT = RLSTT-RLVTT
131      print *, 'Thermodynamic, transition of phase:'      print *, 'Thermodynamic, transition of phase:'
132      print '('' Fusion point = '', E13.7)', RTT      print '('' Fusion point = '', E13.7)', RTT
     print '('' RLvTt = '', E13.7)', RLVTT  
133      print '('' RLsTt = '', E13.7)', RLSTT      print '('' RLsTt = '', E13.7)', RLSTT
134      print '('' RLMlt = '', E13.7)', RLMLT      print '('' RLMlt = '', E13.7)', RLMLT
135      print '('' Normal pressure = '', E13.7)', RATM      print '('' Normal pressure = '', E13.7)', RATM
Legend:
Removed from v.198  
changed lines
  Added in v.266
  ViewVC Help
Powered by ViewVC 1.1.21