Re: few data points in map

If you upgrade to 6.1.0-beta or 6.0.0, this should fix your problem.

--Mary

On Jul 31, 2012, at 11:39 AM, Fan Fang wrote:

> Hi Mary,
>
> We use 5.2.1. Yes the X11 window flashes. If we download in a file the
> file got created but empty. It is also a ' Segmentation fault' so even
> if there is a plot, our system using NCL would quit without showing any
> plots.
>
> -Fan
>
> On 07/31/2012 01:13 PM, Mary Haley wrote:
>> Fan,
>>
>> Using the information you gave me below, I was still unable to reproduce the problem. I do get warnings, but I also get a plot.
>>
>> What version of NCL are you using?
>>
>> Can you send me your script offline so I can take a look?
>>
>> Here's the script I used to plot your data:
>>
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>>
>> a = addfile("test.nc","r")
>> data = a->Water_Vapor_Near_Infrared(0,:,:)
>> printVarSummary(data)
>>
>> wks = gsn_open_wks("x11","test")
>> res = True
>> res@cnFillMode = "RasterFill"
>> res@cnFillOn = True
>> res@mpMinLatF = min(data&latitude)
>> res@mpMaxLatF = max(data&latitude)
>> res@mpMinLonF = min(data&longitude)
>> res@mpMaxLonF = max(data&longitude)
>> res@mpCenterLonF = 26
>> res@gsnAddCyclic = False
>> plot = gsn_csm_contour_map_ce(wks,data,res)
>>
>>
>> --Mary
>>
>> On Jul 31, 2012, at 10:33 AM, Fan Fang wrote:
>>
>>> Hi Mary,
>>>
>>> The error message is
>>>
>>> fatal:A minimum of 3 coordinates points must be provided to
>>> CreateSplineApprox
>>> fatal:A minimum of 3 coordinates points must be provided to
>>> CreateSplineApprox
>>>
>>> The CDL attached below. The error is gone when we increase the grid
>>> points to 3x3.
>>>
>>> netcdf testmap.nc {
>>> dimensions:
>>> time = 1 ;
>>> latitude = 2 ;
>>> longitude = 2 ;
>>>
>>> variables:
>>> float latitude(latitude) ;
>>> latitude:units = "degrees_north" ;
>>> latitude:standard_name = "Latitude" ;
>>> float longitude(longitude) ;
>>> longitude:units = "degrees_east" ;
>>> longitude:standard_name = "Longitude" ;
>>> float Water_Vapor_Near_Infrared(time, latitude, longitude) ;
>>> Water_Vapor_Near_Infrared:long_name = "Total Column
>>> Precipitable Water Vapor - Near IR" ;
>>> Water_Vapor_Near_Infrared:units = "cm" ;
>>> Water_Vapor_Near_Infrared:quantity_type = "water Vapor" ;
>>>
>>> // global attributes:
>>> :DESCRIPTION = "Projected image" ;
>>> :BOUNDINGBOX = -15.f, 25.f, -6.f, 43.f ;
>>> :BBOX_SEQUENCE = "lat0,lon0,lat1,lon1" ;
>>> :LOWER_LEFT = "lat0,lon0" ;
>>> :UPPER_RIGHT = "lat1,lon1" ;
>>>
>>> data:
>>>
>>> latitude = -15.0, -14.0 ;
>>> longitude = 25.0, 27.0 ;
>>> Water_Vapor_Near_Infrared = -9 , -99 , -99 , -9 ;
>>> }
>>>
>>> On 07/31/2012 10:23 AM, Mary Haley wrote:
>>>> Dear Fan,
>>>>
>>>> Can you can provide a short example that illustrates the problem? Also, what is the error that you get?
>>>>
>>>> I was able to get a plot with 2 x 2 points, although I did get some errors:
>>>>
>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>>>>
>>>> data = (/ (/1,2/), (/3,4/) /)
>>>> data!0 = "lat"
>>>> data!1 = "lon"
>>>> data&lat = (/-10,10/)
>>>> data&lon = (/-20,20/)
>>>> data&lat@units = "degrees_north"
>>>> data&lon@units = "degrees_east"
>>>>
>>>> wks = gsn_open_wks("x11","test")
>>>> res = True
>>>> res@mpMinLatF = min(data&lat)
>>>> res@mpMaxLatF = max(data&lat)
>>>> res@mpMinLonF = min(data&lon)
>>>> res@mpMaxLonF = max(data&lon)
>>>> res@gsnAddCyclic = False
>>>> plot = gsn_csm_contour_map_ce(wks,data,res)
>>>>
>>>> --Mary
>>>>
>>>>
>>>> On Jul 31, 2012, at 7:54 AM, Fan Fang wrote:
>>>>
>>>>> Would gsn_csm_contour_map_ce() plot for few data points? It seems to
>>>>> complain without plotting when we feed less than 3x3 contiguous points.
>>>>> Also any limit for scattered points?
>>>>>
>>>>> Thanks.
>>>>>
>>>>> -Fan
>>>>> _______________________________________________
>>>>> 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
>
> _______________________________________________
> 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 Tue Jul 31 12:53:58 2012