Re: few data points in map

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
Received on Tue Jul 31 11:39:39 2012