#from matmod6.volume_frac import G0_func
from volume_frac import G0_func

def grain_size_site_saturation(N=1.6 * 10**14, T = 300, M0 = 1e4, Q = 156000, P_D = 1e5):
    D = (1 / N) ** (1 / 3)  # Grain size
    return D

def grain_size_Johnson_Mehl(N=1.6 * 10**14, T = 300, M0 = 1e4, Q = 156000, P_D = 1e5):
    G0 = G0_func(M0=M0, Q=Q, P_D=P_D, T=T)
    D = (G0 / N)**(1/4)
    return D

def grain_size_variable_growth_rate(N=1.6 * 10**14, T = 300, M0 = 1e4, Q = 156000, P_D = 1e5):
    return grain_size_site_saturation(N=1.6 * 10**14, T = 300, M0 = 1e4, Q = 156000, P_D = 1e5)




#NB! har endret navn på alle funksjonene, lagt til grain_size_