Re: gsn_add_polymarker problem

On Apr 10, 2012, at 6:50 PM, Cecille Villanueva-Birriel wrote:

> Hello,
>
> I am trying to mark my locations on top of my contour map using "gsn_add_polymarker". However, for some reason is only marking a few of them and not all of them. I looked for errors and I can;t figure out why is not plotting all the locations. The only error I get is
> warning:asciiread: End of file reached and only (325) elements were read from the file, filling remaining elements with the default missing value for the requested type
>
> but I am not sure of this error is related to the markers not been plotted. I've copied my code below.

Cecille,

Since you only have one call to gsn_add_polymarker, and since some of these markers are appearing and some aren't, then one of the following must be happening:

     - rlon and/or rlat contain missing values
     - rlon and/or rlat are falling outside the range of the map you're plotting

The warning you're getting may be the issue, because it's telling you that you don't have enough values to fill the 65 x 7 array that you've told it to read in. This means it's going to fill the rest of the "data" array with missing values, which potentially means that rlat and/or rlon contain missing values.

The warning is also telling you that it's only read in 325 values, which would be the equivalent of a 65 x 5 array. My guess is that your ASCII file is being read incorrectly, and that rlat and rlon are getting the wrong values. This would happen if you told it to read 7 columns of data, and you really only have 5.

--Mary

>
> Cecille
>
> ;*************************************************
> ; rdm2grid_2.ncl
> ;*************************************************
> ; interpolates sparse random data to grid for smoother
> ; contours, written by SGLT Dec 2011
> ;***********************************************
>
> 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"
>
> begin
> ;
> ;------------------------------------------------------------------
> ; Read data
> ;------------------------------------------------------------------
> ;
> diri = "./"
> ; fili = "Occurrences.txt"
> ; fili = "CCN.txt"
> ; fili = "CBtemp.txt"
> ; fili = "FLdepth.txt"
> ; fili = "WatFL.txt"
> fili = "timetoFL.txt"
> ; fili = "Alwc_FL.txt"
> ; fili = "PwcFL.txt"
> ; fili = "Pwc12C.txt"
> ; fili = "CCN_sens.txt" ;%changes
> ; fili = "CCN_sens_threshdelta.txt"
> ; fili = "Driz_FL.txt"
> ; fili2 = "Locations.txt"
> ; fili2 = "Locations_old.txt"
> ; data = asciiread(diri+fili,(/65,5/),"float")
> ; data = asciiread(diri+fili,(/65,3/),"float")
> ; data = asciiread(diri+fili,(/65,6/),"float") ;CCN sensitivity file %
> data = asciiread(diri+fili,(/65,7/),"float") ;CCN sensitivity file
> ; data = asciiread(diri+fili,(/50,5/),"float")
>
> ; below not needed if coords in data file above
> ; data2 = asciiread(diri+fili2,(/65,2/),"float")
> ; data2 = asciiread(diri+fili2,(/50,2/),"float")
>
> rlat = data(:,1)
> rlon = data(:,0)
> ; Z = data(:,2) ; # sndgs going into averaged past sndgs
> ; Z1 = data(:,3) ; # sndgs going into averaged future sndgs
> ; Z2 = data(:,4) ; # difference in past and future # sndgs
> ; Z = data(:,2) ; # xxx in past sndgs
> ; Z1 = data(:,3) ; # xxx in future sndgs
> ; Z2 = data(:,4) ; # difference of xxx in past and future
> ; Z = data(:,4) ; # xxx in past sndgs -CCN sens
> ; Z1 = data(:,5) ; # xxx in future sndgs -CCN sens
>
> Z = data(:,2) ; C values
> ; Z = data(:,4) ; FL depth change in m
> ; Z = data(:,4) ; change in w at FL in m/s
> ; Z = data(:,4) ; change in time to FL in s
> ; Z = data(:,4) ; change in alwc at FL in g/m3
> ; lat = data2(:,1)
> ; lon = data2(:,0)
> ;
> rlat@units = "degrees_north"
> rlon@units = "degrees_east"
> ; lat@units = "degrees_north"
> ; lon@units = "degrees_east"
> ;------------------------------------------------------------------
> ; Read grid point locations-- for interpolating data to grid
> ;------------------------------------------------------------------
> diri = "./"
> fil1 = "grid_lat.dat"
> fil2 = "grid_lon.dat"
> glat = asciiread(diri+fil1,30,"float")
> glon = asciiread(diri+fil2,70,"float")
>
> glat@units = "degrees_north"
> glon@units = "degrees_east"
> ;******************************************************************
> ; cssgrid
> ;******************************************************************
> zcss = cssgrid_Wrap (rlat, rlon, Z, glat,glon)
> ; zcss = cssgrid_Wrap (rlat, rlon, Z1, glat,glon)
> ; zcss = cssgrid_Wrap (rlat, rlon, Z2, glat,glon)
> ;------------------------------------------------------------------
> ; plot
> ;------------------------------------------------------------------
>
> ; wks = gsn_open_wks("ncgm","Occurrences_diff")
> ; wks = gsn_open_wks("ncgm","Occurrences_future")
> ; wks = gsn_open_wks("ncgm","Occurrences_past")
> ; wks = gsn_open_wks("ncgm","CCN_Cvalues")
> ; wks = gsn_open_wks("ncgm","TESTold") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","CBtemp_past") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","CBtemp_future") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","CBtemp_diff") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","FLdepth_past") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","FLdepth_future") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","FLdepth_diff") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","WatFL_past") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","WatFL_future") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","WatFL_diff") ; Open a workstation.
> wks = gsn_open_wks("ncgm","timetoFL_past") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","timetoFL_fut") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","timetoFL_diff") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","AlwcFL_past") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","AlwcFL_fut") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","AlwcFL_diff") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","PwcFL_past") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","PwcFL_fut") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","PwcFL_diff") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","Pwc12C_past") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","Pwc12C_fut") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","Pwc12C_diff") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","DrizFL_past") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","DrizFL_fut") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","DrizFL_diff") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","CCNsens_past") ; Open a workstation.
> ; wks = gsn_open_wks("ncgm","CCNsens_fut") ; Open a workstation.
>
> getvalues wks
> "wkColorMap" : cmap ; Retrieve the default color map.
> end getvalues
>
> cmap(0,:) = (/1.,1.,1./) ; Change background to white.
> cmap(1,:) = (/0.,0.,0./) ; and foreground to black.
> gsn_define_colormap(wks, cmap)
>
> ; gsn_define_colormap (wks,"GreenYellow")
> ; gsn_define_colormap (wks,"BlWhRe")
> ; gsn_define_colormap(wks,"amwg")
>
> res = True
> ;
> ; res@tiXAxisString = ":F25:longitude"
> ; res@tiYAxisString = ":F25:latitude"
>
> res@cnLinesOn = False ; turn off contour lines
> res@cnFillOn = True ; turn on color
> res@cnInfoLabelOn = False ; turn off info label
> res@pmLabelBarDisplayMode = "Always"
> res@lbPerimOn = False ;turn off label bar perimeter
>
> res@mpLimitMode = "LatLon"
> res@mpMinLonF = -127.0
> res@mpMaxLonF = -65.0
> res@mpMinLatF = 23.0
> res@mpMaxLatF = 52.0
> ; new
> ; res@sfXCStartV = -127.
> ; res@sfXCEndV = -58.
> ; res@sfYCStartV = 23.
> ; res@sfYCEndV = 52.
> res@gsnAddCyclic = False
>
> res@mpPerimOn = True ; Turn on map perimeter.
>
> ; res@vpXF = 0.1 ; Change the size and location of the
> ; res@vpYF = 0.9 ; plot on the viewport.
> ; res@vpWidthF = 0.7
> ; res@vpHeightF = 0.7
>
> res@gsnScale = True ; Force X/Y axes labels to be the same size.
>
> txres = True ; Set some resources
> txres@txFontHeightF = 0.010 ; for a text string.
> txres@txFontColor = 1
> ; gsn_text_ndc(wks,":F25:C (cm-3)",.45,.25,txres) ; Draw a text string.
> ; gsn_text_ndc(wks,":F25:deg C",1.45,.25,txres) ; Draw a text string.
> ; gsn_text_ndc(wks,":F25:meters ",.45,.25,txres) ; Draw a text string.
> ; gsn_text_ndc(wks,":F25:m/s ",.45,.25,txres) ; Draw a text string.
> gsn_text_ndc(wks,":F25:seconds ",.45,.25,txres) ; Draw a text string.
> ; gsn_text_ndc(wks,":F25:g/m3 ",.45,.25,txres) ; Draw a text string.
> ; on the viewport.
>
> delete(cmap)
> cmap = (/(/1.00, 1.00, 1.00/), (/0.00, 0.00, 0.00/), \
> (/.560, .500, .700/), (/.300, .300, .700/), \
> (/.100, .100, .700/), (/.000, .100, .700/), \
> (/.000, .300, .700/), (/.000, .500, .500/), \
> (/.000, .400, .200/), (/.000, .600, .000/), \
> (/.000, 1.00, 0.00/), (/.550, .550, .000/), \
> (/.570, .420, .000/), (/.700, .285, .000/), \
> (/.700, .180, .000/), (/.870, .050, .000/), \
> (/1.00, .000, .000/), (/.700, .700, .700/)/)
>
> ; gsn_define_colormap(wks, cmap) ; Define a new color map.
> ; gsn_define_colormap(wks, "BlueDarkOrange18") ; Define a new color map.
> ; gsn_reverse_colormap(wks) ; reverse the color map.
> gsn_define_colormap(wks, "BrownBlue12") ; Define a new color map.
> ;gsn_define_colormap(wks, "BlueYellowRed") ; Define a new color map.
> ; gsn_define_colormap(wks, "gsltod") ; Define a new color map.
>
> res@mpFillOn = True ; Turn on map fill.
> res@mpFillAreaSpecifiers = (/"Water","Land","USStatesWater"/)
> ; res@mpSpecifiedFillColors = (/17,17,17/)
> ;res@mpSpecifiedFillColors = (/1,1,1/)
> res@mpSpecifiedFillColors = (/0,0,0/)
>
> res@lbLabelFontHeightF = 0.008 ; label bar font height
> res@lbLabelStride = 4
>
> res@gsnDraw = False ;so we can add poly stuff
> res@gsnFrame = False ; Don't advance frame.
>
> res@mpAreaMaskingOn = True ; Indicate we want to
> res@mpMaskAreaSpecifiers = "USStatesLand" ; mask land.
> res@mpPerimOn = True ; Turn on a perimeter.
> ; res@mpGridMaskMode = "MaskLand" ; Mask grid over land.
> res@cnFillDrawOrder = "Predraw" ; Draw contours first.
>
> ; res@sfXArray = rlon
> ; res@sfYArray = rlat
>
> res@gsnSpreadColors = True
> ; res@gsnSpreadColorStart = 12 ;for Bluedarkorange18
> ; res@gsnSpreadColorEnd = 20 ;for Bluedarkorange18
> ; res@gsnSpreadColorStart = 8 ;for Bluedarkorange18
> ; res@gsnSpreadColorEnd = 13 ;for Bluedarkorange18
> ; res@gsnSpreadColorStart = 8 ;for Bluedarkorange18
> ; res@gsnSpreadColorEnd = 13 ;for Bluedarkorange18
> res@gsnSpreadColorStart = 2 ;for Bluedarkorange18
> res@gsnSpreadColorEnd = 13 ;for Bluedarkorange18
>
>
> res@cnLevelSelectionMode = "ExplicitLevels" ; explicit levels
> res@cnLevels = (/700.,900.,1100.,1300.,1500.,1700.,\
> 1900.,2100.,2300.,2500.,2700.,2900.,3100/) ;set levels
> ; res@cnLevels = (/10.,12.,14.,16.,18.,20.,22./) ;CB temp
> ; res@cnLevels = (/0.5,1.,1.5,2.,2.5,3./) ;CB temp diff
> ; res@cnLevels = (/0.,100.,200.,300.,400.,500.,600.,\
> ; 700,800,900.,1000.,1100.,1200.,1300.,1400.,1500./) ;set levels
> ; res@cnLevels = (/400.,500.,600.,700.,900.,1000.,\
> ; 1200,1500,2000.,2500.,3000.,4000./) ;set levels
> ; res@cnLevels = (/1500.,2000.,2500.,3000.,3500.,4000.,\
> ; 4500.,5000.,5500.,6000./) ;FL depth
> ; res@cnLevels = (/300.,400.,500.,600.,700.,800.,900./) ;FL depth diffs
> ; res@cnLevels = (/8.,9.,10.,11.,12.,13.,14.,15.,16.,17.,18.,19.,20./)
> ; res@cnLevels = (/1.,1.5,2.,2.5,3.,3.5,4./) ;w diffs
> ; res@cnLevels = (/-2.,-1.5,-1.,-0.5,0.,0.5,1.,1.5,2.,2.5,3./)
> ; res@cnLevels = (/0.,1.,2.,3.,4./)
> ; res@cnLevels = (/2.,2.5,3.,3.5,4.,4.5,5.,5.5,6.,6.5,7./) ;alwc
> ; res@cnLevels = (/0.2,0.4,0.6,0.8,1./) ;alwc diff
> ; res@cnLevels = (/0.5,1.,1.5,2.,2.5,3.,3.5,4.,4.5,5.,5.5,6.,6.5/) ;pwc
> ; res@cnLevels = (/0.5,1.,1.5,2.,2.5/) ;pwc diff
> ; res@cnLevels = (/0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.,1.1,1.2/) ;driz
> ; res@cnLevels = (/0.001,0.01,0.03,0.05,0.07,0.09,0.1/) ;driz diffs
> ; res@cnLevels = (/-100.,-10.,-1.,-0.1,0.1/) ;CCN sensitivity%
> ; res@cnLevels = (/-1.,-0.75,-0.5,-0.25,0./) ;CCN sensitivity delta
>
> res@pmLabelBarSide = "Bottom"
> res@lbOrientation = "Horizontal"
> res@lbLabelFontHeightF = 0.01
> ; res@tiMainString = "# Sndgs going into avgd profiles"
> ; res@tiMainString = "C values for CCN"
> ; res@tiMainString = "Temp at CB"
> ; res@tiMainString = "Depth to FL"
> ; res@tiMainString = "Vertical velocity at FL"
> ; res@tiMainString = "Change in parcel ascent time to FL"
> ; res@tiMainString = "ALWC at FL"
> ; res@tiMainString = "PWC at FL"
> ; res@tiMainString = "Driz at FL"
> ; res@tiMainString = "PWC at -12C"
> ; res@tiMainString = "PWC sensitivity to CCN ,% change from C=700 to 2000"
> ; res@tiMainString = "PWC sensitivity to CCN , delta values"
> res@tiMainString = "Time to FL"
> map = gsn_csm_contour_map_ce(wks,zcss,res)
> ; --------------plot sites--------
> res_mark = True
> ; res_mark@gsMarkerIndex = 17 ; Filled circle
> ;res_mark@gsMarkerIndex = 16 ; dots
> res_mark@gsMarkerIndex = 2 ; crosses
> res_mark@gsMarkerSizeF = 15.
> ; res_mark@gsMarkerColor = "White"
> res_mark@gsMarkerColor = "Black"
>
> ; add polymarkers
> dum = gsn_add_polymarker(wks, map, rlon, rlat, res_mark)
> ; dum = gsn_add_polymarker(wks, map, lon, lat, res_mark)
> draw(map)
> frame(wks)
> end
>
>
>
>> ********************************************************
>> Cecille M. Villanueva Birriel
>> Graduate Research Assistant
>> Cloud Microphysics Research Group
>> Purdue University
>> Department of Earth and Atmospheric Sciences
>> 550 Stadium Mall Drive, West Lafayette, IN 47907-2051
>> email: cvillanu@purdue.edu
>> ********************************************************
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Apr 11 08:29:41 2012