pkg/pp81/pp81_init.F

C $Header: /u/gcmpack/MITgcm/pkg/pp81/pp81_init.F,v 1.2 2004/10/17 23:12:59 jmc Exp $
C $Name:  $

#include "PP81_OPTIONS.h"

      SUBROUTINE PP81_INIT( myThid )
C     /==========================================================\
C     | SUBROUTINE PP81_INIT                                     |
C     | o Routine to initialize PP81 parameters and variables.   |
C     |==========================================================|
C     | Initialize PP81 parameters and variables.                |
C     \==========================================================/
      IMPLICIT NONE

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

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

#ifdef ALLOW_PP81

C     === Local variables ===
C     i,j,k,bi,bj - Loop counters

      INTEGER i, j, k, bi, bj

C-----------------------------------------------------------------------
C Initialize PP variables PPviscAr, PPdiffKrT
C-----------------------------------------------------------------------

      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
           PPviscAr (i,j,k,bi,bj) = viscAr
           PPdiffKr (i,j,k,bi,bj) = diffKrNrT(k)
          end do
         end do
        end do
       end do
      end do

#endif /* ALLOW_PP81 */

      return
      end
1	jmc	1.3	C $Header: /u/gcmpack/MITgcm/pkg/pp81/pp81_init.F,v 1.2 2004/10/17 23:12:59 jmc Exp $
2	mlosch	1.1	C $Name: $
3
4			#include "PP81_OPTIONS.h"
5
6			SUBROUTINE PP81_INIT( myThid )
7			C /==========================================================\
8			C \| SUBROUTINE PP81_INIT \|
9			C \| o Routine to initialize PP81 parameters and variables. \|
10			C \|==========================================================\|
11			C \| Initialize PP81 parameters and variables. \|
12			C \==========================================================/
13			IMPLICIT NONE
14
15			C === Global variables ===
16			#include "SIZE.h"
17			#include "EEPARAMS.h"
18			#include "PARAMS.h"
19			#include "GRID.h"
20			#include "PP81.h"
21
22			C === Routine arguments ===
23			C myThid - Number of this instance of PP81_INIT
24			INTEGER myThid
25
26			#ifdef ALLOW_PP81
27
28			C === Local variables ===
29			C i,j,k,bi,bj - Loop counters
30
31			INTEGER i, j, k, bi, bj
32
33			C-----------------------------------------------------------------------
34			C Initialize PP variables PPviscAr, PPdiffKrT
35			C-----------------------------------------------------------------------
36
37			do bj = myByLo(myThid), myByHi(myThid)
38			do bi = myBxLo(myThid), myBxHi(myThid)
39			do k = 1, Nr
40			do j = 1-OLy, sNy+OLy
41			do i = 1-OLx, sNx+OLx
42			PPviscAr (i,j,k,bi,bj) = viscAr
43	jmc	1.2	PPdiffKr (i,j,k,bi,bj) = diffKrNrT(k)
44	mlosch	1.1	end do
45			end do
46			end do
47			end do
48			end do
49
50			#endif /* ALLOW_PP81 */
51
52			return
53			end