#!/bin/sh -evx

# load modules
# . /opt/modules/default/init/sh
module load netcdf nco cdo

BASEDIR=/projects/NS9039K/shared/noresm/cases

mkdir -p /scratch/$USER/micom_interp
cd /scratch/$USER/micom_interp

MEM1=1
MEMN=6
YR1=1870
YRN=2012
MON1=1
MONN=12

for MEM in `seq -w $MEM1 $MEMN`
do
rm -f in_????.??.nc
OFILE=${BASEDIR}/TAUOCN-PAC_ext/N20TREXTAERCN_f19_g16_TAUOCN${MEM}_PAC/ocn/hist/N20TREXTAERCN_f19_g16_TAUOCN${MEM}_PAC.tauy.micom.hm.${YR1}-${YRN}.nc
for YR in `seq -w $YR1 $YRN`
do
for MON in `seq -w $MON1 $MONN`
do

# set paths etc
IFILE=${BASEDIR}/TAUOCN-PAC_ext/N20TREXTAERCN_f19_g16_TAUOCN${MEM}_PAC/ocn/hist/N20TREXTAERCN_f19_g16_TAUOCN${MEM}_PAC.micom.hm.${YR}-${MON}.nc
GFILE=${BASEDIR}/grid.nc

# prepare input file for interpolation
ncpdq -6 -O -U -v tauy $IFILE in_${YR}.${MON}.nc

done
done

# merge files 
ncrcat -h -O in_????.??.nc in_${YR1}-${YRN}.nc

# add coordinates
ncks -6 -A -v vlon,vlat -o in_${YR1}-${YRN}.nc $GFILE 

# interpolate
cdo remapbil,r360x180 in_${YR1}-${YRN}.nc $OFILE

# convert to netcdf4 and pack
#ncpdq -M flt_sht -4 -L 5 -h -O $OFILE $OFILE

# remove temporary directory
done  
rm -rf /scratch/$USER/micom_interp