Model help:TopoFlow-Meteorology: Difference between revisions
No edit summary |
No edit summary |
||
Line 98: | Line 98: | ||
!Parameter!!Description!!Unit | !Parameter!!Description!!Unit | ||
|-valign="top" | |-valign="top" | ||
|width="20%"| | |width="20%"|T_air type | ||
|width="60%"| | |width="60%"|allowed input types (Scalar/Grid/Time_Series/Grid_Sequence) | ||
|width="20%"| | |width="20%"| - | ||
|- | |||
|T_air | |||
|temperature of air | |||
|degree C | |||
|- | |||
|T_surf type | |||
|allowed input types (Scalar/Grid/Time_Series/Grid_Sequence) | |||
| | |||
|- | |||
|T_surf | |||
|temperature of surface | |||
|degree C | |||
|- | |||
|RH type | |||
|allowed input types (Scalar/Grid/Time_Series/Grid_Sequence) | |||
| | |||
|- | |||
|RH | |||
|relative humidity | |||
| | |||
|- | |||
|p0 type | |||
|allowed input types (Scalar/Grid/Time_Series/Grid_Sequence) | |||
| | |||
|- | |||
|p0 | |||
|atmospheric pressure | |||
| mbar | |||
|- | |||
|uz type | |||
|allowed input types (Scalar/Grid/Time_Series/Grid_Sequence) | |||
| | |||
|- | |||
|uz | |||
|wind velocity at height z | |||
| m / s | |||
|- | |||
|z type | |||
|allowed input types (Scalar/Grid/Time_Series/Grid_Sequence) | |||
| | |||
|- | |||
|z | |||
|wind reference height | |||
| m | |||
|- | |||
|z0_air type | |||
|allowed input types (Scalar/Grid/Time_Series/Grid_Sequence) | |||
| | |||
|- | |||
|z0_air | |||
|surface roughness length scale for wind | |||
| m | |||
|- | |||
|} | |} | ||
Line 108: | Line 161: | ||
!Parameter!!Description!!Unit | !Parameter!!Description!!Unit | ||
|-valign="top" | |-valign="top" | ||
|width="20%"| | |width="20%"|albedo type | ||
|width="60%"| | |width="60%"|allowed input types (Scalar/Grid/Time_Series/Grid_Sequence) | ||
|width="20%"| | |width="20%"| | ||
|- | |||
|albedo | |||
|surface albedo in [0,1] | |||
| | |||
|- | |||
|emissivity type | |||
|allowed input types (Scalar/Grid/Time_Series/Grid_Sequence) | |||
| | |||
|- | |||
|emissivity | |||
|surface emissivity in [0,1] | |||
| | |||
|- | |||
|dust atten. type | |||
|allowed input types (Scalar/Grid/Time_Series/Grid_Sequence) | |||
| | |||
|- | |||
|dust attenuation | |||
|dust attenuation factor in [0,0.3] | |||
| | |||
|- | |||
|cloud factor type | |||
|allowed input types (Scalar/Grid/Time_Series/Grid_Sequence) | |||
| | |||
|- | |||
|cloud factor | |||
|cloudiness factor in [0,1] | |||
| | |||
|- | |||
|canopy factor type | |||
|allowed input types (Scalar/Grid/Time_Series/Grid_Sequence) | |||
| | |||
|- | |||
|canopy factor | |||
|canopy coverage factor in [0,1] | |||
| | |||
|- | |||
|slope grid file | |||
|flat binary, row-major file with grid of 4-byte slopes | |||
| | |||
|- | |||
|aspect grid file | |||
|flat binary, row-major file with grid of 4-byte aspects | |||
| | |||
|- | |||
|time zone offset | |||
|time zone offset from GMT | |||
| | |||
|- | |||
|start month | |||
| | |||
| | |||
|- | |||
|start_day | |||
| | |||
| | |||
|- | |||
|start_hour | |||
| | |||
| | |||
|- | |||
|} | |} | ||
Line 118: | Line 232: | ||
!Parameter!!Description!!Unit | !Parameter!!Description!!Unit | ||
|-valign="top" | |-valign="top" | ||
|width="20%"| | |width="20%"|Save grid timestep | ||
|width=" | |width="60%"|time interval between saved grids | ||
| | |width="20%"|sec | ||
|- | |||
|Save e_air grids toggle | |||
| option to save grids of vapor pressure (in air) | |||
| - | |||
|- | |||
|Save e_air grids file | |||
| filename for grid stack of vapor pressure (in air) | |||
| mbar | |||
|- | |||
|Save e_surf grids toggle | |||
| option to save grids of vapor pressure (at surface) | |||
| - | |||
|- | |||
|Save e_surf grids file | |||
| filename for grid stack of vapor pressure (surf) | |||
| mbar | |||
|- | |||
|Save Qn_SW grids toggle | |||
| option to save grids of net short wave radiation | |||
| - | |||
|- | |||
|Save Qn_SW grids file | |||
| filename for grid stack of short wave radiation | |||
| W / m<sup>2</sup> | |||
|- | |||
|Save Qn_LW grids toggle | |||
| option to save grids of net longwave radiation | |||
| - | |||
|- | |||
|Save Qn_LW grids file | |||
| filename for grid stack of net longwave radiation | |||
| W / m<sup>2</sup> | |||
|- | |||
|Save em_air grids toggle | |||
| option to save grids of air emissivity | |||
| | |||
|- | |||
|Save em_air grids file | |||
| filename for grid stack of air emissivity in [0,1] | |||
| | |||
|- | |||
|} | |} | ||
Line 128: | Line 283: | ||
!Parameter!!Description!!Unit | !Parameter!!Description!!Unit | ||
|-valign="top" | |-valign="top" | ||
|width="20%"| | |width="20%"|Save pixels timestep | ||
|width=" | |width="60%"|time interval between time series values | ||
| | |width="20%"|sec | ||
|- | |||
|Save e_air pixels toggle | |||
| option to save time series of sat. vapor pressure (in air) | |||
| - | |||
|- | |||
|Save e_air pixels file | |||
| filename for time series of sat. vapor pressure (in air) | |||
| mbar | |||
|- | |||
|Save e_surf pixels toggle | |||
| option to save time series of sat. vapor pressure (at surface) | |||
| - | |||
|- | |||
|Save e_surf pixels file | |||
| filename for time series of sat. vapor pressure (surf) | |||
| mbar | |||
|- | |||
|Save Qn_SW pixels toggle | |||
| option to save time series of net shortwave radiation | |||
| - | |||
|- | |||
|Save Qn_SW pixels file | |||
| filename for time series of short wave radiation | |||
| W / m<sup>2</sup> | |||
|- | |||
|Save Qn_LW pixels toggle | |||
| option to save time series of net longwave radiation | |||
| - | |||
|- | |||
|Save Qn_LW pixels file | |||
| filename for time series of net longwave radiation | |||
| W / m<sup>2</sup> | |||
|- | |||
|Save em_air pixels toggle | |||
| option to save time series of air emissivity | |||
| | |||
|- | |||
|Save em_air pixels file | |||
| filename for time series of air emissivity in [0,1] | |||
| | |||
|- | |||
|} | |} | ||
<headertabs/> | <headertabs/> | ||
Revision as of 15:03, 4 May 2011
TopoFlow-Meteorology
The module is the meteorology process component for a D8-based, spatial hydrologic model
Model introduction
This component reads a variety of variables for the atmosphere and for the land surface from input files or as simple scalars. It then computes many additional variables, such as vapor pressure, eair, and net shortwave (solar) radiation, QnSW, using built-in [shortwave radiation] and [longwave radiation] calculators that are based on celestial mechanics and widely-used empirical relationships. These additional variables are needed by the Snowmelt → Energy Balance and Evaporation → Energy Balance components. Direct, diffuse and back-scattered radiation fluxes are all modeled. Properties of the atmosphere (e.g. precipitation rate, P, air temperature, Tair, relative humidity, RH, and dust attenuation, γ), are used as well as surface/topographic properties (e.g. slope angle, aspect angle and surface albedo. The approach used here closely follows the one outlined in Appendix E of Dingman (2002)[1]. However, instantaneous vs. day-integrated radiation fluxes are used and the optical air mass is modeled using the widely used method of Kasten and Young (1989)[2].
Model parameters
Uses ports
This will be something that the CSDMS facility will add
Provides ports
This will be something that the CSDMS facility will add
Main equations
- Absolute optical air mass equation
[math]\displaystyle{ m_{abs} \left ( \gamma \right ) = \rho _{0} \int _{0} ^ \infty \left ( \rho / \rho _{0} \right ) \{ 1 - [ 1 + 2 \delta _{0} \left ( 1 - \rho / \rho _{0} \right ) ] \ast [ \cos \gamma / \left ( 1 + h / R \right ) ] ^2 \} ^ {\frac{-1}{2}} dh }[/math] (1)
[math]\displaystyle{ m \left ( \gamma \right ) = m _{abs} \left ( \gamma \right ) / m_{abs} \left ( 90^o \right ) }[/math] (2)
Symbol | Description | Unit |
---|---|---|
ρH2O | density of water | kg / m^3 |
Cpair | heat capacity of air | J / kg / K |
ρair | density of air | kg / m^3 |
precip. rate | precipitation rate | mm / hr |
Tair | air temperature | deg C |
Tsurface | surface temperature | deg C |
RH | relative humidity (in [0,1]) | - |
p0 | atmospheric pressure | mbar |
z | reference height for uz | m |
uz | wind velocity at reference height | m / s |
z0 | surface roughness length scale for wind | m |
albedo | surface albedo in [0,1] | - |
em_air | surface emissivity in [0,1] | - |
dust atten. in [0, 0.3] | dust attenuation factor | - |
cloud factor | cloud factor in [0,1], 0 for no clouds | - |
canopy factor | forest canopy factor in [0,1], 0 for no canopy | - |
slope | topographic slope in [0, infinity] | - |
slope grid file | as flat binary, row-major file with 4-byte floats | - |
aspect angle | aspect angle [radians] in [0,1] | - |
aspect grid file | as flat binary, row-major file with 4-byte floats | - |
time zone offset | offset, in hours, from Greenwich Mean Time (GMT), negative for east of prime meridian, positive otherwise | - |
start month | start month for solar radiation calculations | - |
start day | start day for solar radiation calculations | - |
start hour | start hour for solar radiation calculations [decimal, 24-hour clock] | - |
h | height above mean sea level | m |
ρ | equals to ρ(h), air density at height h | kg / m3 |
ρ0 | air density at h = 0 | kg / m3 |
δ0 | equals to n0 - 1 | - |
n0 | refractive index for air at 0.7μm wavelength at h = 0 | - |
R | mean earth radius | m |
m(γ) | the relative optical air mass at solar elevation γ | - |
mabs(γ) | absolute optical air mass | - |
Output
Symbol | Description | Unit |
---|---|---|
eair | vapor pressure of air | mbar |
esurf | vapor pressure at the surface | mbar |
QnSW | net shortwave radiation | W / m^2 |
QnLW | net longwave radiation, equals to LWin - LWout | W / m^2 |
emair | air emissivity in [0,1] | - |
Notes
- Note on input parameters
For each input variable, you may choose from the droplist of data types. For the "Scalar" data type, enter a numeric value with the units indicated in the dialog. For the other data types, enter a filename. Values in files must also use the indicated units.
Single grids and grid sequences are assumed to be stored as RTG and RTS files, respectively. Time series are assumed to be stored as text files, with one value per line. For a time series or grid sequence, the time between values must coincide with the timestep provided.
For DEMs with pixel geometry and bounding box given in terms of Geographic coordinates, the latitude and longitude of each pixel is used in the calculations. For DEMs with a "fixed-length" pixel geometry (e.g. UTM coordinates), which tend to span smaller areas, the dialog prompts for a single lat/lon pair to be used in the calculations.
- Note on Equations
time zone:Boundaries of time zones can be very irregular and a time zone map should be consulted if you are unsure. The time zone is not simply a function of the longitude. You can select an adjacent time zone to include the effect of Daylight Savings Time. Time zones with non-integer offsets from GMT are not yet supported.
slope:Topographic slopes (not slope angles) are specified as dimensionless numbers [m/m]. A RiverTools grid (RTG file) with extension "_slope.rtg", "_mf-slope.rtg" or "_dinf-slope.rtg" can be used.
aspect:Aspect is specified as an angle measured in radians counter-clockwise from due east (the standard convention). A RiverTools grid (RTG file) with extension "_mf-angle.rtg" or "_dinf-angle.rtg" can be used for the (continuous-angle) aspect grid.
QSW is set to zero between the times of local sunset and local sunrise, so frames in the RTS file that correspond to nighttime hours will contain only zeros.
Examples
An example run with input parameters, BLD files, as well as a figure / movie of the output
Follow the next steps to include images / movies of simulations:
- Upload file: http://csdms.colorado.edu/wiki/Special:Upload
- Create link to the file on your page: [[Image:<file name>]].
See also: Help:Images or Help:Movies
Developer(s)
References
Dingman, S.L (2002) Physical Hydrology, 2nd ed., Prentice Hall, New Jersey. (see Appendix E)
Kasten and Young (1989) Revised optical air mass tables and approximation formula. Applied Optics, 28 (22): 4735~4738. (for the optical air mass equation)
Liston, G. *******
Marks and Dozier (1992) ******* Water Resources Research.
Whitman, A.M. (2003) A simple expression for the equation of time, online document, http://www.sunspot.noao.edu/sunspot/pr/answerbook/expl-5.html