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
Received on Thu Jul 02 2009 - 07:23:49 MDT
This archive was generated by hypermail 2.2.0 : Thu Jul 02 2009 - 11:39:38 MDT