model/src/initialise_varia.F

C $Header: /u/gcmpack/development/heimbach/ecco_env/model/src/initialise_varia.F,v 1.2 2000/08/31 20:47:05 heimbach Exp $

#include "CPP_OPTIONS.h"

CStartOfInterface
      SUBROUTINE INITIALISE_VARIA(myThid)
C     /==========================================================\
C     | SUBROUTINE INITIALISE_VARIA                              |
C     | o Set the initial conditions for dynamics variables      |
C     |   and time dependent arrays                              |
C     |==========================================================|
C     | This routine reads/writes  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     \==========================================================/
      IMPLICIT NONE

C     == Global variables ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"

C     == Routine arguments ==
      INTEGER myThid
CEndOfInterface

C     == Local variables ==
#ifdef ALLOW_OBCS
      INTEGER bi,bj,K
#endif

C--   Initialise model fields.
C     Starting values of U, V, W, temp., salt. and tendency terms
C     are set here. Fields are either set to default or read from
C     stored files.
      CALL INI_FIELDS( myThid )
      _BARRIER

C--   Check dependances between packages
      CALL PACKAGES_INIT( myThid )

C--   Initialise time-stepping (predictor variables are carried in G)
      IF ( startTime .EQ. 0. ) THEN
        CALL INI_PREDICTOR( myThid )
        _BARRIER
      END IF

#ifdef ALLOW_OBCS
C--   Calculate values on open boundaries at iter=nIter0
      IF (openBoundaries) THEN
        DO bj = myByLo(myThid), myByHi(myThid)
         DO bi = myBxLo(myThid), myBxHi(myThid)
          DO K=1,Nr
           CALL SET_OBCS( K, bi, bj, startTime, myThid )
Caja       CALL CYCLE_OBCS( K, bi, bj, myThid )
          ENDDO
         ENDDO
        ENDDO
      ENDIF
#endif

C--   Finally summarise the model state
      CALL STATE_SUMMARY( myThid )

      END
1	C $Header: /u/gcmpack/development/heimbach/ecco_env/model/src/initialise_varia.F,v 1.2 2000/08/31 20:47:05 heimbach Exp $
2
3	#include "CPP_OPTIONS.h"
4
5	CStartOfInterface
6	SUBROUTINE INITIALISE_VARIA(myThid)
7	C /==========================================================\
8	C \| SUBROUTINE INITIALISE_VARIA \|
9	C \| o Set the initial conditions for dynamics variables \|
10	C \| and time dependent arrays \|
11	C \|==========================================================\|
12	C \| This routine reads/writes data from an input file and \|
13	C \| from various binary files. \|
14	C \| Each thread invokes an instance of this routine as does \|
15	C \| each process in a multi-process parallel environment like\|
16	C \| MPI. \|
17	C \==========================================================/
18	IMPLICIT NONE
19
20	C == Global variables ==
21	#include "SIZE.h"
22	#include "EEPARAMS.h"
23	#include "PARAMS.h"
24
25	C == Routine arguments ==
26	INTEGER myThid
27	CEndOfInterface
28
29	C == Local variables ==
30	#ifdef ALLOW_OBCS
31	INTEGER bi,bj,K
32	#endif
33
34	C-- Initialise model fields.
35	C Starting values of U, V, W, temp., salt. and tendency terms
36	C are set here. Fields are either set to default or read from
37	C stored files.
38	CALL INI_FIELDS( myThid )
39	_BARRIER
40
41	C-- Check dependances between packages
42	CALL PACKAGES_INIT( myThid )
43
44	C-- Initialise time-stepping (predictor variables are carried in G)
45	IF ( startTime .EQ. 0. ) THEN
46	CALL INI_PREDICTOR( myThid )
47	_BARRIER
48	END IF
49
50	#ifdef ALLOW_OBCS
51	C-- Calculate values on open boundaries at iter=nIter0
52	IF (openBoundaries) THEN
53	DO bj = myByLo(myThid), myByHi(myThid)
54	DO bi = myBxLo(myThid), myBxHi(myThid)
55	DO K=1,Nr
56	CALL SET_OBCS( K, bi, bj, startTime, myThid )
57	Caja CALL CYCLE_OBCS( K, bi, bj, myThid )
58	ENDDO
59	ENDDO
60	ENDDO
61	ENDIF
62	#endif
63
64	C-- Finally summarise the model state
65	CALL STATE_SUMMARY( myThid )
66
67	END