1 |
C $Header: $ |
2 |
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 |
PPdiffKr (i,j,k,bi,bj) = diffKrT |
44 |
end do |
45 |
end do |
46 |
end do |
47 |
end do |
48 |
end do |
49 |
|
50 |
#endif /* ALLOW_PP81 */ |
51 |
|
52 |
return |
53 |
end |