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 |
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, RCPD, 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:: RMO3 = 47.9942 |
real, parameter:: RMO3 = 47.9942 |
35 |
real, parameter:: RMV = 18.0153 |
real, parameter:: RMV = 18.0153 |
36 |
REAL RKAPPA, RETV |
REAL RKAPPA, RETV |
102 |
|
|
103 |
! 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE. |
! 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE. |
104 |
|
|
|
R = RNAVO * RKBOL |
|
|
RD = 1000. * R / RMD |
|
105 |
RV = 1000.*R/RMV |
RV = 1000.*R/RMV |
106 |
RCPD = 3.5*RD |
RCPD = 3.5*RD |
107 |
RCVD = RCPD-RD |
RCVD = RCPD-RD |
111 |
RETV = RV/RD-1. |
RETV = RV/RD-1. |
112 |
print *, 'Thermodynamic, gas ' |
print *, 'Thermodynamic, gas ' |
113 |
print '('' Perfect gas = '', e13.7)', R |
print '('' Perfect gas = '', e13.7)', R |
|
print '('' Dry air mass = '', e13.7)', RMD |
|
114 |
print '('' Ozone mass = '', e13.7)', RMO3 |
print '('' Ozone mass = '', e13.7)', RMO3 |
115 |
print '('' Vapour mass = '', e13.7)', RMV |
print '('' Vapour mass = '', e13.7)', RMV |
116 |
print '('' Dry air constant = '', e13.7)', RD |
print '('' Dry air constant = '', e13.7)', RD |