def __init__(self, domain, z0=None, k0=None, dk=None, beta=None): """ set up mapping :param domain: domain of the mapping :type domain: ``Domain`` :param z0: depth weighting offset. If not present no depth scaling is applied. :type z0: scalar :param k0: reference density, defaults to 0 :type k0: scalar :param dk: susceptibility scale, defaults to 1 :type dk: scalar :param beta: depth weighting exponent, defaults to 2 :type beta: ``float`` """ if k0 is None: k0 = 0.0 if beta is None: beta = 2.0 if dk is None: dk = 1.0 self.domain = domain if z0 is not None: DIM = self.domain.getDim() l_z = boundingBoxEdgeLengths(domain)[DIM - 1] a = dk * (clip(z0 - domain.getX()[DIM - 1], minval=0) / l_z) ** (beta / 2) else: a = dk super(SusceptibilityMapping, self).__init__(a=a, p0=k0)
def __init__(self, domain, z0=None, k0=None, dk=None, beta=None): """ set up mapping :param domain: domain of the mapping :type domain: ``Domain`` :param z0: depth weighting offset. If not present no depth scaling is applied. :type z0: scalar :param k0: reference density, defaults to 0 :type k0: scalar :param dk: susceptibility scale, defaults to 1 :type dk: scalar :param beta: depth weighting exponent, defaults to 2 :type beta: ``float`` """ if k0 is None: k0 = 0. if beta is None: beta = 2. if dk is None: dk = 1. self.domain = domain if z0 is not None: DIM = self.domain.getDim() l_z = boundingBoxEdgeLengths(domain)[DIM - 1] a = dk * (clip(z0 - domain.getX()[DIM - 1], minval=0) / l_z)**(beta / 2) else: a = dk super(SusceptibilityMapping, self).__init__(a=a, p0=k0)
def __init__(self, domain, z0=None, rho0=None, drho=None, beta=None): """ initializes the mapping :param domain: domain of the mapping :type domain: ``Domain`` :param z0: depth weighting offset. If not present no depth scaling is applied. :type z0: scalar :param rho0: reference density, defaults to 0 :type rho0: scalar :param drho: density scale. By default density of granite = 2750kg/m**3 is used. :type drho: scalar :param beta: depth weighting exponent, defaults to 2 :type beta: ``float`` """ if rho0 is None: rho0 = 0.0 if drho is None: drho = 2750 * U.kg / U.m ** 3 if beta is None: beta = 2.0 self.domain = domain if z0 is not None: DIM = self.domain.getDim() l_z = boundingBoxEdgeLengths(domain)[DIM - 1] a = drho * (clip(z0 - domain.getX()[DIM - 1], minval=0) / l_z) ** (beta / 2) else: a = drho super(DensityMapping, self).__init__(a=a, p0=rho0)
def __init__(self, domain, z0=None, rho0=None, drho=None, beta=None): """ initializes the mapping :param domain: domain of the mapping :type domain: ``Domain`` :param z0: depth weighting offset. If not present no depth scaling is applied. :type z0: scalar :param rho0: reference density, defaults to 0 :type rho0: scalar :param drho: density scale. By default density of granite = 2750kg/m**3 is used. :type drho: scalar :param beta: depth weighting exponent, defaults to 2 :type beta: ``float`` """ if rho0 is None: rho0 = 0. if drho is None: drho = 2750 * U.kg / U.m**3 if beta is None: beta = 2. self.domain = domain if z0 is not None: DIM = self.domain.getDim() l_z = boundingBoxEdgeLengths(domain)[DIM - 1] a = drho * (clip(z0 - domain.getX()[DIM - 1], minval=0) / l_z)**(beta / 2) else: a = drho super(DensityMapping, self).__init__(a=a, p0=rho0)