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

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

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

revision 1.6 by jmc, Sun Dec 16 18:54:49 2001 UTC revision 1.7 by heimbach, Sun Mar 24 02:33:16 2002 UTC
# Line 40  C     Initialize arrays in common blocks Line 40  C     Initialize arrays in common blocks
40          DO k=1,Nr          DO k=1,Nr
41           DO j=1-Oly,sNy+OLy           DO j=1-Oly,sNy+OLy
42            DO i=1-Olx,sNx+Olx            DO i=1-Olx,sNx+Olx
43             Kwx(i,j,k,bi,bj) = 0.             Kwx(i,j,k,bi,bj) = 0. _d 0
44             Kwy(i,j,k,bi,bj) = 0.             Kwy(i,j,k,bi,bj) = 0. _d 0
45             Kwz(i,j,k,bi,bj) = 0.             Kwz(i,j,k,bi,bj) = 0. _d 0
46  #ifdef GM_EXTRA_DIAGONAL  #ifdef GM_EXTRA_DIAGONAL
47             Kuz(i,j,k,bi,bj) = 0.             Kuz(i,j,k,bi,bj) = 0. _d 0
48             Kvz(i,j,k,bi,bj) = 0.             Kvz(i,j,k,bi,bj) = 0. _d 0
49  #endif  #endif
50  #ifdef GM_NON_UNITY_DIAGONAL  #ifdef GM_NON_UNITY_DIAGONAL
51             Kux(i,j,k,bi,bj) = 0.             Kux(i,j,k,bi,bj) = 0. _d 0
52             Kvy(i,j,k,bi,bj) = 0.             Kvy(i,j,k,bi,bj) = 0. _d 0
53  #endif  #endif
54  #ifdef GM_BOLUS_ADVEC  #ifdef GM_BOLUS_ADVEC
55             GM_PsiX(i,j,k,bi,bj) = 0.             GM_PsiX(i,j,k,bi,bj) = 0. _d 0
56             GM_PsiY(i,j,k,bi,bj) = 0.             GM_PsiY(i,j,k,bi,bj) = 0. _d 0
57  #endif  #endif
58  #ifdef GM_VISBECK_VARIABLE_K  #ifdef GM_VISBECK_VARIABLE_K
59            VisbeckK(i,j,bi,bj) = 0.            VisbeckK(i,j,bi,bj) = 0. _d 0
60  #endif  #endif
61            ENDDO            ENDDO
62           ENDDO           ENDDO
# Line 68  C     Initialize averages to zero Line 68  C     Initialize averages to zero
68          CALL TIMEAVE_RESET(GM_Kwy_T,Nr, bi,bj,myThid)          CALL TIMEAVE_RESET(GM_Kwy_T,Nr, bi,bj,myThid)
69          CALL TIMEAVE_RESET(GM_Kwz_T,Nr, bi,bj,myThid)          CALL TIMEAVE_RESET(GM_Kwz_T,Nr, bi,bj,myThid)
70          DO k=1,Nr          DO k=1,Nr
71           GM_TimeAve(k,bi,bj)=0.           GM_TimeAve(k,bi,bj)=0. _d 0
72          ENDDO          ENDDO
73  #ifdef GM_VISBECK_VARIABLE_K  #ifdef GM_VISBECK_VARIABLE_K
74          CALL TIMEAVE_RESET(Visbeck_K_T, 1, bi,bj,myThid)          CALL TIMEAVE_RESET(Visbeck_K_T, 1, bi,bj,myThid)
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22