import arrayUtil as au import numpy as np nrow, ncol = 301, 501 ibnd_l1 = au.loadArrayFromFile(nrow, ncol, 'ref\\icbnd_l1.ref') #ibnd = au.loadArrayFromFile(nrow,ncol,'ref\\icbund.ref') au.plotArray(ibnd_l1, 500, 500) au.plotArray(ibnd, 500, 500)
import arrayUtil as au import numpy as np import pylab #--load hard data harddata_file = '..\\tbl29_pro.dat' title,harddata_names,harddata = gu.loadGslibFile(harddata_file) #print np.shape(harddata) nrow,ncol = 459,615 delr,delc = 500,500 offset = [668350.,580985.] max_elev - 15.0 top = au.loadArrayFromFile(nrow,ncol,'model_top.ref') top[np.where(top>max_elev)] = max_elev au.plotArray(top,delr,delc,gpts=harddata[:,0:2],offset=offset) # #l1_bot = top - au.loadArrayFromFile(nrow,ncol,'H_resample_500_500.ref') #l2_bot = l1_bot - au.loadArrayFromFile(nrow,ncol,'Q5_resample_500_500.ref') #l3_bot = l2_bot - au.loadArrayFromFile(nrow,ncol,'Q4_resample_500_500.ref') #l4_bot = l3_bot - au.loadArrayFromFile(nrow,ncol,'Q3_resample_500_500.ref') #l5_bot = l4_bot - au.loadArrayFromFile(nrow,ncol,'Q2_resample_500_500.ref') #l6_bot = l5_bot - au.loadArrayFromFile(nrow,ncol,'Q1_resample_500_500.ref') #l7_bot = l6_bot - au.loadArrayFromFile(nrow,ncol,'T3_resample_500_500.ref') #l8_bot = l7_bot - au.loadArrayFromFile(nrow,ncol,'T2_resample_500_500.ref') #l9_bot = l8_bot - au.loadArrayFromFile(nrow,ncol,'T1_resample_500_500.ref') #au.writeArrayToFile(top,'model_top_25ft.ref') #au.writeArrayToFile(l1_bot,'l1_bot.ref') #au.writeArrayToFile(l2_bot,'l2_bot.ref')
offset_new = [728600., 577350.]
filter = np.array([[1, 4, 7, 4, 1], [4, 16, 26, 16, 4], [7, 26, 41, 26, 7],
[4, 16, 26, 16, 4], [1, 4, 7, 4, 1]])
#print np.shape(filter),np.cumsum(filter)
min_elev, max_elev = 0.0, 10.0
nrow, ncol = 301, 501
iterations = 40
#--load hard data
print 'loading hard data...'
harddata_file = '..\\tbl29_pro.dat'
title, harddata_names, harddata = gu.loadGslibFile(harddata_file)
#print np.shape(harddata)
print 'loading top array...'
top = au.loadArrayFromFile(nrow, ncol, 'filter_14_edge.ref')
top[np.where(top > max_elev)] = max_elev
top[np.where(top < min_elev)] = min_elev
#print top
org_top = top
#au.plotArray(org_top,500,500,offset=offset,output='save',min=0,gpts=harddata[:,0:2],title='org')
new_top = np.zeros_like(top)
xi, yi = 2, 2
print 'beginning smoothing iterations...'
for iter in range(14, iterations):
print iter
while xi < ncol - 3:
while yi < nrow - 3:
nrow, ncol = 301, 501
sp = 100
fileNew = 'bro.chd'
f_out = open(fileNew, 'w')
offset = [728600., 577350.]
newMXBND = -1e+30
orgMXBND = -1e+30
#write_file_header()
cont = True
count = 0
ibnd_l1 = au.loadArrayFromFile(nrow, ncol, 'ref\\ibound_l1.ref')
ibnd = au.loadArrayFromFile(nrow, ncol, 'ref\\ibound.ref')
l1_idx = (np.argwhere(ibnd_l1 == -1)) + 1
idx = (np.argwhere(ibnd == -1)) + 1
print np.shape(l1_idx), np.shape(idx)
chd_stage = 0.0
newCells = []
for cell in range(0, np.shape(l1_idx)[0]):
#print l1_idx[cell,0],l1_idx[cell,1]
newCells.append(
mf.MODFLOW_HDF('HDF', 1, l1_idx[cell, 0], l1_idx[cell, 1],
[chd_stage, chd_stage]))
offset_new = [728600.,577350.]
filter = np.array([[1,4,7,4,1],[4,16,26,16,4],[7,26,41,26,7],[4,16,26,16,4],[1,4,7,4,1]])
#print np.shape(filter),np.cumsum(filter)
min_elev,max_elev = 0.0,10.0
nrow,ncol = 301,501
iterations = 40
#--load hard data
print 'loading hard data...'
harddata_file = '..\\tbl29_pro.dat'
title,harddata_names,harddata = gu.loadGslibFile(harddata_file)
#print np.shape(harddata)
print 'loading top array...'
top = au.loadArrayFromFile(nrow,ncol,'filter_14_edge.ref')
top[np.where(top>max_elev)] = max_elev
top[np.where(top<min_elev)] = min_elev
#print top
org_top = top
#au.plotArray(org_top,500,500,offset=offset,output='save',min=0,gpts=harddata[:,0:2],title='org')
new_top = np.zeros_like(top)
xi,yi = 2,2
print 'beginning smoothing iterations...'
for iter in range(14,iterations):
print iter
while xi < ncol-3:
while yi < nrow-3:
data.resize(nrow,ncol)
return(data)
d += 1
file_in.close()
data.resize(nrow,ncol)
return(data)
#--run fieldgen with an input file
os.system('fieldgen.exe <fieldgen.in')
#--get and plot the realizations
files = os.listdir('reals')
nlay,ncol = 50,151
for i,f in enumerate(files):
array = au.loadArrayFromFile(nlay,ncol,'reals\\'+f)
fig = pylab.figure(figsize=(15,5))
ax = pylab.axes((0.05,0.05,0.85,0.85))
cax = pylab.axes((0.9,0.05,0.025,0.85))
p = ax.pcolor(np.log10(np.flipud(array)))
#p = ax.pcolor(np.flipud(array/array.max()))
#p = ax.pcolor((np.flipud(poro)))
ax.set_xlim(0,ncol)
ax.set_ylim(0,nlay)
ax.set_xlabel('col')
ax.set_ylabel('row')
pylab.colorbar(p,cax=cax)
ax.set_title('log10 hydraulic conductivity realization '+str(i+1))
np.savetxt('reals\\'+f,array,fmt='%15.3e')
#break
pylab.show()
import arrayUtil as au
import ssmClass
reload(au)
reload(ssmClass)
nrow,ncol = 301,501
ibnd_l1 = au.loadArrayFromFile(nrow,ncol,'ibound_l1.ref')
idx = np.argwhere(ibnd_l1 == -1)
print np.shape(idx)
#f = open('simple.ssm','w')
cells1 = []
map1 = au.mapBndCellsArray(nrow,ncol,cells1,parm='concen')
map2 = au.mapBndCellsArray(nrow,ncol,cells2,parm='concen')
au.plotArray(map1,500,500,output=None)
au.plotArray(map2,500,500,output=None)
#f.close()
import gslibUtil as gu import arrayUtil as au import numpy as np import pylab #--load hard data harddata_file = '..\\tbl29_pro.dat' title, harddata_names, harddata = gu.loadGslibFile(harddata_file) #print np.shape(harddata) nrow, ncol = 301, 501 delr, delc = 500, 500 offset = [728600., 577350.] top = au.loadArrayFromFile(nrow, ncol, 'filter_28_edge.ref') l1_thk = au.loadArrayFromFile(nrow, ncol, 'Layer_thk\\H_resample_500_500.ref') l2_thk = au.loadArrayFromFile(nrow, ncol, 'Layer_thk\\Q5_resample_500_500.ref') l3_thk = au.loadArrayFromFile(nrow, ncol, 'Layer_thk\\Q4_resample_500_500.ref') l4_thk = au.loadArrayFromFile(nrow, ncol, 'Layer_thk\\Q3_resample_500_500.ref') l5_thk = au.loadArrayFromFile(nrow, ncol, 'Layer_thk\\Q2_resample_500_500.ref') l6_thk = au.loadArrayFromFile(nrow, ncol, 'Layer_thk\\Q1_resample_500_500.ref') l7_thk = au.loadArrayFromFile(nrow, ncol, 'Layer_thk\\T3_resample_500_500.ref') l8_thk = au.loadArrayFromFile(nrow, ncol, 'Layer_thk\\T2_resample_500_500.ref') l9_thk = au.loadArrayFromFile(nrow, ncol, 'Layer_thk\\T1_resample_500_500.ref') l1_bot = top - l1_thk l2_bot = l1_bot - l2_thk l3_bot = l2_bot - l3_thk l4_bot = l3_bot - l4_thk l5_bot = l4_bot - l5_thk
sp = 100
fileNew = 'bro.chd'
f_out = open(fileNew,'w')
offset = [728600.,577350.]
newMXBND = -1e+30
orgMXBND = -1e+30
#write_file_header()
cont = True
count = 0
ibnd_l1 = au.loadArrayFromFile(nrow,ncol,'ref\\ibound_l1.ref')
ibnd = au.loadArrayFromFile(nrow,ncol,'ref\\ibound.ref')
l1_idx = (np.argwhere(ibnd_l1==-1))+1
idx = (np.argwhere(ibnd==-1))+1
print np.shape(l1_idx),np.shape(idx)
chd_stage = 0.0
newCells = []
for cell in range(0,np.shape(l1_idx)[0]):
#print l1_idx[cell,0],l1_idx[cell,1]
newCells.append(mf.MODFLOW_HDF('HDF',1,l1_idx[cell,0],l1_idx[cell,1],[chd_stage,chd_stage]))
for ilay in range(2,9):
import arrayUtil as au import numpy as np nrow,ncol = 301,501 ibnd_l1 = au.loadArrayFromFile(nrow,ncol,'ref\\icbnd_l1.ref') #ibnd = au.loadArrayFromFile(nrow,ncol,'ref\\icbund.ref') au.plotArray(ibnd_l1,500,500) au.plotArray(ibnd,500,500)
import os,sys,pylab
import numpy as np
import gslibUtil as gu
import arrayUtil as au
import numpy as np
import pylab
#--load hard data
harddata_file = 'tbl29_pro.dat'
title,harddata_names,harddata = gu.loadGslibFile(harddata_file)
#print np.shape(harddata)
nrow,ncol = 197,116
delr,delc = 2650,2650
offset = [668350.,288415.]
path = 'array_sk\\'
files = os.listdir(path)
print files
for file in files:
this_array = au.loadArrayFromFile(nrow,ncol,path+file)
print np.mean(this_array)
au.plotArray(this_array,delr,delc,gpts=harddata[:,0:2],offset=offset)
import arrayUtil as au import numpy as np import pylab #--load hard data harddata_file = '..\\tbl29_pro.dat' title,harddata_names,harddata = gu.loadGslibFile(harddata_file) #print np.shape(harddata) nrow,ncol = 301,501 delr,delc = 500,500 offset = [728600.,577350.] top = au.loadArrayFromFile(nrow,ncol,'filter_28_edge.ref') l1_thk = au.loadArrayFromFile(nrow,ncol,'H_resample_500_500.ref') l2_thk = au.loadArrayFromFile(nrow,ncol,'Q5_resample_500_500.ref') l3_thk = au.loadArrayFromFile(nrow,ncol,'Q4_resample_500_500.ref') l4_thk = au.loadArrayFromFile(nrow,ncol,'Q3_resample_500_500.ref') l5_thk = au.loadArrayFromFile(nrow,ncol,'Q2_resample_500_500.ref') l6_thk = au.loadArrayFromFile(nrow,ncol,'Q1_resample_500_500.ref') l7_thk = au.loadArrayFromFile(nrow,ncol,'T3_resample_500_500.ref') l8_thk = au.loadArrayFromFile(nrow,ncol,'T2_resample_500_500.ref') l9_thk = au.loadArrayFromFile(nrow,ncol,'T1_resample_500_500.ref') l1_bot = top - l1_thk l2_bot = l1_bot - l2_thk l3_bot = l2_bot - l3_thk l4_bot = l3_bot - l4_thk l5_bot = l4_bot - l5_thk
for y1 in range(0,np.shape(map_y)[0]):
for y2 in range(1,np.shape(halfmile_y)[0]):
node_y = new_y[y1+1]-(new_delc/2)
if node_y <= halfmile_y[y2]:
map_y[y1] = y2-1
break
#print map_y.astype(int)
#--load hard data
harddata_file = 'tbl29_pro.dat'
title,harddata_names,harddata = gu.loadGslibFile(harddata_file)
#print np.shape(harddata)
files = os.listdir('array_sk\\')
for file in files:
print file
new_array = np.zeros((new_nrow,new_ncol),dtype='float')+1.0e+32
org_array = np.flipud(au.loadArrayFromFile(hm_nrow,hm_ncol,'array_sk\\'+file))
for row in range(0,new_nrow):
for col in range(0,new_ncol):
new_array[row,col] = org_array[map_y[row],map_x[col]]
title = file.split('_')[0]+'_org.ref'
au.plotArray(np.flipud(org_array),hm_delr,hm_delc,output='save',offset=offset,title=title,gpts=harddata[:,0:2])
au.plotArray(np.flipud(new_array),new_delr,new_delc,offset=offset_new,output='save',title=title,gpts=harddata[:,0:2])
# au.writeArrayToFile(np.flipud(new_array),title)
# pylab.show()
import arrayUtil as au
import ssmClass
reload(au)
reload(ssmClass)
nrow, ncol = 301, 501
ibnd_l1 = au.loadArrayFromFile(nrow, ncol, 'ibound_l1.ref')
idx = np.argwhere(ibnd_l1 == -1)
print np.shape(idx)
#f = open('simple.ssm','w')
cells1 = []
map1 = au.mapBndCellsArray(nrow, ncol, cells1, parm='concen')
map2 = au.mapBndCellsArray(nrow, ncol, cells2, parm='concen')
au.plotArray(map1, 500, 500, output=None)
au.plotArray(map2, 500, 500, output=None)
#f.close()
import arrayUtil as au import numpy as np import pylab #--load hard data harddata_file = '..\\tbl29_pro.dat' title,harddata_names,harddata = gu.loadGslibFile(harddata_file) #print np.shape(harddata) nrow,ncol = 301,501 delr,delc = 500,500 offset = [728600.,577350.] top = au.loadArrayFromFile(nrow,ncol,'filter_28_edge.ref') l1_thk = au.loadArrayFromFile(nrow,ncol,'Layer_thk\\H_resample_500_500.ref') l2_thk = au.loadArrayFromFile(nrow,ncol,'Layer_thk\\Q5_resample_500_500.ref') l3_thk = au.loadArrayFromFile(nrow,ncol,'Layer_thk\\Q4_resample_500_500.ref') l4_thk = au.loadArrayFromFile(nrow,ncol,'Layer_thk\\Q3_resample_500_500.ref') l5_thk = au.loadArrayFromFile(nrow,ncol,'Layer_thk\\Q2_resample_500_500.ref') l6_thk = au.loadArrayFromFile(nrow,ncol,'Layer_thk\\Q1_resample_500_500.ref') l7_thk = au.loadArrayFromFile(nrow,ncol,'Layer_thk\\T3_resample_500_500.ref') l8_thk = au.loadArrayFromFile(nrow,ncol,'Layer_thk\\T2_resample_500_500.ref') l9_thk = au.loadArrayFromFile(nrow,ncol,'Layer_thk\\T1_resample_500_500.ref') l1_bot = top - l1_thk l2_bot = l1_bot - l2_thk l3_bot = l2_bot - l3_thk l4_bot = l3_bot - l4_thk l5_bot = l4_bot - l5_thk
import gslibUtil as gu import arrayUtil as au import numpy as np import pylab #--load hard data harddata_file = '..\\tbl29_pro.dat' title, harddata_names, harddata = gu.loadGslibFile(harddata_file) #print np.shape(harddata) nrow, ncol = 301, 501 delr, delc = 500, 500 offset = [728600., 577350.] top = au.loadArrayFromFile(nrow, ncol, 'filter_28_edge.ref') l1_thk = au.loadArrayFromFile(nrow, ncol, 'H_resample_500_500.ref') l2_thk = au.loadArrayFromFile(nrow, ncol, 'Q5_resample_500_500.ref') l3_thk = au.loadArrayFromFile(nrow, ncol, 'Q4_resample_500_500.ref') l4_thk = au.loadArrayFromFile(nrow, ncol, 'Q3_resample_500_500.ref') l5_thk = au.loadArrayFromFile(nrow, ncol, 'Q2_resample_500_500.ref') l6_thk = au.loadArrayFromFile(nrow, ncol, 'Q1_resample_500_500.ref') l7_thk = au.loadArrayFromFile(nrow, ncol, 'T3_resample_500_500.ref') l8_thk = au.loadArrayFromFile(nrow, ncol, 'T2_resample_500_500.ref') l9_thk = au.loadArrayFromFile(nrow, ncol, 'T1_resample_500_500.ref') l1_bot = top - l1_thk l2_bot = l1_bot - l2_thk l3_bot = l2_bot - l3_thk l4_bot = l3_bot - l4_thk l5_bot = l4_bot - l5_thk
