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