Re: NCL function "linint2_points_Wrap"

[1]
I am surprised there was no error message.
The linint2_points_Wrap is prototyped as follows

function linint2_points_Wrap (
                xi : numeric,
                yi : numeric,
                fi : numeric,
                fiCyclicX [1] : logical,
                xo [*] : numeric,
                yo [*] : numeric,
                Option [1] : integer
        )

the 'xo' and 'yo' must be one dimensional.

Your code indicates otherwise

  tintpre = linint2_points_Wrap(tdpre&longitude, tdpre&latitude, tdpre,\
                                True, ti@lat2d, ti@lon2d, 0)

Perhaps, the 2D attributes of a variable are passed as
one dimemsional [*] ?

===
[2]
Do the ti@lon2d range from -180 to +180 ?
If so, that is the problem because the NCEP data are 0-360

    ti@lon2d = where (ti@lon2d.lt.0, ti@lon2d + 360, ti@lon2d)

On 5/19/11 1:49 PM, Sha Feng wrote:
>> Hi everyone,
>>
>> I try to use the NCL function "linint2_points_Wrap" to match two datasets.
>>
>> One is from NCEP with dimension like (ntime, nz, nlat, nlon). Another is
>> CPR from CloudSat with dimension like (nbin, nray)and with attributes of
>> @lat2d and @lon2d. Now I need to interpolate the NCEP data to the CPR's
>> dimension in a selected region. Horizontally, I try to use
>> linint2_points_Wrap to accomplish this purpose. The output array, however,
>> are all missing value.
>>
>> I checked all the inputs and no one is full of missing value. Could anyone
>> help me to figure out it, or give me some suggestions how to accomplish my
>> purpose?
>>
>> The code is as follows, and attached is the output of all the "print".
>>
>> print(tdpre&longitude)
>> print("-------------------------------------------")
>> print(tdpre&latitude)
>> print("-------------------------------------------")
>> print(tdpre(1,:,:))
>> print("-------------------------------------------")
>> print(ti@lat2d)
>> print("-------------------------------------------")
>> print(ti@lon2d)
>>
>> tintpre = linint2_points_Wrap(tdpre&longitude, tdpre&latitude, tdpre,
>> True, ti@lat2d, ti@lon2d, 0)
>> print(tintpre)
>>
>> Thanks and regards,
>>
>> S.
>>
>>
>>>
>>>
>>
>>
>>
>> _______________________________________________
>> 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 Thu May 19 17:05:17 2011