#!/bin/sh
# interpolate from Gaussian grid to 2.5x2.5 grid (144x41): lon=0,2.5,...,357.5; lat=-50,-47.5,...,47.5,40

# grid description file
echo "#
# gridID 0
#
gridtype  = lonlat
gridsize  = 10512
xname     = longitude
xlongname = longitude
xunits    = degrees_east
yname     = latitude
ylongname = latitude
yunits    = degrees_north
xsize     = 144
ysize     = 41
xfirst    = 0
xinc      = 2.5
yfirst    = -50
yinc      = 2.5" > grid_144x41.txt

ifile=$1
ofile=$2
cdo remapcon,grid_144x41.txt $ifile $ofile

rm -f grid_144x41.txt