Any requested data manipulation or post-processing is carried out by the MARS client, except in the case of a local Member State's client where data is first processed at ECMWF prior to its transmission over the network. The post-processing is carried out by a set of routines present in the EMOSLIB library. Please, refer to the Field interpolation routines for in-depth details about such post-processing.

grid

specifies the output grid which can be either a Gaussian grid or a Latitude/Longitude grid. MARS requests specifying grid=av will return the archived model grid.

For Latitude/Longitude output, the format is in degrees, i.e. grid = 2.5/2.5, where the first number denotes the east-west resolution (longitude) and the second denotes the north-south resolution (latitude). The grid spacing needs to be an integer fraction of 90 degrees (latitude range from the Equator to the Pole), e.g. 0.225, 0.25, 0.28125, 0.3, 0.4, ...

Output on a Gaussian grid is specified by a letter denoting the type of Gaussian grid followed by an integer (the grid number) representing the number of lines between the Pole and Equator, e.g. 

Not all grid names and grid number combinations are supported. Currently, the Gaussian grid names and numbers supported for MARS retrievals are:

Octahedral reduced Gaussian (O grids)O1280O640O512O400O320O256O200O160O128O96O64O48O32
Original reduced Gaussian  (N grids) N640N512N400N320N256N200N160N128N96N64N48N32
Regular (Full) Gaussian (F grids)F1280F640F512F400F320F256F200F160F128F96F64F48F32

MARS requests specifying an unsupported grid will fail. For example, a retrieval request with grid=O123 will return an error.  If the letter denoting the type of Gaussian grid is omitted, e.g. grid=320, a full (or regular) Gaussian grid with 320 grid lines is returned.

For legacy requests see also the section on keyword gaussian below.

gaussian

This keyword is deprecated and should no longer be used. 

specifies the desired type of Gaussian grid for the output.

Valid Gaussian grids are quasi-regular (reduced) or regular. Keyword gaussian can only be specified together with keyword grid. Gaussian without grid has no effect.

When gaussian = reduced is selected, only certain values are valid for keyword grid, according to the definitions for Reduced Gaussian Grids used at ECMWF.

area

specifies the desired sub-area of data to be extracted.

Some predefined areas exist (Europe, and four quadrants for each hemisphere). Alternatively, you may specify an area of your own, by inputting the area limits as North/West/South/East; if South > North, the values are swapped and a warning issued; southern latitudes and western longitudes must be given as negative numbers. For example, Europe may be defined by area=75/-20/10/60.

Areas can be defined to wrap around the globe, e.g. in order to retrieve all longitudes but Europe one would specify area=75/60/10/-20. Sub-area extraction only works on regular latitude/longitude fields, regular Gaussian fields and wave fields. For observations, it is used as a filter. 

NORTH = 90° WEST = -180°SOUTH = -90°EAST = 179.5°

 

NORTH = 75° WEST = -25°SOUTH = 10°EAST = 60°

 

NORTH = 75° WEST = -15°SOUTH = 30°EAST = 42.5°

NORTH = 70° WEST = -20°SOUTH = 35°EAST = 20°

NORTH = 75° WEST = 10°SOUTH = 30°EAST = 40°

NORTH = 40° WEST = -20°SOUTH = -40°EAST = 60°

NORTH = 50° WEST = -135°SOUTH = -60°EAST = -35°

NORTH = 40° WEST = -110°SOUTH = -10°EAST = -60°

NORTH = 80° WEST = 115°SOUTH = -5°EAST = 180°

 

frame

specifies the number of points to be selected from a sub-area inwards. It works together with the area keyword. For instance, area=europe, grid=2/2, frame=2 will select a frame with outer limits 74/60/10/-20 and inner limits 72/58/8/-18.

resol

specifies the desired triangular truncation of retrieved data, before carrying out any other selected post-processing.

The default is automatic truncation (auto), by which the lowest resolution compatible with the value specified in grid is automatically selected for the retrieval (see Truncation before interpolation).

Users wanting to perform post-processing from full spectral resolution should specify Archived Value (av). The following are examples of existing resolutions found in the archive: 63, 106, 159, 213, 255, 319, 399, 511, 799 or 1279. This keyword has no meaning/effect if the archived data is not in spherical harmonics representation.

For applications requiring VarEPS products starting from or spanning step T+240, users wishing to obtain products interpolated to regular lat/lon grids are advised to have the products from Leg 1 truncated to the N128 grid prior to interpolation to their choice of lat/lon grid. This can be achieved by specifying resol=N128 (or reduced gaussian 128) which will reduce the resolution to N128 prior to further interpolation to regular grids.

rotation

specifies a rotation for the output fields. The values for this keyword consist of the latitude and longitude of the southern pole in degrees, e.g. rotation = -30.0/10.0.

Fields on spherical harmonics or Gaussian grid can be rotated with the directive rotation. The rotation is performed prior to any other conversion. Therefore:

If the input is a spectral field and the output is a rotated spectral field, a file of rotation coefficients is generated and used in the processing. The convergence of the algorithm for rotating the spectral fields is sensitive to the spectral truncation. For some high resolutions and some large angles of rotation, it is necessary to split the rotation into three steps (two forwards and one backwards!). The size of the rotation coefficient files can be very large.

If the input is a spectral field and the output is a rotated grid point field, the spectral field is interpolated to a non-rotated Gaussian grid which is then transformed to the required rotated grid point field using a 12-point interpolation based on the FULL-POS scheme used in the ECMWF forecast model.

A rotated grid-point field is created from an input grid-point field by finding for each rotated grid-point its nearest four neighbours in the input field and carrying out a bilinear interpolation.

interpolation

can be used to force the interpolation method to bilinear or nearest neighbour for interpolations from grids.

accuracy

specifies the number of bits per value to be used in the generated GRIB coded fields.

A positive integer may be given to specify the preferred number of bits per packed value. This must not be greater than the number of bits normally used for a Fortran integer on the processor handling the request (typically 32 or 64 bit).

Within a compute request the accuracy of the original fields can be passed to the result field by specifying accuracy=av.

packing

specifies the packing method of the output fields. Valid values for latitude/longitude output are simple (si) or second order (so). Valid values for spectral fields are complex packing (co) or simple packing (si). In order to preserve the packing of the archived data one would specify archived value (packing = av).

specification

forces an old grid definition to be used (Research Department designated _12). Can only be used with N80 or N160 reduced Gaussian grids, e.g. specification = 12 will force a grid definition which has 12 points in the northernmost and southernmost latitudes.

style

specifies the style of post-processing.

A valid directive is style = dissemination, in order to instruct MARS to perform post-processing in a way as close to the dissemination system as possible. For example:

bitmap

specifies a UNIX filename containing a bitmap definition in the format defined in appendix 1 of the  Dissemination manual.