| 1 | ce107 | 1.2 | C $Id: initialise.F,v 1.1 2006/05/12 21:58:06 ce107 Exp $ | 
| 2 | ce107 | 1.1 | CStartOfInterface | 
| 3 |  |  | SUBROUTINE INITIALISE | 
| 4 |  |  | C     /==========================================================\ | 
| 5 |  |  | C     | SUBROUTINE INITIALISE                                    | | 
| 6 |  |  | C     | o Master controlling routine for setting model initial   | | 
| 7 |  |  | C     |   conditions.                                            | | 
| 8 |  |  | C     |==========================================================| | 
| 9 |  |  | C     | INITIALISE is invoked at the start of the model to set   | | 
| 10 |  |  | C     | the model arrays. It reads data from an input file and   | | 
| 11 |  |  | C     | from various binary files.                               | | 
| 12 |  |  | C     | Each thread invokes an instance of this routine as does  | | 
| 13 |  |  | C     | each process in a multi-process parallel environment like| | 
| 14 |  |  | C     | MPI.                                                     | | 
| 15 |  |  | C     | Notes                                                    | | 
| 16 |  |  | C     | =====                                                    | | 
| 17 |  |  | C     | C*P* comments indicating place holders for which code is | | 
| 18 |  |  | C     |      presently being developed.                          | | 
| 19 |  |  | C     \==========================================================/ | 
| 20 |  |  |  | 
| 21 |  |  | C     == Routine arguments == | 
| 22 |  |  | CEndOfInterface | 
| 23 |  |  |  | 
| 24 |  |  |  | 
| 25 |  |  | C--   Set model parameters. | 
| 26 |  |  | C     Parameters are set to defaults and then updates are read from | 
| 27 |  |  | C     an input file called data. | 
| 28 |  |  | CALL INI_PARMS | 
| 29 |  |  |  | 
| 30 |  |  | C--   Initialise the con. grad solver operators. | 
| 31 |  |  | C     Laplace operator diagonals are set here. | 
| 32 |  |  | CALL INI_CG2D | 
| 33 |  |  |  | 
| 34 |  |  | C--   Initialise the right-hand side. | 
| 35 |  |  | C     Only used during testing to set right-hand side for elliptic equation | 
| 36 |  |  | C     so that  stand alone inversion can be analysed. | 
| 37 |  |  | CALL INI_RHS | 
| 38 |  |  |  | 
| 39 |  |  | END | 
| 40 |  |  |  |