Skip to content

array.py

Arrays dataclass 

Arrays(rsm: RichterSlopeModel, beta_max_frac: Rational | float, n_pts: int = 1001)

XXXX

Parameters

XXX

Attributes

TBD

Returns

XXX

Parameters

XXX

Attributes

TBD

Returns

XXX

Used by:
Source code in erosionfront/theory/arrays.py 
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
def __init__(
        self, 
        rsm: RichterSlopeModel,
        beta_max_frac: Rational | float,
        n_pts: int=1001,
    ) -> None:
    """
    XXXX

    Parameters
    ----------
    XXX

    Attributes
    ----------
    TBD

    Returns
    -------
    XXX
    """
    beta = np.linspace(0.0, (np.pi/2)*float(beta_max_frac), n_pts)
    if (
        rsm.px_lambda is None 
        or rsm.pz_lambda is None
        or rsm.vx_lambda is None
        or rsm.vz_lambda is None
        or rsm.p_lambda is None
        ): 
        return
    px = rsm.px_lambda(beta) 
    pz = rsm.pz_lambda(beta)
    vx = rsm.vx_lambda(beta)
    vz = rsm.vz_lambda(beta)
    p = rsm.p_lambda(beta)
    v = np.sqrt(vx**2+vz**2)

    self.n_pts = n_pts
    self.β = beta[np.isfinite(vx) & np.isfinite(p)]

    self.px = px[np.isfinite(vx) & np.isfinite(p)]
    self.pz = pz[np.isfinite(vx) & np.isfinite(p)]
    self.p = p[np.isfinite(vx) & np.isfinite(p)]

    self.v = v[np.isfinite(vx) & np.isfinite(p)]
    self.vz = vz[np.isfinite(vx) & np.isfinite(p)]
    self.vx = vx[np.isfinite(vx) & np.isfinite(p)]

    self.ξx = self.px/self.p**2
    self.ξz = self.pz/self.p**2
    self.ξrt = 1/self.px
    self.ξup = 1/self.pz