model/src/ini_pnh.F

C $Header: /u/gcmpack/MITgcm/model/src/ini_pnh.F,v 1.6 2006/08/04 00:40:43 jmc Exp $
C $Name:  $

#include "CPP_OPTIONS.h"

CBOP
C     !ROUTINE: INI_PNH
C     !INTERFACE:
      SUBROUTINE INI_PNH( myThid )

C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE INI_PNH
C     | o Set model initial non-hydrostatic potential.
C     *==========================================================*
C     | There are several options for setting the initial
C     | NH-potential (=Pressure/rhoConst) anomaly:
C     |  1. Inline code
C     |  2. Three-dimensional data from a file.
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE
C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "NH_VARS.h"

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine arguments ==
C     myThid  :: My Thread Id number
      INTEGER myThid

#ifdef ALLOW_NONHYDROSTATIC
C     !LOCAL VARIABLES:
C     == Local variables ==
      INTEGER bi,bj,I,J,K
CEOP


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 (phiNHinitFile .NE. ' ') THEN
c      _BARRIER
c      CALL READ_FLD_XY_RL( phiNHinitFile, ' ', phi_nh, 0, myThid )
c      _EXCH_XYZ_RL(phi_nh, myThid)
c     ENDIF

#endif /* ALLOW_NONHYDROSTATIC */
      RETURN
      END
1	jmc	1.7	C $Header: /u/gcmpack/MITgcm/model/src/ini_pnh.F,v 1.6 2006/08/04 00:40:43 jmc Exp $
2	cnh	1.2	C $Name: $
3	adcroft	1.1
4			#include "CPP_OPTIONS.h"
5
6	cnh	1.2	CBOP
7			C !ROUTINE: INI_PNH
8			C !INTERFACE:
9	adcroft	1.1	SUBROUTINE INI_PNH( myThid )
10	cnh	1.2
11			C !DESCRIPTION: \bv
12			C ==========================================================
13	jmc	1.6	C \| SUBROUTINE INI_PNH
14			C \| o Set model initial non-hydrostatic potential.
15	cnh	1.2	C ==========================================================
16	jmc	1.6	C \| There are several options for setting the initial
17	jmc	1.5	C \| NH-potential (=Pressure/rhoConst) anomaly:
18	jmc	1.6	C \| 1. Inline code
19			C \| 2. Three-dimensional data from a file.
20	cnh	1.2	C ==========================================================
21			C \ev
22
23			C !USES:
24	adcroft	1.1	IMPLICIT NONE
25			C === Global variables ===
26			#include "SIZE.h"
27			#include "EEPARAMS.h"
28			#include "PARAMS.h"
29			#include "GRID.h"
30	jmc	1.5	#include "NH_VARS.h"
31	adcroft	1.1
32	cnh	1.2	C !INPUT/OUTPUT PARAMETERS:
33	adcroft	1.1	C == Routine arguments ==
34	jmc	1.5	C myThid :: My Thread Id number
35	adcroft	1.1	INTEGER myThid
36
37	jmc	1.4	#ifdef ALLOW_NONHYDROSTATIC
38	cnh	1.2	C !LOCAL VARIABLES:
39	adcroft	1.1	C == Local variables ==
40			INTEGER bi,bj,I,J,K
41	cnh	1.2	CEOP
42	adcroft	1.1
43
44			C-- Initialise surface position anomaly to zero
45			DO bj = myByLo(myThid), myByHi(myThid)
46			DO bi = myBxLo(myThid), myBxHi(myThid)
47	jmc	1.5	DO K = 1, Nr
48	adcroft	1.1	DO J=1-Oly,sNy+Oly
49			DO I=1-Olx,sNx+Olx
50			phi_nh(I,J,K,bi,bj) = 0. _d 0
51			ENDDO
52			ENDDO
53			ENDDO
54			ENDDO
55			ENDDO
56	jmc	1.6
57			C-- Read an initial state
58	jmc	1.5	c IF (phiNHinitFile .NE. ' ') THEN
59	jmc	1.6	c _BARRIER
60	jmc	1.5	c CALL READ_FLD_XY_RL( phiNHinitFile, ' ', phi_nh, 0, myThid )
61	jmc	1.6	c _EXCH_XYZ_RL(phi_nh, myThid)
62	adcroft	1.1	c ENDIF
63	jmc	1.5
64	adcroft	1.1	#endif /* ALLOW_NONHYDROSTATIC */
65			RETURN
66			END