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