def mastinEVEM(HT, DRE = 2500):
    """Calculating erupted volume based on plume height.
    Inputs: Plume height (HT) in km
    DRE = dense-rock equivalent (kg/m^3)
    Outputs: Eruptive Volume (km^3 DRE) and Eruptive Mass (kg)
    No time component to this??
    Equation from Mastin et al. 2009: H = 25.9 + 6.64log(V) --> log base 10"""
    p = (HT - 25.9)/6.64
    EV = 10**p
    EM = EV * DRE * (1000**3)
    return EV, EM