| | 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 | |
