[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