Page 1 (printed 4/2/99) PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) AAAA FFFFOOOORRRRTTTTRRRRAAAANNNN RRRReeeeffffeeeerrrreeeennnncccceeee MMMMaaaannnnuuuuaaaallll PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) NNNNAAAAMMMMEEEE PPMLIB - A collection of routines for interpolation and CFD based on the Piecewise Parabolic Method (PPM) DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN PPMLIB consists of a set routines, callable from larger programs, suitable for interpolation and computational fluid dynamics. Included with the library package are several working examples. See the man page for ppmlibtutorial. The routines contined within PPMLIB are ddddoooo____ppppppppmmmmddddeeee0000____1111ddddcccc Example of using the PPMLIB library routine to create a routine to perform one directional pass of a single timestep of a 3d computational hydrodynamics calculation in direct Eulerian style assuming a uniform Cartesian grid and a gamma law equation of state, with numerical diffusion. ddddoooo____ppppppppmmmmddddeeee0000____1111ddddcccc____ggggaaaammmmmmmmaaaa Example of using the PPMLIB library routine to create a routine to perform one timestep of a 1d computational hydrodynamics calculation in direct Eulerian style assuming a uniform Cartesian grid and a gamma law equation of state, with numerical diffusion. ddddoooo____ppppppppmmmmddddeeee0000____1111ddddcccc____ppppaaaaqqqq____ggggaaaammmmmmmmaaaa Example of using the PPMLIB library routine to create a routine to perform one timestep of a 1d computational hydrodynamics calculation in direct Eulerian style assuming a uniform Cartesian grid and a gamma law equation of state, with numerical diffusion. Allowance is provided for passively advected quantities. ddddoooo____ppppppppmmmmddddeeee0000____2222ddddcccc____ggggaaaammmmmmmmaaaa Example of using the PPMLIB library routine to create a routine to perform one directional pass of a single timestep of a 2d computational hydrodynamics calculation in direct Eulerian style assuming a uniform Cartesian grid and a gamma law equation of state, with numerical diffusion. ddddoooo____ppppppppmmmmddddeeee0000____3333ddddcccc____ggggaaaammmmmmmmaaaa Example of using the PPMLIB library routine to create a Page 1 (printed 4/2/99) PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) AAAA FFFFOOOORRRRTTTTRRRRAAAANNNN RRRReeeeffffeeeerrrreeeennnncccceeee MMMMaaaannnnuuuuaaaallll PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) routine to perform one directional pass of a single timestep of a 3d computational hydrodynamics calculation in direct Eulerian style assuming a uniform Cartesian grid and a gamma law equation of state, with numerical diffusion. ddddoooo____ppppppppmmmm____llllaaaagggg0000____fffflllluuuuxxxx computes the one-dimenaional hydrodynamical Lagrangian fluxes on a uniform Cartesian grid. It does NOT assume a gamma law equation of state. ddddoooo____ppppppppmmmm____llllaaaagggg0000____fffflllluuuuxxxx____ggggaaaammmmmmmmaaaa computes the one-dimenaional hydrodynamical Lagrangian fluxes on a uniform Cartesian grid. It assumes a gamma law equation of state. ddddoooo____ppppppppmmmm____llllaaaagggg____fffflllluuuuxxxx computes the one-dimenaional hydrodynamical Lagrangian fluxes on a non-uniform Cartesian grid. It does NOT assume a gamma law equation of state. ddddoooo____ppppppppmmmm____llllaaaagggg____fffflllluuuuxxxx____ggggaaaammmmmmmmaaaa computes the one-dimenaional hydrodynamical Lagrangian fluxes on a non-uniform Cartesian grid. It assumes a gamma law equation of state. ddddoooo____ppppppppmmmmllllrrrr0000____1111ddddcccc Example of using the PPMLIB library routine to create a routine to perform one timestep of a 1d computational hydrodynamics calculation in lag + remap style assuming a uniform Cartesian grid, with numerical diffusion. A Generalized equation of state is used withing the PPMLIB routines. ddddoooo____ppppppppmmmmllllrrrr0000____1111ddddcccc____ggggaaaammmmmmmmaaaa Example of using the PPMLIB library routine to create a routine to perform one timestep of a 1d computational hydrodynamics calculation in Lag + remap style assuming a uniform Cartesian grid and a gamma law equation of state, with numerical diffusion. ddddoooo____ppppppppmmmmllllrrrr0000____1111ddddcccc____ppppaaaaqqqq____ggggaaaammmmmmmmaaaa Example of using the PPMLIB library routine to create a routine to perform one timestep of a 1d computational hydrodynamics calculation in Lagrangian + remap style Page 2 (printed 4/2/99) PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) AAAA FFFFOOOORRRRTTTTRRRRAAAANNNN RRRReeeeffffeeeerrrreeeennnncccceeee MMMMaaaannnnuuuuaaaallll PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) assuming a uniform Cartesian grid and a gamma law equation of state, with numerical diffusion. Allowance is provided for passively advected quantities. ddddoooo____ppppppppmmmmllllrrrr0000____2222ddddcccc____ggggaaaammmmmmmmaaaa Example of using the PPMLIB library routine to create a routine to perform one directional pass of a single timestep of a 2d computational hydrodynamics calculation in lag + remap style assuming a uniform Cartesian grid and a gamma law equation of state, with numerical diffusion. ddddoooo____ppppppppmmmmllllrrrr0000____3333ddddcccc____ggggaaaammmmmmmmaaaa Example of using the PPMLIB library routine to create a routine to perform one directional pass of a single timestep of a 3d computational hydrodynamics calculation in direct Eulerian style assuming a uniform Cartesian grid and a gamma law equation of state, with numerical diffusion. ppppppppmmmm99998888____aaaaddddvvvveeeecccctttt____ppppaaaassssssssiiiivvvveeee performs an entire remap of a passively advected quantity on a nonuniform Cartesian grid. ppppppppmmmm99998888____aaaaddddvvvveeeecccctttt____ppppaaaassssssssiiiivvvveeee0000 performs an entire remap of a passively advected quantity on a uniform Cartesian grid. ppppppppmmmm99998888____bbbbiiiinnnncccclllloooosssseeee Closes a file associated with a designated logical unit number. ppppppppmmmm99998888____bbbbiiiinnnnooooppppeeeennnn Opens a file associated with a designated logical unit number. ppppppppmmmm99998888____bbbbiiiinnnnrrrreeeeaaaadddd Reads a specified number of bytes from a binary file without FORTRAN record marks. ppppppppmmmm99998888____bbbbiiiinnnnwwwwrrrriiiitttteeee Writes a specified number of bytes to a binary file without FORTRAN record marks. Page 3 (printed 4/2/99) PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) AAAA FFFFOOOORRRRTTTTRRRRAAAANNNN RRRReeeeffffeeeerrrreeeennnncccceeee MMMMaaaannnnuuuuaaaallll PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) ppppppppmmmm99998888____ddddeeeeuuuullll____uuuuppppddddaaaatttteeee____ccccaaaarrrrtttt updates zone-averaged hydrodynamic quantities given zone interface fluxes. 1D. ppppppppmmmm99998888____ddddiiiifffftttteeeerrrrpppp0000____ddddiiiissssccccoooonnnntttt____ccccoooommmmpppp computes piecewise parabolic interpolations based on data which are differences of zone-averages. With monotonicity constraints and discontinuity detection. Provision for a secondary condition on the discontinuity detection is provided. for uniform grids. ppppppppmmmm99998888____ddddiiiifffftttteeeerrrrpppp0000____ddddiiiissssccccoooonnnnttttiiiinnnnuuuuiiiittttyyyy computes piecewise parabolic interpolations based on data which are differences of zone-averages. With monotonicity constraints and discontinuity detection. for uniform grids. ppppppppmmmm99998888____ddddiiiifffftttteeeerrrrpppp0000____iiiinnnntttteeeerrrrffffaaaacccceeee interpolatesthird-order accurate interface point values based on data which are differences of zone-averages. Assumes a smooth data distribution, with no monotonicity constraints and discontinuity detection. For uniform grids. ppppppppmmmm99998888____ddddiiiifffftttteeeerrrrpppp0000____mmmmoooonnnnoooo computes piecewise parabolic interpolations based on data which are differences of zone-averages, with monotonicity constraints. for uniform grids. ppppppppmmmm99998888____ddddiiiifffftttteeeerrrrpppp____ddddiiiissssccccoooonnnntttt____ccccoooommmmpppp computes piecewise parabolic interpolations based on data which are differences of zone-averages. With monotonicity constraints and discontinuity detection. Provision for a secondary condition on the discontinuity detection is provided. For non-uniform grids. ppppppppmmmm99998888____ddddiiiifffftttteeeerrrrpppp____ddddiiiissssccccoooonnnnttttiiiinnnnuuuuiiiittttyyyy computes piecewise parabolic interpolations based on data which are differences of zone-averages. With monotonicity constraints and discontinuity detection. For non-uniform grids. ppppppppmmmm99998888____ddddiiiifffftttteeeerrrrpppp____iiiinnnntttteeeerrrrffffaaaacccceeee Page 4 (printed 4/2/99) PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) AAAA FFFFOOOORRRRTTTTRRRRAAAANNNN RRRReeeeffffeeeerrrreeeennnncccceeee MMMMaaaannnnuuuuaaaallll PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) interpolatesthird-order accurate interface point values based on data which are differences of zone-averages. Assumes a smooth data distribution, with no monotonicity constraints and discontinuity detection. For non-uniform grids. ppppppppmmmm99998888____ddddiiiifffftttteeeerrrrpppp____mmmmoooonnnnoooo computes piecewise parabolic interpolations based on data which are differences of zone-averages, with monotonicity constraints. For non-niform grids. ppppppppmmmm99998888____ddddiiiiffffuuuusssseeeefffflllluuuuxxxx____hhhhyyyyddddrrrroooo____ccccaaaarrrrtttt constructd zone interface centered fluxes from zone- averaged difusion velocities and zone averaged of the hydrodynamic quantities. ppppppppmmmm99998888____ddddiiiiffffuuuusssseeeefffflllluuuuxxxx____ppppaaaassssssssiiiivvvveeee constructd zone interface centered passive fluxes from zone interface cenetered difusion fluxes. ppppppppmmmm99998888____ddddiiiissssccccoooonnnnttttiiiinnnnuuuuiiiittttyyyy0000 Computes PPM unsmoothness parameter providing a measure of how well is data are resolved. Used zone-averaged data, and the PPM discontinuity parameter providing a measure of how likely a given zone contains a discontinuity. For zone-averged data. For use on uniform grids. ppppppppmmmm99998888____ddddiiiissssccccoooonnnnttttiiiinnnnuuuuiiiittttyyyy Computes PPM unsmoothness parameter providing a measure of how well is data are resolved. Used zone-averaged data, and the PPM discontinuity parameter providing a measure of how likely a given zone contains a discontinuity. For zone-averged data. For use on non-uniform grids. ppppppppmmmm99998888____eeee____rrrriiiieeeemmmmaaaannnnnnnn solves Riemann's shock tube problem to provide time- averaged values of density, pressure, and velocity. A generalized equation of state is used. ppppppppmmmm99998888____eeee____rrrriiiieeeemmmmaaaannnnnnnn____ggggaaaammmmmmmmaaaa solves Riemann's shock tube problem to provide time- averaged values of density, pressure, and velocity. A generalized equation of state is assumed. Page 5 (printed 4/2/99) PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) AAAA FFFFOOOORRRRTTTTRRRRAAAANNNN RRRReeeeffffeeeerrrreeeennnncccceeee MMMMaaaannnnuuuuaaaallll PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) ppppppppmmmm99998888____eeeeuuuullll____ddddiiiiffffuuuusssseeee0000____2222dddd Constructs zone centered Eulerian difusion velocities. Assumes 2 dimensions and a uniform Cartesian grid. No equation of state dependence is assumed. ppppppppmmmm99998888____eeeeuuuullll____ddddiiiiffffuuuusssseeee0000____3333dddd Constructs zone centered Eulerian difusion velocities. Assumes 3 dimensions and a uniform Cartesian grid. No equation of state dependence is assumed. ppppppppmmmm99998888____eeeeuuuullll____ddddiiiiffffuuuusssseeee0000____ggggaaaammmmmmmmaaaa____2222dddd Constructs zone centered Eulerian difusion velocities. Assumes 2 dimensions and a uniform Cartesian grid. A gamma law equation of state is assumed. ppppppppmmmm99998888____eeeeuuuullll____ddddiiiiffffuuuusssseeee0000____ggggaaaammmmmmmmaaaa____3333dddd Constructs zone centered Eulerian difusion velocities. Assumes 3 dimensions and a uniform Cartesian grid. A gamma law equation of state is assumed. ppppppppmmmm99998888____eeeeuuuullll____ddddiiiiffffuuuusssseeee____1111dddd Constructs zone centered Eulerian difusion velocities. Assumes 1 dimension and a Cartesian grid, which need not be uniform. No equation of state dependence is assumed. ppppppppmmmm99998888____eeeeuuuullll____ddddiiiiffffuuuusssseeee____2222dddd Constructs zone centered Eulerian difusion velocities. Assumes 2 dimensions and a Cartesian grid, which need not be uniform. No equation of state dependence is assumed. ppppppppmmmm99998888____eeeeuuuullll____ddddiiiiffffuuuusssseeee____3333dddd Constructs zone centered Eulerian difusion velocities. Assumes 3 dimensions and a Cartesian grid, which need not be uniform. No equation of state dependence is assumed. ppppppppmmmm99998888____eeeeuuuullll____ddddiiiiffffuuuusssseeee____ggggaaaammmmmmmmaaaa____2222dddd Constructs zone centered Eulerian difusion velocities. Assumes 2 dimensions and a Cartesian grid, which need not be uniform. A gamma law equation of state is assumed. Page 6 (printed 4/2/99) PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) AAAA FFFFOOOORRRRTTTTRRRRAAAANNNN RRRReeeeffffeeeerrrreeeennnncccceeee MMMMaaaannnnuuuuaaaallll PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) ppppppppmmmm99998888____eeeeuuuullll____ddddiiiiffffuuuusssseeee____ggggaaaammmmmmmmaaaa____3333dddd Constructs zone centered Eulerian difusion velocities. Assumes 3 dimensions and a Cartesian grid, which need not be uniform. A gamma law equation of state is assumed. ppppppppmmmm99998888____eeeeuuuulllleeeerrrr0000____ddddoooommmmaaaaiiiinnnn____aaaavvvvgggg Interpolates density, pressure and longitudinal velocity. Determines spatial averages of pressure and longitudinal velocity over the + and - Eulerian domains of dependence. Assumes a uniform grid for efficiency. ppppppppmmmm99998888____eeeeuuuulllleeeerrrr____ddddoooommmmaaaaiiiinnnn____aaaavvvvgggg Interpolates density, pressure and longitudinal velocity. Determines spatial averages of pressure and longitudinal velocity over the + and - Eulerian domains of dependence. ppppppppmmmm99998888____iiiinnnntttteeeerrrrffffaaaacccceeee____ffffaaaaccccttttoooorrrrssss Precomputes various geometrical factors for interpolation of interface values on non-uniform grids. ppppppppmmmm99998888____iiiinnnntttteeeerrrrpppp0000____ddddiiiissssccccoooonnnntttt____ccccoooommmmpppp computes piecewise parabolic interpolations based on zone-averaged data. With monotonicity constraints and discontinuity detection. Provision for a secondary condition on the discontinuity detection is provided. For uniform grids. ppppppppmmmm99998888____iiiinnnntttteeeerrrrpppp0000____ddddiiiissssccccoooonnnnttttiiiinnnnuuuuiiiittttyyyy computes piecewise parabolic interpolations based on zone-averaged data. With monotonicity constraints and discontinuity detection. For uniform grids. ppppppppmmmm99998888____iiiinnnntttteeeerrrrpppp0000____iiiinnnntttteeeerrrrffffaaaacccceeee interpolatesthird-order accurate interface point values based on zone-averaged data. Assumes a smooth data distribution, with no monotonicity constraints and discontinuity detection. For uniform grids. ppppppppmmmm99998888____iiiinnnntttteeeerrrrpppp0000____mmmmoooonnnnoooo computes piecewise parabolic interpolations based on zone-averaged data. with monotonicity constraints. For uniform grids. Page 7 (printed 4/2/99) PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) AAAA FFFFOOOORRRRTTTTRRRRAAAANNNN RRRReeeeffffeeeerrrreeeennnncccceeee MMMMaaaannnnuuuuaaaallll PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) ppppppppmmmm99998888____iiiinnnntttteeeerrrrpppp0000____rrrrhhhhoooo Computed piecewise parabolic interpolations of density based on zone-averaged data. With monotonicity constraints and contact-discontinuity detection. For use on uniform grids. ppppppppmmmm99998888____iiiinnnntttteeeerrrrpppp0000____uuuupppp Computed piecewise parabolic interpolations of pressure and velocity based on zone-averaged data. With monotonicity constraints and discontinuity detection. For use on uniform grids. ppppppppmmmm99998888____iiiinnnntttteeeerrrrpppp____ddddiiiissssccccoooonnnntttt____ccccoooommmmpppp computes piecewise parabolic interpolations based on zone-averaged data. With monotonicity constraints and discontinuity detection. Provision for a secondary condition on the discontinuity detection is provided. For non-uniform grids. ppppppppmmmm99998888____iiiinnnntttteeeerrrrpppp____ddddiiiissssccccoooonnnnttttiiiinnnnuuuuiiiittttyyyy computes piecewise parabolic interpolations based on zone-averaged data. Does not impose monotonicity constraints - for smooth data distributions only. For non-uniform grid. ppppppppmmmm99998888____iiiinnnntttteeeerrrrpppp____iiiinnnntttteeeerrrrffffaaaacccceeee interpolatesthird-order accurate interface point values based on zone-averaged data. Assumes a smooth data distribution, with no monotonicity constraints and discontinuity detection. For non-uniform grids. ppppppppmmmm99998888____iiiinnnntttteeeerrrrpppp____mmmmoooonnnnoooo computes piecewise parabolic interpolations based on zone-averaged data. with monotonicity constraints. For non-uniform grids. ppppppppmmmm99998888____iiiinnnntttteeeerrrrpppp____rrrrhhhhoooo Computed piecewise parabolic interpolations of density based on zone-averaged data. With monotonicity constraints and contact-discontinuity detection. For use on non-uniform grids. ppppppppmmmm99998888____iiiinnnntttteeeerrrrpppp____uuuupppp Computed piecewise parabolic interpolations of pressure and velocity based on zone-averaged data. With Page 8 (printed 4/2/99) PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) AAAA FFFFOOOORRRRTTTTRRRRAAAANNNN RRRReeeeffffeeeerrrreeeennnncccceeee MMMMaaaannnnuuuuaaaallll PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) monotonicity constraints and discontinuity detection. For use on non-uniform grids. ppppppppmmmm99998888____iiiinnnntttteeeerrrrppppoooollllaaaattttiiiioooonnnn____ffffaaaaccccttttoooorrrrssss Precomputes various geometrical factors for piecewise parabolic interpolations values on non-uniform grids. ppppppppmmmm99998888____llll____rrrriiiieeeemmmmaaaannnnnnnn solves Riemann's shock tube problem to provide time- averaged values of pressure and velocity. A generalized equation of state is used. ppppppppmmmm99998888____llll____rrrriiiieeeemmmmaaaannnnnnnn____ggggaaaammmmmmmmaaaa solves Riemann's shock tube problem to provide time- averaged values of pressure and velocity. A gamma law equation of state is assumed. ppppppppmmmm99998888____llllaaaagggg0000____ddddoooommmmaaaaiiiinnnn____aaaavvvvgggg Determines piecewise parabolic interpolations for pressure and longitudinal velocity, and finds averages of pressure and velocity in the plus and minus Lagrangian domains of dependence for each lefthand zone interface. For uniform grids. ppppppppmmmm99998888____llllaaaagggg____ddddoooommmmaaaaiiiinnnn____aaaavvvvgggg Determines piecewise parabolic interpolations for pressure and longitudinal velocity, and finds averages of pressure and velocity in the plus and minus Lagrangian domains of dependence for each lefthand zone interface. For non-uniform grids. ppppppppmmmm99998888____llllaaaagggg____ddddiiiiffffuuuusssseeee0000____2222dddd Constructs zone centered Eulerian difusion velocities. Assumes 2 dimensions and a uniform Cartesian grid. No equation of state dependence is assumed. A special Courant (stability) calculation is performed to ensure that no wave traverses more than one Lagrangian zone. ppppppppmmmm99998888____llllaaaagggg____ddddiiiiffffuuuusssseeee0000____3333dddd Constructs zone centered Eulerian difusion velocities. Assumes 3 dimensions and a uniform Cartesian grid. No equation of state dependence is assumed. A special Courant (stability) calculation is performed to ensure Page 9 (printed 4/2/99) PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) AAAA FFFFOOOORRRRTTTTRRRRAAAANNNN RRRReeeeffffeeeerrrreeeennnncccceeee MMMMaaaannnnuuuuaaaallll PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) that no wave traverses more than one Lagrangian zone. ppppppppmmmm99998888____llllaaaagggg____ddddiiiiffffuuuusssseeee0000____ggggaaaammmmmmmmaaaa____2222dddd Constructs zone centered Eulerian difusion velocities. Assumes 2 dimensions and a uniform Cartesian grid. Assumes a gamma law equation of state. A special Courant (stability) calculation is performed to ensure that no wave traverses more than one Lagrangian zone. ppppppppmmmm99998888____llllaaaagggg____ddddiiiiffffuuuusssseeee0000____ggggaaaammmmmmmmaaaa____3333dddd Constructs zone centered Eulerian difusion velocities. Assumes 3 dimensions and a uniform Cartesian grid. Assumes a gamma law equation of state. A special Courant (stability) calculation is performed to ensure that no wave traverses more than one Lagrangian zone. ppppppppmmmm99998888____llllaaaagggg____ddddiiiiffffuuuusssseeee____1111dddd Constructs zone centered Eulerian difusion velocities. Assumes 1 dimensions and a Cartesian grid which need not be uniform. No equation of state dependence is assumed. A special Courant (stability) calculation is performed to ensure that no wave traverses more than one Lagrangian zone. ppppppppmmmm99998888____llllaaaagggg____uuuuppppddddaaaatttteeee____ccccaaaarrrrtttt performs a Lagrangian update of the zone-averaged quantities given Lagrangian zone interface fluxes (pressure and velocity) in Cartesian coordinates. 1D. ppppppppmmmm99998888____mmmmaaaapppp____fffflllluuuuxxxx____1111dddd Computes the Conservative zone interface fluxes (mass, momentum and energy) necessary to map hydrodynamic quantities from one grid to another. based on zone averaged data. Assumes a generalized equation of state. 1D ppppppppmmmm99998888____mmmmaaaapppp____fffflllluuuuxxxx____ggggaaaammmmmmmmaaaa____1111dddd Computes the Conservative zone interface fluxes (mass, momentum and energy) necessary to map hydrodynamic quantities from one grid to another. based on zone averaged data. Assumes a gamma law equation of state. 1D Page 10 (printed 4/2/99) PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) AAAA FFFFOOOORRRRTTTTRRRRAAAANNNN RRRReeeeffffeeeerrrreeeennnncccceeee MMMMaaaannnnuuuuaaaallll PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) ppppppppmmmm99998888____mmmmaaaapppp____fffflllluuuuxxxx____ggggaaaammmmmmmmaaaa____2222dddd Computes the Conservative zone interface fluxes (mass, momentum and energy) necessary to map hydrodynamic quantities from one grid to another. based on zone averaged data. Assumes a gamma law equation of state. 2D ppppppppmmmm99998888____mmmmaaaapppp____fffflllluuuuxxxx____ggggaaaammmmmmmmaaaa____3333dddd Computes the Conservative zone interface fluxes (mass, momentum and energy) necessary to map hydrodynamic quantities from one grid to another. based on zone averaged data. Assumes a gamma law equation of state. 3D ppppppppmmmm99998888____mmmmaaaapppp____uuuuppppddddaaaatttteeee____1111dddd Updates the zone averaged conserved hydrodynamic quantities given the advective fluxes. 1D ppppppppmmmm99998888____mmmmaaaapppp____uuuuppppddddaaaatttteeee____2222dddd Updates the zone averaged conserved hydrodynamic quantities given the advective fluxes. 2D ppppppppmmmm99998888____mmmmaaaapppp____uuuuppppddddaaaatttteeee____3333dddd Updates the zone averaged conserved hydrodynamic quantities given the advective fluxes. 3D ppppppppmmmm99998888____mmmmeeeesssshhhhoooo2222 writes the 2D grid information into a compressed dump file ppppppppmmmm99998888____mmmmeeeesssshhhhoooo3333 writes the 3D grid information into a compressed dump file ppppppppmmmm99998888____mmmmoooonnnnoooo____ppppaaaarrrraaaabbbbllllooooaaaa Determines monotonized piecewise parabolic interpolation coefficients from the zone averages and un-monotonized left and right hand zone interface point values. Monotonized and unmonotonized coefficients are blended using UNSMTH. ppppppppmmmm99998888____ppppaaaassssssssiiiivvvveeee____fffflllluuuuxxxx determines piecewise parabolic interpolations zone- averages of a passive quantity and determines the Page 11 (printed 4/2/99) PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) AAAA FFFFOOOORRRRTTTTRRRRAAAANNNN RRRReeeeffffeeeerrrreeeennnncccceeee MMMMaaaannnnuuuuaaaallll PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) fluxes necessary to map it from one grid to another. For non-unifrom grids. ppppppppmmmm99998888____ppppaaaassssssssiiiivvvveeee____fffflllluuuuxxxx0000 determines piecewise parabolic interpolations zone- averages of a passive quantity and determines the fluxes necessary to map it from one grid to another. For unifrom grids. ppppppppmmmm99998888____pppprrrriiiinnnntttt____ccccoooonnnnssssttttaaaannnnttttssss displays the current values of the PPMLIB internal constants which determine the behaviour of the interpolation and difusion routines. ppppppppmmmm99998888____rrrreeeesssseeeetttt____ccccoooonnnnssssttttaaaannnnttttssss sets all PPMLIB internal constants, which determine the behaviour of the interpolation and difusion routines, to their default values. ppppppppmmmm99998888____sssseeeetttt____ccccoooonnnnssssttttaaaannnnttttssss allows the user to set values of individual PPMLIB internal constants, which determine the behaviour of the interpolation and difusion routines. ppppppppmmmm99998888____ttttiiiimmmmeeeessssttttaaaammmmpppp record the date and time of creation in a PPMLIB compressed dump file or restart dump ppppppppmmmm99998888____uuuunnnnssssmmmmooooooootttthhhh0000 Computes PPM unsmoothness parameter providing a measure of how well is data are resolved. Used zone-averaged data. For use on uniform grids. ppppppppmmmm99998888____uuuunnnnssssmmmmooooooootttthhhh Computes PPM unsmoothness parameter providing a measure of how well is data are resolved. Used zone-averaged data. For use on non-uniform grids. ppppppppmmmm99998888____vvvvaaaarrrr0000dddd writes a scalar real*4 floating point value into a compressed dump file Page 12 (printed 4/2/99) PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) AAAA FFFFOOOORRRRTTTTRRRRAAAANNNN RRRReeeeffffeeeerrrreeeennnncccceeee MMMMaaaannnnuuuuaaaallll PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) ppppppppmmmm99998888____vvvvaaaarrrr2222dddd writes a 2D floating point array into a compressed dump file ppppppppmmmm99998888____vvvvaaaarrrr3333dddd writes a 3D floating point array into a compressed dump file ppppppppmmmm99998888____wwwweeeeiiiigggghhhhtttteeeedddd____uuuuppppddddaaaatttteeee updates zone-averages given weighted interface fluxes and the old and new zone-average weightings ppppppppmmmm99998888____wwwwrrrriiiitttteeee____ccccoooommmmmmmmeeeennnntttt writes an ASCII comment into a PPMLIB compressed dump file or restart dump LLLLIIIINNNNKKKKIIIINNNNGGGG WWWWIIIITTTTHHHH TTTTHHHHEEEE LLLLIIIIBBBBRRRRAAAARRRRYYYY The library is callable from wither FORTRAN90 or FORTRAN77. Examples or each are +o f90 myprog.f $PPMLIBO ppmlib.???? +o f77 myprog.f $PPMLIBO ppmlib.???? -lfortran HHHHIIIINNNNTTTTSSSS +o the names of double precision versions of the above routines are preceeded by "d_". For example dddd____ppppppppmmmm99998888____iiiinnnntttteeeerrrrpppp____uuuupppp is the double precision version of ppppppppmmmm99998888____iiiinnnntttteeeerrrrpppp____uuuupppp. +o In order to access the man pages, the map $PPMLIBPATH/Man must be added to the user's MANPATH environmental variable, where $PPMLIBPATH is the location of the ppmlib man page directory. +o SSSSEEEEEEEE AAAALLLLSSSSOOOO PPMLIBTUTORIAL, DO_PPM, and individual man pages Page 13 (printed 4/2/99) PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) AAAA FFFFOOOORRRRTTTTRRRRAAAANNNN RRRReeeeffffeeeerrrreeeennnncccceeee MMMMaaaannnnuuuuaaaallll PPPPPPPPMMMMLLLLIIIIBBBB((((1111)))) BBBBUUUUGGGG RRRREEEEPPPPOOOORRRRTTTTSSSS TTTTOOOO ppmlib@sapphire.lcse.umn.edu AAAAUUUUTTTTHHHHOOOORRRRSSSS PPPPaaaauuuullll RRRR.... WWWWooooooooddddwwwwaaaarrrrdddd paul@lcse.umn.edu BBBB.... KKKKeeeevvvviiiinnnn EEEEddddggggaaaarrrr bke@lcse.umn.edu Department of Astronomy University of Minnesota Minneapolis, MN 55455 USA Page 14 (printed 4/2/99)