/[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.30 by jmc, Mon Sep 27 14:42:40 2004 UTC revision 1.31 by heimbach, Wed Sep 29 04:53:30 2004 UTC
# Line 244  C--   Make local copy of tracer array an Line 244  C--   Make local copy of tracer array an
244         ENDDO         ENDDO
245        ENDDO        ENDDO
246    
247    #ifndef ALLOW_AUTODIFF_TAMC
248        IF (useCubedSphereExchange) THEN        IF (useCubedSphereExchange) THEN
249          withSigns = .FALSE.          withSigns = .FALSE.
250          CALL FILL_CS_CORNER_UV_RS(          CALL FILL_CS_CORNER_UV_RS(
251       &            withSigns, maskLocW,maskLocS, bi,bj, myThid )       &            withSigns, maskLocW,maskLocS, bi,bj, myThid )
252        ENDIF        ENDIF
253    #endif
254    
255  C--   Multiple passes for different directions on different tiles  C--   Multiple passes for different directions on different tiles
256  C--   For cube need one pass for each of red, green and blue axes.  C--   For cube need one pass for each of red, green and blue axes.
# Line 303  C       a) needed in overlap only Line 305  C       a) needed in overlap only
305  C   and b) the overlap of myTile are not cube-face Edges  C   and b) the overlap of myTile are not cube-face Edges
306         IF ( .NOT.overlapOnly .OR. N_edge .OR. S_edge ) THEN         IF ( .NOT.overlapOnly .OR. N_edge .OR. S_edge ) THEN
307    
308    #ifndef ALLOW_AUTODIFF_TAMC
309  C-     Internal exchange for calculations in X  C-     Internal exchange for calculations in X
310  #ifdef MULTIDIM_OLD_VERSION  #ifdef MULTIDIM_OLD_VERSION
311          IF ( useCubedSphereExchange ) THEN          IF ( useCubedSphereExchange ) THEN
# Line 312  C-     Internal exchange for calculation Line 315  C-     Internal exchange for calculation
315  #endif  #endif
316           CALL FILL_CS_CORNER_TR_RL( .TRUE., localTij, bi,bj, myThid )           CALL FILL_CS_CORNER_TR_RL( .TRUE., localTij, bi,bj, myThid )
317          ENDIF          ENDIF
318    #endif
319    
320  C-     Advective flux in X  C-     Advective flux in X
321          DO j=1-Oly,sNy+Oly          DO j=1-Oly,sNy+Oly
# Line 346  CADJ &     comlev1_bibj_k_gad_pass, key= Line 350  CADJ &     comlev1_bibj_k_gad_pass, key=
350  C-     Advective flux in X : done  C-     Advective flux in X : done
351         ENDIF         ENDIF
352    
353    #ifndef ALLOW_AUTODIFF_TAMC
354  C-     Internal exchange for next calculations in Y  C-     Internal exchange for next calculations in Y
355         IF ( overlapOnly .AND. ipass.EQ.1 ) THEN         IF ( overlapOnly .AND. ipass.EQ.1 ) THEN
356           CALL FILL_CS_CORNER_TR_RL(.FALSE., localTij, bi,bj, myThid )           CALL FILL_CS_CORNER_TR_RL(.FALSE., localTij, bi,bj, myThid )
357         ENDIF         ENDIF
358    #endif
359    
360  C-     Update the local tracer field where needed:  C-     Update the local tracer field where needed:
361    
# Line 436  C       a) needed in overlap only Line 442  C       a) needed in overlap only
442  C   and b) the overlap of myTile are not cube-face edges  C   and b) the overlap of myTile are not cube-face edges
443         IF ( .NOT.overlapOnly .OR. E_edge .OR. W_edge ) THEN         IF ( .NOT.overlapOnly .OR. E_edge .OR. W_edge ) THEN
444    
445    #ifndef ALLOW_AUTODIFF_TAMC
446  C-     Internal exchange for calculations in Y  C-     Internal exchange for calculations in Y
447  #ifdef MULTIDIM_OLD_VERSION  #ifdef MULTIDIM_OLD_VERSION
448          IF ( useCubedSphereExchange ) THEN          IF ( useCubedSphereExchange ) THEN
# Line 445  C-     Internal exchange for calculation Line 452  C-     Internal exchange for calculation
452  #endif  #endif
453           CALL FILL_CS_CORNER_TR_RL(.FALSE., localTij, bi,bj, myThid )           CALL FILL_CS_CORNER_TR_RL(.FALSE., localTij, bi,bj, myThid )
454          ENDIF          ENDIF
455    #endif
456    
457  C-     Advective flux in Y  C-     Advective flux in Y
458          DO j=1-Oly,sNy+Oly          DO j=1-Oly,sNy+Oly
# Line 479  CADJ &     comlev1_bibj_k_gad_pass, key= Line 487  CADJ &     comlev1_bibj_k_gad_pass, key=
487  C-     Advective flux in Y : done  C-     Advective flux in Y : done
488         ENDIF         ENDIF
489    
490    #ifndef ALLOW_AUTODIFF_TAMC
491  C-     Internal exchange for next calculations in X  C-     Internal exchange for next calculations in X
492         IF ( overlapOnly .AND. ipass.EQ.1 ) THEN         IF ( overlapOnly .AND. ipass.EQ.1 ) THEN
493           CALL FILL_CS_CORNER_TR_RL( .TRUE., localTij, bi,bj, myThid )           CALL FILL_CS_CORNER_TR_RL( .TRUE., localTij, bi,bj, myThid )
494         ENDIF         ENDIF
495    #endif
496    
497  C-     Update the local tracer field where needed:  C-     Update the local tracer field where needed:
498    
Legend:
Removed from v.1.30  
changed lines
  Added in v.1.31
  ViewVC Help
Powered by ViewVC 1.1.22