1 |
C $Header: /u/gcmpack/MITgcm/pkg/bulk_force/BULKF_PARAMS.h,v 1.1 2003/11/23 01:36:55 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_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 |
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.. 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 |
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 |
|
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 |
_RL cdrag_1, cdrag_2, cdrag_3 |
66 |
_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 |
#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 |
_RL blk_taveFreq |
106 |
|
107 |
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 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
147 |
|
148 |
#endif /*ALLOW_BULK_FORCE*/ |