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

Diff of /trunk/phylmd/suphec.f

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

revision 171 by guez, Tue Sep 29 19:48:59 2015 UTC revision 198 by guez, Tue May 31 16:17:35 2016 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, 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 52  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 100  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 109  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.171  
changed lines
  Added in v.198
  ViewVC Help
Powered by ViewVC 1.1.21