/[MITgcm]/MITgcm/pkg/seaice/seaice_preconditioner.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/seaice_preconditioner.F

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

revision 1.17 by mlosch, Wed Feb 13 13:26:50 2013 UTC revision 1.18 by mlosch, Thu Feb 14 17:03:34 2013 UTC
# Line 214  C--     prevent tri-diagonal solver from Line 214  C--     prevent tri-diagonal solver from
214         ENDDO         ENDDO
215        ENDDO        ENDDO
216  #endif /* OBCS_UVICE_OLD */  #endif /* OBCS_UVICE_OLD */
 C     this is a hack to avoid divisions by zero  
 C     I do not know why I need BU/BV are zero in some places  
       DO bj=myByLo(myThid),myByHi(myThid)  
        DO bi=myBxLo(myThid),myBxHi(myThid)  
         DO J=jMin,jMax  
          DO I=iMin,iMax  
           IF ( BU(I,J,bi,bj) .EQ. 0. _d 0 ) BU(I,J,bi,bj) = 1. _d 0  
           IF ( BV(I,J,bi,bj) .EQ. 0. _d 0 ) BV(I,J,bi,bj) = 1. _d 0  
          ENDDO  
         ENDDO  
        ENDDO  
       ENDDO  
217    
218  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
219    
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22