MESA

General Overview
- Good textbook by Edward Brown

- Run MESA as a supplementary test since many calculations on various conditions are already made by NuGrid collaboration and stored in CANFAR

- MESA web is quite useful as well. (you can get abundance at certain metalicity, etc).

MESA versions
- Sony VAIO laptop --> mesa-r11554 (2019 Mar. release)

- GCA server --> mesa-r10398 (2018 Mar. release). Updated to r12115 (2019 Sep. release) (2019/10/10).

Installation
Quite straight forward. Follow the official website. Also, when I updated the version on GCA server (because I didn't update "mesasdk"), it was much easier as well. Download the .zip file from MESA page, unzip, move to the directory ($MESA_DIR), and "./install".

Useful test_suits
- 25M_pre_ms_to_core_collapse, since my first goal is studying s-process in massive stars.

- agb, c13_pocket are also very useful.

- Note the test_suite might not work if the version is different (my GCA couldn't run Sony Linux ones).

How to run?
- Mostly, just make by ".mk", and then run ".rn". If you want to restart from some point, then ".re # of model" like ".re 1000". Note "mk" is a file and you may need to change the command like inlist file name.

- You may want to set pgstar, but some test_suites don't set pgstar. Set pgstar_file = .true. in inlist (names are different from suite to suite). You also need to set pgstar details in other inlists too. You can check my past inlists in some mesa test_suite directories.

Plot output data
- py_mesa_reader. I made various .py files for various conditions, by coping from original sample .py files.

MESA FOR NUGRID
I succeeded calculating 7M AGB using inlist files listed in Canfar (7/30/2020). Note keys for success are

- Always check what's written in "mk", "rn", files etc. "./star" basically do calculations, but some needs to use "rn" which is sort of shell script and copy files to inlist and so on.

- You need to change the reaction file to o18_and_ne22.net, which is much simpler than Canfar's agamass***.net.

- When you run very metal poor (Z=2e-5, e.g.), the calculations may take crazily long time.