model/src/ini_mixing.F

C $Header: /u/gcmpack/MITgcm/model/src/ini_mixing.F,v 1.3 2004/10/17 23:03:38 jmc Exp $
C $Name:  $

#include "CPP_OPTIONS.h"
# if (defined (ALLOW_GMREDI) && defined (ALLOW_KAPGM_CONTROL))
# include "GMREDI_OPTIONS.h"
#endif


CBOP
C     !ROUTINE: INI_MIXING
C     !INTERFACE:
      SUBROUTINE INI_MIXING( myThid )
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE INI_MIXING                                     
C     | o Initialise diapycnal diffusivity to default constant val.
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE
C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "DYNVARS.h"
# if (defined (ALLOW_GMREDI) && defined (ALLOW_KAPGM_CONTROL))
# include "GMREDI.h"
#endif

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine arguments ==
C     myThid -  Number of this instance
      INTEGER myThid

C     !LOCAL VARIABLES:
C     == Local variables ==
C     bi,bj  - Loop counters
C     I,J
      INTEGER bi, bj
      INTEGER i, j, k
CEOP

       DO bj = myByLo(myThid), myByHi(myThid)
        DO bi = myBxLo(myThid), myBxHi(myThid)
         DO k=1,Nr
          DO j=1,sNy
           DO i=1,sNx
#if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))
            diffKr(i,j,k,bi,bj) = diffKrNrS(k)
#endif
#if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_KAPGM_CONTROL))
            KapGM(i,j,k,bi,bj) = GM_background_K
#endif
           ENDDO
          ENDDO
         ENDDO
        ENDDO
       ENDDO

#if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))
      _EXCH_XYZ_R8(diffKr , myThid )
#endif
#if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_KAPGM_CONTROL))
      _EXCH_XYZ_R8(KapGM,   myThid )
#endif

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/model/src/ini_mixing.F,v 1.3 2004/10/17 23:03:38 jmc Exp $
2	C $Name: $
3
4	#include "CPP_OPTIONS.h"
5	# if (defined (ALLOW_GMREDI) && defined (ALLOW_KAPGM_CONTROL))
6	# include "GMREDI_OPTIONS.h"
7	#endif
8
9
10	CBOP
11	C !ROUTINE: INI_MIXING
12	C !INTERFACE:
13	SUBROUTINE INI_MIXING( myThid )
14	C !DESCRIPTION: \bv
15	C ==========================================================
16	C \| SUBROUTINE INI_MIXING
17	C \| o Initialise diapycnal diffusivity to default constant val.
18	C ==========================================================
19	C \ev
20
21	C !USES:
22	IMPLICIT NONE
23	C === Global variables ===
24	#include "SIZE.h"
25	#include "EEPARAMS.h"
26	#include "PARAMS.h"
27	#include "GRID.h"
28	#include "DYNVARS.h"
29	# if (defined (ALLOW_GMREDI) && defined (ALLOW_KAPGM_CONTROL))
30	# include "GMREDI.h"
31	#endif
32
33	C !INPUT/OUTPUT PARAMETERS:
34	C == Routine arguments ==
35	C myThid - Number of this instance
36	INTEGER myThid
37
38	C !LOCAL VARIABLES:
39	C == Local variables ==
40	C bi,bj - Loop counters
41	C I,J
42	INTEGER bi, bj
43	INTEGER i, j, k
44	CEOP
45
46	DO bj = myByLo(myThid), myByHi(myThid)
47	DO bi = myBxLo(myThid), myBxHi(myThid)
48	DO k=1,Nr
49	DO j=1,sNy
50	DO i=1,sNx
51	#if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))
52	diffKr(i,j,k,bi,bj) = diffKrNrS(k)
53	#endif
54	#if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_KAPGM_CONTROL))
55	KapGM(i,j,k,bi,bj) = GM_background_K
56	#endif
57	ENDDO
58	ENDDO
59	ENDDO
60	ENDDO
61	ENDDO
62
63	#if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))
64	_EXCH_XYZ_R8(diffKr , myThid )
65	#endif
66	#if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_KAPGM_CONTROL))
67	_EXCH_XYZ_R8(KapGM, myThid )
68	#endif
69
70	RETURN
71	END