/[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.35 by baylor, Thu Mar 10 03:45:11 2005 UTC
# Line 121  C--   (at least in part) Line 121  C--   (at least in part)
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.35
  ViewVC Help
Powered by ViewVC 1.1.22