/[MITgcm]/MITgcm/pkg/mom_vecinv/mom_vecinv.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_vecinv/mom_vecinv.F

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

revision 1.31 by jmc, Wed Nov 10 03:05:04 2004 UTC revision 1.36 by jmc, Wed Apr 6 18:42:12 2005 UTC
# Line 66  C     myThid - Instance number for this Line 66  C     myThid - Instance number for this
66  #ifdef ALLOW_MOM_VECINV  #ifdef ALLOW_MOM_VECINV
67    
68  C     == Functions ==  C     == Functions ==
69        LOGICAL  DIFFERENT_MULTIPLE        LOGICAL  DIFF_BASE_MULTIPLE
70        EXTERNAL DIFFERENT_MULTIPLE        EXTERNAL DIFF_BASE_MULTIPLE
71    
72  C     == Local variables ==  C     == Local variables ==
73        _RL      vF (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL      vF (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
# Line 114  C--   (at least in part) Line 114  C--   (at least in part)
114        fVerV(1,1,kUp) = fVerV(1,1,kUp)        fVerV(1,1,kUp) = fVerV(1,1,kUp)
115  #endif  #endif
116    
117        writeDiag = DIFFERENT_MULTIPLE(diagFreq, myTime,        writeDiag = DIFF_BASE_MULTIPLE(baseTime, diagFreq,
118       &                                         myTime-deltaTClock)       &                               myTime, deltaTClock)
119    
120  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
121        IF (useMNC .AND. snapshot_mnc .AND. writeDiag) THEN        IF (useMNC .AND. snapshot_mnc .AND. writeDiag) THEN
122          IF ((bi .EQ. 1).AND.(bj .EQ. 1).AND.(k .EQ. 1)) THEN          IF ((bi .EQ. 1).AND.(bj .EQ. 1).AND.(k .EQ. 1)) THEN
123            CALL MNC_CW_SET_UDIM('mom_vi', -1, myThid)            CALL MNC_CW_SET_UDIM('mom_vi', -1, myThid)
124            CALL MNC_CW_I_W_S('I','mom_vi',0,0,'iter',myIter,myThid)            CALL MNC_CW_I_W_S('I','mom_vi',0,0,'T',myIter,myThid)
125            CALL MNC_CW_SET_UDIM('mom_vi', 0, myThid)            CALL MNC_CW_SET_UDIM('mom_vi', 0, myThid)
126          ENDIF          ENDIF
127          DO i = 1,9          DO i = 1,9
# Line 210  C      Calculate del^2 u and del^2 v for Line 210  C      Calculate del^2 u and del^2 v for
210       &     .OR. viscA4D.NE.0. .OR. viscA4Z.NE.0.       &     .OR. viscA4D.NE.0. .OR. viscA4Z.NE.0.
211       &     .OR. viscA4Grid.NE.0.       &     .OR. viscA4Grid.NE.0.
212       &     .OR. viscC4leith.NE.0.       &     .OR. viscC4leith.NE.0.
213         &     .OR. viscC4leithD.NE.0.
214       &    ) THEN       &    ) THEN
215           CALL MOM_VI_DEL2UV(bi,bj,k,hDiv,vort3,hFacZ,           CALL MOM_VI_DEL2UV(bi,bj,k,hDiv,vort3,hFacZ,
216       O                      del2u,del2v,       O                      del2u,del2v,
# Line 224  C      in terms of vorticity and diverge Line 225  C      in terms of vorticity and diverge
225       &    .OR. viscA4D.NE.0. .OR. viscA4Z.NE.0.       &    .OR. viscA4D.NE.0. .OR. viscA4Z.NE.0.
226       &    .OR. viscAhGrid.NE.0. .OR. viscA4Grid.NE.0.       &    .OR. viscAhGrid.NE.0. .OR. viscA4Grid.NE.0.
227       &    .OR. viscC2leith.NE.0. .OR. viscC4leith.NE.0.       &    .OR. viscC2leith.NE.0. .OR. viscC4leith.NE.0.
228         &    .OR. viscC2leithD.NE.0. .OR. viscC4leithD.NE.0.
229       &    ) THEN       &    ) THEN
230           CALL MOM_VI_HDISSIP(bi,bj,k,hDiv,vort3,hFacZ,dStar,zStar,           CALL MOM_VI_HDISSIP(bi,bj,k,hDiv,vort3,hFacZ,dStar,zStar,
231       O                       guDiss,gvDiss,       O                       guDiss,gvDiss,
232       &                       myThid)       &                       myThid)
233         ENDIF         ENDIF
234  C      or in terms of tension and strain  C      or in terms of tension and strain
235         IF (viscAstrain.NE.0. .OR. viscAtension.NE.0.) THEN         IF (viscAstrain.NE.0. .OR. viscAtension.NE.0.
236         O      .OR. viscC2smag.ne.0) THEN
237           CALL MOM_CALC_TENSION(bi,bj,k,uFld,vFld,           CALL MOM_CALC_TENSION(bi,bj,k,uFld,vFld,
238       O                         tension,       O                         tension,
239       I                         myThid)       I                         myThid)
# Line 445  c      CALL MOM_VI_V_CORIOLIS_C4(bi,bj,K Line 448  c      CALL MOM_VI_V_CORIOLIS_C4(bi,bj,K
448         ENDIF         ENDIF
449    
450  #ifdef ALLOW_TIMEAVE  #ifdef ALLOW_TIMEAVE
451  #ifndef HRCUBE  #ifndef MINIMAL_TAVE_OUTPUT
452         IF (taveFreq.GT.0.) THEN         IF (taveFreq.GT.0.) THEN
453           CALL TIMEAVE_CUMUL_1K1T(uZetatave,vCf,deltaTClock,           CALL TIMEAVE_CUMUL_1K1T(uZetatave,vCf,deltaTClock,
454       &                           Nr, k, bi, bj, myThid)       &                           Nr, k, bi, bj, myThid)
455           CALL TIMEAVE_CUMUL_1K1T(vZetatave,uCf,deltaTClock,           CALL TIMEAVE_CUMUL_1K1T(vZetatave,uCf,deltaTClock,
456       &                           Nr, k, bi, bj, myThid)       &                           Nr, k, bi, bj, myThid)
457         ENDIF         ENDIF
458  #endif /* ndef HRCUBE */  #endif /* ndef MINIMAL_TAVE_OUTPUT */
459  #endif /* ALLOW_TIMEAVE */  #endif /* ALLOW_TIMEAVE */
460    
461  C--   Vertical shear terms (-w*du/dr & -w*dv/dr)  C--   Vertical shear terms (-w*du/dr & -w*dv/dr)
Legend:
Removed from v.1.31  
changed lines
  Added in v.1.36
  ViewVC Help
Powered by ViewVC 1.1.22