...

Code Block

language	py
theme	Confluence
title	Metview macro for plotting images from simulated satellite data
collapse	true

# Metview Macro - created from Display Window contents

#  **************************** LICENSE START ***********************************
# 
#  Copyright 2015 ECMWF. This software is distributed under the terms
#  of the Apache License version 2.0. In applying this license, ECMWF does not
#  waive the privileges and immunities granted to it by virtue of its status as
#  an Intergovernmental Organization or submit itself to any jurisdiction.
# 
#  ***************************** LICENSE END ************************************
# 
# Date: Thu 17 March 2016
# Author:  ECMWF User Support
# 
  

format  = "png"
satelliteID  = 57  # Meteosat-10
instrument   = 207 # SEVIRI  
channel = 9        # Channel = 5, 6 or 9
clbt    = 260510   # Cloudy Brightness Temperature

outfile = "ssd_channel" & channel

# Device Description 
File_ps = ps_output(  
                OUTPUT_NAME : outfile
        ) 
  
File_png = png_output(  
                OUTPUT_NAME : outfile
        ) 

if format = "png" then
   File = File_png
else
   File = File_ps
end if

# Checks running mode  
mode = runmode()  
if mode = 'batch' or mode = 'execute' then  
   setoutput(File)  
end if  


# Call function to build layout (defined at end of macro)
display_window = build_layout()


# Importing : /usd/EMOSLIB/SSD/get_ssd.mars

ssd_common = (
    stream     : "OPER",
    class      : "od",
    expver     : 0001,
    date       : 0,
    time       : 0,
    type       : "ssd",
    ident      : satelliteID,
    instrument : instrument,
    param      : clbt,
    channel    : channel
    )

ssd_3h = retrieve ( 
    ssd_common,             
    step       : [0,"to",144,"by",3]
    )
    
ssd_6h = retrieve ( 
    ssd_common,             
    step       : [150,"to",240,"by",6]
    )

ssd = ssd_3h & ssd_6h

# Countours for IR channel / clouds

ir_clouds = mcont(
    legend                           : "on",
    contour                          : "off",
    contour_level_selection_type     : "level_list",
    contour_level_list               : [185,187,189,191,193,195,197,199,201,203,205,207,209,211,213,215,217,219,221,223,225,227,229,231,233,235,237,239,241,243,245,247,249,251,253,255,257,259,261,263,265,267,269,271,273,275,277,279,281,283,285,287,289,291,293,295,297,299,301,303,305,307,310,315,320,325],
    contour_label_frequency          : 5,
    contour_shade                    : "on",
    contour_shade_technique          : "cell_shading",
    contour_shade_colour_method      : "list",
    contour_shade_cell_resolution    : 40,
    contour_shade_colour_list        : ["RGB(0.28,0.76,1)","RGB(0.35,0.78,1)","RGB(0.41,0.8,1)","RGB(0.48,0.83,1)","RGB(0.55,0.85,1)","RGB(0.62,0.87,1)","RGB(0.69,0.9,1)","RGB(0.76,0.92,1)","RGB(0.83,0.94,1)","RGB(0.9,0.97,1)","white","RGB(0.99,0.99,0.99)","RGB(0.98,0.98,0.98)","RGB(0.98,0.98,0.98)","RGB(0.97,0.97,0.97)","RGB(0.97,0.97,0.97)","RGB(0.97,0.97,0.97)","RGB(0.97,0.97,0.97)","RGB(0.96,0.96,0.96)","RGB(0.96,0.96,0.96)","RGB(0.96,0.96,0.96)","RGB(0.95,0.95,0.95)","RGB(0.95,0.95,0.95)","RGB(0.95,0.95,0.95)","RGB(0.94,0.94,0.94)","RGB(0.94,0.94,0.94)","RGB(0.93,0.93,0.93)","RGB(0.93,0.93,0.93)","RGB(0.92,0.92,0.92)","RGB(0.91,0.91,0.91)","RGB(0.9,0.9,0.9)","RGB(0.89,0.89,0.89)","RGB(0.87,0.87,0.87)","RGB(0.86,0.86,0.86)","RGB(0.84,0.84,0.84)","RGB(0.82,0.82,0.82)","RGB(0.81,0.81,0.81)","RGB(0.78,0.78,0.78)","RGB(0.76,0.76,0.76)","RGB(0.74,0.74,0.74)","RGB(0.71,0.71,0.71)","RGB(0.69,0.69,0.69)","RGB(0.66,0.66,0.66)","RGB(0.63,0.63,0.63)","RGB(0.6,0.6,0.6)","RGB(0.57,0.57,0.57)","RGB(0.54,0.54,0.54)","RGB(0.51,0.51,0.51)","RGB(0.48,0.48,0.48)","RGB(0.44,0.44,0.44)","RGB(0.41,0.41,0.41)","RGB(0.38,0.38,0.38)","RGB(0.35,0.35,0.35)","RGB(0.32,0.32,0.32)","RGB(0.29,0.29,0.29)","RGB(0.26,0.26,0.26)","RGB(0.23,0.23,0.23)","RGB(0.21,0.21,0.21)"],
    contour_shade_max_level_density  : 25,
    grib_scaling_of_retrieved_fields : "off"
    )

# Contours for Water Vaopour channels

wv_image_common = (
    legend                           : "on",
    contour                          : "off",
    contour_level_selection_type     : "level_list",
    contour_label                    : "off",
    contour_shade                    : "on",
    contour_shade_technique          : "cell_shading",
    contour_shade_colour_method      : "list",
    contour_shade_label_blanking     : "off",
    contour_shade_cell_resolution    : 40,
    contour_shade_colour_list        : ["rgb(1,1,1)","rgb(0.971428571429,0.971428571429,0.971428571429)","rgb(0.942857142857,0.942857142857,0.942857142857)","rgb(0.914285714286,0.914285714286,0.914285714286)","rgb(0.885714285714,0.885714285714,0.885714285714)","rgb(0.857142857143,0.857142857143,0.857142857143)","rgb(0.828571428571,0.828571428571,0.828571428571)","rgb(0.8,0.8,0.8)","rgb(0.771428571429,0.771428571429,0.771428571429)","rgb(0.742857142857,0.742857142857,0.742857142857)","rgb(0.714285714286,0.714285714286,0.714285714286)","rgb(0.685714285714,0.685714285714,0.685714285714)","rgb(0.657142857143,0.657142857143,0.657142857143)","rgb(0.628571428571,0.628571428571,0.628571428571)","rgb(0.6,0.6,0.6)","rgb(0.571428571429,0.571428571429,0.571428571429)","rgb(0.542857142857,0.542857142857,0.542857142857)","rgb(0.514285714286,0.514285714286,0.514285714286)","rgb(0.485714285714,0.485714285714,0.485714285714)","rgb(0.457142857143,0.457142857143,0.457142857143)","rgb(0.428571428571,0.428571428571,0.428571428571)","rgb(0.4,0.4,0.4)","rgb(0.371428571429,0.371428571429,0.371428571429)","rgb(0.342857142857,0.342857142857,0.342857142857)","rgb(0.314285714286,0.314285714286,0.314285714286)","rgb(0.285714285714,0.285714285714,0.285714285714)","rgb(0.257142857143,0.257142857143,0.257142857143)","rgb(0.228571428571,0.228571428571,0.228571428571)","rgb(0.2,0.2,0.2)","rgb(0.171428571429,0.171428571429,0.171428571429)","rgb(0.142857142857,0.142857142857,0.142857142857)","rgb(0.114285714286,0.114285714286,0.114285714286)","rgb(0.0857142857143,0.0857142857143,0.0857142857143)","rgb(0.0571428571429,0.0571428571429,0.0571428571429)","rgb(0.0285714285714,0.0285714285714,0.0285714285714)","rgb(0,0,0)"],
    grib_scaling_of_retrieved_fields : "off"
   )

wv_image_shade5 = mcont(
    wv_image_common,
    contour_level_list               : [221,222,223,224,225,226,227,228,229,230,231,232,233,234,235,236,237,238,239,240,241,242,243,244,245,246,247,248,249,250,251,252,253,254,255,270]
   )
   
wv_image_shade6 = mcont(
    wv_image_common,
    contour_level_list               : [220,222,224,226,228,230,232,234,236,238,240,242,244,246,248,250,252,254,256,258,260,262,264,266,268,270,272,274,276,278,280,282,284,286,288,295]
   )

if channel = 5 then
  image_contour = wv_image_shade5
else if channel = 6 then
  image_contour = wv_image_shade6
else if channel = 9 then
  image_contour = ir_clouds
end if
                 
# Plot command  
plot ( display_window[1], ssd, image_contour) 

# End of main program

# Function to build the layout.
function build_layout( )
    
    full_globe = geoview()
    
    _Page_ = plot_page(  
                TOP        : 0, 
                BOTTOM     : 100, 
                LEFT       : 0, 
                RIGHT      : 100, 
                ROWS       : 1, 
                COLUMNS    : 1, 
                PAGE_X_GAP : 0, 
                PAGE_Y_GAP : 0, 
                VIEW       : full_globe 
        )  
      
    # SuperPageDescription 
    _Display_Window_ = plot_superpage(  
                LAYOUT_SIZE        : "A4", 
                LAYOUT_ORIENTATION : "LANDSCAPE", 
                CUSTOM_WIDTH       : 29.7, 
                CUSTOM_HEIGHT      : 21.0, 
                PAGES                    : [ _Page_ ] 
        ) 
    # plot_superpage returns a list of drop identifiers. 
    # Index 1 is for first page, and so on. 
    # end of superpage definition  
      
    return _Display_Window_ 

end build_layout

Image Modified	Image Modified
48-hour global simulated Meteosat-10 SEVIRI imagery at 10.8 μm (channel 9) from the ECMWF model cycle 41r2 at T_CO1279 / O1280 horizontal resolution with 137 vertical levels run on 1 March 2016 at 00 UTC.	48-hour global simulated Meteosat-10 SEVIRI imagery at 6.30 μm (channel 5)from the ECMWF model cycle 41r2 at T_CO1279 / O1280 horizontal resolution with 137 vertical levels run on 1 March 2016 at 00 UTC.

Image Modified	Image Modified
Animation of 10 day global simulated Meteosat-10 SEVIRI imagery at 10.8 μm (channel 9) from the ECMWF model cycle 41r2 at T_CO1279 / O1280 horizontal resolution with 137 vertical levels run on 1 March 2016 at 00 UTC.	Animation of 10 day global simulated Meteosat-10 SEVIRI imagery at 6.30 μm (channel 5)from the ECMWF model cycle 41r2 at T_CO1279 / O1280 horizontal resolution with 137 vertical levels run on 1 March 2016 at 00 UTC.

Simulated satellite data as ECMWF web charts

...

The simulated satellite data will be available in dissemination. Initially, the products will be made available to ECMWF Member and Co-operating States. Availability in the real-time catalogue will need to be discussed and possibly agreed by ECMWF Committees.

Hocking, J., Rayer, P., Rundle, D., Saunders, R., Matricardi, M., Geer, A., Brunel, P. and Vidot J., 2013: RTTOV v11 Users Guide , NWP SAF report, Met Office, 107 pp.RTTOV v11 web site
Lupu C. and A. J. Geer, 2015: Evaluation and operational implementation of the RTTOV-11 in the IFS, ECMWF Tech. Memo., 748, ECMWF, Reading, UK.

...

. For the simulated satellite data based on Meteosat-10 these are as in the table below:

Satellite	Instrument	Satellite ID	Instrument Type	Channel	Central wavelength (μm)	Description
Meteosat-10	SEVIRI (Spinning Enhanced Visible and InfraRed Imager)	57	207	5	6.30	Water Vapour at ~300 hPa
				6	7.36	Water Vapour at ~500 hPa
				9	10.79	Atmospheric Window channel ~clouds, surface

GRIB encoding

The simulated satellite data fields are encoded in GRIB edition 2 with the parameter set to Cloudy brightness temperature:

paramId	shortName	name	units	GRIB edition
260510	clbt	Cloudy brightness temperature	Kelvin	2

The product description uses productDefinitionTemplateNumber=32 (analysis or forecast at a horizontal level or in a horizontal layer at a point in time for simulated (synthetic) satellite data. The relevant GRIB API keys are encoded as:

GRIB key name
satelliteSeries	satelliteNumber ident	instrumentType instrument	channelNumber channel	NB numberOfContributingSpectralBands	scaleFactorOfCentralWaveNumber	scaledValueOfCentralWaveNumber (units: m^-1)
333	57	207	5	1	0	158800
			6	1	0	135800
			9	1	0	92700

Other GRIB API keys describing the product are:

Key	Value	Description
discipline	3	Space products
parameterCategory	1	Quantitative products
parameterNumber	14	Cloudy brightness temperature

GRIB API version at least 1.14.5 needs to be used in order to decode the simulated satellite data products.

...

Prior to the implementation of IFS cycle 41r2, the simulated satellite image data were stored in MARS with TYPE=SIM and REPRES=SV (Space View) for the IDENT, INSTRUMENT and CHANNEL as specified in the table below.

Sector	Name	Operator	Longitude	Instrument Type	MARS availability		MARS Keywords			Nominal central wavelength (μm)	Dimensions (pixels)	Satellite subpoint (pixels)
Sector	Name	Operator	Longitude	Instrument Type	From	To	IDENT	INSTRUMENT	CHANNEL	Nominal central wavelength (μm)	Dimensions (pixels)	Satellite subpoint (pixels)
Indian Ocean	Meteosat-7	EUMETSAT	54.5° E	MVIRI	2009-03-10	2016-03-08	54	205	1	6.34	2500x2500	1250x1250
Indian Ocean	Meteosat-7	EUMETSAT	54.5° E	MVIRI	2009-03-10	2016-03-08	54	205	2	11.52	2500x2500	1250x1250
East Atlantic	Meteosat-8	EUMETSAT	0°	SEVIRI	2005-06-22	2009-03-09	55	207	5	6.30	1856x1856	928x928
									6	7.36
									9	10.79
	Meteosat-9	EUMETSAT	0°	SEVIRI	2009-03-10	2013-01-23	56	207	5	6.30	1856x1856	928x928
									6	7.36
									9	10.79
	Meteosat-10	EUMETSAT	0°	SEVIRI	2013-01-24	2016-03-08	56	207	5	6.30	1856x1856	928x928
									6	7.36
									9	10.79
East Pacific	GOES-11	NOAA	135.0° W	Imager	2009-03-10	2012-01-17	255	615	3	6.74	1250x1250	625x625
	GOES-11	NOAA	135.0° W	Imager	2009-03-10	2012-01-17	255	615	4	10.72	1250x1250	625x625
	GOES-15	NOAA	135.0° W	Imager	2012-01-18	2016-03-08	259	615	3	6.55	1250x1250	625x625
	GOES-15	NOAA	135.0° W	Imager	2012-01-18	2016-03-08	259	615	4	10.68	1250x1250	625x625
West Atlantic	GOES-12	NOAA	75.0° W	Imager	2009-03-10	2011-05-17	256	615	3	6.48	1250x1250	625x625
	GOES-12	NOAA	75.0° W	Imager	2009-03-10	2011-05-17	256	615	4	10.71	1250x1250	625x625
	GOES-13	NOAA	75.0° W	Imager	2011-05-18	2016-03-08	257	615	3	6.54	1250x1250	625x625
	GOES-13	NOAA	75.0° W	Imager	2011-05-18	2016-03-08	257	615	4	10.67	1250x1250	625x625
West Pacific	MTSAT-1R	JMA	140.2° E	JAMI	2009-03-10	2011-05-17	171	294	2	10.8	1250x1250	625x625
	MTSAT-1R	JMA	140.2° E	JAMI	2009-03-10	2011-05-17	171	294	4	6.75	1250x1250	625x625
	MTSAT-2	JMA	145.0° E	Imager	2011-05-18	2016-03-08	172	294	2	10.8	1250x1250	625x625
	MTSAT-2	JMA	145.0° E	Imager	2011-05-18	2016-03-08	172	294	4	6.77	1250x1250	625x625

These products have been encoded using GRIB edition 1.

...

Space shortcuts

Page tree

Versions Compared

Old Version 24

New Version 25

Key

Simulated satellite data as ECMWF web charts

Further reading

GRIB encoding

Space shortcuts

Page tree

Page History

Versions Compared

Old Version 24

New Version 25

Key

Simulated satellite data as ECMWF web charts

Further reading

GRIB encoding