This repository contains scripts and data for Figures 2-7 in the paper by Ibanez-Mejia, Mac Low, & Klessen (2021, Astrophys. J., submitted, arXiv:2108.04967) entitled "Gravity Versus Magnetic Fields in Forming Molecular Clouds".  The authors may be contacted at mordecai@amnh.org. Problems with data access can be addressed by the Library (libref@amnh.org).

This material is available by ftp at

ftp://ftp.amnh.org/pub/group/publish/mordecai/sd.astro.4

or by Globus at

https://app.globus.org/file-manager?origin_id=d8d218f0-e092-11ea-8bbc-029866a337f1&origin_path=%2Fftp%2Fdata%2Fpub%2Fgroup%2Fpublish%2Fmordecai%2F%2Fsd.astro.4

The scripts for each Figure are in the Scripts subdirectory, while the data is in the Data subdirectory or in sd.astro.3 as described below.

The base data for the plots are data cubes produced from zoom-in simulations of 100 pc cubes centered on newly forming molecular clouds of final mass 3000-8000 solar masses.  These data cubes are available from the repository Chira et al. (2018, https://dx.doi.org/10.5531/sd.astro.3) in the directories M[348]e3 at

ftp://ftp.amnh.org/pub/group/publish/mordecai/sd.astro.3/src_data/ 

That data can also be accessed via Globus in directories M[348]e3 at

https://app.globus.org/file-manager?origin_id=d8d218f0-e092-11ea-8bbc-029866a337f1&origin_path=%2Fftp%2Fdata%2Fpub%2Fgroup%2Fpublish%2Fmordecai%2F%2Fsd.astro.3%2Fsrc_data

The intermediate data products in the Data directory were generated from this base data using New_Follow_a_Cloud.py, which also accesses particle data from the original simulation in the directory Data/Particles.

The scripts rely on the Python visualization package yt (Turk et al. 2010) version 3.2.3 (Changeset = bdea84d95099, 11 May 2016), available from https://yt-project.org. Prior to running the scripts, yt must be activated from bin/activate within the installed package.

Supplemental functions for the scripts can be found in the repository directory Functions, which should be added to the system variable PYTHONPATH.  Although the scripts may well work with more recent versions of yt, they have not been tested and may need work to port them forward.

The individual figures are generated with the following scripts.

Figure 2: Tile_Energy_Velocity.py with energies = True

Figure 3: Tile_Energy_Velocity.py with velocities = True

Figure 4: Tile_ForcePlot.py

Figures 5 and 6: B-n-figures.py with axis = "x" or "total"

Figure 7: Routine_Figure2.py