/[MITgcm]/MITgcm/pkg/kl10/kl10_init_varia.F
ViewVC logotype

Annotation of /MITgcm/pkg/kl10/kl10_init_varia.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.1 - (hide annotations) (download)
Wed Jul 30 03:28:05 2014 UTC (9 years, 9 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint65h, checkpoint65i, checkpoint65b, checkpoint65c, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e
-  add new pkg "kl10" for mixing due to internal wave breaking
  ( http://www.sciencedirect.com/science/article/pii/S1463500310000144 )

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

  ViewVC Help
Powered by ViewVC 1.1.22