;---------------------------------------------------------------------- ; Read and plot a 2D variable read off the fnl_20070703_06_00.grib1 ; file. ; ; Data is rectilinear, 181 x 360 (lat x lon) ;---------------------------------------------------------------------- ; Read data dir = "../Data/" filename = "fnl_20070703_06_00.grib1" a = addfile(dir+filename,"r") print(getfilevarnames(a)) ; Print all variables on the file var_name = "TMP_3_HTGL" ; pick a variable to plot var = a->$var_name$ ;---Print some information about the variable. printVarSummary(var) printMinMax(var,0) wks = gsn_open_wks("png","fnl_plot_" + var_name + "_ncl") ;---Set some plot options res = True res@cnFillOn = True ; turn on contour fill res@cnLinesOn = False ; turn off contour lines res@cnLineLabelsOn = False ; turn off line labels ;res@cnFillPalette = "MPL_jet" ; chose a color map res@tiMainString = filename + " : " + var_name plot = gsn_csm_contour_map(wks,var,res)