#!/bin/sh -evx 

YEAR1=1850
YEARN=1949
CNAM=NHISTfrc2_f09_tn14_20191001
IDIR=/projects/NS9560K/noresm/cases/${CNAM}/ocn/hist
ODIR=/projects/NS2345K/www/diagnostics/esmvaltool/jene/KeyClim/netcdf/
WDIR=/scratch/$USER/tmp/ocn/ref
VAR=sss

mkdir -p $ODIR $WDIR
cd $WDIR 

for YYYY in `seq -w $YEAR1 $YEARN` 
do 
  FOUT2=NorESM2-MM.mon_${YYYY}.nc 
  echo ${FOUT2}
  for MM in `seq -w 01 12`
  do 
    FBGC=${CNAM}.micom.hm.${YYYY}-${MM}.nc 
    ncks -h -O -v ${VAR} -o in_bgc_${MM}.nc ${IDIR}/${FBGC}
    ncpdq -U -O in_bgc_${MM}.nc in_bgc_${MM}.nc
  done 
  ncrcat in_bgc_*.nc ${FOUT2}
  rm in_*.nc
done 

YEAR1=1950
YEARN=2014
CNAM=NHISTfrc2_f09_tn14_20191025
IDIR=/projects/NS9560K/noresm/cases/${CNAM}/ocn/hist

mkdir -p $ODIR $WDIR
cd $WDIR 

for YYYY in `seq -w $YEAR1 $YEARN` 
do 
  FOUT2=NorESM2-MM.mon_${YYYY}.nc 
  echo ${FOUT2}
  for MM in `seq -w 01 12`
  do 
    FBGC=${CNAM}.micom.hm.${YYYY}-${MM}.nc 
    ncks -h -O -v ${VAR} -o in_bgc_${MM}.nc ${IDIR}/${FBGC}
    ncpdq -U -O in_bgc_${MM}.nc in_bgc_${MM}.nc
  done 
  ncrcat in_bgc_*.nc ${FOUT2}
  rm in_*.nc
done 

YEAR1=2015
YEARN=2100
CNAM=NSSP585frc2_f09_tn14_20191105
IDIR=/projects/NS9560K/noresm/cases/${CNAM}/ocn/hist

for YYYY in `seq -w $YEAR1 $YEARN` 
do 
  FOUT2=NorESM2-MM.mon_${YYYY}.nc 
  echo ${FOUT2}
  for MM in `seq -w 01 12`
  do 
    FBGC=${CNAM}.micom.hm.${YYYY}-${MM}.nc 
    ncks -h -O -v ${VAR} -o in_bgc_${MM}.nc ${IDIR}/${FBGC}
    ncpdq -U -O in_bgc_${MM}.nc in_bgc_${MM}.nc
  done 
  ncrcat in_bgc_*.nc ${FOUT2}
  rm in_*.nc
done 

ncrcat NorESM2-MM.mon_*.nc NorESM2-MM_all_months.nc 
ncpdq NorESM2-MM_all_months.nc ${ODIR}/KeyCLIM_ref_${VAR}_1850-2100.nc

rm  NorESM2-MM*.nc