;---------------------------------------------------------------------- ; This script draws a county map of Colorado, to show that Broomfield ; is now a county in NCL V6.4.0. If you have an older version of NCL, ; you will not see Broomfield. ;---------------------------------------------------------------------- wks = gsn_open_wks("png","colorado_map_640") res = True res@gsnMaximize = True res@gsnDraw = False res@gsnFrame = False res@pmTickMarkDisplayMode = "Always" ; nicer tickmarks res@mpLimitMode = "LatLon" res@mpMinLatF = 37 ; zoom in on Colorado res@mpMaxLatF = 41 res@mpMinLonF = -109.05 res@mpMaxLonF = -102.05 res@mpFillOn = False ; turn off map fill res@mpOutlineOn = True res@mpOutlineBoundarySets = "AllBoundaries" res@mpDataBaseVersion = "MediumRes" res@mpDataSetName = "Earth..2" ; U.S. counties map = gsn_csm_map(wks,res) ;--Add box to show location of updated counties lat_min = 39.5 lat_max = 40.2 lon_min = -105.2 lon_max = -104.4 lnres = True lnres@gsLineColor = "NavyBlue" lnres@gsLineThicknessF = 3.0 id = gsn_add_polyline(wks,map,(/lon_min,lon_max,lon_max,lon_min,lon_min/),\ (/lat_min,lat_min,lat_max,lat_max,lat_min/),lnres) draw(map) frame(wks)