Re: few data points in map

Ok thanks.

-Fan

On 07/31/2012 02:53 PM, Mary Haley wrote:
> 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 13:30:33 2012