Difference between revisions of "Coupling Water-Atmosphere User Manual"
From MohidWiki
(→Heat Fluxes) |
(→Heat Fluxes) |
||
| Line 38: | Line 38: | ||
If one wants to simulate heat fluxes on water interface and account its role in water temperature some keywords and properties must be added to data files. | If one wants to simulate heat fluxes on water interface and account its role in water temperature some keywords and properties must be added to data files. | ||
The next image ilustrates heat fluxes in MOHID and atmosphere, interface water-air and water properties that are called to compute them. | The next image ilustrates heat fluxes in MOHID and atmosphere, interface water-air and water properties that are called to compute them. | ||
| − | [[Image:FluxosCalor_MOHID.png|thumb|center| | + | [[Image:FluxosCalor_MOHID.png|thumb|center|400px|Water-Atmosphere heat fluxes and related properties in MOHID]] |
==== WaterProperties Data File ==== | ==== WaterProperties Data File ==== | ||
Revision as of 17:03, 8 November 2010
Contents
Overview
This user manual intends to help MOHID users to construct properties in the atmospheric and interface water-air compartments, in a way that his modeling needs are fulfilled. The document tries to answer to three main processes in what atmosphere takes part: i) wind forcing; ii) heat fluxes between air and water; iii) mass fluxes between air and water.
Momentum fluxes, Oxygen and Carbon fluxes are not yet considered in this user manual.
Heat fluxes between air and water command water warming and cooling processes and thermal stratification formation. Wind can play an important role on surface stress and on mixture depth control, destroying stratification. Thus, when thermal stratification in natural systems occurs, and 3D model is used, both processes should be simulated.
Mass fluxes can be important if level variations due to precipitation and evaporation or water abstraction are important.
Inputs Required
Wind Forcing
If one wants to simulate wind stress on water interface and account its role in hydrodynamic circulation some keywords and properties must be added to data files.
Hydrodynamic Data File
- Keyword “WIND” value has to be equal 1 (one): Ex: WIND : 1.
InterfaceWaterAir Data File
- Property “wind stress X” and Property “wind stress Y” must be defined. These properties can be a constant value, entered as timeserie or computed by the model (REMAIN_CONSTANT : 0).
Atmosphere Data File
- If Properties “wind stress X” and “wind stress Y” in InterfaceWaterAir data file are to be computed, then Property “wind velocity X” and Property “wind velocity Y” must occur in atmosphere data file. These properties can be a constant value, entered as timeserie or computed by the model (REMAIN_CONSTANT : 0).
- If Property “wind velocity X” and Property “wind velocity Y” are to be computed, this is done with wind modulus and wind angle. Then Property “wind modulus” and Property “wind angle” must occur in atmosphere data file. These properties can be a constant value or entered as timeserie.
Heat Fluxes
In MOHID heat fluxes in water column are computed with:
i) a heat source: solar radiation that enters through the water surface (surface radiation - see description in Module InterfaceWaterAir) suffering a decay with depth (light extinction).
ii) a boundary condition for surface: non solar flux (see description in Module InterfaceWaterAir)
iii) a boundary condition for bottom: no flux, all radiation reaching the bottom is transformed in heat
If one wants to simulate heat fluxes on water interface and account its role in water temperature some keywords and properties must be added to data files.
The next image ilustrates heat fluxes in MOHID and atmosphere, interface water-air and water properties that are called to compute them.
WaterProperties Data File
- Property “temperature” must be defined because it is the state variable in water column affected by heat fluxes. Inside this property, block keyword “SURFACE_FLUXES” value has to be equal 1 (one) so that water-air fluxes can be considered.
InterfaceWaterAir Data File
- Property “surface radiation” must be defined in InterfaceWaterAir data file. This property accounts for solar flux trough the water surface and can be a constant value, entered as timeserie or computed by the model (REMAIN_CONSTANT : 0).
- Property “non solar flux” must be defined in InterfaceWaterAir data file. This property accounts for non solar flux (latent heat, sensible heat, etc) trough the water surface and can be a constant value, entered as timeserie or computed by the model (REMAIN_CONSTANT : 0).
- If Property “non solar flux” has to be computed than this is a balance between latent heat, sensible heat and long wave radiation trough the water surface. Thus Property “latent heat”, Property “sensible heat” and Property “net long wave radiation” must be defined in InterfaceWaterAir data file. These properties can be a constant value, entered as timeserie or computed by the model (REMAIN_CONSTANT : 0).
- Remark: latent heat and sensible heat depend on water temperature so they should be computed by the model. Net long wave radiation is a balance between downward long wave radiation (emission from the atmosphere) and upward long wave radiation (emission from water) and should be computed by the model.
- If Property “net long wave radiation” has to be computed than Property "downward long wave radiation" and Property "upward long wave radiation" must be defined in InterfaceWaterAir data file. These properties can be a constant value, entered as timeserie or computed by the model (REMAIN_CONSTANT : 0).
- Remark: upward long wave radiation depends on water temperature so it should be computed by the model.
Atmosphere Data File
- If Property “surface radiation” in InterfaceWaterAir is to be computed by the model, Property “solar radiation” must be defined in Atmosphere data file. This property can be a constant value, entered as timeserie or computed by the model (REMAIN_CONSTANT : 0). If this property is to be computed by the model, in Atmosphere data file keyword RADIATION_METHOD (global variable) must define if it is MOHID method (value 1 - by default) or CE-QUAL based (value 2). See Solar Radiation theory for more details.
- If Property “solar radiation” is to be computed also Property “cloud cover” is needed in Atmosphere data file. If Cloud Cover is computed, Keyword CLOUD_COVER_METHOD (global variable) must define if it is computed from sun hours (value 1), from radiation (value 3) or from random values (value 2 - by default).
- Remark: Solar Radiation is a common property in INAG meteorological stations so, if available, data should be entered in timeserie format instead of computing with the model. Notice that if radiation is computed the model needs cloud cover and if cloud cover is computed (from radiation)it needs radiation.
- If Property “sensible heat” in InterfaceWaterAir is to be computed it is needed air temperature and wind velocity for the formulation . Thus, Property “air temperature” and Property “wind velocity X” and Property “wind velocity Y” must be defined in Atmosphere data file. These properties can be a constant value or entered as timeserie; wind velocity X and wind velocity Y can in addition be computed from Property “wind modulus” and Property “wind angle”, as stated previously (see 2.1.3). These last properties can be a constant value or entered as timeserie.
- If Property “latent heat” in InterfaceWaterAir is to be computed properties needed are the same as in Property “sensible heat” plus humidity. Thus, also Property “humidity” has to be added in Atmosphere data file. This property can be a constant value or entered as timeserie.
- If Property “downward long wave radiation” in InterfaceWaterAir is to be computed, then Property “air temperature” and Property “cloud cover” need to be defined in Atmosphere data file. These properties can be a constant value or entered as timeserie; cloud cover can in addition be computed by the model. If cloud cover is computed, Keyword CLOUD_COVER_METHOD (global variable) must define if it is computed from sun hours (value 1), from radiation (value 3) or from random values (value 2 - by default).
Mass Fluxes
Mass fluxes should be computed if level variations due to precipitation, evaporation or water abstraction are important. Surface water flux is defined as the balance between Precipitation, Evaporation and Irrigation.
Hydrodynamic Data File
- Keyword SURFACEWATERFLUX value has to be 1 (one). Ex: SURFACEWATERFLUX : 1.
InterfaceWaterAir Data File
- Property “surface water flux” must be present in InterfaceWaterAir data file. This property can be a constant value, entered as timeserie or computed by the model.
- If Property “surface water flux” is to be computed by the model then the three properties needed for the balance can appear: Precipitation, Irrigation and Evaporation (this is not mandatory, properties that not appear, are not considered). Thus, Property “evaporation” can be defined. This property can be a constant value, entered as timeserie or can be computed by the model from latent heat. As so, if Property “evaporation” is to be computed by the model, also Property “latent heat” must be defined (see 2.2.2).
AtmosphereDataFile
- If Property “surface water flux” in InterfaceWaterAir data file is to be computed by the model then the three properties needed for the balance can appear: Precipitation, Irrigation and Evaporation (this is not mandatory, properties that not appear, are not considered). Thus, Property “precipitation” and Property “irrigation” can be defined. These two properties can be a constant value or entered as timeserie.
Example Data Files
InterfaceWaterAir Data File
OUTPUT_TIME : 0. 86400 <begin_rugosity> INITIALIZATION_METHOD : CONSTANT DEFAULTVALUE : 0.0025 REMAIN_CONSTANT : 0 OUTPUT_HDF : 1 TIME_SERIE : 0 <end_rugosity> <beginproperty> NAME : wind stress X UNITS : N/m2 DESCRIPTION : Calculated wind stress X DEFAULTVALUE : 0 REMAIN_CONSTANT : 0 TYPE_ZUV : z OUTPUT_HDF : 1 TIME_SERIE : 0 <endproperty> <beginproperty> NAME : wind stress Y UNITS : N/m2 DESCRIPTION : Calculated wind stress Y DEFAULTVALUE : 0 TYPE_ZUV : z REMAIN_CONSTANT : 0 OUTPUT_HDF : 1 TIME_SERIE : 0 <endproperty> <beginproperty> NAME : latent heat UNITS : W/m^2 DESCRIPTION : Calculated latent heat FILE_IN_TIME : NONE REMAIN_CONSTANT : 0 DEFAULTVALUE : 0. TIME_SERIE : 0 OUTPUT_HDF : 1 <endproperty> <beginproperty> NAME : sensible heat UNITS : W/m^2 DESCRIPTION : Calculated sensible heat FILE_IN_TIME : NONE REMAIN_CONSTANT : 0 DEFAULTVALUE : 0. TIME_SERIE : 0 OUTPUT_HDF : 1 <endproperty> <beginproperty> NAME : surface radiation UNITS : W/m^2 DESCRIPTION : Calculated infrared radiation ALBEDO : 0.05 FILE_IN_TIME : NONE REMAIN_CONSTANT : 0 DEFAULTVALUE : 0. TIME_SERIE : 0 OUTPUT_HDF : 1 <endproperty> <beginproperty> NAME : net long wave radiation UNITS : W/m^2 DESCRIPTION : Calculated net long wave radiation FILE_IN_TIME : NONE REMAIN_CONSTANT : 0 DEFAULTVALUE : 0. TIME_SERIE : 0 OUTPUT_HDF : 1 <endproperty> <beginproperty> NAME : downward long wave radiation UNITS : W/m^2 DESCRIPTION : downward long wave radiation data FILE_IN_TIME : HDF FILENAME : ..\..\GeneralData\Atmosphere\MM5_20070606_Portugal.hdf5 DEFAULTVALUE : 0. TIME_SERIE : 0 OUTPUT_HDF : 1 <endproperty> <beginproperty> NAME : upward long wave radiation UNITS : W/m^2 DESCRIPTION : Calculated upward long wave radiation FILE_IN_TIME : NONE REMAIN_CONSTANT : 0 DEFAULTVALUE : 0. TIME_SERIE : 0 OUTPUT_HDF : 1 <endproperty> <beginproperty> NAME : non solar flux UNITS : W/m^2 DESCRIPTION : Calculated infrared radiation FILE_IN_TIME : NONE REMAIN_CONSTANT : 0 DEFAULTVALUE : 0. TIME_SERIE : 0 OUTPUT_HDF : 1 <endproperty>
Atmosphere Data File
RUGOSITY : 0.0025 OUTPUT_TIME : 0. 86400 <beginproperty> NAME : wind velocity X UNITS : m/s DESCRIPTION : wind velocity X DEFAULTVALUE : 0. FILE_IN_TIME : HDF FILENAME : ..\..\GeneralData\Atmosphere\MM5_20070606_Portugal.hdf5 TIME_SERIE : 0 OUTPUT_HDF : 1 <endproperty> <beginproperty> NAME : wind velocity Y UNITS : m/s DESCRIPTION : wind velocity Y DEFAULTVALUE : 0. FILE_IN_TIME : HDF FILENAME : ..\..\GeneralData\Atmosphere\MM5_20070606_Portugal.hdf5 OUTPUT_HDF : 1 <endproperty> <beginproperty> NAME : solar radiation UNITS : W/m^2 DESCRIPTION : meteoIST Solar Radiation FILE_IN_TIME : HDF FILENAME : ..\..\GeneralData\Atmosphere\MM5_20070606_Portugal.hdf5 DEFAULTVALUE : 0.0 TIME_SERIE : 0 OUTPUT_HDF : 1 <endproperty> <beginproperty> NAME : air temperature UNITS : ºC DESCRIPTION : Temperature DEFAULTVALUE : 0. FILE_IN_TIME : HDF FILENAME : ..\..\GeneralData\Atmosphere\MM5_20070606_Portugal.hdf5 TIME_SERIE : 0 OUTPUT_HDF : 1 <endproperty> <beginproperty> NAME : relative humidity UNITS : fraction DESCRIPTION : Humidity DEFAULTVALUE : 0. FILE_IN_TIME : HDF FILENAME : ..\..\GeneralData\Atmosphere\MM5_20070606_Portugal.hdf5 TIME_SERIE : 0 OUTPUT_HDF : 1 <endproperty> <beginproperty> NAME : cloud cover UNITS : fraction DESCRIPTION : cloud cover from random (default) DEFAULTVALUE : 0. CLOUD_COVER_METHOD : 2 FILE_IN_TIME : NONE TIME_SERIE : 0 OUTPUT_HDF : 0 <endproperty>
Links
- Module_Atmosphere
- Module_InterfaceWaterAir
- ConvertToHDF5
- Meteo - Available meteorological data
