#!/bin/bash

export ddir='/home/disk/rocinante/DATA/temp/chico/hyak/data_rbm/'
export odir='/home/disk/rocinante/DATA/temp/chico/hyak/writeup/pub_temp/'
export rdir='/home/disk/rocinante/DATA/temp/chico/hyak/data_ref/air_temp/'
export pdir='/home/disk/rocinante/DATA/temp/chico/hyak/writeup/plots/'
tbl='/home/disk/rocinante/DATA/temp/chico/hyak/data_ref/tbl.csv'

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
	  kitsap_run84_pnnl)

pds=('1990s' '2050s' '2080s', 'full')

function process_site {
    ifile=$1.hour
    site=`basename $1`
    gcm=$2
    prd=$3

    echo $site $gcm $prd
    out=$odir/$site/
    pfx=${site}_${gcm}_${prd}
    mkdir -p $out
    
    if [ $prd == '1990s' ]; then
	styr=1980
	edyr=2009

    elif [ $prd == '2050s' ]; then
	styr=2040
	edyr=2069
    elif [ $prd == '2080s' ]; then
	styr=2070
	edyr=2099
    else
	python $script -dm --name ${site}_${gcm} $ifile $out
	return
    fi
        
    
    python $script -sle --lq 1 --exlyr 2,10 --start $styr --end $edyr --name ${pfx}_1day $ifile $out 
    python $script -sle --lq 7 --exlyr 2,10 --start $styr --end $edyr --name ${pfx}_7day $ifile $out
    python $script -spe --pq 1 --expyr 2,5,10,25 --start $styr --end $edyr --name ${pfx}_1day $ifile $out
    
    #python $script -dm --name ${pfx} $ifile $out
    

}
export -f process_site

for g in ${gcms[@]}; do
    for p in ${pds[@]}; do		
	echo $g $p
	find $ddir/$g/*.hour | parallel -j10 --eta process_site "{=s/.hour//=}" $g $p

	#exit
    done
done

exit

data = sorted(glob("{}/*".format(ddir)))
for d in data:
    gcm = os.path.basename(d)
    print(d)
    print("gcm: ", gcm)

    
for



    
# Get list of sites
# create streamflow data

#low flows 2,10 yr (1-day, 7-day)
#peak stats, 2,5,10,25, 100 yr (1-hour, 1-day)
#timeseries
#monthly flows
#time periods: 1990s, 2050s, 2080s

# creating summary table for streamflow statistics