#!/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/cam_interp
cd /scratch/$USER/cam_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/atm/hist/N20TREXTAERCN_f19_g16_TAUOCN${MEM}_PAC.icefrac.cam2.h0.${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/atm/hist/N20TREXTAERCN_f19_g16_TAUOCN${MEM}_PAC.cam2.h0.${YR}-${MON}.nc

# prepare input file for interpolation
ncpdq -6 -O -U $IFILE in1_${YR}.${MON}.nc
cdo select,name=ICEFRAC in1_${YR}.${MON}.nc in_${YR}.${MON}.nc
done
done

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


# convert to netcdf4 and pack
#ncpdq -M flt_sht -4 -L 5 -h -O in_${YR1}-${YRN}.nc $OFILE
done
# remove temporary directory 
rm -rf /scratch/$USER/cam_interp