#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)