shape_dir = 'shapes\\'
#sys.exit()
#--heads
ibound = np.loadtxt('ref\ibound.ref')
#--loop over each month, using the last day in the month
tot_days = 0
for m in range(1,13):
days = calendar.mdays[m]
tot_days += days
hds_handle = mfb.MODFLOW_Head(nlay,nrow,ncol,results+'bro_6lay.hds')
totim,kstp,kper,h,success = hds_handle.get_record(tot_days)
print 'day',tot_days,' processing for month ',m
for l in range(nlay):
this_head = h[l,:,:]
this_head[np.where(ibound==0)] = -9999.0
au.ref2grd('shapes\\head_layer'+str(l+1)+'_month'+str(m)+'.txt',\
h[l,:,:],nrow,ncol,offset,500.0)
#sys.exit()
#--heads
ibound = np.loadtxt('ref\ibound.ref')
hds_handle = mfb.MODFLOW_Head(nlay,nrow,ncol,results+'bro_6lay.hds')
totim,kstp,kper,h,success = hds_handle.get_record()
conc_handle = mfb.MT3D_Concentration(nlay,nrow,ncol,'MT3D001.UCN')
totim_c,kstp_c,kper_c,c,success = conc_handle.get_record()
#--write concentrations at the end of the model run
for l in range(nlay):
au.ref2grd('shapes\\conc_layer'+str(l+1)+'.txt',c[l,:,:],nrow,ncol,offset,500.0)
#--loop over each month, using the last day in the month
#tot_days = 0
#for m in range(1,13):
#
# days = calendar.mdays[m]
# tot_days += days
#
# hds_handle = mfb.MODFLOW_Head(nlay,nrow,ncol,results+'bro_7lay.hds')
# totim,kstp,kper,h,success = hds_handle.get_record(tot_days)
#
# conc_handle = mfb.MT3D_Concentration(nlay,nrow,ncol,'MT3D001.UCN')
# totim_c,kstp_c,kper_c,c,success = conc_handle.get_record(tot_days)
#
