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