| 1 | jmc | 1.2 | C $Header: /u/gcmpack/MITgcm/pkg/pp81/pp81_init.F,v 1.1 2004/09/02 09:11:54 mlosch 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 |