Re: Decreasing coordinate array in linint2_points_Wrap function

Oops ....

    Use NCL syntax to pass the lat/lon associated with the variable.
    I am using the names 'lat' and 'lon' as placeholders

    xPoints = linint2_points_Wrap(x&lon, x&lat, x, False, lono, lato,0)

Dennis Shea wrote:
> Hi Lonneke
>
> f = addfile ("....grb", "r")
> x = f->FOO ; (;,:,:) N=>S
> printVarSummary(x) ; original order
>
> x = x(:,::-1,:) ; reverse lat CV *and* data array
> printVarSummary(x) ; Note the latitude CV
>
> The solution is simple but first a few comments about NCL variables.
>
> printVarSummary is a procedure I use more than any other.
>
> Good luck
> D
>
>
>
> Lonneke Goddijn-Murphy wrote:
>> Hi all,
>>
>> I want to interpolate ECMWF model data to fit stations locations using
>> the function
>>
>> fo = *linint2_points_Wrap*(xi, yi, fi, False, xo, yo, 0)
>>
>> The variables of the ECMWF data, fi, are of the form fi(ti, lati, loni),
>> so if I understand correctly,
>>
>> fo = *linint2_points_Wrap*(loni, lati, fi, False, lon_station,
>> lat_station,0)
>>
>> should give me the interpolated values at stations locations lon_station
>> and lat_station, at times ti (because the two rightmost dimensions
>> (/nyi/ x /nxi/) are the dimensions to be used in the interpolation)
>> The loni and lati should be increasing arrays, but the ECMWF latitudes
>> are decreasing lati= [60, 59.75, 59.5, ... 50]. Can I then simply use
>> -lati and apply *linint2_points_Wrap*?
>>
>> fo = *linint2_points_Wrap*(loni, -lati, fi, False, lon_station,
>> -lat_station, 0)
>>
>> Hope I made myself clear, and you can give me advise,
>> Thanks
>>
>> Lonneke
>>
>>
>>
>>
>>
>>
>> ------------------------------------------------------------------------
>>
>> _______________________________________________
>> 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 Thu Jul 02 2009 - 07:28:21 MDT