import pandas as pd
from glob import glob
import matplotlib as plt
import seaborn as sns
import os

ddir = '/home/disk/rocinante/DATA/temp/chico/hyak/data_atmos/'
odir = '/home/disk/rocinante/DATA/temp/chico/hyak/writeup/'
rdir = '/home/disk/rocinante/DATA/temp/chico/hyak/data_ref/stream/'
pdir = '/home/disk/rocinante/DATA/temp/chico/hyak/writeup/plots/'
tbl = ''

sites={
    'CHICOMAIN':
    'CHICO':
    'CHICOabvDKRSN':
    'CHICOblwDKRSN':
    'DICKERSON':
    'KITSAP':
    'LOST':
    'OUTLET':
    'WILDCAT':

}

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

python $script --lq 1 --pfreq 2,10 --styr styr --edyr edyr


    
# 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