16 | | 3. |
| 16 | 3. tol = |
| 17 | |
| 18 | II. Main Program Body |
| 19 | 1. Call P* routine - find and record p* |
| 20 | 2. Call Identification routine - find types of waves present, record |
| 21 | 3. Call U*() - reads in vars and returns U* |
| 22 | 4. Call rho* - reads in vars and returns rho* |
| 23 | 2. Print to File - u*, rho*, p* |
| 24 | |
| 25 | III. Functions |
| 26 | 1. FL = () |
| 27 | 2. FR = () |
| 28 | 3. U* = () |
| 29 | 4. rho*L = {} |
| 30 | 5. rho*R = {} |
| 31 | |
| 32 | IV. Subroutines |
| 33 | 1. Find p* using initial estimate |
| 34 | a. type of approximation (1)-(5) for different options (on screen options) |
| 35 | b. let p0=(approximation) |
| 36 | c. Begin iterative process, that ends once | | < tol |
| 37 | |
| 38 | 2. Determine types of waves present |
| 39 | a. comparing p* to pL and pR |
| 40 | b. series of if statements corresponding to flow chart |
| 41 | c. left_wave = 1 for shock, 2 for rarefaction |