Difference between revisions of "WET DRY"

From WikiROMS
Jump to navigationJump to search
(New page: Wetting/drying. The grid file should use land sea masking. Cells with a rho mask = 0 will always be dry. Cells with a rho mask = 1 can be wet (active) or dry (inactive) depending on the ...)   (change visibility)
 
 
(5 intermediate revisions by 2 users not shown)
Line 1: Line 1:
Wetting/drying.
<div class="title">WET_DRY</div>
The grid file should use land sea masking.  
<!-- Import short description from Options page -->
Cells with a rho mask = 0 will always be dry.  
{{#lst:Options|WET_DRY}}
Cells with a rho mask = 1 can be wet (active) or dry (inactive) depending on the depth of water computed on that cell.  
<br />
In the present formulation for wetting/drying:
----
At initialiation the model computes for each cell the total depth of water  
<br />
D(i,j)=h(i,j)+zeta(i,j,1)  
* The Land/Sea [[MASKING|masking]] arrays ([[Variables#mask_rho|mask_rho]], [[Variables#mask_u|mask_u]], [[Variables#mask_v|mask_v]]) are required in the grid NetCDF file. Recall that the Land/Sea masking is set in some of the [[Test Cases|test cases]] with an analytical routine by activating [[ANA_MASKING]].
If the total depth is less than a user specified minimum depth (Dcrit), then the water level is increased to = Dcrit - h(i,j).
** Cells with a [[Variables#mask_rho|mask_rho]] '''= 0''' will always be '''dry'''.  
This 'drapes' the water level over all the land cells.  
** Cells with a [[Variables#mask_rho|mask_rho]] '''= 1''' can be '''wet''' (active) or '''dry''' (inactive) depending on the depth of water computed on that cell.
During computation the water level in each cell is compared to Dcrit. If the water level is less than Dcrit than no flux is allowed out of that cell. Water can always flow into a cell.
: Cells that are permanently dry (''i.e.'' cells with a [[Variables#mask_rho|mask_rho]] '''= 0''') have a water depth equal to [[Variables#Dcrit|Dcrit]].
* In the model, the wet/dry masks are stored in the [[Variables#rmask_wet|rmask_wet]], [[Variables#umask_wet|umask_wet]], and [[Variables#vmask_wet|vmask_wet]] arrays.
* In the present formulation for wetting/drying, at initialization, the model computes for each cell the total depth of water as:<div class="box">D(i,j)=[[Variables#h|h(i,j)]]+[[Variables#zeta|zeta(i,j,1)]]</div> If the total depth is less than a user specified minimum depth ([[Variables#Dcrit|Dcrit]]), then the water level is increased to<div class="box">[[Variables#zeta|zeta(i,j,1)]] = [[Variables#Dcrit|Dcrit]] - [[Variables#h|h(i,j)]]</div> This ''drapes'' the water level over all the land cells. The default value for [[Variables#Dcrit|Dcrit]] is 0.10 m.
* During computation, the water level in each cell is compared to [[Variables#Dcrit|Dcrit]]. If the water level is less than [[Variables#Dcrit|Dcrit]], then no flux is allowed out of that cell. Water can always flow into a cell.
* For output, a variable [[Variables#rmask_wet|wet_dry_masking]] provides the time history of the wet/dry masking at &rho;-points.
* The bathymetry [[Variables#h|h]] is measured positive downwards, with a typical refernce to mean sea level (MSL). Therefore value of [[Variables#h|h]] below MSL are positive, and values of [[Variables#h|h]] above MSL are negative. ROMS internally deals with the positive and negative values of [[Variables#h|h]].
* To activate this feature you need to use:
:<div class="box">  #define WET_DRY</div>

Latest revision as of 17:56, 31 March 2008

WET_DRY
WET_DRY
Option to activate wetting and drying.
parameter = Dcrit
variable = rmask_wet, umask_wet, vmask_wet
routine = step2d.F




  • The Land/Sea masking arrays (mask_rho, mask_u, mask_v) are required in the grid NetCDF file. Recall that the Land/Sea masking is set in some of the test cases with an analytical routine by activating ANA_MASKING.
    • Cells with a mask_rho = 0 will always be dry.
    • Cells with a mask_rho = 1 can be wet (active) or dry (inactive) depending on the depth of water computed on that cell.
Cells that are permanently dry (i.e. cells with a mask_rho = 0) have a water depth equal to Dcrit.
  • In the model, the wet/dry masks are stored in the rmask_wet, umask_wet, and vmask_wet arrays.
  • In the present formulation for wetting/drying, at initialization, the model computes for each cell the total depth of water as: If the total depth is less than a user specified minimum depth (Dcrit), then the water level is increased to This drapes the water level over all the land cells. The default value for Dcrit is 0.10 m.
  • During computation, the water level in each cell is compared to Dcrit. If the water level is less than Dcrit, then no flux is allowed out of that cell. Water can always flow into a cell.
  • For output, a variable wet_dry_masking provides the time history of the wet/dry masking at ρ-points.
  • The bathymetry h is measured positive downwards, with a typical refernce to mean sea level (MSL). Therefore value of h below MSL are positive, and values of h above MSL are negative. ROMS internally deals with the positive and negative values of h.
  • To activate this feature you need to use:
#define WET_DRY