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 |
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 |
|
|
53 |
real, parameter:: RTT = 273.16 |
real, parameter:: RTT = 273.16 |
54 |
real, parameter:: RLVTT = 2.5008E+6 |
real, parameter:: RLVTT = 2.5008E+6 |
55 |
real, parameter:: RLSTT = 2.8345E+6 |
real, parameter:: RLSTT = 2.8345E+6 |
56 |
real, parameter:: RATM = 100000. |
real, parameter:: RATM = 1e5 |
57 |
|
|
58 |
! A1.8 Curve of saturation |
! A1.8 Curve of saturation |
59 |
REAL RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS |
REAL RALPW, RBETW, RGAMW, RALPS, RBETS, RGAMS |
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 |
|
|
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 |
|
|
140 |
! 8. DEFINE THERMODYNAMIC CONSTANTS, TRANSITION OF PHASE. |
! 8. DEFINE THERMODYNAMIC CONSTANTS, TRANSITION OF PHASE. |
141 |
|
|
142 |
RLMLT = RLSTT-RLVTT |
RLMLT = RLSTT-RLVTT |
143 |
print *, 'Thermodynamic, trans. ' |
print *, 'Thermodynamic, transition of phase:' |
144 |
print '('' Fusion point = '', E13.7)', RTT |
print '('' Fusion point = '', E13.7)', RTT |
145 |
print '('' RLvTt = '', E13.7)', RLVTT |
print '('' RLvTt = '', E13.7)', RLVTT |
146 |
print '('' RLsTt = '', E13.7)', RLSTT |
print '('' RLsTt = '', E13.7)', RLSTT |
147 |
print '('' RLMlt = '', E13.7)', RLMLT |
print '('' RLMlt = '', E13.7)', RLMLT |
148 |
print '('' Normal press. = '', E13.7)', RATM |
print '('' Normal pressure = '', E13.7)', RATM |
149 |
|
|
150 |
! 9. SATURATED VAPOUR PRESSURE. |
! 9. SATURATED VAPOUR PRESSURE. |
151 |
|
|