1 |
jmc |
1.2 |
C $Header: /u/gcmpack/MITgcm/pkg/bulk_force/BULKF_PARAMS.h,v 1.1 2003/11/23 01:36:55 jmc Exp $ |
2 |
jmc |
1.1 |
C $Name: $ |
3 |
|
|
|
4 |
|
|
#ifdef ALLOW_BULK_FORCE |
5 |
|
|
C *==========================================================* |
6 |
|
|
C | BULK_PARAMS.h |
7 |
|
|
C | o Header file for BULK_FORC package parameters: |
8 |
|
|
C | - basic parameter ( I/O frequency, etc ...) |
9 |
|
|
C | - physical constants (not used in therm_seaice pkg) |
10 |
|
|
C *==========================================================* |
11 |
|
|
|
12 |
jmc |
1.2 |
C-- COMMON / BULK_PHYS_CONST / physical constants |
13 |
jmc |
1.1 |
C.. densities |
14 |
|
|
C rhoa :: density of air (kg/m^3) |
15 |
|
|
C rhofw :: density of fresh water (kg/m^3) |
16 |
|
|
C.. specific heats |
17 |
|
|
C cpair :: specific heat of air (J/kg/K) |
18 |
|
|
C cpwv :: specific heat of water vapour (J/kg/KC) |
19 |
|
|
C .. latent heat |
20 |
|
|
C Lvap :: latent heat of vaporizn at 0 C (J/kg) |
21 |
|
|
C Lfresh :: latent heat of melting of pure ice (J/kg) |
22 |
|
|
C .. melting |
23 |
|
|
C Tf0kel :: Freezing temp of fresh ice in Kelvin = 273.15 |
24 |
|
|
C .. wind drag |
25 |
|
|
C drag_[n] :: n = 1,2,3 coefficients used to evaluate drag coefficient |
26 |
|
|
C .. constants |
27 |
|
|
C stefan :: Stefan-Boltzmann constant (W/m^2 K^4) |
28 |
|
|
C xkar :: Von Karman constant QQ units? |
29 |
|
|
C Rvap :: gas constant for H2O vapor (J/kg/K) |
30 |
|
|
C.. Miscellaneous |
31 |
|
|
C p0 :: surface pressure (mb) |
32 |
|
|
C.. Combinations used for efficiency |
33 |
|
|
C Qcoef :: another constant for latent heat flux |
34 |
|
|
C .. for bulk formula |
35 |
|
|
C humid_fac :: const. for the evaluation of the virtual temp. |
36 |
|
|
C saltsat :: reduction of sat. vapor pressure over salt water |
37 |
|
|
C gamma_blk :: adiabatic lapse rate |
38 |
|
|
C .. for Vince bulk formula QQ check units |
39 |
|
|
C Lvap_ice :: latent heat from sublimation |
40 |
|
|
C Rgas :: gas constant for dry air |
41 |
|
|
C Sha :: |
42 |
|
|
C .. emissivities |
43 |
|
|
C atm_emissivity :: |
44 |
|
|
C ocean_emissivity :: |
45 |
|
|
C snow_emissivity :: |
46 |
|
|
C ice_emissivity :: |
47 |
jmc |
1.2 |
C .. for BULKF_FORMULA_AIM |
48 |
|
|
C FWIND0 :: ratio of near-sfc wind to lowest-level wind |
49 |
|
|
C CHS :: heat exchange coefficient over sea |
50 |
|
|
C VGUST :: wind speed for sub-grid-scale gusts |
51 |
|
|
C DTHETA :: Potential temp. gradient for stability correction |
52 |
|
|
C dTstab :: potential temp. increment for stability function derivative |
53 |
|
|
C FSTAB :: Amplitude of stability correction (fraction) |
54 |
|
|
C .. Albedo |
55 |
|
|
C ocean_albedo :: ocean surface albedo [0-1] |
56 |
jmc |
1.1 |
|
57 |
|
|
_RL rhoa |
58 |
|
|
_RL rhosw |
59 |
|
|
_RL rhofw |
60 |
|
|
_RL cpair |
61 |
|
|
_RL cpwv |
62 |
|
|
_RL Lvap |
63 |
|
|
_RL Lfresh |
64 |
|
|
_RL Tf0kel |
65 |
jmc |
1.2 |
_RL cdrag_1, cdrag_2, cdrag_3 |
66 |
jmc |
1.1 |
_RL stefan |
67 |
|
|
_RL xkar |
68 |
|
|
_RL Rvap |
69 |
|
|
_RL p0 |
70 |
|
|
_RL Qcoef |
71 |
|
|
_RL humid_fac |
72 |
|
|
_RL saltsat |
73 |
|
|
_RL gamma_blk |
74 |
|
|
_RL Lvap_ice |
75 |
|
|
_RL Rgas |
76 |
|
|
_RL Sha |
77 |
|
|
_RL atm_emissivity |
78 |
|
|
_RL ocean_emissivity |
79 |
|
|
_RL snow_emissivity |
80 |
|
|
_RL ice_emissivity |
81 |
jmc |
1.2 |
#ifdef ALLOW_FORMULA_AIM |
82 |
|
|
_RL FWIND0, CHS, VGUST, DTHETA, dTstab, FSTAB |
83 |
|
|
#endif |
84 |
|
|
_RL ocean_albedo |
85 |
|
|
|
86 |
|
|
COMMON / BULK_PHYS_CONST / |
87 |
|
|
& rhoa, rhosw, rhofw, |
88 |
|
|
& cpair, cpwv, |
89 |
|
|
& Lvap, Lfresh, |
90 |
|
|
& Tf0kel, |
91 |
|
|
& cdrag_1, cdrag_2, cdrag_3, |
92 |
|
|
& stefan, xkar, Rvap, |
93 |
|
|
& p0, Qcoef, |
94 |
|
|
& humid_fac, saltsat, gamma_blk, |
95 |
|
|
& Lvap_ice, Rgas, Sha, |
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_R / real parameter |
104 |
|
|
|
105 |
jmc |
1.1 |
_RL blk_taveFreq |
106 |
|
|
|
107 |
jmc |
1.2 |
COMMON / BULK_PAR_R / |
108 |
|
|
& blk_taveFreq |
109 |
|
|
|
110 |
|
|
C-- COMMON / BULK_PAR_L / Logical parameters |
111 |
|
|
C .. for BULKF_FORMULA_AIM |
112 |
|
|
C calcWindStress :: True to calculate Wind-Stress from surface wind |
113 |
|
|
C useFluxFormula_AIM :: set to T when using AIM flux formula rather |
114 |
|
|
C than the default formula (LANL) |
115 |
|
|
|
116 |
|
|
LOGICAL calcWindStress |
117 |
|
|
LOGICAL useFluxFormula_AIM |
118 |
|
|
|
119 |
|
|
COMMON / BULK_PAR_L / |
120 |
|
|
& calcWindStress, |
121 |
|
|
& useFluxFormula_AIM |
122 |
|
|
|
123 |
|
|
C-- COMMON / BULK_PAR_C / Character string parameters |
124 |
|
|
CHARACTER*(MAX_LEN_FNAM) AirTempFile |
125 |
|
|
CHARACTER*(MAX_LEN_FNAM) AirHumidityFile |
126 |
|
|
CHARACTER*(MAX_LEN_FNAM) RainFile |
127 |
|
|
CHARACTER*(MAX_LEN_FNAM) SolarFile |
128 |
|
|
CHARACTER*(MAX_LEN_FNAM) LongwaveFile |
129 |
|
|
CHARACTER*(MAX_LEN_FNAM) UWindFile |
130 |
|
|
CHARACTER*(MAX_LEN_FNAM) VWindFile |
131 |
|
|
CHARACTER*(MAX_LEN_FNAM) WSpeedFile |
132 |
|
|
CHARACTER*(MAX_LEN_FNAM) RunoffFile |
133 |
|
|
CHARACTER*(MAX_LEN_FNAM) QnetFile |
134 |
|
|
CHARACTER*(MAX_LEN_FNAM) EmPFile |
135 |
|
|
CHARACTER*(MAX_LEN_FNAM) CloudFile |
136 |
|
|
CHARACTER*(MAX_LEN_FNAM) SnowFile |
137 |
|
|
CHARACTER*(MAX_LEN_FNAM) airPotTempFile |
138 |
|
|
|
139 |
|
|
COMMON / BULK_PAR_C / |
140 |
|
|
& AirTempFile, AirHumidityFile, RainFile, |
141 |
|
|
& SolarFile, LongwaveFile, UWindFile, VWindFile, |
142 |
|
|
& WSpeedFile, RunoffFile, |
143 |
|
|
& QnetFile, EmPFile, CloudFile, SnowFile, |
144 |
|
|
& airPotTempFile |
145 |
|
|
|
146 |
jmc |
1.1 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
147 |
|
|
|
148 |
|
|
#endif /*ALLOW_BULK_FORCE*/ |