#fn.save_figures_minerals(carb_rt, max_pqty, time_points[j], path, nn, ptype=m) #fn.save_figures_mols(carb_rt, time_points[j], path, nn, ptype=m) j +=1 carb_rt.advance() results = fn.append_results(carb_rt, results, step = S ) itr += 1 #%% SIMULATION TIME simulation_time = time.time()-it fn.print_time(simulation_time, carb_rt) #%% SAVE if(False): fn.save_obj(results, path + str(nn) +'_results') np.save(path + 'SI', carb_rt.phrqc.selected_output()['SI_calcite'] ) np.save(path + 'pH', carb_rt.phrqc.selected_output()['pH'] ) np.save(path + 'Ca', carb_rt.phrqc.selected_output()['Ca'] ) np.save(path + 'C', carb_rt.phrqc.selected_output()['C'] ) np.save(path + 'De', carb_rt.fluid.Ca.De ) np.save(path + 'poros', carb_rt.fluid.Ca.poros) #%% PLOT if(True): fn.plot_species(results, names=[])#['calcite'] fn.plot_avg(results, names=['avg_poros', 'avg_D_eff']) fn.plot_points(results, names=['calcite', 'portlandite', 'poros', 'Ca', 'C']) fn.plot_fields(carb_rt, names=['calcite', 'portlandite', 'Ca', 'C', 'poros'],fsize=(15,1)) #%% PRINT if(True):
print('C +ss %s' %str(np.array(csh.fluid.C.c[1,:]) + np.array(csh.fluid.C._ss[1,:])/np.array(csh.phrqc.poros[1,:]))) print('H +ss %s' %str(np.array(csh.fluid.H.c[1,:]) + np.array(csh.fluid.H._ss[1,:])/np.array(csh.phrqc.poros[1,:]))) print('O +ss %s' %str(np.array(csh.fluid.O.c[1,:]) + np.array(csh.fluid.O._ss[1,:])/np.array(csh.phrqc.poros[1,:]))) print('CH %s' %str(np.array(csh.solid.portlandite.c[1,:]))) print('dCH %s' %str(np.array(csh.phrqc.dphases['portlandite'][1,:]))) print('CC %s' %str(np.array(csh.solid.calcite.c[1,:]))) print('dCC %s' %str(np.array(csh.phrqc.dphases['calcite'][1,:]))) print('Vol %s' %str(np.array(csh.solid.vol[1,:]))) print('D %s' %str(np.array(csh.fluid.C.De[1,:]))) print('SI %s' %str(np.array(csh.solid.target_SI[1,:]))) print('pH %s' %str(np.array(csh.phrqc.selected_output()['pH'][1,:]))) print('poros %s' %str(np.array(csh.solid.poros[1,:]))) print('phrqc poros %s' %str(np.array(np.array(csh.phrqc.poros[1,:])))) print('CSHQ TobD %s' %str(np.array(csh.solid.CSHQ_TobD.c[1,:]))) print('CSHQ TobH %s' %str(np.array(csh.solid.CSHQ_TobH.c[1,:]))) print('CSHQ JenD %s' %str(np.array(csh.solid.CSHQ_JenD.c[1,:]))) print('CSHQ JenH %s' %str(np.array(csh.solid.CSHQ_JenH.c[1,:]))) print('SiO2 %s' %str(np.array(csh.solid.SiO2am.c[1,:]))) #''' #%% SAVE ''' fn.save_obj(results, path + str(nn) +'_results') np.save(path + 'SI', csh.phrqc.selected_output()['SI_calcite'] ) np.save(path + 'pH', csh.phrqc.selected_output()['pH'] ) np.save(path + 'Ca', csh.phrqc.selected_output()['Ca'] ) np.save(path + 'C', csh.phrqc.selected_output()['C'] ) np.save(path + 'CH', csh.phrqc.selected_output()['portlandite'] ) np.save(path + 'CC', csh.phrqc.selected_output()['calcite'] ) np.save(path + 'De', csh.fluid.Ca.De ) np.save(path + 'poros', csh.fluid.Ca.poros) '''