time   calls  line
             1   14 g           = 9.8;       % acceleration due to gravity [m/s^2] 
             1   15 sigmaSB     = 5.6697e-8; % Stefan-Boltzmann constant [W/m^2/K^4] 
             1   16 eps_air     = 0.62197;   % molecular weight ratio (water/air) 
             1   17 gas_const_R = 287.04;    % gas constant for dry air [J/kg/K] 
             1   18 CtoK        = 273.16;    % conversion factor for [C] to [K] 
                 19 
                 20 
                 21 % ------- meteorological constants
                 22 
             1   23 kappa          = 0.4;    % von Karman's constant 
             1   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. 
             1   29 R_roughness   = 0.11;    % limiting roughness Reynolds # for aerodynamically  
                 30                          % smooth flow         
                 31                          
                 32                         
                 33 % ------ defaults suitable for boundary-layer studies
                 34 
             1   35 cp            = 1004.7;   % heat capacity of air [J/kg/K] 
             1   36 rho_air       = 1.22;     % air density (when required as constant) [kg/m^2] 
             1   37 Ta_default    = 10;       % default air temperature [C] 
             1   38 P_default     = 1020;     % default air pressure for Kinneret [mbars] 
             1   39 psych_default = 'screen'; % default psychmometer type (see relhumid.m) 
             1   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 
             1   47 CVB_depth     = 600; % depth of convective boundary layer in atmosphere [m] 
             1   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)
             1   51 beta_conv     = 1.25;% scaling constant for gustiness 
                 52 
                 53 
                 54 % ------ short-wave flux calculations
                 55 
             1   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 
             1   64 emiss_lw = 0.985;     % long-wave emissivity of ocean from Dickey et al 
                 65                       % (1994), J. Atmos. Oceanic Tech., 11, 1057-1076.
                 66 
             1   67 bulkf_default = 'berliand';  % default bulk formula when downward long-wave 

Other subfunctions in this file are not included in this listing.