#!/bin/bash


#-----------------------------------------------------------------------------------------
# Set run conditions
#-----------------------------------------------------------------------------------------
#rundhsvm=true      # Runs DHSVM
#RBM_script=true    # Prepares and compiles RBM related codes
prepareRBM=true    # Prepares dir, segmap, and net files for RBM run
forcing=true       # Prepare forcing files for RBM run
#remakeForcing=true # Force overwrite existing forcing if set to true. Otherwise use existing
runRBM=true        # Runs RBM
extract=true       # Extracts temp data
tplot=true         # Create time series plots


#extract_all=true   # Extract temp data by rbm id
#getTribIDs=true    # Grabs all DHSVM IDs between given start and end DHSVM ID
#combine=true       # Combine temp data for all rbm id in a tributary
#dplot=true         # Create distance plots
#dplot_ymax=true    # Create distance plots based on annual median
#dplot_dmax=true    # Create distance plots based on daily median
#dplot_wk_ymax=true # Create distance plots based on weekly median per each year
#dplot_wk_dmax=true # Create distance plots based on weekly median per each week
#cplot=true         # Create cold patches plots

#-----------------------------------------------------------------------------------------
# Scenario conditions
#-----------------------------------------------------------------------------------------
basin=kitsap
project=kitsap_40m75k_run84
scen=base

#-----------------------------------------------------------------------------------------
# Directory paths
#-----------------------------------------------------------------------------------------
home_dir=/home/disk/tsuga2/jswon11/workdir/2019_03_kitsap-dhsvm/hyak/data_atmos/
basin_dir=/home/disk/tsuga2/jswon11/workdir/2019_03_kitsap-dhsvm/model/kitsap/
ref_dir=${home_dir}/data_ref/
stream_map=${basin_dir}/input/stream.map_40m-75k
stream_net=${basin_dir}/input/stream.network_40m-75k

#-----------------------------------------------------------------------------------------
# RBM Conditions
#-----------------------------------------------------------------------------------------
rvegfile=NA
origfilename=NA
chgheadwater=NA
chgDmin=NA

# Mohseni Params
f_smooth=0.1 #0.1
alpha=5.8  #5.8
beta=29  #29
gamma=0.49  #0.49
temp_min=1  # 1

# Leopold Params
D_a=0.2 # 0.2
D_b=0.4 # 0.4
D_min=0.5  # 0.5
U_a=0.9  # 0.9
U_b=0.21 # 0.21
U_min=0.5 #0.5


#-----------------------------------------------------------------------------------------
# DHSVM Conditions
#-----------------------------------------------------------------------------------------

dhsvmstyr=NA
dhsvmendyr=NA
dhsvmout=NA
csvstdate=NA
csvstyear=NA
tplotstyear=NA
tplotendyear=NA

#-----------------------------------------------------------------------------------------
# Plot conditions
#-----------------------------------------------------------------------------------------
plotday=NA
stsummer=NA
endsummer=NA
tirdate=NA
obswtemp=NA




sites=183,196,216,266,276,277
#Kitsap_Lk_outlet,183
#Dickr,196
#Wildc,216
#Lost,266
#ChiGC,276
#ChicoMS,277




#-----------------------------------------------------------------------------------------
#-----------------------------------------------------------------------------------------
#-----------------------------------------------------------------------------------------


# RBM - forcing file, Moh param, Leopold Param, convergence, stream connectivity


##-------------------------------------------------------------------
# Prepare RBM
#home_dir -> user supplied
#code_dir = ${home_dir}/code/


# make_stream_connectivity <mapfile> <networkfile> <output directory> <NSeg> <skip>
# code/Scripts/make_stream_connectivity
# <mapfile> - stream.map file : basin_dir/input/stream.map
# <network> - stream.network file: basin_dir/input/stream.network.dat
# <outdir> - output directory location
# <no_segs> - -> user supplied
# <skip> - number of lines in the header of stream map file : 9
# output: convergence.txt

# Convergence.txt
# <segment id> <destination segment id> <length-meters> <depth-meters> <avg azimuth> <upstream segment id>

# build_DHSVM_network.pl
#"${filename}\n ${f_smooth} \n ${alpha}, ${beta}, ${gamma}, ${temp_min} \n ${d_a}, ${d_b}, ${depth_min} \n ${U_a}, ${U_b}, ${U_min}\n"

# write_Moh_params.bsh
# write_Moh_params.bsh <filename> <f_smooth> <alpha> <beta> <gamma> <temp_min> <basin> 
# filename -
# f_smooth
# alpha
# beta
# gamma
# temp_min
# basin


# write_Leo_params.bsh
# write_Leo_params.bsh <filename> <d_a> <d_b> <depth_min> <U_a> <U_b> <U_min> <basin> 
# filename -
# d_a
# d_b
# depth_min
# U_a
# U_b
# U_min
# basin


##-------------------------------------------------------------------
# Create forcing file
# Create_file <Directory with forcing files (*.Only)> <Project Name>
# InDirectory -> Directory with forcing files
# Project -> Project Name
# Required - Seg map
# <project>.segmap --> Made in Prepare RBM
# Required - Forcings
# ATP, NLW, NSW, VP, WND, Inflow, Outflow --> Made from DHSVM


# workdir -> user supplied
# work_dir = ${workdir}/Work_space


##-------------------------------------------------------------------
# Run RBM
# basin -> user supplied
# scen -> user supplied
# tp -> user supplied
# filename = ${basin}_#{scen}_${tp}

# RBM ${work_dir}/${filename}



##-------------------------------------------------------------------
##-------------------------------------------------------------------
##-------------------------------------------------------------------
##-------------------------------------------------------------------
##-------------------------------------------------------------------
##-------------------------------------------------------------------
##-------------------------------------------------------------------
##-------------------------------------------------------------------