{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Run one binary on your laptop for debugging" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "For debugging purposes you might want to run a single binary on your laptop to test a quick modification of the inlists, run_star_extras.f or run_binary_extras.f. Assume we want to run the same binary as in the previous tutorial. To use the POSYDON API to run a single binary, you can use the following command:\n", "\n", "```bash\n", "posydon-setup-grid --grid-type fixed --inifile HMS-HMS.ini --submission-type shell\n", "```\n", "\n", "Similar to before, in our directorty we have:" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "binary\t\t grid_test.csv\t\tlaptop.ipynb star1\n", "grid_command.sh HMS-HMS_yggdrasil.ini\tmk\t star2\n" ] } ], "source": [ "!ls" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The `grid_command.sh` will run all the binaries in `gird_test.csv`. In case you like to run only one specific simulation you can either remove all other initial conditions from the `gird_test.csv` before running `posydon-setup-grid` (You can rerun it whenever you made changes). Or you can modify the loop in `grid_command.sh` by removing all the indexes you do not like to run.\n", "\n", "We can now run the MESA binary simulation in the terminal with the following command:\n", "\n", "```bash\n", "./grid_command.sh\n", "``````\n", "\n", "The terminal will not display the MESA output, the screen output of MESA is saved to an `out.txt` file in the simulation directory." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "binary\n", "column_lists\n", "grid_command.sh\n", "grid_test.csv\n", "HMS-HMS.ini\n", "laptop.ipynb\n", "mk\n", "star1\n", "star2\n", "Zbase_0.0014_m1_30.0000_m2_21.0000_initial_z_1.4200e-03_initial_period_in_days_1.0000e+01_grid_index_0\n" ] } ], "source": [ "!ls" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "binary_history_columns.list inlist_grid_star1_binary_controls\n", "binary_history.data.gz\t inlist_grid_star2_binary_controls\n", "final_star1.mod.gz\t LOGS1\n", "history_columns.list\t LOGS2\n", "initial_star1.mod.gz\t out.txt.gz\n", "initial_star2.mod.gz\t profile_columns.list\n", "inlist\t\t\t tmp.hdf5\n", "inlist_grid_points\n" ] } ], "source": [ "!ls Zbase_0.0014_m1_30.0000_m2_21.0000_initial_z_1.4200e-03_initial_period_in_days_1.0000e+01_grid_index_0" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "\n", " read /home/users/k/kruckow/tutorial_run_MESA_2/binary/inlist_project\n", " read inlist_grid_points\n", " version_number 11701\n", " read inlist_grid_star1_binary_controls\n", " set_eos_PC_parameters\n", " mass_fraction_limit_for_PC 1.0000000000000000D-03\n", " logRho1_PC_limit 2.9990000000000001D+00\n", " logRho2_PC_limit 2.7999999999999998D+00\n", " log_Gamma_all_HELM 1.0000000000000000D+00\n", " log_Gamma_all_PC 1.3010299956000000D+00\n", " PC_Gamma_start_crystal 1.5000000000000000D+02\n", " PC_Gamma_full_crystal 1.7500000000000000D+02\n", " PC_min_Z 9.9900000000000000D-01\n", " change rates preference to 2\n", " set_initial_age 0.0000000000000000D+00\n", " set_initial_model_number 0\n", " change to \"approx21.net\"\n", " number of species 21\n", " new_rotation_flag T\n", " new_surface_rotation_v 0 0.0000000000000000D+00 0.0000000000000000D+00\n", " net name approx21.net\n", " rotation_flag T\n", " species\n", " 1 neut \n", " 2 h1 \n", " 3 prot \n", " 4 he3 \n", " 5 he4 \n", " 6 c12 \n", " 7 n14 \n", " 8 o16 \n", " 9 ne20 \n", " 10 mg24 \n", " 11 si28 \n", " 12 s32 \n", " 13 ar36 \n", " 14 ca40 \n", " 15 ti44 \n", " 16 cr48 \n", " 17 cr56 \n", " 18 fe52 \n", " 19 fe54 \n", " 20 fe56 \n", " 21 ni56 \n", "\n", " kappa_file_prefix gs98\n", " kappa_lowT_prefix lowT_fa05_gs98\n", " OMP_NUM_THREADS 4\n", "\n", "\n", "\n", "\n", " read inlist_grid_star2_binary_controls\n", " set_eos_PC_parameters\n", " mass_fraction_limit_for_PC 1.0000000000000000D-03\n", " logRho1_PC_limit 2.9990000000000001D+00\n", " logRho2_PC_limit 2.7999999999999998D+00\n", " log_Gamma_all_HELM 1.0000000000000000D+00\n", " log_Gamma_all_PC 1.3010299956000000D+00\n", " PC_Gamma_start_crystal 1.5000000000000000D+02\n", " PC_Gamma_full_crystal 1.7500000000000000D+02\n", " PC_min_Z 9.9900000000000000D-01\n", " change rates preference to 2\n", " set_initial_age 0.0000000000000000D+00\n", " set_initial_model_number 0\n", " change to \"approx21.net\"\n", " number of species 21\n", " new_rotation_flag T\n", " new_surface_rotation_v 0 0.0000000000000000D+00 0.0000000000000000D+00\n", " net name approx21.net\n", " rotation_flag T\n", " species\n", " 1 neut \n", " 2 h1 \n", " 3 prot \n", " 4 he3 \n", " 5 he4 \n", " 6 c12 \n", " 7 n14 \n", " 8 o16 \n", " 9 ne20 \n", " 10 mg24 \n", " 11 si28 \n", " 12 s32 \n", " 13 ar36 \n", " 14 ca40 \n", " 15 ti44 \n", " 16 cr48 \n", " 17 cr56 \n", " 18 fe52 \n", " 19 fe54 \n", " 20 fe56 \n", " 21 ni56 \n", "\n", " kappa_file_prefix gs98\n", " kappa_lowT_prefix lowT_fa05_gs98\n", " OMP_NUM_THREADS 4\n", "\n", "\n", "\n", "\n", "\n", " m2 2.1000000000000000D+01\n", " m1 3.0000000000000000D+01\n", " initial_period_in_days 1.0000000000000000D+01\n", " initial_separation_in_Rsun 7.2417496333000003D+01\n", " jdot_multiplier 1.0000000000000000D+00\n", " fr 1.0000000000000000D-02\n", "\n", "\n", "\n", " The binary terminal output contains the following information\n", "\n", " 'step' is the number of steps since the start of the run,\n", " 'lg_dt' is log10 timestep in years,\n", " 'age_yr' is the simulated years since the start run,\n", " 'M1+M2' is the total mass of the system (Msun),\n", " 'M1' is the mass of the primary (Msun)\n", " 'M2' is the mass of the secondary (Msun)\n", " 'separ' is the semi-major axis of the orbit (Rsun),\n", " 'R1' is the radius of the primary (Rsun)\n", " 'R2' is the radius of the secondary (Rsun)\n", " 'Porb' is the orbital period (days),\n", " 'P1' is the rotation period of star 1 (days, zero if not modeling rotation),\n", " 'P2' is the rotation period of star 2 (days, zero if not modeling rotation),\n", " 'e' orbital eccentricity,\n", " 'dot_e' time derivative of e (1/yr),\n", " 'Eorb' orbital energy G*M1*M2/2*separation (ergs),\n", " 'M2/M1' mass ratio,\n", " 'vorb1' orbital velocity of star 1 (km/s),\n", " 'vorb2' orbital velocity of star 2 (km/s),\n", " 'pm_i' index of star evolved as point mass, zero if both stars are modeled,\n", " 'RL1' Roche lobe radius of star 1 (Rsun),\n", " 'Rl2' Roche lobe radius of star 2 (Rsun),\n", " 'donor_i' index of star taken as donor,\n", " 'RL_gap1' (R1-Rl1)/Rl1,\n", " 'RL_gap2' (R2-Rl2)/Rl2,\n", " 'dot_Mmt', mass transfer rate (Msun/yr),\n", " 'dot_M1', time derivative for the mass of star 1 (Msun/yr),\n", " 'dot_M2', time derivative for the mass of star 2 (Msun/yr),\n", " 'eff', mass transfer efficiency, computed as -dot_M2/dot_M1 (zero if dot_M1=0),\n", " 'dot_Medd', Eddington accretion rate (Msun/yr),\n", " 'L_acc', accretion luminosity when accreting to a point mass (ergs/s),\n", " 'Jorb', orbital angular momentum (g*cm^2/s)\n", " 'spin1', spin angular momentum of star 1 (g*cm^2/s),\n", " 'spin2', spin angular momentum of star 2 (g*cm^2/s),\n", " 'dot_J', time derivative of Jorb (g*cm^2/s^2),\n", " 'dot_Jgr', time derivative of Jorb due to gravitational waves (g*cm^2/s^2),\n", " 'dot_Jml', time derivative of Jorb due to mass loss (g*cm^2/s^2),\n", " 'dot_Jmb', time derivative of Jorb due to magnetic braking (g*cm^2/s^2),\n", " 'dot_Jls', time derivative of Jorb due to spin-orbit coupling (g*cm^2/s^2),\n", " 'rlo_iters', number of iterations for implicit calculation of mass transfer,\n", "\n", " All this and more can be saved in binary_history.data during the run.\n", " num_steps_to_relax_rotation 50\n", " relax to omega: wanted-current, current, wanted 1 7.2722052166430393D-06 0.0000000000000000D+00 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 2 7.1239370299804958D-06 1.4826818666254330D-07 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 3 6.9775031634253124D-06 2.9470205321772722D-07 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 4 6.8325227781707742D-06 4.3968243847226474D-07 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 5 6.6880336129568242D-06 5.8417160368621492D-07 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 6 6.5429721100078805D-06 7.2923310663515874D-07 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 7 6.3960414375365654D-06 8.7616377910647371D-07 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 8 6.2481184946104557D-06 1.0240867220325836D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 9 6.0988143280937228D-06 1.1733908885493166D-06 7.2722052166430393D-06\n", " 10 7.649681 4.541E+04 5.083637 5.083787 30.000000 30.000000 0.746765 0.000508 0.251100 0.747141 -6.159624 2598 0\n", " 3.021833 0.768039 0.749892 -23.011310 3.935137 -99.000000 0.000000 0.251795 0.000441 0.001420 0.251427 0.009963 3 0\n", " 2.3543E+03 16.667399 5.081908 1.622938 4.361600 -8.707435 0.000000 0.000004 0.000133 0.000071 1.432E-03 0.000E+00 max increase\n", "\n", " relax to omega: wanted-current, current, wanted 10 5.9494920708119283D-06 1.3227131458311106D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 11 5.8047898071657684D-06 1.4674154094772709D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 12 5.6656471713199783D-06 1.6065580453230609D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 13 5.5253851216154814D-06 1.7468200950275575D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 14 5.3814004626389407D-06 1.8908047540040986D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 15 5.2360336852655910D-06 2.0361715313774479D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 16 5.0905570727440271D-06 2.1816481438990118D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 17 4.9451037265468788D-06 2.3271014900961605D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 18 4.7996593714595307D-06 2.4725458451835086D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 19 4.6542152930004046D-06 2.6179899236426347D-06 7.2722052166430393D-06\n", " 20 7.648797 4.525E+04 5.069089 5.069239 30.000000 30.000000 0.746765 0.000508 0.251100 0.747141 -6.164975 2613 0\n", " 6.032133 0.764405 0.746596 -23.057889 3.920548 -99.000000 0.000000 0.251795 0.000441 0.001420 0.251427 0.009955 2 0\n", " 2.1538E+06 16.663052 5.069239 1.607921 4.359365 -8.704639 0.000000 0.000004 0.000133 0.000071 1.432E-03 0.000E+00 max increase\n", "\n", " relax to omega: wanted-current, current, wanted 20 4.5087711966205017D-06 2.7634340200225373D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 21 4.3633268036162395D-06 2.9088784130267999D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 22 4.2178823416292592D-06 3.0543228750137802D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 23 4.0724377954490792D-06 3.1997674211939601D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 24 3.9269931389606084D-06 3.3452120776824310D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 25 3.7815483392570430D-06 3.4906568773859964D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 26 3.6361034724003977D-06 3.6361017442426417D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 27 3.4906587782149381D-06 3.7815464384281013D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 28 3.3452144143481404D-06 3.9269908022948989D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 29 3.1997702486009012D-06 4.0724349680421382D-06 7.2722052166430393D-06\n", " 30 7.648796 4.525E+04 5.069078 5.069228 30.000000 30.000000 0.746765 0.000508 0.251100 0.747141 -6.164960 2610 0\n", " 9.042433 0.764409 0.746606 -23.057923 3.920536 -99.000000 0.000000 0.251795 0.000441 0.001420 0.251427 0.009955 2 0\n", " 2.2053E+09 16.663054 5.069228 1.607910 4.359335 -8.704658 0.000000 0.000004 0.000133 0.000071 1.432E-03 0.000E+00 max increase\n", "\n", " relax to omega: wanted-current, current, wanted 30 3.0543261358872070D-06 4.2178790807558323D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 31 2.9088820315467191D-06 4.3633231850963203D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 32 2.7634379268261852D-06 4.5087672898168542D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 33 2.6179938218649762D-06 4.6542113947780632D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 34 2.4725497163505801D-06 4.7996555002924592D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 35 2.3271056116451454D-06 4.9450996049978939D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 36 2.1816615088007590D-06 5.0905437078422803D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 37 2.0362174007229542D-06 5.2359878159200852D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 38 1.8907732940142874D-06 5.3814319226287519D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 39 1.7453291891721467D-06 5.5268760274708927D-06 7.2722052166430393D-06\n", " 40 7.648795 4.525E+04 5.069062 5.069212 30.000000 30.000000 0.746765 0.000508 0.251100 0.747141 -6.164940 2607 0\n", " 12.052733 0.764417 0.746620 -23.057971 3.920520 -99.000000 0.000000 0.251795 0.000441 0.001420 0.251427 0.009955 2 0\n", " 2.2582E+12 16.663057 5.069212 1.607894 4.359294 -8.704686 0.000000 0.000004 0.000133 0.000071 1.432E-03 0.000E+00 max increase\n", "\n", " relax to omega: wanted-current, current, wanted 40 1.5998850831994643D-06 5.6723201334435750D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 41 1.4544409770744422D-06 5.8177642395685971D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 42 1.3089968705342554D-06 5.9632083461087839D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 43 1.1635527639419701D-06 6.1086524527010692D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 44 1.0181086591042797D-06 6.2540965575387596D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 45 8.7266455123154744D-07 6.3995406654114919D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 46 7.2722044498914098D-07 6.5449847716538984D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 47 5.8177633970541403D-07 6.6904288769376253D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 48 4.3633223233433363D-07 6.8358729843087057D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 49 2.9088812460997866D-07 6.9813170920330607D-06 7.2722052166430393D-06\n", " 50 7.648794 4.525E+04 5.069042 5.069192 30.000000 30.000000 0.746765 0.000508 0.251100 0.747141 -6.164913 2607 0\n", " 15.063033 0.764426 0.746639 -23.058033 3.920500 -99.000000 0.000000 0.251795 0.000441 0.001420 0.251427 0.009955 2 0\n", " 2.3124E+15 16.663062 5.069192 1.607873 4.359240 -8.704722 0.000000 0.000004 0.000133 0.000071 1.432E-03 0.000E+00 max increase\n", "\n", " relax to omega: wanted-current, current, wanted 50 1.4544401957845492D-07 7.1267611970645844D-06 7.2722052166430393D-06\n", " final step: wanted-current, current, wanted 51 -8.8872453572253797D-14 7.2722053055154929D-06 7.2722052166430393D-06\n", " num_steps_to_relax_rotation 50\n", " relax to omega: wanted-current, current, wanted 1 7.2722052166430393D-06 0.0000000000000000D+00 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 2 7.1233476044913972D-06 1.4885761215164192D-07 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 3 6.9765158364664393D-06 2.9568938017660038D-07 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 4 6.8315785028308242D-06 4.4062671381221501D-07 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 5 6.6872664451027811D-06 5.8493877154025819D-07 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 6 6.5422553499448946D-06 7.2994986669814492D-07 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 7 6.3946688756986521D-06 8.7753634094438770D-07 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 8 6.2457320323138947D-06 1.0264731843291449D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 9 6.0961527996288285D-06 1.1760524170142110D-06 7.2722052166430393D-06\n", " 2 10 7.626746 4.039E+04 4.702123 4.702274 21.000000 21.000000 0.746746 0.000504 0.251100 0.747186 -5.802804 2094 0\n", " 3.397968 0.888883 0.660278 -24.009890 3.553286 -99.000000 0.000000 0.251813 0.000447 0.001420 0.251381 0.011524 3 0\n", " 5.6097E+03 16.729885 4.698975 1.243311 4.233606 -8.722599 0.000000 0.000004 0.000133 0.000071 1.433E-03 0.000E+00 max increase\n", "\n", " relax to omega: wanted-current, current, wanted 10 5.9497664059555575D-06 1.3224388106874818D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 11 5.8096280128646104D-06 1.4625772037784294D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 12 5.6702121385033817D-06 1.6019930781396575D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 13 5.5268762195040200D-06 1.7453289971390194D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 14 5.3814893418133806D-06 1.8907158748296590D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 15 5.2359962060341824D-06 2.0362090106088570D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 16 5.0905469533439396D-06 2.1816582632991002D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 17 4.9451025348753598D-06 2.3271026817676795D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 18 4.7996593159162278D-06 2.4725459007268115D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 19 4.6542137048352694D-06 2.6179915118077704D-06 7.2722052166430393D-06\n", " 2 20 7.626057 4.027E+04 4.690568 4.690719 21.000000 21.000000 0.746746 0.000504 0.251100 0.747186 -5.806944 2142 0\n", " 6.408268 0.886095 0.658647 -24.047363 3.541695 -99.000000 0.000000 0.251813 0.000447 0.001420 0.251381 0.011517 2 0\n", " 5.1209E+06 16.726495 4.690719 1.231378 4.231839 -8.720788 0.000000 0.000004 0.000133 0.000071 1.433E-03 0.000E+00 max increase\n", "\n", " relax to omega: wanted-current, current, wanted 20 4.5087702091765951D-06 2.7634350074664442D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 21 4.3633258925594751D-06 2.9088793240835646D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 22 4.2178814610967289D-06 3.0543237555463109D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 23 4.0724368772888120D-06 3.1997683393542269D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 24 3.9269922854427176D-06 3.3452129312003218D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 25 3.7815476480289336D-06 3.4906575686141058D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 26 3.6361029746011353D-06 3.6361022420419041D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 27 3.4906584885813882D-06 3.7815467280616511D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 28 3.3452142083025420D-06 3.9269910083404973D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 29 3.1997700535865694D-06 4.0724351630564699D-06 7.2722052166430393D-06\n", " 2 30 7.626056 4.027E+04 4.690559 4.690709 21.000000 21.000000 0.746746 0.000504 0.251100 0.747186 -5.806932 2161 0\n", " 9.418568 0.886099 0.658656 -24.047394 3.541686 -99.000000 0.000000 0.251813 0.000447 0.001420 0.251381 0.011517 2 0\n", " 5.2432E+09 16.726497 4.690709 1.231368 4.231817 -8.720799 0.000000 0.000004 0.000133 0.000071 1.433E-03 0.000E+00 max increase\n", "\n", " relax to omega: wanted-current, current, wanted 30 3.0543259269715730D-06 4.2178792896714663D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 31 2.9088818039901709D-06 4.3633234126528685D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 32 2.7634376805293879D-06 4.5087675361136514D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 33 2.6179935564114802D-06 4.6542116602315591D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 34 2.4725494315840545D-06 4.7996557850589848D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 35 2.3271053084760396D-06 4.9450999081669997D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 36 2.1816611779797231D-06 5.0905440386633162D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 37 2.0362170531092207D-06 5.2359881635338187D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 38 1.8907729258279075D-06 5.3814322908151319D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 39 1.7453287980464179D-06 5.5268764185966214D-06 7.2722052166430393D-06\n", " 2 40 7.626055 4.027E+04 4.690545 4.690696 21.000000 21.000000 0.746746 0.000504 0.251100 0.747186 -5.806914 2153 0\n", " 12.428868 0.886105 0.658667 -24.047437 3.541672 -99.000000 0.000000 0.251813 0.000447 0.001420 0.251381 0.011517 2 0\n", " 5.3691E+12 16.726501 4.690696 1.231354 4.231787 -8.720815 0.000000 0.000004 0.000133 0.000071 1.433E-03 0.000E+00 max increase\n", "\n", " relax to omega: wanted-current, current, wanted 40 1.5998846692471412D-06 5.6723205473958981D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 41 1.4544405402861448D-06 5.8177646763568946D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 42 1.3089964117734896D-06 5.9632088048695497D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 43 1.1635522807878633D-06 6.1086529358551760D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 44 1.0181081448464963D-06 6.2540970717965430D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 45 8.7266401721243792D-07 6.3995411994306014D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 46 7.2721988481735897D-07 6.5449853318256804D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 47 5.8177575181154132D-07 6.6904294648314980D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 48 4.3633161713273364D-07 6.8358735995103057D-06 7.2722052166430393D-06\n", " relax to omega: wanted-current, current, wanted 49 2.9088748859537223D-07 6.9813177280476671D-06 7.2722052166430393D-06\n", " 2 50 7.626054 4.027E+04 4.690527 4.690678 21.000000 21.000000 0.746746 0.000504 0.251100 0.747186 -5.806891 2153 0\n", " 15.439168 0.886113 0.658682 -24.047492 3.541654 -99.000000 0.000000 0.251813 0.000447 0.001420 0.251381 0.011517 2 0\n", " 5.4979E+15 16.726505 4.690678 1.231337 4.231748 -8.720836 0.000000 0.000004 0.000133 0.000071 1.433E-03 0.000E+00 max increase\n", "\n", " relax to omega: wanted-current, current, wanted 50 1.4544334521563968D-07 7.1267618714273997D-06 7.2722052166430393D-06\n", " final step: wanted-current, current, wanted 51 -7.8765733440910391D-13 7.2722060043003737D-06 7.2722052166430393D-06\n", " 1 7.648422 4.523E+04 5.343064 5.343138 29.999999 29.999999 0.746869 0.000476 0.251100 0.747139 -6.166319 2606 0\n", " 1.307197 0.763236 0.746854 -23.077820 4.293714 -7.350350 0.000000 0.251689 0.000467 0.001420 0.251428 0.009954 4 0\n", " 2.0286E+01 16.661525 5.068843 1.573581 4.358537 -8.705162 0.000000 0.000013 0.000133 0.000071 1.433E-03 0.000E+00 max increase\n", "\n", " 2 1 7.625915 4.026E+04 4.849915 4.850013 21.000000 21.000000 0.746868 0.000475 0.251100 0.747186 -5.807440 2125 0\n", " 1.307197 0.885645 0.658746 -24.055524 3.785636 -7.994602 0.000000 0.251688 0.000475 0.001420 0.251381 0.011516 5 0\n", " 2.0286E+01 16.725936 4.690535 1.203320 4.231541 -8.720928 0.000000 0.000008 0.000133 0.000071 1.433E-03 0.000E+00 b_companion\n", "\n", " saving initial models\n", " 10 7.647466 4.517E+04 5.017842 5.017992 29.999977 29.999977 0.746916 0.000476 0.251100 0.747118 -6.169780 2531 0\n", " 2.019828 0.760299 0.747480 -23.127544 3.892503 -7.353486 0.000000 0.251639 0.000479 0.001420 0.251446 0.009954 4 0\n", " 5.2660E+02 16.657613 5.067885 1.556907 4.356530 -8.706430 0.000000 0.000004 0.000133 0.000071 1.436E-03 0.000E+00 max increase\n", "\n", " 2 10 7.625061 4.021E+04 4.645284 4.645433 20.999995 20.999995 0.746937 0.000465 0.251100 0.747172 -5.811002 2175 0\n", " 2.019828 0.882850 0.659232 -24.101626 3.538730 -7.998093 0.000000 0.251616 0.000491 0.001420 0.251393 0.011514 3 0\n", " 5.2660E+02 16.722315 4.689495 1.180081 4.229992 -8.721628 0.000000 0.000004 0.000133 0.000071 1.435E-03 0.000E+00 max increase\n", "\n", "bin 10 50.999971 72.417537 10.000011 0.000E+00 0.700000 0 1 0.000E+00 -2.267E-01 4.539E+54 -1.281E+38 0.000E+00\n", " 2.019828 29.999977 5.590882 10.055301 0.000E+00 150.922990 29.701262 -8.118E-01 -4.431E-08 1.000E+99 6.342E+51 -1.398E+35 -7.998E+33\n", " 5.2660E+02 20.999995 4.562810 10.036142 -1.651E+49 215.604158 25.237109 -8.192E-01 -1.004E-08 0.000E+00 2.880E+51 -1.280E+38 1\n", "\n", "gzip: stdout: Broken pipe\n" ] } ], "source": [ "!gzip -cd Zbase_0.0014_m1_30.0000_m2_21.0000_initial_z_1.4200e-03_initial_period_in_days_1.0000e+01_grid_index_0/out.txt.gz | head -n 320" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Congratulation, you now know how to harvast the power of POSYDON to run a grid of simulations on a supercomputer!" ] } ], "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 }