/[MITgcm]/MITgcm/pkg/thsice/thsice_advection.F
ViewVC logotype

Diff of /MITgcm/pkg/thsice/thsice_advection.F

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

revision 1.13 by gforget, Fri Dec 17 04:00:14 2010 UTC revision 1.14 by heimbach, Wed Jun 1 22:38:25 2011 UTC
# Line 222  C--   set mask West & South Line 222  C--   set mask West & South
222         ENDDO         ENDDO
223        ENDDO        ENDDO
224    
225  #ifndef ALLOW_AUTODIFF_TAMC  cph-exch2#ifndef ALLOW_AUTODIFF_TAMC
226        IF (useCubedSphereExchange) THEN        IF (useCubedSphereExchange) THEN
227         withSigns = .FALSE.         withSigns = .FALSE.
228         CALL FILL_CS_CORNER_UV_RS(         CALL FILL_CS_CORNER_UV_RS(
229       &      withSigns, maskLocW,maskLocS, bi,bj, myThid )       &      withSigns, maskLocW,maskLocS, bi,bj, myThid )
230        ENDIF        ENDIF
231  #endif  cph-exch2#endif
232    
233  C--   Multiple passes for different directions on different tiles  C--   Multiple passes for different directions on different tiles
234  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 297  C-     Advective flux in X Line 297  C-     Advective flux in X
297            ENDDO            ENDDO
298           ENDDO           ENDDO
299    
300  #ifndef ALLOW_AUTODIFF_TAMC  cph-exch2#ifndef ALLOW_AUTODIFF_TAMC
301  C-     Internal exchange for calculations in X  C-     Internal exchange for calculations in X
302           IF ( overlapOnly ) THEN           IF ( overlapOnly ) THEN
303            CALL FILL_CS_CORNER_TR_RL( 1, .FALSE.,            CALL FILL_CS_CORNER_TR_RL( 1, .FALSE.,
# Line 306  C-     Internal exchange for calculation Line 306  C-     Internal exchange for calculation
306       &    CALL FILL_CS_CORNER_TR_RL( 1, .FALSE.,       &    CALL FILL_CS_CORNER_TR_RL( 1, .FALSE.,
307       &                               iceVol, bi,bj, myThid )       &                               iceVol, bi,bj, myThid )
308           ENDIF           ENDIF
309  #endif  cph-exch2#endif
310    
311  C-     Compute CFL number  C-     Compute CFL number
312           IF ( useGridVolume ) THEN           IF ( useGridVolume ) THEN
# Line 354  C-     Compute CFL number Line 354  C-     Compute CFL number
354       & 'THSICE_ADVECTION: adv. scheme incompatibale with multi-dim'       & 'THSICE_ADVECTION: adv. scheme incompatibale with multi-dim'
355           ENDIF           ENDIF
356    
357  #ifndef ALLOW_AUTODIFF_TAMC  cph-exch2#ifndef ALLOW_AUTODIFF_TAMC
358  C--   Internal exchange for next calculations in Y  C--   Internal exchange for next calculations in Y
359           IF ( overlapOnly .AND. ipass.EQ.1 ) THEN           IF ( overlapOnly .AND. ipass.EQ.1 ) THEN
360            CALL FILL_CS_CORNER_TR_RL( 2, .FALSE.,            CALL FILL_CS_CORNER_TR_RL( 2, .FALSE.,
# Line 363  C--   Internal exchange for next calcula Line 363  C--   Internal exchange for next calcula
363       &    CALL FILL_CS_CORNER_TR_RL( 2, .FALSE.,       &    CALL FILL_CS_CORNER_TR_RL( 2, .FALSE.,
364       &                               iceVol, bi,bj, myThid )       &                               iceVol, bi,bj, myThid )
365           ENDIF           ENDIF
366  #endif  cph-exch2#endif
367    
368  C--   Advective flux in X : done  C--   Advective flux in X : done
369          ENDIF          ENDIF
# Line 512  C-     Advective flux in Y Line 512  C-     Advective flux in Y
512            ENDDO            ENDDO
513           ENDDO           ENDDO
514    
515  #ifndef ALLOW_AUTODIFF_TAMC  cph-exch2#ifndef ALLOW_AUTODIFF_TAMC
516  C-     Internal exchange for calculations in Y  C-     Internal exchange for calculations in Y
517           IF ( overlapOnly ) THEN           IF ( overlapOnly ) THEN
518            CALL FILL_CS_CORNER_TR_RL( 2, .FALSE.,            CALL FILL_CS_CORNER_TR_RL( 2, .FALSE.,
# Line 521  C-     Internal exchange for calculation Line 521  C-     Internal exchange for calculation
521       &    CALL FILL_CS_CORNER_TR_RL( 2, .FALSE.,       &    CALL FILL_CS_CORNER_TR_RL( 2, .FALSE.,
522       &                               iceVol, bi,bj, myThid )       &                               iceVol, bi,bj, myThid )
523           ENDIF           ENDIF
524  #endif  cph-exch2#endif
525    
526  C-     Compute CFL number  C-     Compute CFL number
527           IF ( useGridVolume ) THEN           IF ( useGridVolume ) THEN
# Line 569  C-     Compute CFL number Line 569  C-     Compute CFL number
569       &  'THSICE_ADVECTION: adv. scheme incompatibale with mutli-dim'       &  'THSICE_ADVECTION: adv. scheme incompatibale with mutli-dim'
570           ENDIF           ENDIF
571    
572  #ifndef ALLOW_AUTODIFF_TAMC  cph-exch2#ifndef ALLOW_AUTODIFF_TAMC
573           IF ( overlapOnly .AND. ipass.EQ.1 ) THEN           IF ( overlapOnly .AND. ipass.EQ.1 ) THEN
574            CALL FILL_CS_CORNER_TR_RL( 1, .FALSE.,            CALL FILL_CS_CORNER_TR_RL( 1, .FALSE.,
575       &                               iceFld, bi,bj, myThid )       &                               iceFld, bi,bj, myThid )
# Line 577  C-     Compute CFL number Line 577  C-     Compute CFL number
577       &    CALL FILL_CS_CORNER_TR_RL( 1, .FALSE.,       &    CALL FILL_CS_CORNER_TR_RL( 1, .FALSE.,
578       &                               iceVol, bi,bj, myThid )       &                               iceVol, bi,bj, myThid )
579           ENDIF           ENDIF
580  #endif  cph-exch2#endif
581    
582  C-     Advective flux in Y : done  C-     Advective flux in Y : done
583          ENDIF          ENDIF
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22