New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
2016WP/NocWad (diff) – NEMO

Changes between Version 1 and Version 2 of 2016WP/NocWad


Ignore:
Timestamp:
2016-04-07T12:41:07+02:00 (8 years ago)
Author:
acc
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • 2016WP/NocWad

    v1 v2  
    2626 
    2727== Abstract == 
     28 
     29=== A Wetting and Drying (W/D) Method for NEMO-shelf Based on Flux Limiter === 
     30 
     31==== Hedong Liu and Jason Holt (NOC, UK) ==== 
     32 
     33NEMO discretize its computational domain with the structured C-grid, hence its W/D 
     34algorithm needs to deal with two types of grid cells: velocity grid cells and tracer grid 
     35cells. The essences of current method are to keep the positivity of watermass and avoid 
     36unphysical overshooting of velocity.At the same time, this method must respect the 
     37properties of local and global mass conservation. There are quite many existing methods 
     38for W/D,(see Balzano,A, 1998 and Medeiros,S.C. and Hagen,S.C, 2013). But there are still some 
     39robustness issues, such as model stability, computational efficiency.  In the current 
     40method, a flux limiter has been developed to take care of the mass positivity while a new 
     41barotropic pressure gradient filter is developed to prevent the velocity overshooting at 
     42the wetting and drying cells on a steep slope. 
     43 
     44In the following parts of this note, we will first derive the flux limiter and demonstrate 
     45its applications to the continuity equation, tracer equations. As the momentum equations 
     46is also linked with the W/D processes, we then detail the design of the gravity forcing 
     47filter which aims to remove the errors generated by the normal pressure gradient 
     48calculations in the W/D area. The implementation of the method in NEMO v3.6_STABLE source 
     49code will be illustrated at the last. 
     50 
     51==== Flux limiter ==== 
     52 
     53In NEMO, the water column can be considered as a control volume V. The horizontal 
     54advection is the main process causing the wetting and drying on the grid cells (with the 
     55evaporation and precipitation as the other two reasons which will be dealt separately). 
     56For such a control volume V, we can define a flux limiter as follows: 
     57 
     58[[Image(math figs/fig01.png)]] 
    2859 
    2960This section should be completed before starting to develop the code, in order to find agreement on the method beforehand.