/[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 53 by guez, Fri Oct 7 13:11:58 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, 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      ! (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 50  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 96  contains Line 106  contains
106    
107      ! 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.      ! 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.
108    
     R = RNAVO * RKBOL  
     RD = 1000. * R / RMD  
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
     print '('' Dry air mass = '', e13.7)', RMD  
116      print '('' Ozone mass = '', e13.7)', RMO3      print '('' Ozone mass = '', e13.7)', RMO3
117      print '('' Vapour mass = '', e13.7)', RMV      print '('' Vapour mass = '', e13.7)', RMV
118      print '('' Dry air constant = '', e13.7)', RD      print '('' Dry air constant = '', e13.7)', RD
Legend:
Removed from v.40  
changed lines
  Added in v.53
  ViewVC Help
Powered by ViewVC 1.1.21