Re: Add overlay of ascii data from CSV file with missing values?

Hi Katrina,
Try this:

; untested, check that the coding below does what it is supposed to
dataj@_FillValue = -9999
datadjf = dataj(:,0:2)
datadjfgood = ind(.not.ismissing(datadjf(:,2)))

datadjf_good = dataj(datadjfgood,2)
lats_djfgood = dataj(datadjfgood,0)
lons_djfgood = dataj(datadjfgood,1)

plot(0)=gsn_csm_contour_map_polar(wks,precDJF,res)
dum1=gsn_add_polymarker(wks, plot(0), lons_djfgood, lats_djfgood, pmres)

Also, until you get this working, I would reduce dataj down to say size
(100,3). (Things might slow down a bit if you pass a large number of
points to gsn_add_polymarker. This way you'll know if things are working
sooner rather than later.)

Hope that helps. If not, please respond to ncl-talk.
Adam

On 04/07/2014 07:41 PM, Katrina Bennett wrote:
> Hello NCL Talk,
>
> I made a bit of progress on this but I'm still having issue with what
> I think is the subseting of the lat/lon from the missing values in the
> other columns. So, for this example below I'm working from the second
> ascii file, and then trying to subset it just for the first 3 columns
> (lat, lon, djf sig value with -9999s).
> However, when I ran this in my code I got a segmentation fault.
>
> I think the problem might be that I am not subseting the properly to
> screen out the data...?
>
> dataj@_FillValue = -9999
> datadjf = dataj(:,0:2)
> datadjfgood = ind(.not.ismissing(datadjf(:,2)))
>
> datadjf_good = dataj(datadjfgood)
> lats_djfgood = datadjf_good(:,0)
> lons_djfgood = datadjf_good(:,1)
>
> plot(0)=gsn_csm_contour_map_polar(wks,precDJF,res)
> dum1=gsn_add_polymarker(wks, plot(0), lons_djfgood, lats_djfgood,
> pmres)
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
______________________________________________________________
Adam Phillips                                asphilli@ucar.edu
NCAR/Climate and Global Dynamics Division       (303) 497-1726
P.O. Box 3000				
Boulder, CO 80307-3000    http://www.cgd.ucar.edu/staff/asphilli

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Apr 08 16:46:50 2014