1 |
C $Header: /u/gcmpack/MITgcm/pkg/bulk_force/BULKF_PARAMS.h,v 1.3 2006/03/14 15:57:31 jmc Exp $ |
2 |
C $Name: $ |
3 |
|
4 |
#ifdef ALLOW_BULK_FORCE |
5 |
C *==========================================================* |
6 |
C | BULK_PARAMS.h |
7 |
C | o Header file for BULK_FORCE package parameters: |
8 |
C | - basic parameter ( I/O frequency, etc ...) |
9 |
C | - physical constants |
10 |
C *==========================================================* |
11 |
|
12 |
C-- COMMON / BULK_PHYS_CONST / physical constants |
13 |
C.. densities |
14 |
C rhoA :: density of air [kg/m^3] |
15 |
C rhoFW :: density of fresh water [kg/m^3] |
16 |
C rhoSW :: density of sea water [kg/m^3] |
17 |
C.. specific heats |
18 |
C cpAir :: specific heat of air [J/kg/K] |
19 |
C cpwv :: specific heat of water vapour [J/kg/K] |
20 |
C.. latent heat |
21 |
C Lvap :: latent heat of vaporization at 0.oC [J/kg] |
22 |
C Lfresh :: latent heat of melting of pure ice [J/kg] |
23 |
C Lvap_ice :: latent heat from sublimation [J/kg] |
24 |
C.. constants |
25 |
C Tf0kel :: Freezing temp of fresh water in Kelvin = 273.15 |
26 |
C Rgas :: gas constant for dry air [J/kg/K] |
27 |
C Rvap :: gas constant for H2O vapor [J/kg/K] |
28 |
C xkar :: von Karman constant [-] |
29 |
C stefan :: Stefan-Boltzmann constant [W/m^2/K^4] |
30 |
C.. for transfer coefficient |
31 |
C zref :: reference height [m] for transfer coefficient |
32 |
C zwd :: height [m] of near-surface wind-speed input data |
33 |
C zth :: height [m] of near-surface air-temp. & air-humid. input |
34 |
C cDrag_[n] :: n = 1,2,3 coefficients used to evaluate drag coefficient |
35 |
C cStantonS,U :: coefficients used to evaluate Stanton number (for |
36 |
C sensib. Heat Flx), under Stable / Unstable stratification |
37 |
C cDalton :: coefficient used to evaluate Dalton number (for Evap) |
38 |
C.. for bulk formula |
39 |
C umin :: minimum wind speed used in bulk-formulae [m/s] |
40 |
C humid_fac :: dry-air - water-vapor molecular mass ratio (minus one) |
41 |
C (used to calculate virtual temp.) |
42 |
C saltQsFac :: reduction of sat. vapor pressure over salty water |
43 |
C gamma_blk :: adiabatic lapse rate |
44 |
C.. for Long-Wave Radiation |
45 |
C atm_emissivity :: |
46 |
C ocean_emissivity :: |
47 |
C snow_emissivity :: |
48 |
C ice_emissivity :: |
49 |
C.. for BULKF_FORMULA_AIM |
50 |
C FWIND0 :: ratio of near-sfc wind to lowest-level wind |
51 |
C CHS :: heat exchange coefficient over sea |
52 |
C VGUST :: wind speed for sub-grid-scale gusts |
53 |
C DTHETA :: Potential temp. gradient for stability correction |
54 |
C dTstab :: potential temp. increment for stability function derivative |
55 |
C FSTAB :: Amplitude of stability correction (fraction) |
56 |
C.. Albedo |
57 |
C ocean_albedo :: ocean surface albedo [0-1] |
58 |
|
59 |
_RL rhoA |
60 |
c _RL rhoSW |
61 |
_RL rhoFW |
62 |
_RL cpAir |
63 |
c _RL cpwv |
64 |
_RL Lvap |
65 |
_RL Lfresh |
66 |
c _RL Lvap_ice |
67 |
_RL Tf0kel |
68 |
_RL Rgas |
69 |
c _RL Rvap |
70 |
_RL xkar |
71 |
_RL stefan |
72 |
_RL zref, zwd, zth |
73 |
_RL cDrag_1, cDrag_2, cDrag_3 |
74 |
_RL cStantonS, cStantonU |
75 |
_RL cDalton |
76 |
_RL umin |
77 |
_RL humid_fac |
78 |
_RL saltQsFac |
79 |
_RL gamma_blk |
80 |
_RL atm_emissivity |
81 |
_RL ocean_emissivity |
82 |
_RL snow_emissivity |
83 |
_RL ice_emissivity |
84 |
#ifdef ALLOW_FORMULA_AIM |
85 |
_RL FWIND0, CHS, VGUST, DTHETA, dTstab, FSTAB |
86 |
#endif |
87 |
_RL ocean_albedo |
88 |
COMMON / BULK_PHYS_CONST / |
89 |
& rhoA, rhoFW, |
90 |
& cpAir, Lvap, Lfresh, |
91 |
& Tf0kel, Rgas, |
92 |
& xkar, stefan, |
93 |
& zref, zwd, zth, |
94 |
& cDrag_1, cDrag_2, cDrag_3, |
95 |
& cStantonS, cStantonU, cDalton, |
96 |
& umin, humid_fac, saltQsFac, gamma_blk, |
97 |
& atm_emissivity, ocean_emissivity, |
98 |
& snow_emissivity, ice_emissivity, |
99 |
#ifdef ALLOW_FORMULA_AIM |
100 |
& FWIND0, CHS, VGUST, DTHETA, dTstab, FSTAB, |
101 |
#endif |
102 |
& ocean_albedo |
103 |
|
104 |
C-- COMMON / BULK_PAR_I / Integer parameters |
105 |
C blk_nIter :: Number of iterations to find turbulent transfer coeff. |
106 |
|
107 |
INTEGER blk_nIter |
108 |
COMMON / BULK_PAR_I / |
109 |
& blk_nIter |
110 |
|
111 |
C-- COMMON / BULK_PAR_R / real parameter |
112 |
C blk_taveFreq :: time-average output frequency [s] |
113 |
|
114 |
_RL blk_taveFreq |
115 |
COMMON / BULK_PAR_R / |
116 |
& blk_taveFreq |
117 |
|
118 |
C-- COMMON / BULK_PAR_L / Logical parameters |
119 |
C .. for BULKF_FORMULA_AIM |
120 |
C calcWindStress :: True to calculate Wind-Stress from surface wind |
121 |
C useFluxFormula_AIM :: set to T when using AIM flux formula rather |
122 |
C than the default formula (LANL) |
123 |
|
124 |
LOGICAL calcWindStress |
125 |
LOGICAL useFluxFormula_AIM |
126 |
LOGICAL useQnetch, useEmPch |
127 |
COMMON / BULK_PAR_L / |
128 |
& calcWindStress, |
129 |
& useFluxFormula_AIM, |
130 |
& useQnetch, useEmPch |
131 |
|
132 |
C-- COMMON / BULK_PAR_C / Character string parameters |
133 |
|
134 |
CHARACTER*(MAX_LEN_FNAM) AirTempFile |
135 |
CHARACTER*(MAX_LEN_FNAM) AirHumidityFile |
136 |
CHARACTER*(MAX_LEN_FNAM) RainFile |
137 |
CHARACTER*(MAX_LEN_FNAM) SolarFile |
138 |
CHARACTER*(MAX_LEN_FNAM) LongwaveFile |
139 |
CHARACTER*(MAX_LEN_FNAM) UWindFile |
140 |
CHARACTER*(MAX_LEN_FNAM) VWindFile |
141 |
CHARACTER*(MAX_LEN_FNAM) WSpeedFile |
142 |
CHARACTER*(MAX_LEN_FNAM) RunoffFile |
143 |
CHARACTER*(MAX_LEN_FNAM) QnetFile |
144 |
CHARACTER*(MAX_LEN_FNAM) EmPFile |
145 |
CHARACTER*(MAX_LEN_FNAM) CloudFile |
146 |
CHARACTER*(MAX_LEN_FNAM) SnowFile |
147 |
CHARACTER*(MAX_LEN_FNAM) airPotTempFile |
148 |
COMMON / BULK_PAR_C / |
149 |
& AirTempFile, AirHumidityFile, RainFile, |
150 |
& SolarFile, LongwaveFile, UWindFile, VWindFile, |
151 |
& WSpeedFile, RunoffFile, |
152 |
& QnetFile, EmPFile, CloudFile, SnowFile, |
153 |
& airPotTempFile |
154 |
|
155 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
156 |
|
157 |
#endif /*ALLOW_BULK_FORCE*/ |