/[MITgcm]/MITgcm/pkg/generic_advdiff/gad_advection.F
ViewVC logotype

Diff of /MITgcm/pkg/generic_advdiff/gad_advection.F

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

revision 1.38 by jmc, Sat Nov 5 01:01:51 2005 UTC revision 1.40 by heimbach, Tue Feb 21 17:20:12 2006 UTC
# Line 193  C     uninitialised but inert locations. Line 193  C     uninitialised but inert locations.
193          fVerT(i,j,1) = 0. _d 0          fVerT(i,j,1) = 0. _d 0
194          fVerT(i,j,2) = 0. _d 0          fVerT(i,j,2) = 0. _d 0
195          rTransKp1(i,j)= 0. _d 0          rTransKp1(i,j)= 0. _d 0
196    #ifdef ALLOW_AUTODIFF_TAMC
197            localTij(i,j) = 0. _d 0
198    #endif
199         ENDDO         ENDDO
200        ENDDO        ENDDO
201    
# Line 314  C-    not CubedSphere Line 317  C-    not CubedSphere
317    
318  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
319  C--   X direction  C--   X direction
320    C-     Advective flux in X
321            DO j=1-Oly,sNy+Oly
322             DO i=1-Olx,sNx+Olx
323              af(i,j) = 0.
324             ENDDO
325            ENDDO
326    C
327    #ifdef ALLOW_AUTODIFF_TAMC
328    # ifndef DISABLE_MULTIDIM_ADVECTION
329    CADJ STORE localTij(:,:)  =
330    CADJ &     comlev1_bibj_k_gad_pass, key=passkey, byte=isbyte
331    CADJ STORE af(:,:)  =
332    CADJ &     comlev1_bibj_k_gad_pass, key=passkey, byte=isbyte
333    # endif
334    #endif /* ALLOW_AUTODIFF_TAMC */
335    C
336        IF (calc_fluxes_X) THEN        IF (calc_fluxes_X) THEN
337    
338  C-     Do not compute fluxes if  C-     Do not compute fluxes if
# Line 333  C-     Internal exchange for calculation Line 352  C-     Internal exchange for calculation
352          ENDIF          ENDIF
353  #endif  #endif
354    
 C-     Advective flux in X  
         DO j=1-Oly,sNy+Oly  
          DO i=1-Olx,sNx+Olx  
           af(i,j) = 0.  
          ENDDO  
         ENDDO  
   
355  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
356  #ifndef DISABLE_MULTIDIM_ADVECTION  # ifndef DISABLE_MULTIDIM_ADVECTION
357  CADJ STORE localTij(:,:)  =  CADJ STORE localTij(:,:)  =
358  CADJ &     comlev1_bibj_k_gad_pass, key=passkey, byte=isbyte  CADJ &     comlev1_bibj_k_gad_pass, key=passkey, byte=isbyte
359  #endif  # endif
360  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
361    
362          IF ( advectionScheme.EQ.ENUM_UPWIND_1RST          IF ( advectionScheme.EQ.ENUM_UPWIND_1RST
# Line 461  C--   End of X direction Line 473  C--   End of X direction
473    
474  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
475  C--   Y direction  C--   Y direction
476    cph-test
477    C-     Advective flux in Y
478            DO j=1-Oly,sNy+Oly
479             DO i=1-Olx,sNx+Olx
480              af(i,j) = 0.
481             ENDDO
482            ENDDO
483    C
484    #ifdef ALLOW_AUTODIFF_TAMC
485    # ifndef DISABLE_MULTIDIM_ADVECTION
486    CADJ STORE localTij(:,:)  =
487    CADJ &     comlev1_bibj_k_gad_pass, key=passkey, byte=isbyte
488    CADJ STORE af(:,:)  =
489    CADJ &     comlev1_bibj_k_gad_pass, key=passkey, byte=isbyte
490    # endif
491    #endif /* ALLOW_AUTODIFF_TAMC */
492    C
493        IF (calc_fluxes_Y) THEN        IF (calc_fluxes_Y) THEN
494    
495  C-     Do not compute fluxes if  C-     Do not compute fluxes if
Legend:
Removed from v.1.38  
changed lines
  Added in v.1.40
  ViewVC Help
Powered by ViewVC 1.1.22