| 1 |
C $Id: initialise.F,v 1.1 2006/05/12 21:58:06 ce107 Exp $ |
| 2 |
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 |
|