This is a simpler version of the code. It just gives a map fully
covered in black.
For some reason, that set of lats/lons results in a problem.
Mike
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/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
wks = gsn_open_wks ("ps", "test")
plot=new(1,graphic)
setvalues wks
"wkColorMap" : "BlWhRe"
end setvalues
res = True
res_at_gsnFrame=False
res_at_tiMainFontHeightF = 0.033
res_at_tiMainString = "Plant Distribution"
res_at_tmYROn = False
res_at_tmXTOn = False
res_at_mpOceanFillColor = -1
res_at_mpLandFillColor = -1
res_at_mpMinLatF = 10.
res_at_mpMaxLatF = 75.
res_at_mpMinLonF = -180.00
res_at_mpMaxLonF = -50.00
res_at_mpCenterLonF = -115.
res_at_mpOutlineBoundarySets = "AllBoundaries"
res_at_mpGeophysicalLineThicknessF = 2.5
res_at_mpOutlineOn = True
plot = gsn_csm_map(wks,res)
lat=new((/8/),double)
lon=new((/8/),double)
lat=(/
54.4126815796,54.4171485901,54.4163818359,54.4119338989,54.4063796997,54
.4074935913,54.4108200073,54.4126815796/)
lon=
(/-130.6879730225,-130.674041748,-130.6725158691,-130.6683349609,-130.66
88842773,-130.6794433594,-130.686126709,-130.6879730225/)
lat!0="lat"
lat_at_units="degrees_north"
lat&lat=lat
lon!0="lon"
lon_at_units="degrees_east"
lon&lon=lon
res1=True
res1_at_gsFillScaleF = 1.0
res1_at_gsFillColor = "forestgreen"
gsn_polygon(wks,plot,lon,lat,res1)
frame(wks)
end
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jan 21 2009 - 11:59:54 MST
This archive was generated by hypermail 2.2.0 : Mon Feb 02 2009 - 07:53:06 MST