/[MITgcm]/MITgcm_contrib/cg2d_bench/ini_cg2d.F
ViewVC logotype

Diff of /MITgcm_contrib/cg2d_bench/ini_cg2d.F

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

revision 1.1 by ce107, Fri May 12 21:58:05 2006 UTC revision 1.2 by ce107, Fri May 12 22:23:42 2006 UTC
# Line 1  Line 1 
1    C       $Id$    
2  CStartOfInterface  CStartOfInterface
3        SUBROUTINE INI_CG2D        SUBROUTINE INI_CG2D
4  C     /==========================================================\  C     /==========================================================\
# Line 42  C     aS2d: integral Ay/dY Line 43  C     aS2d: integral Ay/dY
43          aW2d(I,J) = 1. _d 0          aW2d(I,J) = 1. _d 0
44          aS2d(I,J) = 1. _d 0          aS2d(I,J) = 1. _d 0
45          IF ( jG .EQ. 1 .OR. jG .EQ. nY+1 ) THEN          IF ( jG .EQ. 1 .OR. jG .EQ. nY+1 ) THEN
46           aS2d(I,J) = 0.           aS2d(I,J) = 0. _d 0
47          ENDIF          ENDIF
48         ENDDO         ENDDO
49        ENDDO        ENDDO
# Line 56  C--   Initialise preconditioner Line 57  C--   Initialise preconditioner
57          IF (          IF (
58       &   aW2d(I,J) + aW2d(I+1,J)       &   aW2d(I,J) + aW2d(I+1,J)
59       &  +aS2d(I,J) + aS2D(I,J+1)       &  +aS2d(I,J) + aS2D(I,J+1)
60       &   .EQ. 0.       &   .EQ. 0. _d 0
61       &     )  pC(I,J) = 0. _d 0       &     )  pC(I,J) = 0. _d 0
62          pW(I,J) = 0.          pW(I,J) = 0. _d 0
63          pS(I,J) = 0.          pS(I,J) = 0. _d 0
64         ENDDO         ENDDO
65        ENDDO        ENDDO
66  C--   Update overlap regions  C--   Update overlap regions
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22