Skip to content

Commit

Permalink
predict conductance functionality
Browse files Browse the repository at this point in the history
  • Loading branch information
DSeiferth committed Apr 9, 2024
1 parent 3ca6f54 commit f9c7087
Show file tree
Hide file tree
Showing 2 changed files with 3 additions and 3 deletions.
2 changes: 1 addition & 1 deletion PoreAnalyser/conductance.py
Original file line number Diff line number Diff line change
Expand Up @@ -48,7 +48,7 @@ def bullk_conduct(z,a,b, conduct=bulk_conductivity):
conduct = 1/R # * bulk_conductivity
return conduct, R_vec

def no_bulk_conduct(z,a,b, popt, conduct=bulk_conductivity, plot=True, shift_bond=1.85):
def no_bulk_conduct(z,a,b, popt, conduct=bulk_conductivity, plot=False, shift_bond=1.85):
'''
Calculate the conductance of a pore using a conductivity model
'''
Expand Down
4 changes: 2 additions & 2 deletions PoreAnalyser/poreanalyser.py
Original file line number Diff line number Diff line change
Expand Up @@ -311,11 +311,11 @@ def conductance_estimation(self, index_model=0, f_size=15):
ax[0].set_title('Conductivity along pore axis', fontsize=f_size, loc='center')
ax[1].set_title('Resistance along pore axis', fontsize=f_size, loc='center')
ax[0].set_xlabel("z ($\AA$)", fontsize=f_size)
ax[0].set_ylabel(r"$\kappa$(a,b)/$\kappa_{bulk}$", fontsize=f_size)
ax[0].set_ylabel(r"$\kappa$(a,b)/\kappa_{bulk}$", fontsize=f_size)
ax[1].set_xlabel("z ($\AA$)", fontsize=f_size)
ax[1].set_ylabel(r"Resistance ($G\Omega$)", fontsize=f_size)
ax[0].tick_params(axis='both', which='major', labelsize=f_size)
ax[1].tick_params(axis='both', which='major', labelsize=f_size)
fig.tight_layout()
#fig.tight_layout()
plt.show()
return hole1/pS, pf1/pS, hole_c/pS, pf1_c/pS, fig

0 comments on commit f9c7087

Please sign in to comment.