1 |
jmc |
1.1 |
C $Header: /u/gcmpack/MITgcm/pkg/pp81/pp81_init_varia.F,v 1.2 2009/10/08 20:08:20 jmc Exp $ |
2 |
|
|
C $Name: $ |
3 |
|
|
|
4 |
|
|
#include "KL10_OPTIONS.h" |
5 |
|
|
|
6 |
|
|
SUBROUTINE KL10_INIT_VARIA( myThid ) |
7 |
|
|
C *==========================================================* |
8 |
|
|
C | SUBROUTINE KL10_INIT_VARIA |
9 |
|
|
C | o Routine to initialize KL10 parameters and variables. |
10 |
|
|
C *==========================================================* |
11 |
|
|
C | Initialize KL10 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 |
|
|
c#include "GRID.h" |
20 |
|
|
#include "KL10.h" |
21 |
|
|
|
22 |
|
|
C === Routine arguments === |
23 |
|
|
C myThid - Number of this instance of KL10_INIT |
24 |
|
|
INTEGER myThid |
25 |
|
|
|
26 |
|
|
#ifdef ALLOW_KL10 |
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 KL variables KLviscAr, KLdiffKrT |
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 |
|
|
KLviscAr (i,j,k,bi,bj) = viscArNr(k) |
43 |
|
|
KLdiffKr (i,j,k,bi,bj) = diffKrNrT(k) |
44 |
|
|
ENDDO |
45 |
|
|
ENDDO |
46 |
|
|
ENDDO |
47 |
|
|
ENDDO |
48 |
|
|
ENDDO |
49 |
|
|
|
50 |
|
|
#endif /* ALLOW_KL10 */ |
51 |
|
|
|
52 |
|
|
RETURN |
53 |
|
|
END |