#!/bin/sh
set -e
# interpolate from a global Gaussian grid to uniform unshifted lat/lon grid 144x73 (2.5x2.5 degree)

inp=$1
out=$2

# interpolate
echo "C* GGALL    144   73    5" | ccc /home/rsk/bin_linux/ggall $inp ${inp}_145x73

# strip off cyclic meridian
echo "C*WINDOW      1  144    1   73" | ccc window ${inp}_145x73 $out

release .ccc_cards ${inp}_145x73