#nstp = [] for td in flow.sp_len: perlen.append(td.days) # nstp.append(td.days) dis = mfdis(mf,flow.nrow,flow.ncol,flow.nlay,flow.nper,nstp=nstp,delr=500,delc=500,top=flow.ref_dir+'top_mod.ref',botm=botm,perlen=perlen,steady=False,laycbd=0) bas = mfbas(mf,ibound=ibound,strt=strt,hnoflo=1.0e+10) #upw = mfupw(mf,laytyp=0,hk=4500.0,vka=4500.0) #nwt = mfnwt(mf,headtol=0.15,fluxtol=100000.0,maxiterout=500) lpf = mflpf(mf,laytyp=0,hk=4500.0,vka=4500.0) gmg = mfgmg(mf,mxiter=100,hclose=0.3,rclose=3000.0) oc = mfoc(mf,words=['head','budget'],save_head_every=1) #--externally generated... ghb = mfaddoutsidefile(mf,'GHB','ghb',125) ets = mfaddoutsidefile(mf,'ETS','ets',126) rch = mfaddoutsidefile(mf,'RCH','rch',127) #mnw = mfaddoutsidefile(mf,'MNW2','mnw',128) wel = mfaddoutsidefile(mf,'WEL','wel',128) #--swi swi = mfaddoutsidefile(mf,'SWI','swi',129) mf.add_external(modelname+'.zta',130,binflag=True) #--add additional CBC units mf.add_external(flow.root+'_wel.cbc',flow.well_unit,binflag=True) mf.add_external(flow.root+'_rch.cbc',flow.rch_unit,binflag=True) mf.add_external(flow.root+'_ets.cbc',flow.ets_unit,binflag=True) mf.add_external(flow.root+'_ghb.cbc',flow.ghb_unit,binflag=True)
bas = mfbas(mf, ibound, init_heads) #--set ghbs along the coast ghb_list = [] ghb_stage = 0.0 ghb_cond = 10000000000.0 for i in range(nrow): for j in range(ncol): if ibound[i, j] == 34: for k in range(nlay): ghb_list.append([k + 1, i + 1, j + 1, ghb_stage, ghb_cond]) lpf = mflpf(mf, hk=[10.0, 0.1, 500.], vka=[10.0, 0.1, 500.], laytyp=1) gmg = mfgmg(mf, mxiter=1000, hclose=1e-2, rclose=1e-2) oc = mfoc(mf, words=['head', 'budget'], save_head_every=1) ghb = mfghb(mf, layer_row_column_head_cond=[ghb_list]) rch = mfrch(mf, rech=0.000114, external=False) swr = mfaddoutsidefile(mf, 'SWR', 'swr', 25) mf.add_external('tsala.fls', 101, binflag=True) mf.add_external('tsala.stg', 102, binflag=True) mf.add_external('tsala.aqx', 103, binflag=True) mf.add_external('tsala.pqm', 104, binflag=True) mf.add_external('tsala.str', 105) mf.add_external('tsala.riv', 106) mf.write_input() os.system('MF2005-SWR_x64.exe tsala.nam')
bas = mfbas(mf,ibound,init_heads) #--set ghbs along the coast ghb_list = [] ghb_stage = 0.0 ghb_cond = 10000000000.0 for i in range(nrow): for j in range(ncol): if ibound[i,j] == 34: for k in range(nlay): ghb_list.append([k+1,i+1,j+1,ghb_stage,ghb_cond]) lpf = mflpf(mf,hk=[10.0,0.1,500.],vka=[10.0,0.1,500.],laytyp=1) gmg = mfgmg(mf,mxiter=1000,hclose=1e-2,rclose=1e-2) oc = mfoc(mf,words=['head','budget'],save_head_every=1) ghb = mfghb(mf,layer_row_column_head_cond=[ghb_list]) rch = mfrch(mf,rech=0.000114,external=False) swr = mfaddoutsidefile(mf,'SWR','swr',25) mf.add_external('tsala.fls',101,binflag=True) mf.add_external('tsala.stg',102,binflag=True) mf.add_external('tsala.aqx',103,binflag=True) mf.add_external('tsala.pqm',104,binflag=True) mf.add_external('tsala.str',105) mf.add_external('tsala.riv',106) mf.write_input() os.system('MF2005-SWR_x64.exe tsala.nam')
delc=500, top=seawat.ref_dir + top, botm=botm, perlen=perlen, steady=False, laycbd=0) bas = mfbas(mf, ibound=ibound, strt=strt, hnoflo=1.0e+10) #upw = mfupw(mf,laytyp=0,hk=4500.0,vka=4500.0) #nwt = mfnwt(mf,headtol=0.15,fluxtol=100000.0,maxiterout=500) lpf = mflpf(mf, laytyp=0, hk=hk, vka=hk, ss=ss, ilpfcb=0) gmg = mfgmg(mf, mxiter=100, hclose=0.75, rclose=100000.0) oc = mfoc(mf, words=['head', 'budget'], save_head_every=1, compact=True) #--externally generated... ghb = mfaddoutsidefile(mf, 'GHB', 'ghb', 125) shutil.copy(flow.root + '.ets', seawat.root + '.ets') ets = mfaddoutsidefile(mf, 'ETS', 'ets', 126) shutil.copy(flow.root + '.rch', seawat.root + '.rch') rch = mfaddoutsidefile(mf, 'RCH', 'rch', 127) wel = mfaddoutsidefile(mf, 'WEL', 'wel', 128) shutil.copy(flow.root + '.riv', seawat.root + '.riv') riv = mfaddoutsidefile(mf, 'RIV', 'riv', 124) #--add additional CBC units mf.add_external(seawat.root + '_wel.cbc', seawat.well_unit, binflag=True) mf.add_external(seawat.root + '_rch.cbc', seawat.rch_unit, binflag=True) mf.add_external(seawat.root + '_ets.cbc', seawat.ets_unit, binflag=True) mf.add_external(seawat.root + '_ghb.cbc', seawat.ghb_unit, binflag=True)
delc=500, top=flow.ref_dir + 'top_mod.ref', botm=botm, perlen=perlen, steady=False, laycbd=0) bas = mfbas(mf, ibound=ibound, strt=strt, hnoflo=1.0e+10) #upw = mfupw(mf,laytyp=0,hk=4500.0,vka=4500.0) #nwt = mfnwt(mf,headtol=0.15,fluxtol=100000.0,maxiterout=500) lpf = mflpf(mf, laytyp=0, hk=hk, vka=hk) gmg = mfgmg(mf, mxiter=100, hclose=0.25, rclose=1000.0) oc = mfoc(mf, words=['head', 'budget'], save_head_every=1) #--externally generated... ghb = mfaddoutsidefile(mf, 'GHB', 'ghb', 125) ets = mfaddoutsidefile(mf, 'ETS', 'ets', 126) rch = mfaddoutsidefile(mf, 'RCH', 'rch', 127) #mnw = mfaddoutsidefile(mf,'MNW2','mnw',128) wel = mfaddoutsidefile(mf, 'WEL', 'wel', 128) #--swi shutil.copy('..\\..\\_model\\bro.03\\flow.swi', flow.root + '.swi') swi = mfaddoutsidefile(mf, 'SWI', 'swi', 129) mf.add_external(modelname + '.zta', 130, binflag=True) #--add additional CBC units mf.add_external(flow.root + '_wel.cbc', flow.well_unit, binflag=True) mf.add_external(flow.root + '_rch.cbc', flow.rch_unit, binflag=True) mf.add_external(flow.root + '_ets.cbc', flow.ets_unit, binflag=True) mf.add_external(flow.root + '_ghb.cbc', flow.ghb_unit, binflag=True)
dis = mfdis(mf,seawat.nrow,seawat.ncol,seawat.nlay,seawat.nper,nstp=nstp,delr=500,delc=500,top=seawat.ref_dir+'top_mod.ref',botm=botm,perlen=perlen,steady=False,laycbd=0) bas = mfbas(mf,ibound=ibound,strt=strt,hnoflo=1.0e+10) #upw = mfupw(mf,laytyp=0,hk=4500.0,vka=4500.0) #nwt = mfnwt(mf,headtol=0.15,fluxtol=100000.0,maxiterout=500) #hk = [] #for lay in seawat.layer_botm_names: # hk.append(seawat.ref_dir+lay+'_k.ref') lpf = mflpf(mf,laytyp=0,hk=seawat.layer_k,vka=seawat.layer_k,ilpfcb=0) gmg = mfgmg(mf,mxiter=100,hclose=0.25,rclose=10000.0) oc = mfoc(mf,words=['head','budget'],save_head_every=1,compact=True) #--externally generated... ghb = mfaddoutsidefile(mf,'GHB','ghb',125) shutil.copy(flow.root+'.ets',seawat.root+'.ets') ets = mfaddoutsidefile(mf,'ETS','ets',126) shutil.copy(flow.root+'.rch',seawat.root+'.rch') rch = mfaddoutsidefile(mf,'RCH','rch',127) wel = mfaddoutsidefile(mf,'WEL','wel',128) shutil.copy(flow.root+'.riv',seawat.root+'.riv') riv = mfaddoutsidefile(mf,'RIV','riv',124) #--add additional CBC units mf.add_external(seawat.root+'_wel.cbc',seawat.well_unit,binflag=True) mf.add_external(seawat.root+'_rch.cbc',seawat.rch_unit,binflag=True) mf.add_external(seawat.root+'_ets.cbc',seawat.ets_unit,binflag=True) mf.add_external(seawat.root+'_ghb.cbc',seawat.ghb_unit,binflag=True)