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