#!/bin/bash

gcms=(
    #access1.0_RCP45
    #access1.0_RCP85
    #access1.3_RCP85
    #bcc-csm1.1_RCP85
    #canesm2_RCP85
    #ccsm4_RCP85
    #csiro-mk3.6.0_RCP85
    #fgoals-g2_RCP85
    #gfdl-cm3_RCP85
    #giss-e2-h_RCP85
    #miroc5_RCP85
    #mri-cgcm3_RCP85
    #noresm1-m_RCP85
    pnnl
)

merge=~/workdir/script_repo/data/mergeOnly.py
reform=~/workdir/script_repo/data/reformat_streamflow.py
stats=mk_stats.py
ftype='Streamflow.Only'
#ftype='Aggregated.Values'

run_merge() {
    for gcm in ${gcms[@]}; do
	echo $gcm
	python $merge $ftype ../data/${gcm}/ ../merge/
    done
}

run_extract() {
    for gcm in ${gcms[@]}; do
	echo $gcm
	p=`find ../merge/${gcm}*/Streamflow.Only`
	python $reform $p --steps 1h --units cfs
    done
}

run_stats() {
    for gcm in ${gcms[@]}; do
	echo $gcm
	p=`find /home/disk/rocinante/DATA/temp/chico/hyak/data_atmos/${gcm}*`
	python $stats $p
    done
}

run_extract_all() {
    find ../merge/*RCP*/Streamflow.Only -maxdepth 0 | parallel -j13 --eta python $reform {} --steps 1h --units cfs
}

run_stats_all() {    
    find ../merge/*RCP* -maxdepth 0 | parallel -j13 --eta python $stats {}
}


#run_merge
#run_extract
run_stats

#run_extract_all
#run_stats_all