Improve 3rd order dynamic advection scheme
Context
Upstream biased flux form advection scheme for momentum implementation has :
- a 4th order divergence operator implemented along horizontal directions
- a 2nd order operator implemented along vertical direction
-
Reference configuration/test case with ln_dynadv_ubs -
Modifications of versioned files: Fortran routines (dynadv.F90), namelists ( namelist_ref/cfg
)
Proposal
Upstream biased flux form can be changed into a proper 3rd order upstream scheme by :
- setting
gamma2 = 0._wp
or removing it from the code It turns the 4th order divergence operator implemented along horizontal directions into a 2nd order divergence operator. This 2nd order divergence operator is now compatible with the divergence from the continuity equation. - implementing a 3rd order operator along vertical direction