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