;*********************** ; project.ncl ;*********************** load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" ;*********************** ;in_grib1 = addfile("batch-mars-web240-20130319140652-TnMhMt.grib","r") ;in_grib2 = addfile("batch-mars-web240-20130319145852-7dOvP2.grib","r") ;in_grib3 = addfile("batch-mars-web240-20130319174030-8SXgmI.grib","r") function create_map(wks,BBS) local a, res begin res = True res@gsnDraw = False ; don't draw yet res@gsnFrame = False ; don't advance frame yet res@gsnMaximize = True ; maximize plot in frame res@mpProjection = "LambertConformal" ; choose projection res@mpLambertParallel1F = 33 ; first parallel res@mpLambertParallel2F = 45 ; second parallel res@mpLambertMeridianF = -98 ; meridian res@mpLimitMode = "Corners" ; corner method of zoom res@mpLeftCornerLatF = 30 ; left corner res@mpLeftCornerLonF = -125 ; left corner res@mpRightCornerLatF = 65 ; right corner res@mpRightCornerLonF = -35 ; right corner res@pmTickMarkDisplayMode = "Always" ; turn on tickmarks res@tiMainString = "BBS data" ;---Create map. map = gsn_csm_map(wks,res) return(map) end ;************************************************* ; main code ;************************************************* begin f = addfile("s_06se12_Project.shp", "r") ; Open shapefile wks = gsn_open_wks("pdf","project2") ;---Create the map map = create_map(wks,"2009") ;---Attach the filled polygons pres = True poly = gsn_add_shapefile_polygons(wks,map,f,pres) ;---Drawing the map will also draw the attached polygons. draw(map) frame(wks) end ; ; Read data off shapefile ; ;segments = f->segments ;geometry = f->geometry ;segsDims = dimsizes(segments) ;geomDims = dimsizes(geometry) ; ; Read global attributes ; ;geom_segIndex = f@geom_segIndex ;geom_numSegs = f@geom_numSegs ;segs_xyzIndex = f@segs_xyzIndex ;segs_numPnts = f@segs_numPnts ;lines = new(segsDims(0),graphic) ; Array to hold polygons ;numFeatures = geomDims(0) ;x = f->LON ;y = f->LAT ;pgres = True ;pgres@gsEdgesOn = True ;pgres@gsLineColor = "red" ;pgres@gsFillOpacityF = 0.0 ;pgres@gsFillIndex = 0 ;---Resource list for shapefile outlines ;lnres = True ;lnres@gsLineThicknessF = 2. ; Default is 1 ;lnres@gsLineColor = "Red" ;pres = True ;dum = gsn_add_polygons(wks, plot,x,y, pgres) ;draw(plot) ;frame(wks) ;maximize_output(wks,False) ;end