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

Hello Adam, thanks so much for you help. That works in my code! Thanks
again.

Cheers,
Katrina

~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Katrina E. Bennett
PhD Candidate
University of Alaska Fairbanks
International Arctic Research Center
907-474-1939 kebennett@alaska.edu

On Tue, Apr 8, 2014 at 2:46 PM, Adam Phillips <asphilli@ucar.edu> wrote:

> 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 20:35:56 2014