1 | % AS_CONSTS: returns values of many constants used in AIR_SEA TOOLBOX
|
---|
2 | % AS_CONSTS: returns values of many constants used in the AIR_SEA
|
---|
3 | % TOOLBOX. At end of this file are values of constants from COARE
|
---|
4 | % to be used when running the test program t_hfbulktc.m
|
---|
5 |
|
---|
6 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
---|
7 | % 8/19/98: version 1.1 (contributed by RP)
|
---|
8 | % 4/7/99: version 1.2 (revised to include COARE test values by AA)
|
---|
9 | % 8/5/99: version 2.0
|
---|
10 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
---|
11 |
|
---|
12 | % ------- physical constants
|
---|
13 |
|
---|
14 | g = 9.8; % acceleration due to gravity [m/s^2]
|
---|
15 | sigmaSB = 5.6697e-8; % Stefan-Boltzmann constant [W/m^2/K^4]
|
---|
16 | eps_air = 0.62197; % molecular weight ratio (water/air)
|
---|
17 | gas_const_R = 287.04; % gas constant for dry air [J/kg/K]
|
---|
18 | CtoK = 273.16; % conversion factor for [C] to [K]
|
---|
19 |
|
---|
20 |
|
---|
21 | % ------- meteorological constants
|
---|
22 |
|
---|
23 | kappa = 0.4; % von Karman's constant
|
---|
24 | Charnock_alpha = 0.011; % Charnock constant (for determining roughness length
|
---|
25 | % at sea given friction velocity), used in Smith
|
---|
26 | % formulas for drag coefficient and also in Fairall
|
---|
27 | % and Edson. use alpha=0.011 for open-ocean and
|
---|
28 | % alpha=0.018 for fetch-limited (coastal) regions.
|
---|
29 | R_roughness = 0.11; % limiting roughness Reynolds # for aerodynamically
|
---|
30 | % smooth flow
|
---|
31 |
|
---|
32 |
|
---|
33 | % ------ defaults suitable for boundary-layer studies
|
---|
34 |
|
---|
35 | cp = 1004.7; % heat capacity of air [J/kg/K]
|
---|
36 | rho_air = 1.22; % air density (when required as constant) [kg/m^2]
|
---|
37 | Ta_default = 10; % default air temperature [C]
|
---|
38 | P_default = 1020; % default air pressure for Kinneret [mbars]
|
---|
39 | psych_default = 'screen'; % default psychmometer type (see relhumid.m)
|
---|
40 | Qsat_coeff = 0.98; % satur. specific humidity coefficient reduced
|
---|
41 | % by 2% over salt water
|
---|
42 |
|
---|
43 |
|
---|
44 | % the following are useful in hfbulktc.m
|
---|
45 | % (and are the default values used in Fairall et al, 1996)
|
---|
46 |
|
---|
47 | CVB_depth = 600; % depth of convective boundary layer in atmosphere [m]
|
---|
48 | min_gustiness = 0.5; % min. "gustiness" (i.e., unresolved fluctuations) [m/s]
|
---|
49 | % should keep this strictly >0, otherwise bad stuff
|
---|
50 | % might happen (divide by zero errors)
|
---|
51 | beta_conv = 1.25;% scaling constant for gustiness
|
---|
52 |
|
---|
53 |
|
---|
54 | % ------ short-wave flux calculations
|
---|
55 |
|
---|
56 | Solar_const = 1368.0; % the solar constant [W/m^2] represents a
|
---|
57 | % mean of satellite measurements made over the
|
---|
58 | % last sunspot cycle (1979-1995) taken from
|
---|
59 | % Coffey et al (1995), Earth System Monitor, 6, 6-10.
|
---|
60 |
|
---|
61 |
|
---|
62 | % ------ long-wave flux calculations
|
---|
63 |
|
---|
64 | emiss_lw = 0.985; % long-wave emissivity of ocean from Dickey et al
|
---|
65 | % (1994), J. Atmos. Oceanic Tech., 11, 1057-1076.
|
---|
66 |
|
---|
67 | bulkf_default = 'berliand'; % default bulk formula when downward long-wave
|
---|
68 | % measurements are not made.
|
---|
69 |
|
---|
70 |
|
---|
71 | % ------ constants used for COARE; to use simply delete the %
|
---|
72 | % g = 9.7803; % acceleration due to gravity [m/s^2]
|
---|
73 | % sigmaSB = 5.67e-8; % Stefan-Boltzmann constant [m^2/K^4]
|
---|
74 | % gas_const_R = 287.1; % gas constant for dry air [J/kg/K]
|
---|
75 | % cp = 1004.67; % heat capacity of air [J/kg/K]
|
---|
76 | % beta_conv = 1.20; % scaling constant for gustiness
|
---|
77 | % emiss_lw = 0.97; % long-wave emissivity
|
---|