...
The examples below show the variables to change from normal forecast settings.
| Warning |
|---|
If using the Euler method, make sure a smaller suitable choice of timestep compared to the semi-Lagrangian setup is used. Also, ensure the Asselin time-filter is on. |
| Code Block |
|---|
| title | semi-Lagrangian advection namelist |
|---|
| collapse | true |
|---|
|
! Shallow-water model with semi-Lagrangian
&NAMDYN
HDIRVOR=1.e15, ! horizontal diffusion parameters
HDIRDIV=1.e15, ! see code for further description
HDIRSP=1.e15,
LHDIFFM=true,
NDIFFACT=6,
TSTEP=1800.0, ! model timestep (resolution dependent)
LSETTLS=true, ! ensure SL scheme is enabled
LSETTLST=true,
SITR=350., ! reference temperature
REFGEO=78452.0, ! reference geopotential for SW model
/
&NAMCT0
LRFRIC=false, ! turn of Rayleigh friction
LFPOS=false, ! disable fullpos post-processing
N2DINI=1, ! 2D initialisation switch
NCONF=201, ! select shallow-water configuration
LTWOTL=true, ! enable two-time level SL scheme
LSLAG=true, ! enable SL scheme
LSLPHY=false, ! turn off Sl physics option
LSPRT=false, ! turn off 'virtual temperature' as spectral variable
/
&NAEPHY
LEPHYS=false, ! turn off ECMWF Physics package (master switch)
LERADI=false, ! turn off radiation scheme
/
&NAMDYNA
LGRADSP=.false., ! disable de-aliasing the pressure gradient term
/
&NAMDIM
NUNDEFLD=0, ! make sure uninitialised variables are all set to zero
/ |
| Code Block |
|---|
| title | Euler advection namelist |
|---|
| collapse | true |
|---|
|
!
! Shallow-water model with Euler stepping
!
&NAMDYN
TSTEP=300.0, ! reduce timestep for Euler stepping
REPS1=0.01, ! turn on Asselin time-filtering coefficient
REPS2=0.01,
LHDIFFM=.false., ! horizontal diffusion on / off
LSETTLS=.false., ! extrapolations in SL scheme
LSETTLST=.false., ! ditto
/
&NAMCT0
NCONF=201, ! model configuration: 201 = shallow-water (see yomct0.F90)
LSLAG=.false., ! turn off semi-lagrangian scheme
LTWOTL=.false., ! disable two-time-level SL scheme.
LRFRIC=.false., ! turn off Rayleigh friction
LSLPHY=.false., ! turn off split time-step physics
LVERTFE=.false., ! turn off vertical finite element scheme
N2DINI=1, ! initialise 1 = Haurwitz wave, 2 = real fields
LSPRT=.false., ! if T temperature is 'virtual temperature'
LFPOS=false, ! turn off fullpos diagnostics, does not work with SW
N3DINI=0, ! no 3D initialisation
NSTOP=600, ! no. of steps to run
/
&NAEPHY
LEPHYS=false, ! turn off ECMWF Physics package (master switch)
LERADI=false, ! turn off radiation scheme
/
&NAMDYNA
LGRADSP=.false., ! disable de-aliasing the pressure gradient term
/
&NAMDIM
NUNDEFLD=0, ! make sure uninitialised variables are all set to zero
/ |
...