group8_code
Differences
This shows you the differences between two versions of the page.
|
group8_code [2014/06/06 11:06] bartl created |
group8_code [2014/06/06 11:15] (current) bartl added link to code and a few details |
||
|---|---|---|---|
| Line 1: | Line 1: | ||
| ====== Battle of the Winds: Inside the Code ====== | ====== Battle of the Winds: Inside the Code ====== | ||
| - | Other than the indat file, we modified/added four files in VH1/src/Serial . These can be [[|downloaded here]] and are explained on this page. | + | Other than the indat file, we modified/added four files in VH1/src/Serial . These can be {{::battleofthewinds_code.zip|downloaded here}} and are explained on this page. |
| === init.f90 === | === init.f90 === | ||
| Line 53: | Line 53: | ||
| rshock = 2.24e4 !density of shock material | rshock = 2.24e4 !density of shock material | ||
| </code> | </code> | ||
| + | See our [[group8|logbook]] for details on why we chose these units. | ||
| We then output all parameters to the hst file (code not shown here) and initialize our grid. As we need to stabilize the wind before introducing the shock, we only set up the solar wind source and the ambient medium here. As the solar wind is supposed to radiate spherically, we calculate the angle of the velocity of the solar wind for a given point on the injection radius. | We then output all parameters to the hst file (code not shown here) and initialize our grid. As we need to stabilize the wind before introducing the shock, we only set up the solar wind source and the ambient medium here. As the solar wind is supposed to radiate spherically, we calculate the angle of the velocity of the solar wind for a given point on the injection radius. | ||
| Line 151: | Line 152: | ||
| end module solarwind | end module solarwind | ||
| </code> | </code> | ||
| + | |||
| + | === zonemod.f90 === | ||
| + | As usual, the number of zones should be set according to the chosen geometry of the problem. If the number of zones is too small, the solar wind will be nowhere close to spherical. We ended up using 750x250 zones for our 12x4 AU problem. | ||
group8_code.1402067199.txt.gz · Last modified: 2014/06/06 11:06 by bartl
Page Tools
Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 4.0 International
