#!/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

YR1=1870
YRN=2012
MON1=1
MONN=12

rm -f in_????.??.nc
rm -f in1_????.??.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_TAUOCN1_PAC/atm/hist/N20TREXTAERCN_f19_g16_TAUOCN1_PAC.cam2.h0.${YR}-${MON}.nc
OFILE=${BASEDIR}/TAUOCN-PAC_ext/N20TREXTAERCN_f19_g16_TAUOCN1_PAC/atm/hist/N20TREXTAERCN_f19_g16_TAUOCN1_PAC.sst.cam2.h0.${YR1}-${YRN}.nc

# prepare input file for interpolation
ncpdq -6 -O -U $IFILE in1_${YR}.${MON}.nc
#cdo select,name=V,PS in1_${YR}.${MON}.nc in_${YR}.${MON}.nc
cdo select,name=SST 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

# interpolate onto pressure levels
#cdo ml2pl,100000,92500,85000,70000,50000,30000,20000 -chname,PS,sp in_${YR1}-${YRN}.nc $OFILE

# convert to netcdf4 and pack
#ncpdq -M flt_sht -4 -L 5 -h -O $OFILE $OFILE
#ncpdq -M flt_sht -4 -L 5 -h -O in_${YR1}-${YRN}.nc $OFILE

# remove temporary directory 
rm -rf /scratch/$USER/cam_interp