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