4 |
#ifdef ALLOW_BULK_FORCE |
#ifdef ALLOW_BULK_FORCE |
5 |
C *==========================================================* |
C *==========================================================* |
6 |
C | BULK_PARAMS.h |
C | BULK_PARAMS.h |
7 |
C | o Header file for BULK_FORC package parameters: |
C | o Header file for BULK_FORCE package parameters: |
8 |
C | - basic parameter ( I/O frequency, etc ...) |
C | - basic parameter ( I/O frequency, etc ...) |
9 |
C | - physical constants (not used in therm_seaice pkg) |
C | - physical constants |
10 |
C *==========================================================* |
C *==========================================================* |
11 |
|
|
12 |
C-- COMMON / BULK_PHYS_CONST / physical constants |
C-- COMMON / BULK_PHYS_CONST / physical constants |
13 |
C.. densities |
C.. densities |
14 |
C rhoa :: density of air (kg/m^3) |
C rhoA :: density of air [kg/m^3] |
15 |
C rhofw :: density of fresh water (kg/m^3) |
C rhoFW :: density of fresh water [kg/m^3] |
16 |
|
C rhoSW :: density of sea water [kg/m^3] |
17 |
C.. specific heats |
C.. specific heats |
18 |
C cpair :: specific heat of air (J/kg/K) |
C cpAir :: specific heat of air [J/kg/K] |
19 |
C cpwv :: specific heat of water vapour (J/kg/KC) |
C cpwv :: specific heat of water vapour [J/kg/K] |
20 |
C .. latent heat |
C.. latent heat |
21 |
C Lvap :: latent heat of vaporizn at 0 C (J/kg) |
C Lvap :: latent heat of vaporization at 0.oC [J/kg] |
22 |
C Lfresh :: latent heat of melting of pure ice (J/kg) |
C Lfresh :: latent heat of melting of pure ice [J/kg] |
23 |
C .. melting |
C Lvap_ice :: latent heat from sublimation [J/kg] |
24 |
C Tf0kel :: Freezing temp of fresh ice in Kelvin = 273.15 |
C.. constants |
25 |
C .. wind drag |
C Tf0kel :: Freezing temp of fresh water in Kelvin = 273.15 |
26 |
C drag_[n] :: n = 1,2,3 coefficients used to evaluate drag coefficient |
C Rgas :: gas constant for dry air [J/kg/K] |
27 |
C .. constants |
C Rvap :: gas constant for H2O vapor [J/kg/K] |
28 |
C stefan :: Stefan-Boltzmann constant (W/m^2 K^4) |
C xkar :: von Karman constant [-] |
29 |
C xkar :: Von Karman constant QQ units? |
C stefan :: Stefan-Boltzmann constant [W/m^2/K^4] |
30 |
C Rvap :: gas constant for H2O vapor (J/kg/K) |
C.. for transfer coefficient |
31 |
C.. Miscellaneous |
C zref :: reference height [m] for transfer coefficient |
32 |
C p0 :: surface pressure (mb) |
C zwd :: height [m] of near-surface wind-speed input data |
33 |
C.. Combinations used for efficiency |
C zth :: height [m] of near-surface air-temp. & air-humid. input |
34 |
C Qcoef :: another constant for latent heat flux |
C cDrag_[n] :: n = 1,2,3 coefficients used to evaluate drag coefficient |
35 |
C .. for bulk formula |
C cStantonS,U :: coefficients used to evaluate Stanton number (for |
36 |
C humid_fac :: const. for the evaluation of the virtual temp. |
C sensib. Heat Flx), under Stable / Unstable stratification |
37 |
C saltsat :: reduction of sat. vapor pressure over salt water |
C cDalton :: coefficient used to evaluate Dalton number (for Evap) |
38 |
C gamma_blk :: adiabatic lapse rate |
C.. for bulk formula |
39 |
C .. for Vince bulk formula QQ check units |
C umin :: minimum wind speed used in bulk-formulae [m/s] |
40 |
C Lvap_ice :: latent heat from sublimation |
C humid_fac :: dry-air - water-vapor molecular mass ratio (minus one) |
41 |
C Rgas :: gas constant for dry air |
C (used to calculate virtual temp.) |
42 |
C Sha :: |
C saltQsFac :: reduction of sat. vapor pressure over salty water |
43 |
C .. emissivities |
C gamma_blk :: adiabatic lapse rate |
44 |
C atm_emissivity :: |
C.. for Long-Wave Radiation |
45 |
C ocean_emissivity :: |
C atm_emissivity :: |
46 |
C snow_emissivity :: |
C ocean_emissivity :: |
47 |
C ice_emissivity :: |
C snow_emissivity :: |
48 |
C .. for BULKF_FORMULA_AIM |
C ice_emissivity :: |
49 |
|
C.. for BULKF_FORMULA_AIM |
50 |
C FWIND0 :: ratio of near-sfc wind to lowest-level wind |
C FWIND0 :: ratio of near-sfc wind to lowest-level wind |
51 |
C CHS :: heat exchange coefficient over sea |
C CHS :: heat exchange coefficient over sea |
52 |
C VGUST :: wind speed for sub-grid-scale gusts |
C VGUST :: wind speed for sub-grid-scale gusts |
53 |
C DTHETA :: Potential temp. gradient for stability correction |
C DTHETA :: Potential temp. gradient for stability correction |
54 |
C dTstab :: potential temp. increment for stability function derivative |
C dTstab :: potential temp. increment for stability function derivative |
55 |
C FSTAB :: Amplitude of stability correction (fraction) |
C FSTAB :: Amplitude of stability correction (fraction) |
56 |
C .. Albedo |
C.. Albedo |
57 |
C ocean_albedo :: ocean surface albedo [0-1] |
C ocean_albedo :: ocean surface albedo [0-1] |
58 |
|
|
59 |
_RL rhoa |
_RL rhoA |
60 |
_RL rhosw |
c _RL rhoSW |
61 |
_RL rhofw |
_RL rhoFW |
62 |
_RL cpair |
_RL cpAir |
63 |
_RL cpwv |
c _RL cpwv |
64 |
_RL Lvap |
_RL Lvap |
65 |
_RL Lfresh |
_RL Lfresh |
66 |
|
c _RL Lvap_ice |
67 |
_RL Tf0kel |
_RL Tf0kel |
68 |
_RL cdrag_1, cdrag_2, cdrag_3 |
_RL Rgas |
69 |
_RL stefan |
c _RL Rvap |
70 |
_RL xkar |
_RL xkar |
71 |
_RL Rvap |
_RL stefan |
72 |
_RL p0 |
_RL zref, zwd, zth |
73 |
_RL Qcoef |
_RL cDrag_1, cDrag_2, cDrag_3 |
74 |
|
_RL cStantonS, cStantonU |
75 |
|
_RL cDalton |
76 |
|
_RL umin |
77 |
_RL humid_fac |
_RL humid_fac |
78 |
_RL saltsat |
_RL saltQsFac |
79 |
_RL gamma_blk |
_RL gamma_blk |
|
_RL Lvap_ice |
|
|
_RL Rgas |
|
|
_RL Sha |
|
80 |
_RL atm_emissivity |
_RL atm_emissivity |
81 |
_RL ocean_emissivity |
_RL ocean_emissivity |
82 |
_RL snow_emissivity |
_RL snow_emissivity |
85 |
_RL FWIND0, CHS, VGUST, DTHETA, dTstab, FSTAB |
_RL FWIND0, CHS, VGUST, DTHETA, dTstab, FSTAB |
86 |
#endif |
#endif |
87 |
_RL ocean_albedo |
_RL ocean_albedo |
|
|
|
88 |
COMMON / BULK_PHYS_CONST / |
COMMON / BULK_PHYS_CONST / |
89 |
& rhoa, rhosw, rhofw, |
& rhoA, rhoFW, |
90 |
& cpair, cpwv, |
& cpAir, Lvap, Lfresh, |
91 |
& Lvap, Lfresh, |
& Tf0kel, Rgas, |
92 |
& Tf0kel, |
& xkar, stefan, |
93 |
& cdrag_1, cdrag_2, cdrag_3, |
& zref, zwd, zth, |
94 |
& stefan, xkar, Rvap, |
& cDrag_1, cDrag_2, cDrag_3, |
95 |
& p0, Qcoef, |
& cStantonS, cStantonU, cDalton, |
96 |
& humid_fac, saltsat, gamma_blk, |
& umin, humid_fac, saltQsFac, gamma_blk, |
|
& Lvap_ice, Rgas, Sha, |
|
97 |
& atm_emissivity, ocean_emissivity, |
& atm_emissivity, ocean_emissivity, |
98 |
& snow_emissivity, ice_emissivity, |
& snow_emissivity, ice_emissivity, |
99 |
#ifdef ALLOW_FORMULA_AIM |
#ifdef ALLOW_FORMULA_AIM |
101 |
#endif |
#endif |
102 |
& ocean_albedo |
& 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 |
C-- COMMON / BULK_PAR_R / real parameter |
112 |
|
C blk_taveFreq :: time-average output frequency [s] |
113 |
|
|
114 |
_RL blk_taveFreq |
_RL blk_taveFreq |
|
|
|
115 |
COMMON / BULK_PAR_R / |
COMMON / BULK_PAR_R / |
116 |
& blk_taveFreq |
& blk_taveFreq |
117 |
|
|
124 |
LOGICAL calcWindStress |
LOGICAL calcWindStress |
125 |
LOGICAL useFluxFormula_AIM |
LOGICAL useFluxFormula_AIM |
126 |
LOGICAL useQnetch, useEmPch |
LOGICAL useQnetch, useEmPch |
|
|
|
127 |
COMMON / BULK_PAR_L / |
COMMON / BULK_PAR_L / |
128 |
& calcWindStress, |
& calcWindStress, |
129 |
& useFluxFormula_AIM, |
& useFluxFormula_AIM, |
130 |
& useQnetch, useEmPch |
& useQnetch, useEmPch |
131 |
|
|
132 |
C-- COMMON / BULK_PAR_C / Character string parameters |
C-- COMMON / BULK_PAR_C / Character string parameters |
133 |
|
|
134 |
CHARACTER*(MAX_LEN_FNAM) AirTempFile |
CHARACTER*(MAX_LEN_FNAM) AirTempFile |
135 |
CHARACTER*(MAX_LEN_FNAM) AirHumidityFile |
CHARACTER*(MAX_LEN_FNAM) AirHumidityFile |
136 |
CHARACTER*(MAX_LEN_FNAM) RainFile |
CHARACTER*(MAX_LEN_FNAM) RainFile |
145 |
CHARACTER*(MAX_LEN_FNAM) CloudFile |
CHARACTER*(MAX_LEN_FNAM) CloudFile |
146 |
CHARACTER*(MAX_LEN_FNAM) SnowFile |
CHARACTER*(MAX_LEN_FNAM) SnowFile |
147 |
CHARACTER*(MAX_LEN_FNAM) airPotTempFile |
CHARACTER*(MAX_LEN_FNAM) airPotTempFile |
|
|
|
148 |
COMMON / BULK_PAR_C / |
COMMON / BULK_PAR_C / |
149 |
& AirTempFile, AirHumidityFile, RainFile, |
& AirTempFile, AirHumidityFile, RainFile, |
150 |
& SolarFile, LongwaveFile, UWindFile, VWindFile, |
& SolarFile, LongwaveFile, UWindFile, VWindFile, |