;---------------------------------------------------------------------- ; This script demonstrates the timing difference between using the ; default "AreaFill" for color contour plots, versus "RasterFill". ; ; This graphic is used in the Graphics workshop lecture slides. ;---------------------------------------------------------------------- begin start_time = get_cpu_time() g = addfile("arpege_grd.nc","r") f = addfile("IBF_1m_000101_002012_SUTOPRSU.nc","r") data = f->SUTOPRSU(0,0,:) wks = gsn_open_wks("png","arpege_raster") res = True res@gsnMaximize = True res@sfXArray = g->bt42_lon(0,:) res@sfYArray = g->bt42_lat(0,:) res@mpPerimOn = False res@cnFillOn = True ; Turn on contour fill res@cnFillMode = "RasterFill" ; "RasterFill","CellFill" res@cnLinesOn = False ; Turn off contour lines res@cnFillPalette = "BlAqGrYeOrReVi200" res@cnLevelSelectionMode = "ExplicitLevels" res@cnLevels = fspan(0,25,201) ; 201 contour levels res@lbBoxLinesOn = False ; Turn off box lines res@mpProjection = "LambertEqualArea" ; Map projection contour = gsn_csm_contour_map(wks,data,res) ;---Calculate total elapsed time for this script. diff_time = get_cpu_time() - start_time print("=====> CPU Elapsed Time (" + res@cnFillMode + "): " + \ diff_time + " seconds <=====") end