1 |
C $Header: /u/gcmpack/MITgcm/pkg/kl10/kl10_init_varia.F,v 1.1 2014/07/30 03:28:05 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) = diffKrNrS(k) |
44 |
ENDDO |
45 |
ENDDO |
46 |
ENDDO |
47 |
ENDDO |
48 |
ENDDO |
49 |
|
50 |
#endif /* ALLOW_KL10 */ |
51 |
|
52 |
RETURN |
53 |
END |