#!/bin/sh
# interpolate from Gaussian grid to 2.5x2.5 grid (144x73): lon=0,2.5,...,357.5; lat=-90,-87.5,...,87.5,90

# 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     = 73
xfirst    = 0
xinc      = 2.5
yfirst    = -90
yinc      = 2.5" > grid_144x73.txt

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

rm -f grid_144x73.txt