/[MITgcm]/MITgcm/pkg/my82/my82_init_varia.F
ViewVC logotype

Annotation of /MITgcm/pkg/my82/my82_init_varia.F

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

Revision 1.3 - (hide annotations) (download)
Thu Oct 8 20:07:53 2009 UTC (14 years, 8 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint64, checkpoint65, checkpoint62, checkpoint63, checkpoint65h, checkpoint65i, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint61w, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.2: +17 -17 lines 
modif for vertical profile of background viscosity
1 jmc 1.3 C $Header: /u/gcmpack/MITgcm/pkg/my82/my82_init_varia.F,v 1.2 2009/01/20 00:26:04 jmc Exp $
2 mlosch 1.1 C $Name: $
3    
4     #include "MY82_OPTIONS.h"
5    
6     SUBROUTINE MY82_INIT_VARIA( myThid )
7 jmc 1.3 C *==========================================================*
8 mlosch 1.1 C | SUBROUTINE MY82_INIT_VARIA |
9     C | o Routine to initialize MY82 parameters and variables. |
10 jmc 1.3 C *==========================================================*
11 mlosch 1.1 C | Initialize MY92 parameters and variables. |
12 jmc 1.3 C *==========================================================*
13 mlosch 1.1 IMPLICIT NONE
14    
15     C === Global variables ===
16     #include "SIZE.h"
17     #include "EEPARAMS.h"
18     #include "PARAMS.h"
19     #include "GRID.h"
20     #include "MY82.h"
21    
22     C === Routine arguments ===
23     C myThid - Number of this instance of MY_INIT
24     INTEGER myThid
25    
26     #ifdef ALLOW_MY82
27    
28     C === Local variables ===
29     C i,j,k,bi,bj - Loop counters
30    
31     INTEGER i, j, k, bi, bj
32     _RL gam1, gam2
33    
34     C magic parameters( M. Satoh, p. 314)
35 jmc 1.2 gam1 = 1. _d 0/3. _d 0 - 2. _d 0 *A1/B1
36     gam2 = (B2+6. _d 0*A1)/B1
37     alpha1 = 3. _d 0 *A2*gam1
38     alpha2 = 3. _d 0 *A2*(gam1+gam2)
39 mlosch 1.1 beta1 = A1*B1*(gam1-C1)
40 jmc 1.2 beta2 = A1*(B1*(gam1-C1) + 6. _d 0 *A1 + 3. _d 0 *A2)
41 mlosch 1.1 beta3 = A2*B1*gam1
42 jmc 1.2 beta4 = A2*(B1*(gam1+gam2) - 3. _d 0 *A1)
43 mlosch 1.1
44     C-----------------------------------------------------------------------
45     C Initialize MY82 variables MYviscAr, MYdiffKrT
46     C-----------------------------------------------------------------------
47    
48 jmc 1.3 DO bj = myByLo(myThid), myByHi(myThid)
49     DO bi = myBxLo(myThid), myBxHi(myThid)
50     DO k = 1, Nr
51     DO j = 1-OLy, sNy+OLy
52     DO i = 1-OLx, sNx+OLx
53 mlosch 1.1 MYhbl (i,j ,bi,bj) = 0. _d 0
54 jmc 1.3 MYviscAr (i,j,k,bi,bj) = viscArNr(k)
55 mlosch 1.1 MYdiffKr (i,j,k,bi,bj) = diffKrNrT(k)
56 jmc 1.3 ENDDO
57     ENDDO
58     ENDDO
59     ENDDO
60     ENDDO
61 mlosch 1.1
62     #endif /* ALLOW_MY82 */
63    
64 jmc 1.3 RETURN
65     END
  ViewVC Help
Powered by ViewVC 1.1.22