Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

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

Post-processing

Retrieve the 1000,500 and 300 hPa temperatures for tomorrow at noon from the 12 UTC forecast of two days ago on an N48 Gaussian Grid  temperature on model levels 1, 2 and 3 for yesterday's 12 hour forecast. > see in Catalogue

Expand
titlesolution
retrieve, 
    type     = forecastfc, 
    paramdate    = t-1, 
    date time    = -212,
    step     = 4812, 
    gridlevtype = ml,
   = 48, 
   level  levelist = 10001/5002/3003,
    param levtype  = plt,
    target   = "t.gg.grb"

...

_ml_1_2_3"

Time and step

Retrieve the 1000 hPa geopotential from yesterday's 12 UTC forecast (> see in Catalogue) for:

  • the 12 hour forecast

    Expand
    titlesolution
    retrieve,

...

  • 
        type    

...

  • = 

...

  • fc,
        

...

  • date    = 

...

  • -1,
        

...

  • time    = 

...

  • 12,
        step   

...

  •  = 

...

  • 12,
        levtype

...

  •  = pl,
        

...

  • level   = 

...

  • 1000,
        

...

  • param   = 

...

  • z,
        target  

...

  • = "

...

  • z1000_a"
    
  • the 6, 12, 24 and 48 hour forecasts

    Expand
    titlesolution
    retrieve,
        type

...

  •     = fc,
        

...

  • date    = 

...

  • -1,
        

...

  • time    = 

...

  • 12,
        

...

  • step    = 

...

  • 6/

...

  • 12/24/

...

  • 48,
        levtype

...

  •  = pl,
        

...

  • level   = 

...

  • 1000,
        

...

  • param   

...

  • = 

...

  • z,
        

...

  • target  

...

  • = "z1000_b"
  • the 12 hour to 240 hour forecasts at 12 hour intervals

    Expand
    titlesolution
    retrieve,
        type 

...

  •    = 

...

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

  • leaving the data in its original form

    Expand
    titlesolution
    retrieve,fc,
        typedate    = an-1,
        datetime    = -112,
        timestep    = 12/to/240/by/12,
        param  levtype = zpl,
     
       level  levtype = pressure level1000,
        levelparam   = 1000/500z,
        target  = "myfile1z1000_c"
    
  • all forecast time-stepsreducing the truncation to T106

    Expand
    titlesolution
    retrieve,
        type    = anfc,
        date    = -1,
        time    = 12,
    #   step param   = zall, 
    # The value "all" levtypeis =convenient to pressureuse level,
    but MARS can only levelreport on missing = 1000/500,
     steps if 
    resol = 106, target = "myfile1"

    converting the data to a global N48 Gaussian grid

    Expand
    titlesolution
    retrieve,
        type # the full =list an,
    of steps is specified:
    date = -1, step param = z, levtype = pressure level, level = 1000/500, grid = 48, target = "myfile1"= 0/1/2/3/4/5/6/7/8/9/10/11/12/13/14/15/16/17/18/19/20/21/22/23/24/25/26/27/28/29/30/31/32/33/34/35/36/37/38/39/40/41/42/43/44/45/46/47/48/49/50/51/52/53/54/55/56/57/58/59/60/61/62/63/64/65/66/67/68/69/70/71/72/73/74/75/76/77/78/79/80/81/82/83/84/85/86/87/88/89/90/93/96/99/102/105/108/111/114/117/120/123/126/129/132/135/138/141/144/150/156/162/168/174/180/186/192/198/204/210/216/222/228/234/240,
    levtype = pl, level = 1000, param = z, target = "z1000_d"

Post-processing

Retrieve the 1000,500 and 300 hPa temperatures for tomorrow at noon from the 12 UTC forecast of two days ago on an N48 Gaussian Grid  > see in Catalogue

Expand
titlesolution
retrieve, 
    type     = forecast, 
    param    = t, 
    date     = -2,
    step     = 48, 
    grid     = 48, 
    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 Catalogue

Expand
titlesolution
retrieve, 
    type     = forecast,
    param    = t,
    date     = -2,
    levelist = 1000/500/300,
    levtype  = pl,
    step     = 48,
    grid     = 1.5/1.5,
    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 Catalogue

Expand
titlesolution
retrieve,
    type     = fc,
    param    = t,
    date     = -2,
    levelist = 1000/500/300,
    levtype  = pl,
    step     = 48,
    grid     = 1.5/1.5,
    area     = e,
    target   = "t.europe"

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

  • leaving the data in its original form

    Expand
    titlesolution
    retrieve,
        type    = an,
        date    = -1,
        time    = 12,
        param   = z, 
        levtype = pressure level,
        level   = 1000/500,
        target  = "myfile1"
  • reducing the truncation to T106

    Expand
    titlesolution
    retrieve,
        type    = an,
        date    = -1,
        time    = 12,
        param   = z, 
        levtype = pressure level,
        level   = 1000/500,
        resol   = 106,
        target  = "myfile1"
  • converting the data to a global N48 Gaussian grid

    Expand
    titlesolution
    retrieve,
        type    = an,
        date    = -1,
        param   = z, 
        levtype = pressure level,
        level   = 1000/500,
        grid    = 48,
        target  = "myfile1"
  • converting the data to a regular 5.0 by 5.0 degree global latitude/longitude grid

    Expand
    titlesolution
    retrieve

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

    Expand
    titlesolution
    retrieve,
        type    = an,
        date    = -1,
        param   = z, 
        levtype = pressure level,
        level   = 1000/500,
        grid    = 5.0/5.0,
        target  = "myfile1"
  • converting the data to an N80 Gaussian grid for the European area

    Expand
    titlesolution
    retrieve,
        type    = an,
        date    = -1,
        param   = z, 
        levtype = pressure level,
        level   = 1000/500,
        grid    = 80,
        area    = europe,
        target  = "myfile1"
  • 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
    titlesolution
    retrieve,
        type    = an,
        date    = -1,
        param   = z, 
        levtype = pressure level,
        level   = 1000/500,
        grid    = 80,
        area    = 50/-20/10/20,
        target  = "myfile1"
  • converting the data to the same area but convert it to a 2.5 degree by 2.5 degree latitude/longitude grid

    Expand
    titlesolution
    retrieve,
        type    = an,
        date    = -1,
        param   = z, 
        levtype = pressure level,
        level   = 1000/500,
        grid    = 2.5/2.5,
        area    = 50/-20/10/20,
        target  = "myfile1"
    

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

Expand
titlesolution
retrieve,
    type    = fc,
    date    = -1,
    time    = 12,
    step    = 24,
    levtype = sfc,
    param   = 10u/10v,
    area    = europe,       # As the area not compatible with the grid,
    grid    = 0.57/0.57,    # MARS will expand the area to make = 0.57/0.57,it fit.
    target  = "10u_10v"

 

Observations

...

Expand
titlesolution
retrieve,
    class   = e4,
    stream  = moda,
    type    = an,
    levtype = sfc,
    date    = 19640101,
    param   = 69.162/70.162,
    area    = europe    = an,
    levtype = sfc,
    griddate    = 1.125/1.12519640101,
    target  = "e4_moda_fluxes.grb"

 

 

...

  • the 12 hour forecast (solution)
  • the 6, 12, 24 and 48 hour forecasts (solution)
  • the 12 hour to 240 hour forecasts at 12 hour intervals (solution)
  • all forecast time-steps (solution)

...

param   = 69.162/70.162,
    area    = europe,
    grid    = 1.125/1.125,
    target  = "e4_moda_fluxes.grb"

 

 

Retrieve the temperature from model levels 1 and 2 for 1800 on 18th October 1991 from the uninitialised analysis. Try some or all of the above conversions. (solution)

Retrieve the 2 metre temperature and dewpoint for yesterday's four analysis times. Try some or all of the above conversions. (Conversion to spherical harmonics is not possible.)

...

(solution)