;NCL script to read and plot Landsaf data load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" begin DIDSSF_file = addfile("../HDF5_LSASAF_MSG_DIDSSF_Euro_201212050000.nc","r") LAT_file = addfile("../HDF5_LSASAF_MSG_LAT_Euro_4bytesPrecision.h5","r") LON_file = addfile("../HDF5_LSASAF_MSG_LON_Euro_4bytesPrecision.h5","r") lat = LAT_file->LAT lon = LON_file->LON dssf = DIDSSF_file->DSSF lat1d = lat(:,0)*0.01 lon1d = lon(0,:)*0.01 lat1d@units = "degrees_north" lon1d@units = "degrees_east" dssf!0 = "lat" dssf!1 = "lon" dssf&lat = lat1d dssf&lon = lon1d wks = gsn_open_wks("x11", "lsasaf") gsn_define_colormap(wks,"BlGrYeOrReVi200") res = True res@cnFillOn = True ; res@tfDoNDCOverlay = False ; transform to standard lat/lon res@sfXArray = lon res@sfYArray = lat plot = gsn_csm_contour_map(wks,dssf,res) end