C $Id: initialise.F,v 1.2 2006/05/12 22:34:55 ce107 Exp $
CStartOfInterface
      SUBROUTINE INITIALISE
C     /==========================================================\
C     | SUBROUTINE INITIALISE                                    |
C     | o Master controlling routine for setting model initial   |
C     |   conditions.                                            |
C     |==========================================================|
C     | INITIALISE is invoked at the start of the model to set   |
C     | the model arrays. It reads data from an input file and   |
C     | from various binary files.                               |
C     | Each thread invokes an instance of this routine as does  |
C     | each process in a multi-process parallel environment like|
C     | MPI.                                                     |
C     | Notes                                                    |
C     | =====                                                    |
C     | C*P* comments indicating place holders for which code is |
C     |      presently being developed.                          |
C     \==========================================================/

C     == Routine arguments ==
CEndOfInterface


C--   Set model parameters.
C     Parameters are set to defaults and then updates are read from
C     an input file called data.
      CALL INI_PARMS

C--   Initialise the con. grad solver operators.
C     Laplace operator diagonals are set here.
      CALL INI_CG2D

C--   Initialise the right-hand side.
C     Only used during testing to set right-hand side for elliptic equation
C     so that  stand alone inversion can be analysed.
      CALL INI_RHS

      END