#!/bin/sh -e

IDIR=/work/shared/nersc/soutten/Volcanos4Future
ODIR=/work-common/shared/bjerknes/ingo/VOLCANOS4FUTURE/noresm-ver1_r112-r169/inputdata
PREFIX=volc4fut_volcnic_prototype_

mkdir -p $ODIR

for COUNT in `seq -w 60`
do
  echo ${PREFIX}${COUNT}.nc
  IPATH=${IDIR}/${PREFIX}${COUNT}.nc 
  OPATH=${ODIR}/${PREFIX}${COUNT}.nc
  ncdump $IPATH | sed -e 's/date(date)/date(time)/' -e '/date = 1130/d' -e 's/FillValue/_FillValue/g' | ncgen -o $OPATH -
  #
  ncks -O -h -F -d time,1,12 ${OPATH} ${OPATH}_year2005 
  ncrename -h -v date,dateold -v time,timeold -v colmass,colmassold -v MMRVOLC,MMRVOLCOLD ${OPATH}_year2005 
  ncap2 -O -h -s 'date=dateold-10000;time=timeold-1;colmass=colmassold*0;MMRVOLC=MMRVOLCOLD*0;' ${OPATH}_year2005 ${OPATH}_year2005  
  ncap2 -O -h -s 'MMRVOLC@units=MMRVOLCOLD@units;MMRVOLC@_FillValue=MMRVOLCOLD@_FillValue;MMRVOLC@long_name=MMRVOLCOLD@long_name;colmass@_FillValue=colmassold@_FillValue' ${OPATH}_year2005 ${OPATH}_year2005  
  ncrcat -h -O -v time,date,datesec,lat,lev,colmass,MMRVOLC ${OPATH}_year2005 ${OPATH} ${OPATH} 
done