Page History

...

If you have access to the MARS archive at ECMWF then an example of how to retrieve single level fields is:

retrieve,
  class="ei",
  type=an,
  stream=oper,
  date=20100212,
  time=12,
  step=0,
  expver=1,
  param=VOR/D/Z,
  levelist=1,
  levtype=ml,
  target="sh_ml"

 retrieve,
  class="ei",
  type=an,
  stream=oper,
  date=20091227,
  time=12,
  step=0,
  expver=1,
  param=139,
  levelist=off,
  levtype=sfc,
  target="sfc_gg"

Namelists

It is extremely important to correctly set the NAMELIST to configure the shallow-water model correctly. Two examples are given here, one for semi-Lagrangian advection, the other for Euler advection. Euler advection must use a regular gaussian grid. It will not work with a reduced grid.

code

Gotchas

...

Info
title GRIB_GRIBEX_MODE_ON
The model will fail unless the environment variable: export GRIB_GRIBEX_MODE_ON=1

...

is set in the job script. This is because the shallow-water code does not go through the usual code path (if NCONF=1) and uses a older code to output. Otherwise the model will crash at the first output

...

attempt. Users will usually see a failure

...

from the grib_api library (

...

typically grib_set_real8_array).

Info
title Make sure FULLPOS is off
Related to the above, the shallow-water model does not use the FULLPOS part of the model for output (namelist NAMFPC) normally used when NCONF=1 (3D primitive equation model). Make sure FULLPOS is off by setting: &NAMCT0    LFPOS=false, in the fort.4 namelist file.

 

 

• make sure &NAMDIM, NUNDEFLD=0 is set. This ensure all uninitialized variables are set to zero. The SW model can output fields not used.
• FPO FPOS must be off.

Plotting

data is spectral

...