/[MITgcm]/MITgcm/pkg/generic_advdiff/gad_advection.F
ViewVC logotype

Diff of /MITgcm/pkg/generic_advdiff/gad_advection.F

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

revision 1.25 by heimbach, Wed Jun 30 23:45:35 2004 UTC revision 1.26 by cnh, Wed Jul 7 20:09:42 2004 UTC
# Line 306  C      ----------------|--||------- Line 306  C      ----------------|--||-------
306  C                      |  \/  C                      |  \/
307  C                      |   x T(sNx+1,    0)  C                      |   x T(sNx+1,    0)
308         IF ( southWestCorner ) THEN         IF ( southWestCorner ) THEN
309          localTij(0    ,0    )= localTij(0    ,1  )          DO J=1,OLy
310             DO I=1,OLx
311              localTij(1-I, 1-J   )= localTij(1-J  ,1  )
312             ENDDO
313            ENDDO
314         ENDIF         ENDIF
315         IF ( southEastCorner ) THEN         IF ( southEastCorner ) THEN
316          localTij(sNx+1,0    )= localTij(sNx+1,1  )          DO J=1,OLy
317             DO I=1,OLx
318              localTij(sNx+I, 1-J )=localTij(sNx+J, I  )
319             ENDDO
320            ENDDO
321         ENDIF         ENDIF
322         IF ( northWestCorner ) THEN         IF ( northWestCorner ) THEN
323          localTij(0    ,sNy+1)= localTij(0    ,sNy)          DO J=1,OLy
324             DO I=1,OLx
325              localTij( 1-I ,sNy+J)=localTij( 1-J , sNy+1-I )
326             ENDDO
327            ENDDO
328         ENDIF         ENDIF
329         IF ( northEastCorner ) THEN         IF ( northEastCorner ) THEN
330          localTij(sNx+1,sNy+1)= localTij(sNx+1,sNy)          DO J=1,OLy
331             DO I=1,OLx
332              localTij(sNx+I,sNy+J)=localTij(sNx+J, sNy+1-I )
333             ENDDO
334            ENDDO
335         ENDIF         ENDIF
336        ENDIF        ENDIF
337    
# Line 413  C      ----------------|----------- Line 429  C      ----------------|-----------
429  C                      |      C                      |    
430  C         x T(sNx,0) =====>x T(sNx+1,    0)  C         x T(sNx,0) =====>x T(sNx+1,    0)
431         IF ( southWestCorner ) THEN         IF ( southWestCorner ) THEN
432           localTij(    0,0    ) = localTij(  1,0    )          DO J=1,Oly
433             DO I=1,Olx
434              localTij( 1-i , 1-j ) = localTij(j   , 1-i )
435             ENDDO
436            ENDDO
437         ENDIF         ENDIF
438         IF ( southEastCorner ) THEN         IF ( southEastCorner ) THEN
439           localTij(sNx+1,0    ) = localTij(sNx,0    )          DO J=1,Oly
440             DO I=1,Olx
441              localTij(sNx+i, 1-j ) = localTij(sNx+1-j, 1-i )
442             ENDDO
443            ENDDO
444         ENDIF         ENDIF
445         IF ( northWestCorner ) THEN         IF ( northWestCorner ) THEN
446           localTij(0    ,sNy+1) = localTij(  1,sNy+1)          DO J=1,Oly
447             DO I=1,Olx
448              localTij( 1-i ,sNy+j) = localTij(j   ,sNy+i)
449             ENDDO
450            ENDDO
451         ENDIF         ENDIF
452         IF ( northEastCorner ) THEN         IF ( northEastCorner ) THEN
453           localTij(sNx+1,sNy+1) = localTij(sNx,sNy+1)          DO J=1,Oly
454             DO I=1,Olx
455              localTij(sNx+i,sNy+j) = localTij(sNx+1-j,sNy+i)
456             ENDDO
457            ENDDO
458         ENDIF         ENDIF
459        ENDIF        ENDIF
460    
Legend:
Removed from v.1.25  
changed lines
  Added in v.1.26
  ViewVC Help
Powered by ViewVC 1.1.22