WET DRY

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:
  D(i,j)=h(i,j)+zeta(i,j,1)
  If the total depth is less than a user specified minimum depth (Dcrit), then the water level is increased to
  zeta(i,j,1) = Dcrit - h(i,j)
  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