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
Received on Mon Apr 07 19:41:15 2014
This archive was generated by hypermail 2.1.8 : Tue Apr 15 2014 - 10:45:19 MDT