1 |
C $Header: /u/gcmpack/MITgcm/pkg/bulk_force/BULKF_PARAMS.h,v 1.4 2006/05/25 17:30:54 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 |
_RL rhoA |
59 |
c _RL rhoSW |
60 |
_RL rhoFW |
61 |
_RL cpAir |
62 |
c _RL cpwv |
63 |
_RL Lvap |
64 |
_RL Lfresh |
65 |
c _RL Lvap_ice |
66 |
_RL Tf0kel |
67 |
_RL Rgas |
68 |
c _RL Rvap |
69 |
_RL xkar |
70 |
_RL stefan |
71 |
_RL zref, zwd, zth |
72 |
_RL cDrag_1, cDrag_2, cDrag_3 |
73 |
_RL cStantonS, cStantonU |
74 |
_RL cDalton |
75 |
_RL umin |
76 |
_RL humid_fac |
77 |
_RL saltQsFac |
78 |
_RL gamma_blk |
79 |
_RL atm_emissivity |
80 |
_RL ocean_emissivity |
81 |
_RL snow_emissivity |
82 |
_RL ice_emissivity |
83 |
#ifdef ALLOW_FORMULA_AIM |
84 |
_RL FWIND0, CHS, VGUST, DTHETA, dTstab, FSTAB |
85 |
#endif |
86 |
_RL ocean_albedo |
87 |
COMMON / BULK_PHYS_CONST / |
88 |
& rhoA, rhoFW, |
89 |
& cpAir, Lvap, Lfresh, |
90 |
& Tf0kel, Rgas, |
91 |
& xkar, stefan, |
92 |
& zref, zwd, zth, |
93 |
& cDrag_1, cDrag_2, cDrag_3, |
94 |
& cStantonS, cStantonU, cDalton, |
95 |
& umin, humid_fac, saltQsFac, gamma_blk, |
96 |
& atm_emissivity, ocean_emissivity, |
97 |
& snow_emissivity, ice_emissivity, |
98 |
#ifdef ALLOW_FORMULA_AIM |
99 |
& FWIND0, CHS, VGUST, DTHETA, dTstab, FSTAB, |
100 |
#endif |
101 |
& ocean_albedo |
102 |
|
103 |
C-- COMMON / BULK_PAR_I / Integer parameters |
104 |
C blk_nIter :: Number of iterations to find turbulent transfer coeff. |
105 |
|
106 |
INTEGER blk_nIter |
107 |
COMMON / BULK_PAR_I / |
108 |
& blk_nIter |
109 |
|
110 |
C-- COMMON / BULK_PAR_R / real parameter |
111 |
C blk_taveFreq :: time-average output frequency [s] |
112 |
|
113 |
_RL blk_taveFreq |
114 |
COMMON / BULK_PAR_R / |
115 |
& blk_taveFreq |
116 |
|
117 |
C-- COMMON / BULK_PAR_L / Logical parameters |
118 |
C .. for BULKF_FORMULA_AIM |
119 |
C calcWindStress :: True to calculate Wind-Stress from surface wind |
120 |
C useFluxFormula_AIM :: set to T when using AIM flux formula rather |
121 |
C than the default formula (LANL) |
122 |
|
123 |
LOGICAL calcWindStress |
124 |
LOGICAL useFluxFormula_AIM |
125 |
LOGICAL useQnetch, useEmPch |
126 |
COMMON / BULK_PAR_L / |
127 |
& calcWindStress, |
128 |
& useFluxFormula_AIM, |
129 |
& useQnetch, useEmPch |
130 |
|
131 |
C-- COMMON / BULK_PAR_C / Character string parameters |
132 |
|
133 |
CHARACTER*(MAX_LEN_FNAM) AirTempFile |
134 |
CHARACTER*(MAX_LEN_FNAM) AirHumidityFile |
135 |
CHARACTER*(MAX_LEN_FNAM) RainFile |
136 |
CHARACTER*(MAX_LEN_FNAM) SolarFile |
137 |
CHARACTER*(MAX_LEN_FNAM) LongwaveFile |
138 |
CHARACTER*(MAX_LEN_FNAM) UWindFile |
139 |
CHARACTER*(MAX_LEN_FNAM) VWindFile |
140 |
CHARACTER*(MAX_LEN_FNAM) WSpeedFile |
141 |
CHARACTER*(MAX_LEN_FNAM) RunoffFile |
142 |
CHARACTER*(MAX_LEN_FNAM) QnetFile |
143 |
CHARACTER*(MAX_LEN_FNAM) EmPFile |
144 |
CHARACTER*(MAX_LEN_FNAM) CloudFile |
145 |
CHARACTER*(MAX_LEN_FNAM) SnowFile |
146 |
CHARACTER*(MAX_LEN_FNAM) airPotTempFile |
147 |
COMMON / BULK_PAR_C / |
148 |
& AirTempFile, AirHumidityFile, RainFile, |
149 |
& SolarFile, LongwaveFile, UWindFile, VWindFile, |
150 |
& WSpeedFile, RunoffFile, |
151 |
& QnetFile, EmPFile, CloudFile, SnowFile, |
152 |
& airPotTempFile |
153 |
|
154 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
155 |
|
156 |
#endif /*ALLOW_BULK_FORCE*/ |