[RegCNET] interpolation script
moetasim
mashfaq at purdue.edu
Fri Feb 17 03:04:33 CET 2006
Hi Jonathan,
There are regridding functions available in NCL, which can easily do this job
for you. Read in your guassian grid T-62 data file and also the HEAD_OUT.nc
file, which you get from post-processing, to grab lats/lons of RegCM grid, and
then you can easily convert gaussian grid to RegCM grid using the regridding
functions of NCL.
Other then this, you can write a simple program to read in your data file and
also the DOMAIN.INFO file and then use BILINX interpolation subroutine (can be
found in CRU code, I guess) to interpolate gaussian grid to RegCM grid. ECWCP
data is also at T42- gaussian grid so you can use its subroutines from ICBC to
get an idea how you may use ICBC subroutines to regrid your data.
Third option of regridding is to use Ferret. You may use something like
gxy=regcm_variable[d=?] to regrid your gaussian grid to RegCM grid. This would
require a single line in your ferret script.
And if your data is in grib format then you can look for ECMWF EMOSLIB
standard transformation subroutines which can be easily modified to do the
required job.
And I guess, you can also use CDO (Climate Data Operators) to interpolate
gaussain grid to the RegCM gird.
Moet
Quoting Jonathan Winter <jwinter at MIT.EDU>:
> I'm in the process of comparing some of my model output to
> observations/reanalysis data and I'd like to use some NCEP/NCAR data.
> Unfortunately, it uses a messy T62 Gaussian grid with 192x94 points
> (88.542N-88.542S, 0E-358.125E). I have a script that interpolates any
> regular grid to the RegCM grid, but in this dataset the resolution of
> the lat. varies, ie., the spacing between the first two points is 1.878,
> and the spacing between the 2nd and 3rd points is 1.898.
>
> To interpolate this data down to my RegCM grid, I'm going to need a
> script. If anyone has written such a script, has tackled this problem
> in the past, or has ideas, it would be a tremendous help.
>
> thanks,
>
> jonathan
>
> _______________________________________________
> RegCNET mailing list
> RegCNET at lists.ictp.it
> https://lists.ictp.it/mailman/listinfo/regcnet
>
More information about the RegCNET
mailing list