Often in computational astrophysics, we have initial conditions given to us by a one dimensional (perhaps stellar evolution) code. These codes are usually - for convenience - formulated in a Lagrangian coordinate system, where each coordinate (e.g. mass) moves along with the fluid. Multidimensional hydrodynamics codes typically use an Eulerian coordinate system in which the coordinates (e.g. position) are fixed and the fluid moves through the them. Mapping a unequally-spaced Lagrangian model into an Eulerian code with a uniform mesh while preserving hydrostatic equilibrium can sometimes be challenging. Here I describe a variant of the procedure we often use in the CASTRO and MAESTRO codes; this is based on the paper by Mike Zingale et. al.




