Re: few data points in map

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