| 108 | === Program Structure === |
| 109 | |
| 110 | 1. Call init -- Sets up initial conditions. This subroutine (SR) loops through i=1,mx and checks whether the cell center ((i-0.5)*dx) is to left or right of discontinuity (x0), and then assigns appropriate value to that cell. (All initial conditions for the tests in Toro consist of a left and right state separated by a discontinuity. |
| 111 | 2. Call output -- begins the output file -- records the following info - position, density, velocity, pressure, and time and timestep. This output file is named output.out. |
| 112 | 3. Sets time = 0.0 -- starts the clock for iteration scheme |
| 113 | 4. Enters the update cycle, where from n=1 to n=maxitr (n is the counter for the time iteration scheme) it performs the following commands - |
| 114 | |
| 115 | i. Call Bcond -- Boundary condition SR that sets left and right ghost zones to adjacent physical cells (recall boundary conditions are set to transmissive. |
| 116 | ii. Call CFLcond -- |
| 117 | |
| 118 | |