recon1d_mplm_wa module reference
Monotonized Piecewise Linear Method 1D reconstruction.
Data Types
Limited Monotonic PLM reconstruction following White and Adcroft, 2008. |
Functions/Subroutines
Calculate a 1D PLM reconstructions based on h(:) and u(:) |
|
Returns a limited PLM slope following White and Adcroft, 2008, in the same arbitrary units [A] as the input values. |
|
Returns a limited PLM slope following Colella and Woodward 1984, in the same arbitrary units as the input values [A]. |
|
Checks the MPLM_WA reconstruction for consistency. |
|
Runs PLM reconstruction unit tests and returns True for any fails, False otherwise. |
Detailed Description
This implementation of PLM follows White and Adcroft, 2008 [87]. The PLM slopes are first limited following Colella and Woodward, 1984, but are then further limited to ensure the edge values moving across cell boundaries are monotone. The first and last cells are always limited to PCM.
This differs from recon1d_mplm_wa_poly() in the internally not polynomial representations are referred to. in the internally not polynomial representations are referred to.
Type Documentation
- type recon1d_mplm_wa/mplm_wa
Limited Monotonic PLM reconstruction following White and Adcroft, 2008.
- Type fields:
% reconstruct ::
procedure, privateImplementation of the MPLM_WA reconstruction.% check_reconstruction ::
procedure, privateImplementation of check reconstruction for the MPLM_WA reconstruction.% unit_tests ::
procedure, privateImplementation of unit tests for the MPLM_WA reconstruction.% reconstruct_parent ::
procedure, privateDuplicate interface to
Function/Subroutine Documentation
- subroutine recon1d_mplm_wa/reconstruct(this, h, u)
Calculate a 1D PLM reconstructions based on h(:) and u(:)
- Parameters:
this ::
this[inout] This reconstructionh ::
h[in] Grid spacing (thickness) [typically H]u ::
u[in] Cell mean values [A]
- Call to:
- function recon1d_mplm_wa/plm_slope_wa(h_l, h_c, h_r, h_neglect, u_l, u_c, u_r)
Returns a limited PLM slope following White and Adcroft, 2008, in the same arbitrary units [A] as the input values. Note that this is not the same as the Colella and Woodward method.
- Parameters:
h_l ::
h_l[in] Thickness of left cell in arbitrary grid thickness units [H]h_c ::
h_c[in] Thickness of center cell in arbitrary grid thickness units [H]h_r ::
h_r[in] Thickness of right cell in arbitrary grid thickness units [H]h_neglect ::
h_neglect[in] A negligible thickness [H]u_l ::
u_l[in] Value of left cell in arbitrary units [A]u_c ::
u_c[in] Value of center cell in arbitrary units [A]u_r ::
u_r[in] Value of right cell in arbitrary units [A]
- Called from:
- function recon1d_mplm_wa/plm_monotonized_slope(u_l, u_c, u_r, s_l, s_c, s_r)
Returns a limited PLM slope following Colella and Woodward 1984, in the same arbitrary units as the input values [A].
- Parameters:
u_l ::
u_l[in] Value of left cell in arbitrary units [A]u_c ::
u_c[in] Value of center cell in arbitrary units [A]u_r ::
u_r[in] Value of right cell in arbitrary units [A]s_l ::
s_l[in] PLM slope of left cell [A]s_c ::
s_c[in] PLM slope of center cell [A]s_r ::
s_r[in] PLM slope of right cell [A]
- Called from:
- function recon1d_mplm_wa/check_reconstruction(this, h, u)
Checks the MPLM_WA reconstruction for consistency.
- Parameters:
this ::
this[in] This reconstructionh ::
h[in] Grid spacing (thickness) [typically H]u ::
u[in] Cell mean values [A]
- function recon1d_mplm_wa/unit_tests(this, verbose, stdout, stderr)
Runs PLM reconstruction unit tests and returns True for any fails, False otherwise.
- Parameters:
this ::
this[inout] This reconstructionverbose ::
verbose[in] True, if verbosestdout ::
stdout[in] I/O channel for stdoutstderr ::
stderr[in] I/O channel for stderr