Page History

tThe The following MARS requests can be used as a tutorial. The links into the MARS Catalogue might be helpful to create a valid retrieve script. However, the dates might need to be adapted.

...

retrieve,
    class    = od,
    stream   = mnth,
    type     = analysis,
    time     = 06,
    date     = 19970200,
    levtype  = model level,
    param    = 130.128,
    levelist = 1/to/31137, 
    target   = "monthly"

Retrieve the temperature on model levels 1, 2 and 3 for a week old 12 hour forecast. > see in MARS Catalogue

...

Retrieve Total Precipitation, perturbed forecast overlap for 26th June 2023 (> see in MARS Catalogue)

retrieve,

    class = od,

    date = 2023-06-26,

    expver = 1,

    levtype = sfc,

    number = 1/to/50,

    param = 228.128,

    step = 336,

    stream = efov,

    time = 00:00:00,

    type = pf,

    target = [stream]_[date].grb

See also FAQ: "How can I compute the accumulation using data from the overlap stream?"

...

Retrieve the 1000,500 and 300 hPa temperatures valid for today day before yesterday at noon from the 12 UTC forecast of two four days ago on an N48 Gaussian Grid  > see in MARS Catalogue

retrieve, 

    type class    = forecastod,
 
   stream param    = 130.128oper,
 
   expver date     = -21,
    timetype     = 12forecast, 
    stepparam     = 48130.128, 
    griddate     = -4,
    time     = 12, 
    step     = 48, 
    grid     = N48, 
    levelist = 1000/500/300, 
    levtype  = pl, 
    target   = "t.gg.grb"

Retrieve the same on a latitude/longitude grid, 1.5 by 1.5 degrees. > see in MARS Catalogue

retrieve, 
  
      class = od,
  type    stream = forecastoper,
 
   param   expver = 130.1281,
 
   date    type = -2forecast,

    time  param = 130.128,
 = 12,
    levelistdate = 1000/500/300,
-4,
     levtype time = pl12,

    step   step  = 48,

    grid    grid = 1.5/1.5,

    target  levelist = "t.ll.grb"1000/500/300,
      levtype = pl,
      target   = "t.ll.grb"

Retrieve Europe in the same way as in the previous retrieval. Compare the size of the resulting file with the previous one. > see in MARS Catalogue

retrieve,

    class = typeod,
    stream = fcoper,

    paramexpver    = 130.1281,
 
   date  type   = -2forecast,

    timeparam     = 12130.128,

    levelistdate = 1000/500/300,
-4,
    levtypetime  = pl12,

    step     = 48,

    grid     = 1.5/1.5,
   ,
    levelist = 1000/500/300,
    levtype = pl, 
    area     = e,
    target   = "t.europe"

Retrieve the 1000 and 500 hPa geopotential for 1200 yesterday from the analysis UTC analysis from two days ago in the formats below (one job for each form): > see in MARS Catalogue

• leaving the data in its original form

  Expand
  title solution
  retrieve, type class = anod, date = -2, expver = -1, timelevelist = 500/1000, levtype = pressure 12level, param = 129.128, levtypestream = pressure leveloper, leveltime = 12, type= 1000/500, an,   target = "myfile1myfile_d"

  
  

• reducing the truncation to T106

  Expand
  title solution
  retrieve, type class = anod, date = -12, timeexpver = 121, param levelist = 129.128, 500/1000, levtype = pressure level, levelparam = 1000/500, 129.128, resol stream = 106oper, target time = "myfile112, type=an,    truncation = 106, target = "myfile_106"

  
  

• converting the data to a global N48 Gaussian grid

  Expand
  title solution
  retrieve, typeclass = anod, date = -2, expver = -1, param levelist = 129.128, 500/1000, levtype = pressure level, level = 1000/500, grid param = 129.128, stream = oper, time = 12, type=an, grid = N48, target = "myfile1myfile_N48"

  
  

• converting the data to a regular 5.0 by 5.0 degree global latitude/longitude grid

  Expand
  title solution
  retrieve, type class = anod, date = -12, param expver = 129.128, 1, levelist = 500/1000, levtype = pressure level, level = 1000/500, param = 129.128, stream = oper, time = 12, gridtype=an, grid = 0.5/0.0/5.0, target = "myfile1myfile_ll"

  
  

• converting the data to an N80 Gaussian grid for the European area

  Expand
  title solution
  retrieve, type class = an, od, date = -2, expver = -1, param levelist = 500/1000, levtype = pressure level, param = 129.128, stream = oper, time levtype= 12, type=an, grid = pressure level, level = 1000/500, grid = N80, area = europe, target N80, target = "myfile1myfile_N80"

  
  

• converting the data to an N80 Gaussian grid for the area from 50 degrees North to 10 degrees North, and 20 degrees West to 20 degrees East

  Expand
  title solution
  retrieve, type class = anod, date = -1, param = 129.128, 2, expver = 1, levelist = 500/1000, levtype = pressure level, param = 129.128, level stream = 1000/500, oper, time grid= 12, type = N80an, grid = areaN80, area = 50/-20/10/20, target = "myfile1myfile_N80_reg"

  
  

• converting the data to the same area but convert it to a 2.5 degree by 2.5 degree latitude/longitude grid

  Expand
  title solution
  retrieve, typeclass = od, = an, date = -2, date expver = -1, param levelist = 129.128, 500/1000, levtype = pressure level, levelparam = 1000/500, 129.128, gridstream = 2.5/2.5, oper, areatime = 12, type = an, grid = 0.25/0.25, area = 50/-20/10/20, target = "myfile1myfile_ll_reg"

  
  

Retrieve surface winds for yesterdayprevious week's 24 hour forecast on the European area converting to a 0.57 by 0.57 degrees latitude/longitude grid. > see in MARS Catalogue
What happens to the conversion?

retrieve,
    type    = fcclass=od,
    date=-7,
    expver= -1,
    timelevtype=sfc,
    = 12param=165.128/166.128,
    step=24,
    stream= 24oper,
    levtype = sfctime=12:00:00,
    param   = 10u/10v,
    area    type=fc,    
    area = europe,         # As the area is not compatible with the grid,
 
    grid    = 0.57/0.57,      # MARS will expand the area to make it fit.

    target  = "10u_10v"

Observations

Retrieve all the sea surface observations for the 1st February 2015.  > see in MARS Catalogue

retrieve, 
    type    = ob, 
    obsgroup= con,
    obstype = ssd,
#   obstype = 9/11/12/13/14/19/21/22/23/180/181/182,

    date    = 20150201, 
    time    = 00,
    range   = 1439,
    target  = "obs.bufr"

...

Retrieve the 2 metre temperature field of the first 10 members of yesterday's a week old 12 UTC ensemble prediction (type: perturbed forecast) for the next 2, 3 days at noon. and 4 days in advance > see in MARS Catalogue

retrieve,

    class    = od,

    stream   = ef,
date=-7,
    type     = perturbed forecast,
expver=1,
    levtype  = sfc,

    represnumber=1/to/10,
    param= gg,
167.128,
    date     = -1,
    time     = 12,
    ensemble = 1/to/10,
    step     = step=48/to/96/by/24,
    param    = 167.128, 
    target = "ensemble"

Retrieve the same for the control forecast. > see in MARS Catalogue

retrieve,
    class   = od,
    stream  = ef,
    typeenfo,
    time= control forecast,
12,
    levtype type= sfcpf,

    represtarget  = gg,
"ensemble"

Retrieve the same for the control forecast. > see in MARS Catalogue

retrieve,
    dateclass=od,
    date= -17,

    timeexpver=1,
    levtype= 12sfc,

    stepparam=167.128,
    step= 48/to/96/by/24,

    paramstream=enfo,
    time= 167.128,
12,
    type=cf,
    target  = "ensemble.cf"

ERA 5

Retrieve the vertical integral of eastward and northward heat flux as monthly means of daily means on a 1.125 by 1.125 degree latitude/longitude grid over Europe from the ERA 5 archive for January 1964 2021 > see in MARS Catalogue

retrieve,
    class   = ea,
    stream  = oper,
    type    = an,
    levtypedate = sfc20210101,
    datedecade    = 1964-01-01,
2020,
    expver  = 1,

    timelevtype     = 00/to/23sfc,
    param   = 69.162/70.162,
    stream = moda,
 area   type = european,

    gridarea = europe,
    grid = 1.125/1.125,
 
    target  = "ea_moda_fluxes.grb"

Retrieval of snow depth from the ERA 5 archive for November 1993, for all analysis base times > see in MARS Catalogue. How many fields are requested?

retrieve,
  class   = ea,
  stream  = oper,
  expver  = 1,
  date    = 19931101/to/19931130,
  time    = 00/to/23,
  type    = an,
  levtype = sfc, 
  param   = 141.128,
  repres  = gg,
  target  = "era5.199311.sd"

#                                  It retrieves 30 x 24 = 720 fields.

...