#!/bin/sh -evx

EXPS='volcanoes4future_cont volcanoes4future_sens4' 
EXPS='volcanoes4future_sens4' 
VARS='AODVIS DAYFOC'
IDIR=/nird/projects/NS2345K/ingo/Projects/V4F/Volcanoes4future 
ODIR=`pwd`
TMPDIR=/scratch/$USER/extract_bethkeetal2017 
NTHREAD=10 
export OPENBLAS_NUM_THREADS=1 

mkdir -p $TMPDIR
cd $TMPDIR
for EXP in $EXPS
do
  case $EXP in 
    volcanoes4future_sens4)
      EXPID=rcp45Volc
      YEAR1=2005
      YEARN=2099 
      MEM1=01
      MEMN=60  
    ;;
    volcanoes4future_cont)
      EXPID=historicalVolc
      YEAR1=1980
      YEARN=2005 
      MEM1=01
      MEMN=30  
    ;;
  esac 

  for VAR in $VARS 
  do
    THREAD=1 
    rm -f filelist_ncecat 
    for MEM in `seq -w $MEM1 $MEMN` 
    do
      rm -f filelist_ncrcat_$MEM
      for YYYY in `seq $YEAR1 $YEARN`
      do
        for MM in `seq -w 1 12`
        do
          echo ${EXP}_${MEM}.cam2.h0.${YYYY}-${MM}.nc >> filelist_ncrcat_$MEM 
        done
      done 
      OFILE=${VAR}_${EXPID}_mem${MEM}_${YEAR1}01-${YEARN}12.nc 
      if [ ! -e $OFILE ] 
      then 
        cat filelist_ncrcat_$MEM | ncrcat -O -t 1 -v $VAR -p $IDIR/$EXP/${EXP}_${MEM}/atm/hist -o $OFILE & 
	THREAD=$((THREAD+1)) 
	if [ $THREAD -gt $NTHREAD ]
	then 
	  THREAD=1 
	  wait
	fi
      fi 
      echo $OFILE >> filelist_ncecat
    done
    OPATH=$ODIR/${VAR}_${EXPID}_mem${MEM1}-${MEMN}_${YEAR1}01-${YEARN}12.nc
    cat filelist_ncecat | ncecat -O -t 1 -L 5 -u member -o $OPATH 
    ncpdq -O -t 1 -a time,member,lat,lon $OPATH $OPATH 
  done 

done