/[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 49 by guez, Wed Aug 24 11:43:14 2011 UTC revision 53 by guez, Fri Oct 7 13:11:58 2011 UTC
# Line 24  module suphec_m Line 24  module suphec_m
24    
25    ! A1.4 Thermodynamic gas phase    ! A1.4 Thermodynamic gas phase
26    REAL, parameter:: R = RNAVO * RKBOL ! ideal gas constant, in J K-1 mol-1    REAL, parameter:: R = RNAVO * RKBOL ! ideal gas constant, in J K-1 mol-1
27    real RV, RCPD, RCPV, RCVD, RCVV    real RV, RCPV, RCVD, RCVV
28    real, parameter:: RMD = 28.9644 ! molar mass of dry air, in g mol-1    real, parameter:: RMD = 28.9644 ! molar mass of dry air, in g mol-1
29    
30    real, parameter:: RD = 1000. * R / RMD    real, parameter:: RD = 1000. * R / RMD
31    ! 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
32    ! (factor 1000: conversion from g to kg)    ! (factor 1000: conversion from g to kg)
33    
34      real, parameter:: RCPD = 7. / 2 * RD
35      ! specific heat capacity for dry air, in J K-1 kg-1
36    
37    real, parameter:: RMO3 = 47.9942    real, parameter:: RMO3 = 47.9942
38    real, parameter:: RMV = 18.0153    real, parameter:: RMV = 18.0153
39    REAL RKAPPA, RETV    REAL, parameter:: RKAPPA = RD/RCPD
40      real RETV
41    
42    ! A1.5, 6 Thermodynamic liquid, solid phases    ! A1.5, 6 Thermodynamic liquid, solid phases
43    REAL RCW, RCS    REAL RCW, RCS
# Line 56  contains Line 60  contains
60    
61    SUBROUTINE suphec    SUBROUTINE suphec
62    
63      ! From phylmd/suphec.F, v 1.2 2005/06/06 13:16:33      ! From phylmd/suphec.F, version 1.2 2005/06/06 13:16:33
64      ! Initialise certaines constantes et parametres physiques.      ! Initialise certaines constantes et certains paramètres physiques.
65    
66      !------------------------------------------      !------------------------------------------
67    
# Line 103  contains Line 107  contains
107      ! 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.      ! 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.
108    
109      RV = 1000.*R/RMV      RV = 1000.*R/RMV
     RCPD = 3.5*RD  
110      RCVD = RCPD-RD      RCVD = RCPD-RD
111      RCPV = 4. *RV      RCPV = 4. *RV
112      RCVV = RCPV-RV      RCVV = RCPV-RV
     RKAPPA = RD/RCPD  
113      RETV = RV/RD-1.      RETV = RV/RD-1.
114      print *, 'Thermodynamic, gas '      print *, 'Thermodynamic, gas '
115      print '('' Perfect gas = '', e13.7)', R      print '('' Perfect gas = '', e13.7)', R
Legend:
Removed from v.49  
changed lines
  Added in v.53
  ViewVC Help
Powered by ViewVC 1.1.21