model/src/ini_pnh.F

C $Header: $
C $Name: $

#include "CPP_OPTIONS.h"

      SUBROUTINE INI_PNH( myThid )
C     /==========================================================\
C     | SUBROUTINE INI_PNH                                       |
C     | o Set model initial non-hydrostatic pressure.            |
C     |==========================================================|
C     | There are several options for setting the initial        |
C     | surface displacement (r unit) field.                     |
C     |  1. Inline code                                          |
C     |  2. Two-dimensional data from a file.                    |
C     \==========================================================/
      IMPLICIT NONE

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

C     == Routine arguments ==
C     myThid -  Number of this instance of INI_PSURF 
      INTEGER myThid

C     == Local variables ==
      INTEGER bi,bj,I,J,K

#ifdef ALLOW_NONHYDROSTATIC

C--   Initialise surface position anomaly to zero
      DO bj = myByLo(myThid), myByHi(myThid)
       DO bi = myBxLo(myThid), myBxHi(myThid)
        DO K = 1, NR
         DO J=1-Oly,sNy+Oly
          DO I=1-Olx,sNx+Olx
           phi_nh(I,J,K,bi,bj) = 0. _d 0
          ENDDO
         ENDDO
        ENDDO
       ENDDO
      ENDDO
C     Read an initial state
c     IF (pnhInitFile .NE. ' ') THEN
c      _BEGIN_MASTER( myThid )
c      CALL READ_FLD_XY_RL( pnhInitFile, ' ', phi_nh, 0, myThid )
c      _END_MASTER(myThid)
c     ENDIF
C
      _EXCH_XYZ_R8(phi_nh, myThid)

#endif /* ALLOW_NONHYDROSTATIC */
      RETURN
      END
1	adcroft	1.1	C $Header: $
2			C $Name: $
3
4			#include "CPP_OPTIONS.h"
5
6			SUBROUTINE INI_PNH( myThid )
7			C /==========================================================\
8			C \| SUBROUTINE INI_PNH \|
9			C \| o Set model initial non-hydrostatic pressure. \|
10			C \|==========================================================\|
11			C \| There are several options for setting the initial \|
12			C \| surface displacement (r unit) field. \|
13			C \| 1. Inline code \|
14			C \| 2. Two-dimensional data from a file. \|
15			C \==========================================================/
16			IMPLICIT NONE
17
18			C === Global variables ===
19			#include "SIZE.h"
20			#include "EEPARAMS.h"
21			#include "PARAMS.h"
22			#include "GRID.h"
23			#include "DYNVARS.h"
24
25			C == Routine arguments ==
26			C myThid - Number of this instance of INI_PSURF
27			INTEGER myThid
28
29			C == Local variables ==
30			INTEGER bi,bj,I,J,K
31
32			#ifdef ALLOW_NONHYDROSTATIC
33
34			C-- Initialise surface position anomaly to zero
35			DO bj = myByLo(myThid), myByHi(myThid)
36			DO bi = myBxLo(myThid), myBxHi(myThid)
37			DO K = 1, NR
38			DO J=1-Oly,sNy+Oly
39			DO I=1-Olx,sNx+Olx
40			phi_nh(I,J,K,bi,bj) = 0. _d 0
41			ENDDO
42			ENDDO
43			ENDDO
44			ENDDO
45			ENDDO
46			C Read an initial state
47			c IF (pnhInitFile .NE. ' ') THEN
48			c _BEGIN_MASTER( myThid )
49			c CALL READ_FLD_XY_RL( pnhInitFile, ' ', phi_nh, 0, myThid )
50			c _END_MASTER(myThid)
51			c ENDIF
52			C
53			_EXCH_XYZ_R8(phi_nh, myThid)
54
55			#endif /* ALLOW_NONHYDROSTATIC */
56			RETURN
57			END