{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Running a real POSYDON MESA HMS-HMS gird" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "It is now time to become a POSYDON MESA architect. In this notebook, we will run a real HMS-HMS grid. We will run a downsampled version (100 models) of the HMS-HMS, 0.1Zsun, q=0.7 grid slice, using the same `HMS-HMS.ini` file as in the [first notebook](https://posydon.org/POSYDON/tutorials-examples/MESA-grids/1_hms_hms.html), where we run one binary." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Creating the Initialization File to Submit the MESA Grid" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Let's copy some code used in `$PATH_TO_POSYDON/grid_params/POSYDON-MESA-INLISTS/r11701/running_scripts/parameter_space_v2/create_csv.ipynb` to create the `grid_test.csv` file and edit it to run a downsampled version of the grid." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "total resolution Z_n * m1_n * q_n * p_n= 100\n" ] } ], "source": [ "import os\n", "import csv\n", "import numpy as np\n", "\n", "def log_range(x_min,x_max,x_n):\n", " return 10**np.linspace(np.log10(x_min),np.log10(x_max), x_n)\n", "\n", "# digit rounding pick 10 to be sure we resolve 10^-4Zsun\n", "NDIG = 10\n", "\n", "Zsun = 0.0142\n", "m1_min = 5.5\n", "m1_max = 300\n", "m1_n = 10\n", "m1 = log_range(m1_min,m1_max,m1_n)\n", "q_n = 1\n", "q = [0.7]\n", "p_min = 1.0e-1\n", "p_max = 1.0e+4\n", "p_n = 10\n", "p = log_range(p_min,p_max,p_n)\n", "Z_n = 1\n", "met = [0.1*Zsun]\n", "print('total resolution Z_n * m1_n * q_n * p_n=', Z_n * m1_n * q_n * p_n)\n", "\n", "# save entire grid in a single file\n", "with open('./grid_test.csv', 'w', newline='') as file:\n", " writer = csv.writer(file)\n", " writer.writerow(['initial_z','Zbase','m1','m2','initial_period_in_days'])\n", " for Z in met:\n", " for i in range(m1_n):\n", " for j in range(q_n):\n", " for k in range(p_n):\n", " writer.writerow([round(Z,NDIG),round(Z,NDIG),round(m1[i],NDIG),round(m1[i]*q[j],NDIG),round(p[k],NDIG)])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Run the MESA Grid with the POSYDON Submission Script" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We are now ready to run the simulation, with the following commands in your terminal (Please remember to have the POSYDON and MESA variables being set before):\n", "\n", "```bash\n", "posydon-setup-grid --grid-type fixed --inifile HMS-HMS.ini --submission-type slurm\n", "./run_grid.sh\n", "```\n", "\n", "Sit back relax, and wait for the grid to finish." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "binary\n", "cleanup.slurm\n", "column_lists\n", "grid_test.csv\n", "HMS-HMS.ini\n", "job_array_grid_submit.slurm\n", "mesa_cleanup.out\n", "mesa_grid.38800304_0.out\n", "mesa_grid.38800304_10.out\n", "mesa_grid.38800304_11.out\n", "mesa_grid.38800304_12.out\n", "mesa_grid.38800304_13.out\n", "mesa_grid.38800304_14.out\n", "mesa_grid.38800304_15.out\n", "mesa_grid.38800304_16.out\n", "mesa_grid.38800304_17.out\n", "mesa_grid.38800304_18.out\n", "mesa_grid.38800304_19.out\n", "mesa_grid.38800304_1.out\n", "mesa_grid.38800304_20.out\n", "mesa_grid.38800304_21.out\n", "mesa_grid.38800304_22.out\n", "mesa_grid.38800304_23.out\n", "mesa_grid.38800304_24.out\n", "mesa_grid.38800304_25.out\n", "mesa_grid.38800304_26.out\n", "mesa_grid.38800304_27.out\n", "mesa_grid.38800304_28.out\n", "mesa_grid.38800304_29.out\n", "mesa_grid.38800304_2.out\n", "mesa_grid.38800304_30.out\n", "mesa_grid.38800304_31.out\n", "mesa_grid.38800304_32.out\n", "mesa_grid.38800304_33.out\n", "mesa_grid.38800304_34.out\n", "mesa_grid.38800304_35.out\n", "mesa_grid.38800304_36.out\n", "mesa_grid.38800304_37.out\n", "mesa_grid.38800304_38.out\n", "mesa_grid.38800304_39.out\n", "mesa_grid.38800304_3.out\n", "mesa_grid.38800304_40.out\n", "mesa_grid.38800304_41.out\n", "mesa_grid.38800304_42.out\n", "mesa_grid.38800304_43.out\n", "mesa_grid.38800304_44.out\n", "mesa_grid.38800304_45.out\n", "mesa_grid.38800304_46.out\n", "mesa_grid.38800304_47.out\n", "mesa_grid.38800304_48.out\n", "mesa_grid.38800304_49.out\n", "mesa_grid.38800304_4.out\n", "mesa_grid.38800304_50.out\n", "mesa_grid.38800304_51.out\n", "mesa_grid.38800304_52.out\n", "mesa_grid.38800304_53.out\n", "mesa_grid.38800304_54.out\n", "mesa_grid.38800304_55.out\n", "mesa_grid.38800304_56.out\n", "mesa_grid.38800304_57.out\n", "mesa_grid.38800304_58.out\n", "mesa_grid.38800304_59.out\n", "mesa_grid.38800304_5.out\n", "mesa_grid.38800304_60.out\n", "mesa_grid.38800304_61.out\n", "mesa_grid.38800304_62.out\n", "mesa_grid.38800304_63.out\n", "mesa_grid.38800304_64.out\n", "mesa_grid.38800304_65.out\n", "mesa_grid.38800304_66.out\n", "mesa_grid.38800304_67.out\n", "mesa_grid.38800304_68.out\n", "mesa_grid.38800304_69.out\n", "mesa_grid.38800304_6.out\n", "mesa_grid.38800304_70.out\n", "mesa_grid.38800304_71.out\n", "mesa_grid.38800304_72.out\n", "mesa_grid.38800304_73.out\n", "mesa_grid.38800304_74.out\n", "mesa_grid.38800304_75.out\n", "mesa_grid.38800304_76.out\n", "mesa_grid.38800304_77.out\n", "mesa_grid.38800304_78.out\n", "mesa_grid.38800304_79.out\n", "mesa_grid.38800304_7.out\n", "mesa_grid.38800304_80.out\n", "mesa_grid.38800304_81.out\n", "mesa_grid.38800304_82.out\n", "mesa_grid.38800304_83.out\n", "mesa_grid.38800304_84.out\n", "mesa_grid.38800304_85.out\n", "mesa_grid.38800304_86.out\n", "mesa_grid.38800304_87.out\n", "mesa_grid.38800304_88.out\n", "mesa_grid.38800304_89.out\n", "mesa_grid.38800304_8.out\n", "mesa_grid.38800304_90.out\n", "mesa_grid.38800304_91.out\n", "mesa_grid.38800304_92.out\n", "mesa_grid.38800304_93.out\n", "mesa_grid.38800304_94.out\n", "mesa_grid.38800304_95.out\n", "mesa_grid.38800304_96.out\n", "mesa_grid.38800304_97.out\n", "mesa_grid.38800304_98.out\n", "mesa_grid.38800304_99.out\n", "mesa_grid.38800304_9.out\n", "mk\n", "run_grid.sh\n", "running_a_grid.ipynb\n", "star1\n", "star2\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_1.0000e-01_grid_index_70\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_1.0000e+04_grid_index_79\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_1.2915e+00_grid_index_72\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_1.6681e+01_grid_index_74\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_2.1544e+02_grid_index_76\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_2.7826e+03_grid_index_78\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_3.5938e-01_grid_index_71\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_4.6416e+00_grid_index_73\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_5.9948e+01_grid_index_75\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_7.7426e+02_grid_index_77\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_1.0000e-01_grid_index_20\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_1.0000e+04_grid_index_29\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_1.2915e+00_grid_index_22\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_1.6681e+01_grid_index_24\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_2.1544e+02_grid_index_26\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_2.7826e+03_grid_index_28\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_3.5938e-01_grid_index_21\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_4.6416e+00_grid_index_23\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_5.9948e+01_grid_index_25\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_7.7426e+02_grid_index_27\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_1.0000e-01_grid_index_80\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_1.0000e+04_grid_index_89\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_1.2915e+00_grid_index_82\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_1.6681e+01_grid_index_84\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_2.1544e+02_grid_index_86\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_2.7826e+03_grid_index_88\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_3.5938e-01_grid_index_81\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_4.6416e+00_grid_index_83\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_5.9948e+01_grid_index_85\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_7.7426e+02_grid_index_87\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_1.0000e-01_grid_index_30\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_1.0000e+04_grid_index_39\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_1.2915e+00_grid_index_32\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_1.6681e+01_grid_index_34\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_2.1544e+02_grid_index_36\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_2.7826e+03_grid_index_38\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_3.5938e-01_grid_index_31\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_4.6416e+00_grid_index_33\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_5.9948e+01_grid_index_35\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_7.7426e+02_grid_index_37\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_1.0000e-01_grid_index_90\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_1.0000e+04_grid_index_99\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_1.2915e+00_grid_index_92\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_1.6681e+01_grid_index_94\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_2.1544e+02_grid_index_96\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_2.7826e+03_grid_index_98\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_3.5938e-01_grid_index_91\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_4.6416e+00_grid_index_93\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_5.9948e+01_grid_index_95\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_7.7426e+02_grid_index_97\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_1.0000e-01_grid_index_40\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_1.0000e+04_grid_index_49\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_1.2915e+00_grid_index_42\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_1.6681e+01_grid_index_44\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_2.1544e+02_grid_index_46\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_2.7826e+03_grid_index_48\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_3.5938e-01_grid_index_41\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_4.6416e+00_grid_index_43\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_5.9948e+01_grid_index_45\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_7.7426e+02_grid_index_47\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_1.0000e-01_grid_index_50\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_1.0000e+04_grid_index_59\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_1.2915e+00_grid_index_52\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_1.6681e+01_grid_index_54\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_2.1544e+02_grid_index_56\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_2.7826e+03_grid_index_58\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_3.5938e-01_grid_index_51\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_4.6416e+00_grid_index_53\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_5.9948e+01_grid_index_55\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_7.7426e+02_grid_index_57\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_1.0000e-01_grid_index_0\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_1.0000e+04_grid_index_9\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_1.2915e+00_grid_index_2\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_1.6681e+01_grid_index_4\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_2.1544e+02_grid_index_6\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_2.7826e+03_grid_index_8\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_3.5938e-01_grid_index_1\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_4.6416e+00_grid_index_3\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_5.9948e+01_grid_index_5\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_7.7426e+02_grid_index_7\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_1.0000e-01_grid_index_60\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_1.0000e+04_grid_index_69\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_1.2915e+00_grid_index_62\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_1.6681e+01_grid_index_64\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_2.1544e+02_grid_index_66\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_2.7826e+03_grid_index_68\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_3.5938e-01_grid_index_61\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_4.6416e+00_grid_index_63\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_5.9948e+01_grid_index_65\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_7.7426e+02_grid_index_67\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_1.0000e-01_grid_index_10\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_1.0000e+04_grid_index_19\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_1.2915e+00_grid_index_12\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_1.6681e+01_grid_index_14\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_2.1544e+02_grid_index_16\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_2.7826e+03_grid_index_18\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_3.5938e-01_grid_index_11\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_4.6416e+00_grid_index_13\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_5.9948e+01_grid_index_15\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_7.7426e+02_grid_index_17\n" ] } ], "source": [ "!ls" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Congratulations, you have run a MESA grid using POSYDON! You can now process the data using POSYDON's built-in tools, or you can use your own tools to process the data." ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.7" } }, "nbformat": 4, "nbformat_minor": 4 }