/[MITgcm]/MITgcm/model/src/ini_cg2d.F
ViewVC logotype

Diff of /MITgcm/model/src/ini_cg2d.F

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

revision 1.46 by jmc, Tue Apr 28 18:01:14 2009 UTC revision 1.47 by jmc, Sat Jul 11 22:00:40 2009 UTC
# Line 11  C     !INTERFACE: Line 11  C     !INTERFACE:
11    
12  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
13  C     *==========================================================*  C     *==========================================================*
14  C     | SUBROUTINE INI_CG2D                                        C     | SUBROUTINE INI_CG2D
15  C     | o Initialise 2d conjugate gradient solver operators.        C     | o Initialise 2d conjugate gradient solver operators.
16  C     *==========================================================*  C     *==========================================================*
17  C     | These arrays are purely a function of the basin geom.      C     | These arrays are purely a function of the basin geom.
18  C     | We set then here once and them use then repeatedly.        C     | We set then here once and them use then repeatedly.
19  C     *==========================================================*  C     *==========================================================*
20  C     \ev  C     \ev
21    
# Line 64  C     but safer when EXCH do not fill al Line 64  C     but safer when EXCH do not fill al
64            pW(I,J,bi,bj) = 0. _d 0            pW(I,J,bi,bj) = 0. _d 0
65            pS(I,J,bi,bj) = 0. _d 0            pS(I,J,bi,bj) = 0. _d 0
66            pC(I,J,bi,bj) = 0. _d 0            pC(I,J,bi,bj) = 0. _d 0
           cg2d_q(I,J,bi,bj) = 0. _d 0  
67           ENDDO           ENDDO
68          ENDDO          ENDDO
69          DO J=1-1,sNy+1          DO J=1-1,sNy+1
70           DO I=1-1,sNx+1           DO I=1-1,sNx+1
71              cg2d_q(I,J,bi,bj) = 0. _d 0
72            cg2d_r(I,J,bi,bj) = 0. _d 0            cg2d_r(I,J,bi,bj) = 0. _d 0
73            cg2d_s(I,J,bi,bj) = 0. _d 0            cg2d_s(I,J,bi,bj) = 0. _d 0
74    #ifdef ALLOW_CG2D_NSA
75              cg2d_z(I,J,bi,bj) = 0. _d 0
76    #endif /* ALLOW_CG2D_NSA */
77           ENDDO           ENDDO
78          ENDDO          ENDDO
79         ENDDO         ENDDO
# Line 151  CcnhDebugStarts Line 154  CcnhDebugStarts
154  C     CALL PLOT_FIELD_XYRS( aW2d, 'AW2D INI_CG2D.1' , 1, myThid )  C     CALL PLOT_FIELD_XYRS( aW2d, 'AW2D INI_CG2D.1' , 1, myThid )
155  C     CALL PLOT_FIELD_XYRS( aS2d, 'AS2D INI_CG2D.1' , 1, myThid )  C     CALL PLOT_FIELD_XYRS( aS2d, 'AS2D INI_CG2D.1' , 1, myThid )
156  CcnhDebugEnds  CcnhDebugEnds
157  c     _EXCH_XY_RS(aW2d, myThid)        CALL EXCH_UV_XY_RS( aW2d, aS2d, .FALSE., myThid )
 c     _EXCH_XY_RS(aS2d, myThid)  
       CALL EXCH_UV_XY_RS(aW2d,aS2d,.FALSE.,myThid)  
158  CcnhDebugStarts  CcnhDebugStarts
159  C     CALL PLOT_FIELD_XYRS( aW2d, 'AW2D INI_CG2D.2' , 1, myThid )  C     CALL PLOT_FIELD_XYRS( aW2d, 'AW2D INI_CG2D.2' , 1, myThid )
160  C     CALL PLOT_FIELD_XYRS( aS2d, 'AS2D INI_CG2D.2' , 1, myThid )  C     CALL PLOT_FIELD_XYRS( aS2d, 'AS2D INI_CG2D.2' , 1, myThid )
# Line 254  C         pS(I,J,bi,bj) = 0. Line 255  C         pS(I,J,bi,bj) = 0.
255         ENDDO         ENDDO
256        ENDDO        ENDDO
257  C--   Update overlap regions  C--   Update overlap regions
258        _EXCH_XY_RS(pC, myThid)        CALL EXCH_XY_RS( pC, myThid )
259  c     _EXCH_XY_RS(pW, myThid)        CALL EXCH_UV_XY_RS( pW, pS, .FALSE., myThid )
 c     _EXCH_XY_RS(pS, myThid)  
       CALL EXCH_UV_XY_RS(pW,pS,.FALSE.,myThid)  
260  CcnhDebugStarts  CcnhDebugStarts
261  C     CALL PLOT_FIELD_XYRS( pC, 'pC   INI_CG2D.2' , 1, myThid )  C     CALL PLOT_FIELD_XYRS( pC, 'pC   INI_CG2D.2' , 1, myThid )
262  C     CALL PLOT_FIELD_XYRS( pW, 'pW   INI_CG2D.2' , 1, myThid )  C     CALL PLOT_FIELD_XYRS( pW, 'pW   INI_CG2D.2' , 1, myThid )
Legend:
Removed from v.1.46  
changed lines
  Added in v.1.47
  ViewVC Help
Powered by ViewVC 1.1.22