WET DRY

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 depth of water computed on that cell. In the present formulation for wetting/drying: At initialiation the model computes for each cell the total depth of water 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 = Dcrit - h(i,j). This 'drapes' the water level over all the land cells. 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.