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

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

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

revision 1.48 by jmc, Tue Sep 4 14:54:57 2007 UTC revision 1.49 by mlosch, Tue Dec 4 10:22:15 2007 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    #ifdef TARGET_NEC_SX
6    C     set a sensible default for the outer loop unrolling parameter that can
7    C     be overriden in the Makefile with the DEFINES macro or in CPP_OPTIONS.h
8    #ifndef CG2D_OUTERLOOPITERS
9    # define CG2D_OUTERLOOPITERS 10
10    #endif
11    #endif /* TARGET_NEC_SX */
12    
13  CBOP  CBOP
14  C     !ROUTINE: CG2D  C     !ROUTINE: CG2D
# Line 164  C--   Initial residual calculation Line 171  C--   Initial residual calculation
171         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
172          sumRHStile(bi,bj) = 0. _d 0          sumRHStile(bi,bj) = 0. _d 0
173          errTile(bi,bj)    = 0. _d 0          errTile(bi,bj)    = 0. _d 0
174    #ifdef TARGET_NEC_SX
175    !CDIR OUTERUNROLL=CG2D_OUTERLOOPITERS
176    #endif /* TARGET_NEC_SX */
177          DO J=1,sNy          DO J=1,sNy
178           DO I=1,sNx           DO I=1,sNx
179  c         ks = ksurfC(I,J,bi,bj)  c         ks = ksurfC(I,J,bi,bj)
# Line 239  C--    conjugate direction vector "s". Line 249  C--    conjugate direction vector "s".
249         DO bj=myByLo(myThid),myByHi(myThid)         DO bj=myByLo(myThid),myByHi(myThid)
250          DO bi=myBxLo(myThid),myBxHi(myThid)          DO bi=myBxLo(myThid),myBxHi(myThid)
251           eta_qrNtile(bi,bj) = 0. _d 0           eta_qrNtile(bi,bj) = 0. _d 0
252    #ifdef TARGET_NEC_SX
253    !CDIR OUTERUNROLL=CG2D_OUTERLOOPITERS
254    #endif /* TARGET_NEC_SX */
255           DO J=1,sNy           DO J=1,sNy
256            DO I=1,sNx            DO I=1,sNx
257             cg2d_q(I,J,bi,bj) =             cg2d_q(I,J,bi,bj) =
# Line 295  C==    q = A.s Line 308  C==    q = A.s
308         DO bj=myByLo(myThid),myByHi(myThid)         DO bj=myByLo(myThid),myByHi(myThid)
309          DO bi=myBxLo(myThid),myBxHi(myThid)          DO bi=myBxLo(myThid),myBxHi(myThid)
310           alphaTile(bi,bj) = 0. _d 0           alphaTile(bi,bj) = 0. _d 0
311    #ifdef TARGET_NEC_SX
312    !CDIR OUTERUNROLL=CG2D_OUTERLOOPITERS
313    #endif /* TARGET_NEC_SX */
314           DO J=1,sNy           DO J=1,sNy
315            DO I=1,sNx            DO I=1,sNx
316  c          ks = ksurfC(I,J,bi,bj)  c          ks = ksurfC(I,J,bi,bj)
Legend:
Removed from v.1.48  
changed lines
  Added in v.1.49
  ViewVC Help
Powered by ViewVC 1.1.22