Parameters adjusted to the height specified in z_wanted. Adjustment assumes the modified log-profile is valid at both the observation (ref_ht_wind) and the wanted height (z_wanted).
Input Parameters:
List of input parameters. Argument numbers begin at 1.
| Parameter | Type | Description | Units | Argument Number |
| astab | int | Atmospheric stability option | none | 18 |
| CONVECT | float | Convective parameter. Recommended value between 0.7 and 1.25. For details see TOGA NOTES #4 | none | 4 |
| CONV_CRIT | float | Convergence criterion | fraction | 5 |
| air_moist_prm | int | Atmospheric moisture parameter index | none | 7 |
| air_moist_val | float | Value of the parameter corresponding to the above index | see below | 8 |
| dyn_in_prm | int | Dynamic input parameter index | none | 1 |
| dyn_in_val | float | Dynamic input value. Usually the mean wind speed at the height (zref) of the anemometer. Other input options are friction velocity (magnitude), wind stress (magnitude), and equivalent neutral wind speed (scatterometer wind speed). | see below | 2 |
| eqv_neut | int | option for type of output winds | none | 32 |
| pressure | float | Atmospheric surface pressure | Pa | 6 |
| ref_ht_wind | float | Height of the wind observations | m | 16 |
| ref_ht_tq | float | Height of the temperature observations. Note: in the current version of the code this must equal to height of the humidity observations. | m | 17 |
| rel_wind_ang | float | Angle of the mean wind relative to the direction that the dominant waves propagate. The positive direction is counter-clockwise when looking down on the surface. | degrees | 3 |
| salinity | float | Salinity | none | 11 |
| sfc_moist_prm | int | Surface moisture parameter index | none | 9 |
| sfc_moist_val | float | Value of the parameter corresponding to the above index | see below | 10 |
| ss_prm | int | Seastate parameter index | none | 12 |
| ss_val | float | Value of the parameter corresponding to the above index | see below | 13 |
| t_air | float | Air temperature at the reference height of the thermometer and humidity sensor | C | 14 |
| t_skin | float | Skin temperature of the water. | C | 15 |
| warn | int | Warning level: 0 warnings, 1 no warnings. | none | 19 |
| z_wanted | float | height to which the values are adjusted | m | 33 |
Options for type of wind output (eqv_neut)
Equivalent neutral wind speed (perhaps better called 'equivalent only when neutral') are used in applications such as scatterometry. The winds are height adjusted, with friction velocity (and momentum roughness length) determined form non-neutral conditions. However, winds are height adjusted with the atmospheric stability term set to zero. Consequently, friction velocity and roughness length are determined using the atmospheric stability parameterization, but the hight adjustment is made for a neutral atmosphere. If this option is set to '0', then it has no impact on the height adjustment. If this option is set to '1', the potential temperature and specific humidity as well as the winds are output as equivalent neutral values. WARNING: equivalent neutral winds are not equivalent to neutral winds [Liu and Tang, 1996; Vershell et al., 1998].
| eqv_neut | Description | Units |
| 0 | winds | none |
| 1 | equivalent neutral winds | none |
Options for dynamic input:
Typically wind speed is used as an input to boundary-layer models. However, scatterometers are now
producing 'observations' of friction velocity and equivelent neutral wind speed.
| dyn_in | Description | Units |
| 0 | Wind speed, relative to surface current | m/s |
| 1 | Friction velocity (magnitude) | m/s |
| 2 | Surface wind stress (magnitude) | N/m^2 |
| 3 | Equivalent neutral wind speed (relative to the surface current) | m/s |
Options for atmospheric stability condition:
The atmospheric stability in the boundary-layer can be assumed
to neutral, or it can be calculated input parameters.
| astab | Description | Units |
| 0 | Atmospheric stability is assumed to be neutral | none |
| 1 | Stability is calculated | none |
Options for seastate parameterizations:
There are six possible seastate assumptions: any one of the
following can be treated as known: wind-wave stability parameter
(set to 1.0 for local equilibrium), phase speed, or wave age, significant wave height, significant slope, or period of the dominant waves.
Caution: in many cases, these wave characteristics will correspond to swell rather than the phase speed of locally wind induced waves.
| ss_prm | Parameter treated as known (ss_val) | Units |
| 0 | Wind-wave stability parameter | none |
| 1 | Phase speed of the dominant waves. Note: in many cases, this phase speed will correspond to the swell rather than the phase speed of locally wind induced waves. Use of the wrong phase speed can lead to large overestimations of fluxes. | m/s |
| 2 | Wave age the dominant waves (cp/u*) | none |
| 3 | Significant wave height (Hs) | m |
| 4 | Significant slope (Hs/l) | none |
| 5 | Period of the dominant waves (Tp) | s |
Options for atmospheric moisture input:
Choose the moisture parameter that is easiest for you to deal
with:
| air_moist_prm | Parameter for moisture of air (air_moist_val) | Units |
| 0 | Specific humidity at the reference height of the thermometer and humidity sensor | g vapor / g air |
| 1 | Relative humidity | fraction |
| 2 | Dew point temperature | C |
| 3 | Wet bulb temperature | C |
Options for surface moisture input:
Choose the moisture parameter that is easiest for you to deal
with:
| sfc_moist_prm | Parameter for moisture of air (sfc_moist_val) | Units |
| 0 | Specific humidity 'at' (near) the surface | g vapor / g air |
| 1 | Relative humidity | fraction |
| 2 | Dew point temperature | C |
| 3 | Wet bulb temperature | C |
Model Output:
Vector components are calculated parallel and perpendicular to
the direction in which the dominant waves are propagating. The
first component is parallel the direction of wave
propagation, and the second component is perpendicular to the
first (while looking down it is 90 counter-clockwise from the
first component; i.e., in a right handed coordinate system with
the positive vertical axis pointing upward). For most applications
there will be insufficient wave information, requiring the assumption
of local wind-wave equilibrium. This assumption implies that
the wind and the waves are moving in the same direction; which
results in the first component of the vectors being parallel to
the wind direction, and the second component being zero.
All output is single precision floating point.
The routine returns a integer value (i.e., a warning flag). Positive values indicate a lack of specific problems.
If there are problems with missing input, non-convergence within the algorithm, or if the modeled physics obviously fails to apply, then the output is set to -1.
For example, if the thickness of the boundary layer is too small (i.e., the absolute value of the Obhukov scale length less than or equal to 1 m) then the warning flag is set at -1.
| Parameter | Type | Description | Units | Argument Number |
| shf | float | sensible heat flux | W m-2 | 20 |
| lhf | float | latent heat flux | W m-2 | 21 |
| tau | vector float | stress vector. There are more details on the conversion to zonal and meridional components. | N m-2 | 22 |
| u_star | vector float | friction velocity (u*) | m s-1 | 23 |
| t_star | float | scaling term for potential temperature (T*) | C | 24 |
| q_star | float | scaling parameter for moisture (q*) | none | 25 |
| z_over_L | float | dimensionless Monin-Obhukov scale length | none | 26 |
| wave_age | float | wave age, cp/u* | none | 27 |
| dom_phs_spd | float | dominant phase speed of gravity waves | m s-1 | 28 |
| h_sig | float | significant wave height | m | 29 |
| ww_stab | float | wind-wave evolution parameter | none | 30 |
| zo_m | vector float | momentum roughness length | m | 31 |
| u_at_z | float | wind speed at the specified height | m s-1 | 33 |
| t_at_z | float | potential temperature at the specified height | oC | 34 |
| q_at_z | float | specific humidity at the specified height | kg kg-1 | 35 |
Last update: July. 28, 1998