/[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.39 by heimbach, Thu Dec 8 15:44:34 2005 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    cph-test
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 /* ALLOW_AUTODIFF_TAMC */
334    C
335        IF (calc_fluxes_X) THEN        IF (calc_fluxes_X) THEN
336    
337  C-     Do not compute fluxes if  C-     Do not compute fluxes if
# Line 333  C-     Internal exchange for calculation Line 351  C-     Internal exchange for calculation
351          ENDIF          ENDIF
352  #endif  #endif
353    
 C-     Advective flux in X  
         DO j=1-Oly,sNy+Oly  
          DO i=1-Olx,sNx+Olx  
           af(i,j) = 0.  
          ENDDO  
         ENDDO  
   
354  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
355  #ifndef DISABLE_MULTIDIM_ADVECTION  # ifndef DISABLE_MULTIDIM_ADVECTION
356  CADJ STORE localTij(:,:)  =  CADJ STORE localTij(:,:)  =
357  CADJ &     comlev1_bibj_k_gad_pass, key=passkey, byte=isbyte  CADJ &     comlev1_bibj_k_gad_pass, key=passkey, byte=isbyte
358  #endif  # endif
359  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
360    
361          IF ( advectionScheme.EQ.ENUM_UPWIND_1RST          IF ( advectionScheme.EQ.ENUM_UPWIND_1RST
# Line 461  C--   End of X direction Line 472  C--   End of X direction
472    
473  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
474  C--   Y direction  C--   Y direction
475    cph-test
476    C-     Advective flux in Y
477            DO j=1-Oly,sNy+Oly
478             DO i=1-Olx,sNx+Olx
479              af(i,j) = 0.
480             ENDDO
481            ENDDO
482    C
483    #ifdef ALLOW_AUTODIFF_TAMC
484    cph-test
485    CADJ STORE localTij(:,:)  =
486    CADJ &     comlev1_bibj_k_gad_pass, key=passkey, byte=isbyte
487    CADJ STORE af(:,:)  =
488    CADJ &     comlev1_bibj_k_gad_pass, key=passkey, byte=isbyte
489    #endif /* ALLOW_AUTODIFF_TAMC */
490    C
491        IF (calc_fluxes_Y) THEN        IF (calc_fluxes_Y) THEN
492    
493  C-     Do not compute fluxes if  C-     Do not compute fluxes if
Legend:
Removed from v.1.38  
changed lines
  Added in v.1.39
  ViewVC Help
Powered by ViewVC 1.1.22