/[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.136 by jmc, Fri Jul 7 20:10:35 2006 UTC revision 1.139 by jmc, Mon Nov 5 18:52:21 2007 UTC
# Line 93  C     == Global variables === Line 93  C     == Global variables ===
93  # endif  # endif
94  # ifdef ALLOW_PTRACERS  # ifdef ALLOW_PTRACERS
95  #  include "PTRACERS_SIZE.h"  #  include "PTRACERS_SIZE.h"
96  #  include "PTRACERS.h"  #  include "PTRACERS_FIELDS.h"
97  # endif  # endif
98  # ifdef ALLOW_OBCS  # ifdef ALLOW_OBCS
99  #  include "OBCS.h"  #  include "OBCS.h"
# Line 317  cph) Line 317  cph)
317            guDissip(i,j)  = 0. _d 0            guDissip(i,j)  = 0. _d 0
318            gvDissip(i,j)  = 0. _d 0            gvDissip(i,j)  = 0. _d 0
319  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
 cph(  
 c--   need some re-initialisation here to break dependencies  
 cph)  
320  # ifdef NONLIN_FRSURF  # ifdef NONLIN_FRSURF
321  #  ifndef DISABLE_RSTAR_CODE  #  ifndef DISABLE_RSTAR_CODE
322            dWtransC(i,j,bi,bj)  = 0. _d 0            dWtransC(i,j,bi,bj) = 0. _d 0
323            dWtransU(i,j,bi,bj)  = 0. _d 0            dWtransU(i,j,bi,bj) = 0. _d 0
324            dWtransV(i,j,bi,bj)  = 0. _d 0            dWtransV(i,j,bi,bj) = 0. _d 0
325  #  endif  #  endif
326  # endif /* NONLIN_FRSURF */  # endif
327  #endif /* ALLOW_AUTODIFF_TAMC */  #endif
328           ENDDO           ENDDO
329          ENDDO          ENDDO
330    
# Line 425  CADJ STORE gunm1(:,:,k,bi,bj) Line 422  CADJ STORE gunm1(:,:,k,bi,bj)
422  CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte
423  CADJ STORE gvnm1(:,:,k,bi,bj)  CADJ STORE gvnm1(:,:,k,bi,bj)
424  CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte
 #  ifndef DISABLE_RSTAR_CODE  
 CADJ STORE dwtransc(:,:,bi,bj)  
 CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte  
 CADJ STORE dwtransu(:,:,bi,bj)  
 CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte  
 CADJ STORE dwtransv(:,:,bi,bj)  
 CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte  
 #  endif  
425  #  ifdef ALLOW_CD_CODE  #  ifdef ALLOW_CD_CODE
426  CADJ STORE unm1(:,:,k,bi,bj)  CADJ STORE unm1(:,:,k,bi,bj)
427  CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte
# Line 473  C        phiHyd(z=0)=0 Line 462  C        phiHyd(z=0)=0
462  C--      Calculate accelerations in the momentum equations (gU, gV, ...)  C--      Calculate accelerations in the momentum equations (gU, gV, ...)
463  C        and step forward storing the result in gU, gV, etc...  C        and step forward storing the result in gU, gV, etc...
464           IF ( momStepping ) THEN           IF ( momStepping ) THEN
465             IF (.NOT. vectorInvariantMomentum) THEN  #ifdef ALLOW_AUTODIFF_TAMC
466  #ifdef ALLOW_MOM_FLUXFORM  # ifdef NONLIN_FRSURF
467  C  #  ifndef DISABLE_RSTAR_CODE
468  # ifdef ALLOW_AUTODIFF_TAMC  CADJ STORE dWtransC(:,:,bi,bj)
 #  ifdef NONLIN_FRSURF  
 #   ifndef DISABLE_RSTAR_CODE  
 CADJ STORE dwtransc(:,:,bi,bj)  
469  CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte
470  CADJ STORE dwtransu(:,:,bi,bj)  CADJ STORE dWtransU(:,:,bi,bj)
471  CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte
472  CADJ STORE dwtransv(:,:,bi,bj)  CADJ STORE dWtransV(:,:,bi,bj)
473  CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte
 #   endif  
474  #  endif  #  endif
475  # endif /* ALLOW_AUTODIFF_TAMC */  # endif
476    #endif
477               IF (.NOT. vectorInvariantMomentum) THEN
478    #ifdef ALLOW_MOM_FLUXFORM
479  C  C
480                CALL MOM_FLUXFORM(                CALL MOM_FLUXFORM(
481       I         bi,bj,iMin,iMax,jMin,jMax,k,kup,kDown,       I         bi,bj,iMin,iMax,jMin,jMax,k,kup,kDown,

Legend:
Removed from v.1.136  
changed lines
  Added in v.1.139

  ViewVC Help
Powered by ViewVC 1.1.22