/[MITgcm]/MITgcm/pkg/gmredi/gmredi_init.F
ViewVC logotype

Contents of /MITgcm/pkg/gmredi/gmredi_init.F

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

Revision 1.10 - (show annotations) (download)
Tue Jan 4 00:20:35 2005 UTC (19 years, 5 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint57m_post, checkpoint57g_pre, checkpoint57s_post, checkpoint58b_post, checkpoint57g_post, checkpoint57y_post, checkpoint57r_post, checkpoint57d_post, checkpoint57i_post, checkpoint59, checkpoint58, checkpoint58f_post, checkpoint57n_post, checkpoint58d_post, checkpoint58a_post, checkpoint57z_post, checkpoint58y_post, checkpoint58t_post, checkpoint57l_post, checkpoint57h_post, checkpoint57t_post, checkpoint57v_post, checkpoint57f_post, checkpoint60, checkpoint61, checkpoint62, checkpoint57h_pre, checkpoint57x_post, checkpoint58w_post, checkpoint57y_pre, checkpoint58o_post, checkpoint57c_post, checkpoint58p_post, checkpoint58q_post, checkpoint58e_post, mitgcm_mapl_00, checkpoint58m_post, checkpoint57c_pre, checkpoint58r_post, checkpoint58n_post, checkpoint57e_post, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint57p_post, checkpint57u_post, checkpoint57q_post, eckpoint57e_pre, checkpoint58k_post, checkpoint58v_post, checkpoint58l_post, checkpoint57h_done, checkpoint57j_post, checkpoint57f_pre, checkpoint61f, checkpoint58g_post, checkpoint58x_post, checkpoint61n, checkpoint59j, checkpoint58h_post, checkpoint58j_post, checkpoint57o_post, checkpoint61q, checkpoint57k_post, checkpoint61z, checkpoint57w_post, checkpoint61e, checkpoint58i_post, checkpoint58c_post, checkpoint58u_post, checkpoint58s_post, checkpoint61g, checkpoint61d, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61x, checkpoint61y
Changes since 1.9: +19 -11 lines 
rename GMREDI_DIAGS.h to GMREDI_TAVE.h (since it only contains TimeAve vars)
1 C $Header: /u/gcmpack/MITgcm/pkg/gmredi/gmredi_init.F,v 1.9 2004/03/24 15:21:15 adcroft Exp $
2 C $Name: $
3
4 #include "GMREDI_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: GMREDI_INIT
8 C !INTERFACE:
9 SUBROUTINE GMREDI_INIT( myThid )
10
11 C !DESCRIPTION: \bv
12 C *==========================================================*
13 C | SUBROUTINE GMREDI_INIT
14 C | o Routine to initialize GM/Redi variables
15 C *==========================================================*
16 C \ev
17 C !USES:
18 IMPLICIT NONE
19
20 C === Global variables ===
21 #include "SIZE.h"
22 #include "EEPARAMS.h"
23 #include "PARAMS.h"
24 #include "GRID.h"
25 #include "GMREDI.h"
26 #include "GMREDI_TAVE.h"
27
28 C !INPUT/OUTPUT PARAMETERS:
29 C === Routine arguments ===
30 C myThid :: my Thread Id number
31 INTEGER myThid
32 CEOP
33
34 #ifdef ALLOW_GMREDI
35
36 C !LOCAL VARIABLES:
37 C === Local variables ===
38 INTEGER i,j,k,bi,bj
39
40 DO bj = myByLo(myThid), myByHi(myThid)
41 DO bi = myBxLo(myThid), myBxHi(myThid)
42
43 C Initialize arrays in common blocks :
44 DO k=1,Nr
45 DO j=1-Oly,sNy+OLy
46 DO i=1-Olx,sNx+Olx
47 Kwx(i,j,k,bi,bj) = 0. _d 0
48 Kwy(i,j,k,bi,bj) = 0. _d 0
49 Kwz(i,j,k,bi,bj) = 0. _d 0
50 #ifdef GM_EXTRA_DIAGONAL
51 Kuz(i,j,k,bi,bj) = 0. _d 0
52 Kvz(i,j,k,bi,bj) = 0. _d 0
53 #endif
54 #ifdef GM_NON_UNITY_DIAGONAL
55 Kux(i,j,k,bi,bj) = 0. _d 0
56 Kvy(i,j,k,bi,bj) = 0. _d 0
57 #endif
58 #ifdef GM_BOLUS_ADVEC
59 GM_PsiX(i,j,k,bi,bj) = 0. _d 0
60 GM_PsiY(i,j,k,bi,bj) = 0. _d 0
61 #endif
62 #ifdef GM_VISBECK_VARIABLE_K
63 VisbeckK(i,j,bi,bj) = 0. _d 0
64 #endif
65 ENDDO
66 ENDDO
67 ENDDO
68
69 #ifdef ALLOW_TIMEAVE
70 C Initialize averages to zero
71 CALL TIMEAVE_RESET(GM_Kwx_T,Nr, bi,bj,myThid)
72 CALL TIMEAVE_RESET(GM_Kwy_T,Nr, bi,bj,myThid)
73 CALL TIMEAVE_RESET(GM_Kwz_T,Nr, bi,bj,myThid)
74 DO k=1,Nr
75 GM_TimeAve(k,bi,bj)=0. _d 0
76 ENDDO
77 #ifdef GM_VISBECK_VARIABLE_K
78 CALL TIMEAVE_RESET(Visbeck_K_T, 1, bi,bj,myThid)
79 #endif
80 #ifdef GM_BOLUS_ADVEC
81 CALL TIMEAVE_RESET(GM_PsiXtave,Nr, bi,bj,myThid)
82 CALL TIMEAVE_RESET(GM_PsiYtave,Nr, bi,bj,myThid)
83 #endif
84 #endif /* ALLOW_TIMEAVE */
85
86 C- end bi,bj loops
87 ENDDO
88 ENDDO
89 #endif /* ALLOW_GMREDI */
90
91 RETURN
92 END
  ViewVC Help
Powered by ViewVC 1.1.22