#!/bin/sh -evx 

YEAR1=1850
YEARN=2014
CNAM=NHISTfrc2_f09_tn14_cloudxpice_220629
IDIR=/projects/NS9252K/noresm/cases/${CNAM}/ocn/hist
ODIR=/projects/NS2345K/www/diagnostics/esmvaltool/jene/KeyClim/pH/
WDIR=/scratch/$USER/tmp/ocn/cloud

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}.blom.hbgcm.${YYYY}-${MM}.nc 
    ncks -h -O -v phlvl -d depth,0 -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_cloudxpice_220819
IDIR=/projects/NS9252K/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}.blom.hbgcm.${YYYY}-${MM}.nc 
    ncks -h -O -v phlvl -d depth,0 -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_cloud2_Surface_pH_1850-2100.nc

rm  NorESM2-MM*.nc


#echo Preparing pH iceSheet3_hist
#ncrcat -v phlvl /projects/NS9252K/users/heig/NorESM/archive/NHISTfrc2G_f09_tn14_gl4_SMB1_celo/ocn/hist/*.hbgcm.*.nc /projects/NS2345K/www/diagnostics/esmvaltool/jene/KeyClim/tmp.nc
#ncks -h -O -v phlvl -d depth,0 -o /projects/NS2345K/www/diagnostics/esmvaltool/jene/KeyClim/ph_iceSheet3_hist.nc /projects/NS2345K/www/diagnostics/esmvaltool/jene/KeyClim/tmp.nc
#rm -rf /projects/NS2345K/www/diagnostics/esmvaltool/jene/KeyClim/tmp.nc

#echo Preparing pH iceSheet3_ssp585
#ncrcat -v phlvl /projects/NS9252K/users/heig/NorESM/archive/NSSP585frc2G_f09_tn14_gl4_SMB1_celo/ocn/hist/*.hbgcm.*.nc /projects/NS2345K/www/diagnostics/esmvaltool/jene/KeyClim/tmp.nc
#ncks -h -O -v phlvl -d depth,0 -o /projects/NS2345K/www/diagnostics/esmvaltool/jene/KeyClim/ph_iceSheet3_ssp585.nc /projects/NS2345K/www/diagnostics/esmvaltool/jene/KeyClim/tmp.nc
#rm -rf /projects/NS2345K/www/diagnostics/esmvaltool/jene/KeyClim/tmp.nc