/[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.24 by edhill, Thu Oct 7 21:52:29 2004 UTC revision 1.25 by edhill, Fri Oct 8 17:03:21 2004 UTC
# Line 126  C     xxxFac - On-off tracer parameters Line 126  C     xxxFac - On-off tracer parameters
126        _RL vort3(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL vort3(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
127        _RL hDiv(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL hDiv(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
128    
129    #ifdef ALLOW_MNC
130          INTEGER offsets(9)
131    #endif
132    
133  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
134  C--   only the kDown part of fverU/V is set in this subroutine  C--   only the kDown part of fverU/V is set in this subroutine
135  C--   the kUp is still required  C--   the kUp is still required
# Line 144  C--   (at least in part) Line 148  C--   (at least in part)
148    
149  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
150        IF (useMNC .AND. snapshot_mnc .AND. writeDiag) THEN        IF (useMNC .AND. snapshot_mnc .AND. writeDiag) THEN
151          CALL MNC_CW_SET_UDIM('mom_vi', -1, myThid)          IF ((bi .EQ. 1).AND.(bj .EQ. 1).AND.(k .EQ. 1)) THEN
152          CALL MNC_CW_I_W_S('I','mom_vi',0,0,'iter',myIter,myThid)            CALL MNC_CW_SET_UDIM('mom_vi', -1, myThid)
153          CALL MNC_CW_SET_UDIM('mom_vi', 0, myThid)            CALL MNC_CW_I_W_S('I','mom_vi',0,0,'iter',myIter,myThid)
154              CALL MNC_CW_SET_UDIM('mom_vi', 0, myThid)
155            ENDIF
156            DO i = 1,9
157              offsets(i) = 0
158            ENDDO
159            offsets(3) = k
160    C       write(*,*) 'offsets = ',(offsets(i),i=1,9)
161        ENDIF        ENDIF
162  #endif /*  ALLOW_MNC  */  #endif /*  ALLOW_MNC  */
163    
# Line 430  C--   Horizontal Coriolis terms Line 441  C--   Horizontal Coriolis terms
441           ENDIF           ENDIF
442  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
443           IF (useMNC .AND. snapshot_mnc) THEN           IF (useMNC .AND. snapshot_mnc) THEN
444             CALL MNC_CW_RL_W('D','mom_vi',0,0, 'fV', uCf, myThid)             CALL MNC_CW_RL_W_OFFSET('D','mom_vi',bi,bj, 'fV', uCf,
445             CALL MNC_CW_RL_W('D','mom_vi',0,0, 'fU', vCf, myThid)       &          offsets, myThid)
446               CALL MNC_CW_RL_W_OFFSET('D','mom_vi',bi,bj, 'fU', vCf,
447         &          offsets, myThid)
448           ENDIF           ENDIF
449  #endif /*  ALLOW_MNC  */  #endif /*  ALLOW_MNC  */
450         ENDIF         ENDIF
# Line 473  c      CALL MOM_VI_V_CORIOLIS_C4(bi,bj,K Line 486  c      CALL MOM_VI_V_CORIOLIS_C4(bi,bj,K
486           ENDIF           ENDIF
487  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
488           IF (useMNC .AND. snapshot_mnc) THEN           IF (useMNC .AND. snapshot_mnc) THEN
489             CALL MNC_CW_RL_W('D','mom_vi',0,0, 'zV', uCf, myThid)             CALL MNC_CW_RL_W_OFFSET('D','mom_vi',bi,bj, 'zV', uCf,
490             CALL MNC_CW_RL_W('D','mom_vi',0,0, 'zU', vCf, myThid)       &          offsets, myThid)
491               CALL MNC_CW_RL_W_OFFSET('D','mom_vi',bi,bj, 'zU', vCf,
492         &          offsets, myThid)
493           ENDIF           ENDIF
494  #endif /*  ALLOW_MNC  */  #endif /*  ALLOW_MNC  */
495         ENDIF         ENDIF
# Line 526  C--   Bernoulli term Line 541  C--   Bernoulli term
541           ENDIF           ENDIF
542  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
543           IF (useMNC .AND. snapshot_mnc) THEN           IF (useMNC .AND. snapshot_mnc) THEN
544             CALL MNC_CW_RL_W('D','mom_vi',0,0, 'KEx', uCf, myThid)             CALL MNC_CW_RL_W_OFFSET('D','mom_vi',bi,bj, 'KEx', uCf,
545             CALL MNC_CW_RL_W('D','mom_vi',0,0, 'KEy', vCf, myThid)       &          offsets, myThid)
546           ENDIF             CALL MNC_CW_RL_W_OFFSET('D','mom_vi',bi,bj, 'KEy', vCf,
547         &          offsets, myThid)
548            ENDIF
549  #endif /*  ALLOW_MNC  */  #endif /*  ALLOW_MNC  */
550         ENDIF         ENDIF
551    
# Line 567  C--   Set du/dt & dv/dt on boundaries to Line 584  C--   Set du/dt & dv/dt on boundaries to
584          ENDIF          ENDIF
585  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
586          IF (useMNC .AND. snapshot_mnc) THEN          IF (useMNC .AND. snapshot_mnc) THEN
587            CALL MNC_CW_RL_W('D','mom_vi',0,0,'Ds',strain, myThid)            CALL MNC_CW_RL_W_OFFSET('D','mom_vi',bi,bj,'Ds',strain,
588            CALL MNC_CW_RL_W('D','mom_vi',0,0,'Dt',tension, myThid)       &          offsets, myThid)
589            CALL MNC_CW_RL_W('D','mom_vi',0,0,'Du',uDiss, myThid)            CALL MNC_CW_RL_W_OFFSET('D','mom_vi',bi,bj,'Dt',tension,
590            CALL MNC_CW_RL_W('D','mom_vi',0,0,'Dv',vDiss, myThid)       &          offsets, myThid)
591            CALL MNC_CW_RL_W('D','mom_vi',0,0,'Z3',vort3, myThid)            CALL MNC_CW_RL_W_OFFSET('D','mom_vi',bi,bj,'Du',uDiss,
592            CALL MNC_CW_RL_W('D','mom_vi',0,0,'W3',omega3, myThid)       &          offsets, myThid)
593            CALL MNC_CW_RL_W('D','mom_vi',0,0,'KE',KE, myThid)            CALL MNC_CW_RL_W_OFFSET('D','mom_vi',bi,bj,'Dv',vDiss,
594            CALL MNC_CW_RL_W('D','mom_vi',0,0,'D', hdiv, myThid)       &          offsets, myThid)
595              CALL MNC_CW_RL_W_OFFSET('D','mom_vi',bi,bj,'Z3',vort3,
596         &          offsets, myThid)
597              CALL MNC_CW_RL_W_OFFSET('D','mom_vi',bi,bj,'W3',omega3,
598         &          offsets, myThid)
599              CALL MNC_CW_RL_W_OFFSET('D','mom_vi',bi,bj,'KE',KE,
600         &          offsets, myThid)
601              CALL MNC_CW_RL_W_OFFSET('D','mom_vi',bi,bj,'D', hdiv,
602         &          offsets, myThid)
603          ENDIF          ENDIF
604  #endif /*  ALLOW_MNC  */  #endif /*  ALLOW_MNC  */
605        ENDIF        ENDIF
Legend:
Removed from v.1.24  
changed lines
  Added in v.1.25
  ViewVC Help
Powered by ViewVC 1.1.22