<div dir="ltr">Hello,<div><br></div><div>I am having trouble plotting polygon data from a shapefile of a burn scar perimeter over a map of Colorado with counties.  I am using the example shapefiles_15.ncl from the <a href="https://www.ncl.ucar.edu/Applications/shapefiles.shtml" target="_blank">ncl page</a>. I am able to successfully plot the state of Colorado and its counties, but not the burn scar perimeter.  Burn scar shapefile can be found <a href="https://ftp.wildfire.gov/public/incident_specific_data/rocky_mtn/2018/SpringCreek/IR/20180715/PerimeterFileProvidedByIncident/" target="_blank">here</a>.   Can you please help? Thanks in advance!</div><div><br></div><div>Below is a copy of my ncl script:</div><div><br></div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"<br>load "./shapefile_utils.ncl"<br><br>;----------------------------------------------------------------------<br>; This function creates a cylindrical equidistant map of Colorado.<br>;----------------------------------------------------------------------<br><br>function create_colorado_map(wks,res,draw_ncl_outlines)<br>local a, mpres<br>begin<br>  mpres                       = res<br><br>  mpres@gsnMaximize           = True<br>  mpres@gsnPaperOrientation   = "portrait"<br><br>  mpres@gsnDraw               = False<br>  mpres@gsnFrame              = False<br><br>  mpres@mpFillOn              = False  <br><br>;---Turn on fancier tickmark labels.<br>  mpres@pmTickMarkDisplayMode = "Always"<br>  mpres@tmXBLabelFontHeightF  = 0.008      ; smaller tickmark labels<br><br>;---Zoom in on area of interest<br>  mpres@mpLimitMode           = "LatLon"<br>  mpres@mpMinLatF             =   37<br>  mpres@mpMaxLatF             =   41<br>  mpres@mpMinLonF             = -109.05<br>  mpres@mpMaxLonF             = -102.05<br>  mpres@mpFillOn              = False<br><br>  if(draw_ncl_outlines) then<br>    mpres@mpOutlineOn           = True<br>    mpres@mpOutlineBoundarySets    = "AllBoundaries"<br>    mpres@mpDataBaseVersion        = "MediumRes"<br>    mpres@mpDataSetName            = "Earth..4"      ; U.S. counties<br>  else<br>    mpres@mpOutlineOn              = False<br>  end if<br><br>;---Create map.<br>  map = gsn_csm_map(wks,mpres)<br><br>  return(map)<br>end<br><br>;--------------------------------------------------<br>; Main code<br>;--------------------------------------------------<br>begin<br>  wtype          = "png"              ; send graphics to PNG file<br>  wtype@wkWidth  = 2000<br>  wtype@wkHeight = 2000<br>  wks = gsn_open_wks(wtype,"shapefiles")<br><br>  ncl_version = get_ncl_version()<br><br>;---Create two maps of Colorado<br>  res                   = True<br>  res@tiMainFontHeightF = 0.015<br> <br>  res@tiMainString  = "Colorado counties - shapefile"<br>  map_shp = create_colorado_map(wks,res,False)<br><br>  res@tiMainString = "Colorado counties - NCL (version " + ncl_version + ")"<br>  map_ncl  = create_colorado_map(wks,res,True)<br><br>;---Add shapefiles to one of the maps<br>  lnres             = True<br>  lnres@gsLineColor = "red"<br><br>  plot_line_shp = gsn_add_shapefile_polylines(wks,map_shp,"./usa/gadm40_USA_2.shp",lnres)<br><br>;---Add a slightly transparent marker to both maps to show location of county updates<br>  mkres                  = True<br>  mkres@gsMarkerIndex    = 16        ; filled dot<br>  mkres@gsMarkerOpacityF = 0.5       ; make the marker half transparent<br>  mkres@gsMarkerColor    = "red"<br>  mkres@gsMarkerSizeF    = 50.<br><br>  counties_lat_center =   39.8<br>  counties_lon_center = -104.9<br> <br>  plot_marker_ncl = gsn_add_polymarker(wks,map_ncl,counties_lon_center,counties_lat_center,mkres)<br>  plot_marker_shp = gsn_add_polymarker(wks,map_shp,counties_lon_center,counties_lat_center,mkres)<br><br>print("Adding polygons to ncl map and shp map...")<br><br>filename = "Perimeter_20180715_SpringCreek_COHUX001313.shp"<br> <br> ;---Set some options for the polygons<br>   bscar = True<br>   bscar@gsLineColor = "Blue"<br>   bscar@gsLineThicknessF = 3.0  ; 3x thickness<br>   <br>id1 = gsn_add_shapefile_polygons(wks,map_ncl,filename,bscar)<br>id2 = gsn_add_shapefile_polygons(wks,map_shp,filename,bscar)<br><br>;---Panel both plots. Markers, lines, and polygons will be drawn.<br>  pres             = True<br>  pres@gsnMaximize = True<br>  gsn_panel(wks,(/map_ncl,map_shp/),(/1,2/),pres)<br> <br>end<div class="gmail-yj6qo"></div><div class="gmail-adL"><br></div></div><br class="gmail-Apple-interchange-newline"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font color="#000000">Vanessa Vincente<br>Associate Scientist<br>The COMET® Program<br>University Corporation for Atmospheric Research<br></font><br>*My working hours may differ from your own.  Please do not feel obligated to reply outside of standard business hours.*<div dir="ltr"></div></div></div></div></div></div></div></div></div></div></div>