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