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 |
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 |
|
|
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 |