Ooops!
I wan't clear enough. lat(0,:) .ne. lat(1,:) or for that
matter, few, if any, latitude and/or longitude "vectors" are
the same for the case of the original data.
After reading the documentation I thought this was OK for linint2.
Maybe not and I'm not seeing that in what I read?
Bob
On Tue, 28 Feb 2006, Dennis Shea wrote:
>>
>> I have been looking at the web page documenting linint, linint2 and
> linint2_points
>> and can't figure out which, if any of these might suit my needs.
>>
>> I have data on a polar sterographic grid with latitude(89,89) and
> longitude(89,89)
>> and data(89,89). All latitude(n,:) and/or longitude(n,:) are monotonically
>> increasing. I want to interpolate the data to:
>>
>> data2(72,144) grid with lat(72) and lon(144)
>>
>> I've managed to generate a number of (different) error messages with my
> attempts.
>>
>>
>> Can someone provide some guidance?
> =========================================
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>
> LAT = latitude(0,:)
> LON = longitude(0,:)
> LAT_at_units = "degrees_north"
> LON_at_units = "degrees_east"
>
> print(LAT) ; make sure they are monotonically increasing
> print(LON)
>
> data!0 = "LAT"
> data!1 = "LON"
>
> lon = lonGlobeFo (144, "lon", "longitude", "degrees_east")
> lat = latGlobeFo ( 72, "lat", "latitude" , "degrees_north")
>
> data2 = linint2_Warp(LON,LAT,data,True,lon,lat,0)
> printVarSummary(data2)
>
>
> HTH
> D
>
>
-- ------------------------------------------------------------- Robert A. Tomas National Center for Atmospheric Research Climate and Global Dynamics Division - Paleoclimatology 1850 Table Mesa Drive Boulder, CO 80305-3000 e-mail: tomas_at_ucar.edu tel: 303-497-1342 fax: 303-497-1348 _______________________________________________ ncl-talk mailing list ncl-talk_at_ucar.edu http://mailman.ucar.edu/mailman/listinfo/ncl-talkReceived on Tue Feb 28 2006 - 17:14:07 MST
This archive was generated by hypermail 2.2.0 : Wed Mar 01 2006 - 10:20:50 MST