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

Contents of /MITgcm/pkg/my82/my82_init.F

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

Revision 1.3 - (show annotations) (download)
Mon May 30 07:41:45 2005 UTC (19 years, 7 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint57t_post, checkpoint58l_post, checkpoint57m_post, checkpoint57i_post, checkpoint58e_post, checkpoint57v_post, checkpoint57s_post, checkpoint57j_post, checkpoint58b_post, checkpoint58m_post, checkpoint57k_post, checkpoint58r_post, checkpoint57y_post, checkpoint58g_post, checkpoint57x_post, checkpoint58n_post, checkpoint58x_post, checkpoint58h_post, checkpoint58w_post, checkpoint58j_post, checkpoint57y_pre, checkpoint57o_post, checkpoint59p, checkpoint57r_post, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint59j, checkpoint59, checkpoint58, checkpoint58f_post, checkpoint57n_post, checkpoint58d_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint58a_post, checkpoint58i_post, checkpoint57q_post, checkpoint58o_post, checkpoint57z_post, checkpoint58c_post, checkpoint58k_post, checkpoint58u_post, checkpoint58y_post, checkpoint58v_post, checkpoint58s_post, checkpoint58p_post, checkpoint58t_post, checkpoint58q_post, checkpoint57l_post
Changes since 1.2: +2 -2 lines 
added a few "_d 0" and D0 (in PARAMETER statements) in a desperate
effort to make vermix.my82 pass on a Sun, unfortunately no success;
but for cleaner code I check it in anyway.
1 C $Header: /u/gcmpack/MITgcm/pkg/my82/my82_init.F,v 1.2 2004/10/17 23:12:31 jmc Exp $
2 C $Name: $
3
4 #include "MY82_OPTIONS.h"
5
6 SUBROUTINE MY82_INIT( myThid )
7 C /==========================================================\
8 C | SUBROUTINE MY82_INIT |
9 C | o Routine to initialize MY82 parameters and variables. |
10 C |==========================================================|
11 C | Initialize MY92 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 #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 gam1 = 1./3.-2.*A1/B1
36 gam2 = (B2+6.*A1)/B1
37 alpha1 = 3.*A2*gam1
38 alpha2 = 3.*A2*(gam1+gam2)
39 beta1 = A1*B1*(gam1-C1)
40 beta2 = A1*(B1*(gam1-C1) + 6.*A1 + 3.*A2)
41 beta3 = A2*B1*gam1
42 beta4 = A2*(B1*(gam1+gam2) - 3.*A1)
43
44 C-----------------------------------------------------------------------
45 C Initialize MY82 variables MYviscAr, MYdiffKrT
46 C-----------------------------------------------------------------------
47
48 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 MYhbl (i,j ,bi,bj) = 0. _d 0
54 MYviscAr (i,j,k,bi,bj) = viscAr
55 MYdiffKr (i,j,k,bi,bj) = diffKrNrT(k)
56 end do
57 end do
58 end do
59 end do
60 end do
61
62 #endif /* ALLOW_MY82 */
63
64 return
65 end
  ViewVC Help
Powered by ViewVC 1.1.22