/[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 53 by guez, Fri Oct 7 13:11:58 2011 UTC revision 70 by guez, Mon Jun 24 15:39:52 2013 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 ! Boltzmann constant, in J K-1    real, parameter:: KBOL = 1.380658E-23 ! Boltzmann constant, in J K-1
10    real, parameter:: RNAVO = 6.0221367E+23 ! Avogadro number, in mol-1    real, parameter:: NAVO = 6.0221367E23 ! 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 23  module suphec_m Line 23  module suphec_m
23    REAL RSIGMA    REAL RSIGMA
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 = NAVO * KBOL ! ideal gas constant, in J K-1 mol-1
27    real RV, RCPV, RCVD, RCVV    real, parameter:: MV = 18.0153 ! molar mass of water, in g mol-1
   real, parameter:: RMD = 28.9644 ! molar mass of dry air, in g mol-1  
28    
29    real, parameter:: RD = 1000. * R / RMD    real, parameter:: RV = 1e3 * R / MV
30      ! specific ideal gas constant for water vapor, in J K-1 kg-1
31      ! (factor 1e3: conversion from g to kg)
32    
33      real, parameter:: MD = 28.9644 ! molar mass of dry air, in g mol-1
34    
35      real, parameter:: RD = 1e3 * R / MD
36    ! 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
37    ! (factor 1000: conversion from g to kg)    ! (factor 1e3: conversion from g to kg)
38    
39      real RCPV, RCVD, RCVV
40    
41    real, parameter:: RCPD = 7. / 2 * RD    real, parameter:: RCPD = 7. / 2 * RD
42    ! specific heat capacity for dry air, in J K-1 kg-1    ! specific heat capacity for dry air, in J K-1 kg-1
43    
44    real, parameter:: RMO3 = 47.9942    real, parameter:: RMO3 = 47.9942
   real, parameter:: RMV = 18.0153  
45    REAL, parameter:: RKAPPA = RD/RCPD    REAL, parameter:: RKAPPA = RD/RCPD
46    real RETV    real RETV
47    
# Line 75  contains Line 81  contains
81      print '('' PI = '', E13.7, '' -'')', RPI      print '('' PI = '', E13.7, '' -'')', RPI
82      print '('' c = '', E13.7, ''m s-1'')', RCLUM      print '('' c = '', E13.7, ''m s-1'')', RCLUM
83      print '('' h = '', E13.7, ''J s'')', RHPLA      print '('' h = '', E13.7, ''J s'')', RHPLA
     print '('' K = '', E13.7, ''J K-1'')', RKBOL  
     print '('' N = '', E13.7, ''mol-1'')', RNAVO  
84    
85      ! 2. DEFINE ASTRONOMICAL CONSTANTS      ! 2. DEFINE ASTRONOMICAL CONSTANTS
86    
# Line 100  contains Line 104  contains
104    
105      ! 4. DEFINE RADIATION CONSTANTS.      ! 4. DEFINE RADIATION CONSTANTS.
106    
107      rsigma = 2.*rpi**5 * (rkbol/rhpla)**3 * rkbol/rclum/rclum/15.      rsigma = 2.*rpi**5 * (kbol/rhpla)**3 * kbol/rclum/rclum/15.
108      print *, ' Radiation '      print *, ' Radiation '
109      print '('' Stefan-Bol. = '', E13.7, '' W m-2 K-4'')', RSIGMA      print '('' Stefan-Bol. = '', E13.7, '' W m-2 K-4'')', RSIGMA
110    
111      ! 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.      ! 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.
112    
     RV = 1000.*R/RMV  
113      RCVD = RCPD-RD      RCVD = RCPD-RD
114      RCPV = 4. *RV      RCPV = 4. * RV
115      RCVV = RCPV-RV      RCVV = RCPV-RV
116      RETV = RV/RD-1.      RETV = RV / RD - 1.
117      print *, 'Thermodynamic, gas '      print *, 'Thermodynamics, gas'
     print '('' Perfect gas = '', e13.7)', R  
118      print '('' Ozone mass = '', e13.7)', RMO3      print '('' Ozone mass = '', e13.7)', RMO3
119      print '('' Vapour mass = '', e13.7)', RMV      print *, "rd = ", RD, "J K-1 kg-1"
120      print '('' Dry air constant = '', e13.7)', RD      print *, "rv = ", RV, "J K-1 kg-1"
     print '('' Vapour constant = '', e13.7)', RV  
121      print '('' Cpd = '', e13.7)', RCPD      print '('' Cpd = '', e13.7)', RCPD
122      print '('' Cvd = '', e13.7)', RCVD      print '('' Cvd = '', e13.7)', RCVD
123      print '('' Cpv = '', e13.7)', RCPV      print '('' Cpv = '', e13.7)', RCPV
124      print '('' Cvv = '', e13.7)', RCVV      print '('' Cvv = '', e13.7)', RCVV
125      print '('' Rd/Cpd = '', e13.7)', RKAPPA      print '('' Rd/Cpd = '', e13.7)', RKAPPA
126      print '('' Rv/Rd-1 = '', e13.7)', RETV      print '('' Rv / Rd - 1 = '', e13.7)', RETV
127    
128      ! 6. DEFINE THERMODYNAMIC CONSTANTS, LIQUID PHASE.      ! 6. DEFINE THERMODYNAMIC CONSTANTS, LIQUID PHASE.
129    
Legend:
Removed from v.53  
changed lines
  Added in v.70
  ViewVC Help
Powered by ViewVC 1.1.21