# settings for hydraulic conductivity based on distance (in pixels) from stream
# Shullcas example (glacial)
npix_stream_buffer = 1 # number of different buffers, based on pixels from stream
buffer_dist_pix = np.arange(1,npix_stream_buffer+1) # up to buffer_dist_pix pixels from stream
buffer_hydcond = np.array([0.4, 0.2]) # will use the first npix_stream_buffer values
strm_hydcond = 0.6 # for stream pixels
# *****************************************************************************
# %%
# Only run this script if using spatially distributed K
fl_runscript = 0
if len(Settings.hydcond0) == 1:
if type(Settings.hydcond0[0]) is str:
Settings.hydcond0 = Settings.hydcond0[0]
fl_runscript = 1
if fl_runscript == 1:
surfz_fil = Settings.GISinput_dir + slashstr + 'DEM.asc'
NLAY = Settings.NLAY
f = open(surfz_fil, 'r')
sdata = {}
for i in range(6):
line = f.readline()
line = line.rstrip() # remove newline characters
key, value = line.split(': ')
try:
value = int(value)
except:
value = float(value)
# settings for hydraulic conductivity based on distance (in pixels) from stream
# Shullcas example (glacial)
npix_stream_buffer = 1 # number of different buffers, based on pixels from stream
buffer_dist_pix = np.arange(1,npix_stream_buffer+1) # up to buffer_dist_pix pixels from stream
buffer_hydcond = np.array([0.4, 0.2]) # will use the first npix_stream_buffer values
strm_hydcond = 0.6 # for stream pixels
# *****************************************************************************
# %%
# Only run this script if using spatially distributed K
fl_runscript = 0
if len(Settings.hydcond0) == 1:
if type(Settings.hydcond0[0]) is str:
Settings.hydcond0 = settings.hydcond0[0]
fl_runscript = 1
if fl_runscript == 1:
surfz_fil = Settings.GISinput_dir + slashstr + 'DEM.asc'
NLAY = Settings.NLAY
f = open(surfz_fil, 'r')
sdata = {}
for i in range(6):
line = f.readline()
line = line.rstrip() # remove newline characters
key, value = line.split(': ')
try:
value = int(value)
except:
value = float(value)