User Tools

Site Tools


group8_code

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

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.txt ยท Last modified: 2014/06/06 11:15 by bartl