else:
slashstr = '\\'
# add path containing readSettings.py
sys.path.append('..' + slashstr + 'Run')
# Read in user-specified settings
from readSettings import Settings
# Set input file
if len(sys.argv) < 2:
settings_input_file = 'settings.ini'
print 'Using default input file: ' + settings_input_file
else:
settings_input_file = sys.argv[1]
print 'Using specified input file: ' + settings_input_file
Settings = Settings(settings_input_file)
#%% *** SET THE FOLLOWING *****************************************************
# *** Enter list of HRU numbers to identify byhighlighting in plot
hru_to_ID = np.array([1, 23])
#%% *** CHANGE FILE NAMES AS NEEDED *******************************************
# (default is to use entries from Settings File)
projdir_GIS = Settings.GISinput_dir
HRUshp_fil = projdir_GIS + slashstr + "shapefiles/HRUs/HRUs.shp"
segshp_fil = projdir_GIS + slashstr + "shapefiles/segments/segments.shp"
print '\n******************************************'
print 'Plotting from HRU shapefile: ' + HRUshp_fil
print '******************************************\n'
else:
slashstr = '\\'
# add path containing readSettings.py
sys.path.append('..' + slashstr + 'Run')
# Read in user-specified settings
from readSettings import Settings
# Set input file
if len(sys.argv) < 2:
settings_input_file = 'settings.ini'
print 'Using default input file: ' + settings_input_file
else:
settings_input_file = sys.argv[1]
print 'Using specified input file: ' + settings_input_file
Settings = Settings(settings_input_file)
#%% ***SET THE FOLLOWING ******************************************************
sw_scheme = 2 # 1: based on elev, 2: based on streams, 3: based on both
hydcond_default = 0.1
if sw_scheme == 1:
elev_thresh = [4500, 4200, 4000] # elev thresholds for different K, high to low
hydcond_by_elev = [0.1, 0.25, 0.4, 0.5] # K for the different elevation intervals
if sw_scheme >= 2:
# settings for hydraulic conductivity based on distance (in pixels) from stream
# Shullcas example (glacial)