{ "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 ini file as in the previous notebook." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Creating the Initialization File to Submit the MESA Grid" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Let's copy the code use in `$PATH_TO_POSYDON/grid_params/POSYDON-MESA-INLISTS/r11701/running_scripts/parameter_space_v2/create_csv.ipynb` to create the csv file and edit it to run a downsampled version of the grid." ] }, { "cell_type": "code", "execution_count": 3, "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", "import pandas as pd\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", "def lin_range(x_min,x_max,x_n):\n", " return np.linspace(x_min,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 = 10**(-1.)\n", "p_max = 6105\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", "for Z in met:\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 i in range(m1_n):\n", " for j in range(q_n):\n", " for k in range(p_n):\n", " if m1[i]*q[j] >= 0.5:\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:\n", "\n", "```bash\n", "posydon-setup-grid --grid-type fixed --inifile HMS-HMS_yggdrasil.ini --submission-type slurm\n", "sbatch slurm_job_array_grid_submit.sh\n", "```\n", "\n", "Sit back relax, and wait for the grid to finish." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "binary\n", "grid.csv\n", "HMS-HMS_yggdrasil.ini\n", "mesa_grid.28454093_0.out\n", "mesa_grid.28454093_10.out\n", "mesa_grid.28454093_11.out\n", "mesa_grid.28454093_12.out\n", "mesa_grid.28454093_13.out\n", "mesa_grid.28454093_14.out\n", "mesa_grid.28454093_15.out\n", "mesa_grid.28454093_16.out\n", "mesa_grid.28454093_17.out\n", "mesa_grid.28454093_18.out\n", "mesa_grid.28454093_19.out\n", "mesa_grid.28454093_1.out\n", "mesa_grid.28454093_20.out\n", "mesa_grid.28454093_21.out\n", "mesa_grid.28454093_22.out\n", "mesa_grid.28454093_23.out\n", "mesa_grid.28454093_24.out\n", "mesa_grid.28454093_25.out\n", "mesa_grid.28454093_26.out\n", "mesa_grid.28454093_27.out\n", "mesa_grid.28454093_28.out\n", "mesa_grid.28454093_29.out\n", "mesa_grid.28454093_2.out\n", "mesa_grid.28454093_30.out\n", "mesa_grid.28454093_31.out\n", "mesa_grid.28454093_32.out\n", "mesa_grid.28454093_33.out\n", "mesa_grid.28454093_34.out\n", "mesa_grid.28454093_35.out\n", "mesa_grid.28454093_36.out\n", "mesa_grid.28454093_37.out\n", "mesa_grid.28454093_38.out\n", "mesa_grid.28454093_39.out\n", "mesa_grid.28454093_3.out\n", "mesa_grid.28454093_40.out\n", "mesa_grid.28454093_41.out\n", "mesa_grid.28454093_42.out\n", "mesa_grid.28454093_43.out\n", "mesa_grid.28454093_44.out\n", "mesa_grid.28454093_45.out\n", "mesa_grid.28454093_46.out\n", "mesa_grid.28454093_47.out\n", "mesa_grid.28454093_48.out\n", "mesa_grid.28454093_49.out\n", "mesa_grid.28454093_4.out\n", "mesa_grid.28454093_50.out\n", "mesa_grid.28454093_51.out\n", "mesa_grid.28454093_52.out\n", "mesa_grid.28454093_53.out\n", "mesa_grid.28454093_54.out\n", "mesa_grid.28454093_55.out\n", "mesa_grid.28454093_56.out\n", "mesa_grid.28454093_57.out\n", "mesa_grid.28454093_58.out\n", "mesa_grid.28454093_59.out\n", "mesa_grid.28454093_5.out\n", "mesa_grid.28454093_60.out\n", "mesa_grid.28454093_61.out\n", "mesa_grid.28454093_62.out\n", "mesa_grid.28454093_63.out\n", "mesa_grid.28454093_64.out\n", "mesa_grid.28454093_65.out\n", "mesa_grid.28454093_66.out\n", "mesa_grid.28454093_67.out\n", "mesa_grid.28454093_68.out\n", "mesa_grid.28454093_69.out\n", "mesa_grid.28454093_6.out\n", "mesa_grid.28454093_70.out\n", "mesa_grid.28454093_71.out\n", "mesa_grid.28454093_72.out\n", "mesa_grid.28454093_73.out\n", "mesa_grid.28454093_74.out\n", "mesa_grid.28454093_75.out\n", "mesa_grid.28454093_76.out\n", "mesa_grid.28454093_77.out\n", "mesa_grid.28454093_78.out\n", "mesa_grid.28454093_79.out\n", "mesa_grid.28454093_7.out\n", "mesa_grid.28454093_80.out\n", "mesa_grid.28454093_81.out\n", "mesa_grid.28454093_82.out\n", "mesa_grid.28454093_83.out\n", "mesa_grid.28454093_84.out\n", "mesa_grid.28454093_85.out\n", "mesa_grid.28454093_86.out\n", "mesa_grid.28454093_87.out\n", "mesa_grid.28454093_88.out\n", "mesa_grid.28454093_89.out\n", "mesa_grid.28454093_8.out\n", "mesa_grid.28454093_90.out\n", "mesa_grid.28454093_91.out\n", "mesa_grid.28454093_92.out\n", "mesa_grid.28454093_93.out\n", "mesa_grid.28454093_94.out\n", "mesa_grid.28454093_95.out\n", "mesa_grid.28454093_96.out\n", "mesa_grid.28454093_97.out\n", "mesa_grid.28454093_98.out\n", "mesa_grid.28454093_99.out\n", "mesa_grid.28454093_9.out\n", "mk\n", "running_a_grid.ipynb\n", "slurm_job_array_grid_submit.sh\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.1574e+00_grid_index_72\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_1.3396e+01_grid_index_74\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_1.5504e+02_grid_index_76\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_1.7945e+03_grid_index_78\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_3.4021e-01_grid_index_71\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_3.9376e+00_grid_index_73\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_4.5574e+01_grid_index_75\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_5.2747e+02_grid_index_77\n", "Zbase_0.0014_m1_123.3602_m2_86.3521_initial_z_1.4200e-03_initial_period_in_days_6.1050e+03_grid_index_79\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.1574e+00_grid_index_22\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_1.3396e+01_grid_index_24\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_1.5504e+02_grid_index_26\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_1.7945e+03_grid_index_28\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_3.4021e-01_grid_index_21\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_3.9376e+00_grid_index_23\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_4.5574e+01_grid_index_25\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_5.2747e+02_grid_index_27\n", "Zbase_0.0014_m1_13.3755_m2_9.3628_initial_z_1.4200e-03_initial_period_in_days_6.1050e+03_grid_index_29\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.1574e+00_grid_index_82\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_1.3396e+01_grid_index_84\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_1.5504e+02_grid_index_86\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_1.7945e+03_grid_index_88\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_3.4021e-01_grid_index_81\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_3.9376e+00_grid_index_83\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_4.5574e+01_grid_index_85\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_5.2747e+02_grid_index_87\n", "Zbase_0.0014_m1_192.3748_m2_134.6623_initial_z_1.4200e-03_initial_period_in_days_6.1050e+03_grid_index_89\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.1574e+00_grid_index_32\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_1.3396e+01_grid_index_34\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_1.5504e+02_grid_index_36\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_1.7945e+03_grid_index_38\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_3.4021e-01_grid_index_31\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_3.9376e+00_grid_index_33\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_4.5574e+01_grid_index_35\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_5.2747e+02_grid_index_37\n", "Zbase_0.0014_m1_20.8585_m2_14.6009_initial_z_1.4200e-03_initial_period_in_days_6.1050e+03_grid_index_39\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.1574e+00_grid_index_92\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_1.3396e+01_grid_index_94\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_1.5504e+02_grid_index_96\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_1.7945e+03_grid_index_98\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_3.4021e-01_grid_index_91\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_3.9376e+00_grid_index_93\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_4.5574e+01_grid_index_95\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_5.2747e+02_grid_index_97\n", "Zbase_0.0014_m1_300.0000_m2_210.0000_initial_z_1.4200e-03_initial_period_in_days_6.1050e+03_grid_index_99\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.1574e+00_grid_index_42\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_1.3396e+01_grid_index_44\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_1.5504e+02_grid_index_46\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_1.7945e+03_grid_index_48\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_3.4021e-01_grid_index_41\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_3.9376e+00_grid_index_43\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_4.5574e+01_grid_index_45\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_5.2747e+02_grid_index_47\n", "Zbase_0.0014_m1_32.5279_m2_22.7695_initial_z_1.4200e-03_initial_period_in_days_6.1050e+03_grid_index_49\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.1574e+00_grid_index_52\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_1.3396e+01_grid_index_54\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_1.5504e+02_grid_index_56\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_1.7945e+03_grid_index_58\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_3.4021e-01_grid_index_51\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_3.9376e+00_grid_index_53\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_4.5574e+01_grid_index_55\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_5.2747e+02_grid_index_57\n", "Zbase_0.0014_m1_50.7258_m2_35.5080_initial_z_1.4200e-03_initial_period_in_days_6.1050e+03_grid_index_59\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.1574e+00_grid_index_2\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_1.3396e+01_grid_index_4\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_1.5504e+02_grid_index_6\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_1.7945e+03_grid_index_8\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_3.4021e-01_grid_index_1\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_3.9376e+00_grid_index_3\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_4.5574e+01_grid_index_5\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_5.2747e+02_grid_index_7\n", "Zbase_0.0014_m1_5.5000_m2_3.8500_initial_z_1.4200e-03_initial_period_in_days_6.1050e+03_grid_index_9\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.1574e+00_grid_index_62\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_1.3396e+01_grid_index_64\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_1.5504e+02_grid_index_66\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_1.7945e+03_grid_index_68\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_3.4021e-01_grid_index_61\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_3.9376e+00_grid_index_63\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_4.5574e+01_grid_index_65\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_5.2747e+02_grid_index_67\n", "Zbase_0.0014_m1_79.1046_m2_55.3732_initial_z_1.4200e-03_initial_period_in_days_6.1050e+03_grid_index_69\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.1574e+00_grid_index_12\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_1.3396e+01_grid_index_14\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_1.5504e+02_grid_index_16\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_1.7945e+03_grid_index_18\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_3.4021e-01_grid_index_11\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_3.9376e+00_grid_index_13\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_4.5574e+01_grid_index_15\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_5.2747e+02_grid_index_17\n", "Zbase_0.0014_m1_8.5770_m2_6.0039_initial_z_1.4200e-03_initial_period_in_days_6.1050e+03_grid_index_19\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": "posydon_env", "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.5" } }, "nbformat": 4, "nbformat_minor": 2 }