TLUSTY_CLOUDY-Auto_Scripts

Structure and scripts that allow for quick, multithreaded tlusty and cloudy runs.

Set-up

1. Download and unzip the file from github. Name the top level directory "tlusty+cloudy"
2. Download TLUSTY version 205 from a repository. Here's one https://github.com/gabraganca/tlusty
3. Follow the set up instructions in the Tlusty manual. Name the top level directory "tlusty205"
4. Download Cloudy version C25. Here's the link https://gitlab.nublado.org/cloudy/cloudy
5. Follow the set up instructions on the Cloudy github wiki. Name the top level directory "c25.00"
6. Move the "tlusty_runs" folder from the "tlusty+cloudy" directory into the "tlusty205" directory
7. Move the "cloudy_runs" folder from the "tlusty+cloudy" directory into the "c25.00" directory
8. There are some "DELETEME" files in those folders. Make sure you go through each folder and delete them
9. Ensure that the three folders ("tlusty205" , "c25.00" , "tlusty+cloudy") are in your home directory
10. Ensure your python environment has all the modules used by the scripts. They are all at the top of each script, so just open them up and check

Running the Code

This code can automatically handle any Dark Star in a mass range from 1e5 to 1e7 solar masses, both Capture and AC. TLUSTY creates a SED for the star, and Cloudy then simulates a 100 pc hydrogen nebula around it. (To alter that stromgren radius, go to line 208 in "cloudy_auto.py" and change the 100 to your desired number of parsecs)

If you want to run both tlusty and cloudy in sequence:

1. Find the input file in the "tlusty+cloudy" folder
2. On the first line write the log10 solar masses' of the stars you want to simulate, separated by spaces e.g. 5 5.5 6 6.2
3. On the second line write the corresponding Dark Star type for the mass above, separated by spaces e.g. AC AC Capture AC
4. An example of how this should look will be in the downloaded files
5. Run "tlusty+cloudy.py"
6. Outputs will be in the corresponding "auto_outputs" folder in each directory, and graphs will be in the corresponding "plots" folder

If you just want to only run tlusty by itself, run the "tlusty_only.py" script instead of the "tlusty+cloudy.py" script

If you want to re-graph and look at individual output SEDS from tlusty and cloudy, you can use the corresponding "tlusty_viewer.ipynb" and "cloudy_viewer.ipynb" scripts found in the "tlusty+cloudy" folder (TO BE DONE)

If you have any questions or encounter any errors, email me at ggasper@colgate.edu