/[MITgcm]/MITgcm/model/src/ini_mixing.F
ViewVC logotype

Annotation of /MITgcm/model/src/ini_mixing.F

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


Revision 1.1 - (hide annotations) (download)
Mon Aug 13 18:07:35 2001 UTC (22 years, 9 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint40pre7, checkpoint40pre6, checkpoint40pre9, checkpoint40pre8, checkpoint40
Added diapycnal (diffkr) and isopycanl (kapgm) diffusion
to set of control variables.

1 heimbach 1.1
2     #include "CPP_OPTIONS.h"
3    
4     CStartOfInterface
5     SUBROUTINE INI_MIXING( myThid )
6     C /==========================================================\
7     C | SUBROUTINE INI_MIXING |
8     C | o Initialise diapycnal diffusivity |
9     C \==========================================================/
10     IMPLICIT NONE
11    
12     C === Global variables ===
13     #include "SIZE.h"
14     #include "EEPARAMS.h"
15     #include "PARAMS.h"
16     #include "GRID.h"
17     #include "DYNVARS.h"
18    
19     C == Routine arguments ==
20     C myThid - Number of this instance
21     INTEGER myThid
22     CEndOfInterface
23    
24     C == Local variables ==
25     C bi,bj - Loop counters
26     C I,J
27    
28     INTEGER bi, bj
29     INTEGER i, j, k
30    
31     DO bj = myByLo(myThid), myByHi(myThid)
32     DO bi = myBxLo(myThid), myBxHi(myThid)
33     DO j=1,sNy
34     DO i=1,sNx
35     DO k=1,Nr
36     #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))
37     diffKr(i,j,k,bi,bj) = diffKrS
38     #endif
39     #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_KAPGM_CONTROL))
40     KapGM(i,j,k,bi,bj) = GMkbackground
41     #endif
42     ENDDO
43     ENDDO
44     ENDDO
45     ENDDO
46     ENDDO
47    
48     #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))
49     _EXCH_XYZ_R8(diffKr , myThid )
50     #endif
51     #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_KAPGM_CONTROL))
52     _EXCH_XYZ_R8(KapGM, myThid )
53     #endif
54    
55     RETURN
56     END

  ViewVC Help
Powered by ViewVC 1.1.22