/[MITgcm]/MITgcm/model/src/dynamics.F
ViewVC logotype

Diff of /MITgcm/model/src/dynamics.F

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

revision 1.159 by mlosch, Tue Oct 25 15:09:49 2011 UTC revision 1.160 by jmc, Thu Dec 1 14:22:27 2011 UTC
# Line 406  CADJ &     = comlev1_bibj, key=idynkey, Line 406  CADJ &     = comlev1_bibj, key=idynkey,
406  C--   For Stevens boundary conditions velocities need to be extrapolated  C--   For Stevens boundary conditions velocities need to be extrapolated
407  C     (copied) to a narrow strip outside the domain  C     (copied) to a narrow strip outside the domain
408           IF ( useOBCS ) THEN           IF ( useOBCS ) THEN
409            CALL OBCS_COPY_UV_N(            CALL OBCS_COPY_UV_N(
410       U         uVel(1-Olx,1-Oly,1,bi,bj),       U         uVel(1-Olx,1-Oly,1,bi,bj),
411       U         vVel(1-Olx,1-Oly,1,bi,bj),       U         vVel(1-Olx,1-Oly,1,bi,bj),
412       I         Nr, bi, bj, myThid )       I         Nr, bi, bj, myThid )
# Line 594  CADJ STORE gU(:,:,:,bi,bj) = comlev1_bib Line 594  CADJ STORE gU(:,:,:,bi,bj) = comlev1_bib
594  #endif    /* ALLOW_AUTODIFF_TAMC */  #endif    /* ALLOW_AUTODIFF_TAMC */
595            CALL IMPLDIFF(            CALL IMPLDIFF(
596       I         bi, bj, iMin, iMax, jMin, jMax,       I         bi, bj, iMin, iMax, jMin, jMax,
597       I         -1, KappaRU,recip_HFacW,       I         -1, KappaRU, recip_hFacW(1-OLx,1-OLy,1,bi,bj),
598       U         gU,       U         gU,
599       I         myThid )       I         myThid )
600  #ifdef    ALLOW_AUTODIFF_TAMC  #ifdef    ALLOW_AUTODIFF_TAMC
# Line 603  CADJ STORE gV(:,:,:,bi,bj) = comlev1_bib Line 603  CADJ STORE gV(:,:,:,bi,bj) = comlev1_bib
603  #endif    /* ALLOW_AUTODIFF_TAMC */  #endif    /* ALLOW_AUTODIFF_TAMC */
604            CALL IMPLDIFF(            CALL IMPLDIFF(
605       I         bi, bj, iMin, iMax, jMin, jMax,       I         bi, bj, iMin, iMax, jMin, jMax,
606       I         -2, KappaRV,recip_HFacS,       I         -2, KappaRV, recip_hFacS(1-OLx,1-OLy,1,bi,bj),
607       U         gV,       U         gV,
608       I         myThid )       I         myThid )
609          ENDIF          ENDIF
# Line 611  CADJ STORE gV(:,:,:,bi,bj) = comlev1_bib Line 611  CADJ STORE gV(:,:,:,bi,bj) = comlev1_bib
611  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
612  C--      Apply open boundary conditions  C--      Apply open boundary conditions
613          IF ( useOBCS ) THEN          IF ( useOBCS ) THEN
614  C--      but first save intermediate velocities to be used in the  C--      but first save intermediate velocities to be used in the
615  C        next time step for the Stevens boundary conditions  C        next time step for the Stevens boundary conditions
616            CALL OBCS_SAVE_UV_N(            CALL OBCS_SAVE_UV_N(
617       I        bi, bj, iMin, iMax, jMin, jMax, 0,       I        bi, bj, iMin, iMax, jMin, jMax, 0,
618       I        gU, gV, myThid )       I        gU, gV, myThid )
619            CALL OBCS_APPLY_UV( bi, bj, 0, gU, gV, myThid )            CALL OBCS_APPLY_UV( bi, bj, 0, gU, gV, myThid )
620          ENDIF          ENDIF
# Line 627  CADJ STORE vVelD(:,:,:,bi,bj) = comlev1_ Line 627  CADJ STORE vVelD(:,:,:,bi,bj) = comlev1_
627  #endif    /* ALLOW_AUTODIFF_TAMC */  #endif    /* ALLOW_AUTODIFF_TAMC */
628            CALL IMPLDIFF(            CALL IMPLDIFF(
629       I         bi, bj, iMin, iMax, jMin, jMax,       I         bi, bj, iMin, iMax, jMin, jMax,
630       I         0, KappaRU,recip_HFacW,       I         0, KappaRU, recip_hFacW(1-OLx,1-OLy,1,bi,bj),
631       U         vVelD,       U         vVelD,
632       I         myThid )       I         myThid )
633  #ifdef    ALLOW_AUTODIFF_TAMC  #ifdef    ALLOW_AUTODIFF_TAMC
# Line 635  CADJ STORE uVelD(:,:,:,bi,bj) = comlev1_ Line 635  CADJ STORE uVelD(:,:,:,bi,bj) = comlev1_
635  #endif    /* ALLOW_AUTODIFF_TAMC */  #endif    /* ALLOW_AUTODIFF_TAMC */
636            CALL IMPLDIFF(            CALL IMPLDIFF(
637       I         bi, bj, iMin, iMax, jMin, jMax,       I         bi, bj, iMin, iMax, jMin, jMax,
638       I         0, KappaRV,recip_HFacS,       I         0, KappaRV, recip_hFacS(1-OLx,1-OLy,1,bi,bj),
639       U         uVelD,       U         uVelD,
640       I         myThid )       I         myThid )
641          ENDIF          ENDIF
Legend:
Removed from v.1.159  
changed lines
  Added in v.1.160
  ViewVC Help
Powered by ViewVC 1.1.22