star_destroyers
Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
|
star_destroyers [2014/06/06 13:40] richers |
star_destroyers [2014/06/06 13:47] richers |
||
|---|---|---|---|
| Line 30: | Line 30: | ||
| * E_blast = 10^49 erg (all placed within innermost cell) | * E_blast = 10^49 erg (all placed within innermost cell) | ||
| + | |||
| + | **Making VH1 Starter Lagrangian** | ||
| + | |||
| + | ppmlr.f90 needs to be modified so that once the Lagrange step was taken, the values needed to not be put back onto the Eulerian grid. | ||
| + | * remove the "wiggle" part of the file that adds dissipation for carbuncular noise. Since this is only really a problem in higher dimension, we completely removed this for simplicity. | ||
| + | * Instead of defining xf=xf0 define xf0=xf (i.e. make the Eulerian coordinates move to match the Lagrangian ones) | ||
| + | * If the final "remap" command is commented out, the code misbehaves, though we don't know why a remap should be needed. Leave the final remap in. | ||
| **Units** | **Units** | ||
| Line 42: | Line 49: | ||
| * E_cgs = M_sun*R_sun^2 * E_code | * E_cgs = M_sun*R_sun^2 * E_code | ||
| - | **XNet Parameters** | + | **XNet Parameters and Inputs** |
| Three separate initial conditions: | Three separate initial conditions: | ||
| Line 48: | Line 55: | ||
| * 50% N, 50% P | * 50% N, 50% P | ||
| * pure 56Fe | * pure 56Fe | ||
| + | |||
| + | To modify XNet to output 150 nuclei in the plaintext file. (Instead of trying to read in the fortran binary files in python, which is in principle also possible.) You need to change: | ||
| + | |||
| + | * The ''control'' file, you must add in the list of all the nuclei you wish to output. Recall each should use 5 characters. | ||
| + | * the ''inout'' variable in ''control.f90''. | ||
| + | * ''output_nuc'' and | ||
| + | * ''Read(lun_control,"(14a5)") output_nuc'' to ''150a5'', | ||
| + | * The index in the loop on line 198, | ||
| + | * The write statement for ''lun_ev''on lines 352, 465, 473 in the file ''net.f90''. (Do a search for 14 in the files). | ||
| + | |||
| + | As usual, all the code is in the github. | ||
| + | |||
| ------Logbook starts here ------- | ------Logbook starts here ------- | ||
| Line 65: | Line 84: | ||
| {{::screen_shot_2014-06-03_at_6.20.49_pm.png?200|}} | {{::screen_shot_2014-06-03_at_6.20.49_pm.png?200|}} | ||
| - | KJC: Modified XNet to output 150 nuclei in the plaintext file. (Instead of trying to read in the fortran binary files in python, which is in principle also possible.) You need to change: | ||
| - | |||
| - | The ''control'' file, you must add in the list of all the nuclei you wish to output. Recall each should use 5 characters. | ||
| - | |||
| - | Also, change the ''inout'' variable in ''control.f90''. | ||
| - | |||
| - | You also need to change: | ||
| - | |||
| - | ''output_nuc'' and | ||
| - | |||
| - | ''Read(lun_control,"(14a5)") output_nuc'' to ''150a5'', | ||
| - | |||
| - | The index in the loop on line 198, | ||
| - | |||
| - | The write statement for ''lun_ev''on lines 352, 465, 473. | ||
| - | |||
| - | in the file ''net.f90''. (Do a search for 14 in the files). | ||
| - | |||
| - | As usual, all the code is in the github. | ||
| Modified ''plotchart.py'' to allow as many isotopes to be plotted as are in the network you're interested in. The script now outputs an ''.mp4'' file for display. | Modified ''plotchart.py'' to allow as many isotopes to be plotted as are in the network you're interested in. The script now outputs an ''.mp4'' file for display. | ||
star_destroyers.txt ยท Last modified: 2014/06/06 14:04 by cook
Page Tools
Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 4.0 International
