C $Header: /u/gcmpack/MITgcm/eesupp/src/fill_cs_corner_ag_rl.F,v 1.3 2009/06/28 01:02:17 jmc Exp $ C $Name: $ #include "PACKAGES_CONFIG.h" #include "CPP_EEOPTIONS.h" C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| CBOP C !ROUTINE: FILL_CS_CORNER_AG_RL C !INTERFACE: SUBROUTINE FILL_CS_CORNER_AG_RL( I fill4dirX, withSigns, U uFld, vFld, I bi,bj, myThid) C !DESCRIPTION: C *==========================================================* C | SUBROUTINE FILL_CS_CORNER_AG_RL C | o Fill the corner-halo region of CS-grid, C | for a 2 components vector field on A-grid (center of grid cell) C *==========================================================* C | o the corner halo region is filled with valid values C | in order to compute (later on) gradient in X or Y C | direction, on a wide stencil. C *==========================================================* C !USES: IMPLICIT NONE C == Global variables == #include "SIZE.h" #include "EEPARAMS.h" #ifdef ALLOW_EXCH2 #include "W2_EXCH2_SIZE.h" #include "W2_EXCH2_TOPOLOGY.h" #endif /* ALLOW_EXCH2 */ C !INPUT/OUTPUT PARAMETERS: C == Routine arguments == C C fill4dirX :: True = prepare for X direction calculations C otherwise, prepare for Y direction C withSigns :: True = account for sign of vector component C uFld :: u-component of vector field array with empty corners to fill C vFld :: v-component of vector field array with empty corners to fill C bi,bj :: tile indices C myThid :: thread number LOGICAL fill4dirX LOGICAL withSigns _RL uFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy) _RL vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy) INTEGER bi,bj INTEGER myThid C !LOCAL VARIABLES: C == Local variables == C i,j :: loop indices C myTile :: tile number INTEGER i,j LOGICAL southWestCorner LOGICAL southEastCorner LOGICAL northWestCorner LOGICAL northEastCorner _RL negOne #ifdef ALLOW_EXCH2 INTEGER myTile #endif CEOP negOne = 1. IF (withSigns) negOne = -1. IF (useCubedSphereExchange) THEN #ifdef ALLOW_EXCH2 myTile = W2_myTileList(bi,bj) southWestCorner = exch2_isWedge(myTile).EQ.1 & .AND. exch2_isSedge(myTile).EQ.1 southEastCorner = exch2_isEedge(myTile).EQ.1 & .AND. exch2_isSedge(myTile).EQ.1 northEastCorner = exch2_isEedge(myTile).EQ.1 & .AND. exch2_isNedge(myTile).EQ.1 northWestCorner = exch2_isWedge(myTile).EQ.1 & .AND. exch2_isNedge(myTile).EQ.1 #else southWestCorner = .TRUE. southEastCorner = .TRUE. northWestCorner = .TRUE. northEastCorner = .TRUE. #endif IF ( fill4dirX ) THEN C-- Internal exchange for calculations in X C- For cube face corners we need to duplicate the C- i-1 and i+1 values into the null space as follows: C C C o NW corner: copy T( 0,sNy ) into T( 0,sNy+1) e.g. C | C x T(0,sNy+1) | C /\ | C --||------------|----------- C || | C x T(0,sNy) | x T(1,sNy) C | C C o SW corner: copy T(0,1) into T(0,0) e.g. C | C x T(0,1) | x T(1,1) C || | C --||------------|----------- C \/ | C x T(0,0) | C | C C o NE corner: copy T(sNx+1,sNy ) into T(sNx+1,sNy+1) e.g. C | C | x T(sNx+1,sNy+1) C | /\ C ----------------|--||------- C | || C x T(sNx,sNy) | x T(sNx+1,sNy ) C | C o SE corner: copy T(sNx+1,1 ) into T(sNx+1,0 ) e.g. C | C x T(sNx,1) | x T(sNx+1, 1) C | || C ----------------|--||------- C | \/ C | x T(sNx+1, 0) IF ( southWestCorner ) THEN DO j=1,OLy DO i=1,OLx uFld( 1-i , 1-j ) = vFld( 1-j , i )*negOne vFld( 1-i , 1-j ) = uFld( 1-j , i ) ENDDO ENDDO ENDIF IF ( southEastCorner ) THEN DO j=1,OLy DO i=1,OLx uFld(sNx+i, 1-j ) = vFld(sNx+j, i ) vFld(sNx+i, 1-j ) = uFld(sNx+j, i )*negOne ENDDO ENDDO ENDIF IF ( northWestCorner ) THEN DO j=1,OLy DO i=1,OLx uFld( 1-i ,sNy+j) = vFld( 1-j , sNy+1-i ) vFld( 1-i ,sNy+j) = uFld( 1-j , sNy+1-i )*negOne ENDDO ENDDO ENDIF IF ( northEastCorner ) THEN DO j=1,OLy DO i=1,OLx uFld(sNx+i,sNy+j) = vFld(sNx+j, sNy+1-i )*negOne vFld(sNx+i,sNy+j) = uFld(sNx+j, sNy+1-i ) ENDDO ENDDO ENDIF C-- End of X direction ; start Y direction case. ELSE C-- Internal exchange for calculations in Y C- For cube face corners we need to duplicate the C- j-1 and j+1 values into the null space as follows: C C o SW corner: copy T(0,1) into T(0,0) e.g. C | C | x T(1,1) C | C ----------------|----------- C | C x T(0,0)<====== x T(1,0) C | C C o NW corner: copy T( 0,sNy ) into T( 0,sNy+1) e.g. C | C x T(0,sNy+1)<=== x T(1,sNy+1) C | C ----------------|----------- C | C | x T(1,sNy) C | C C o NE corner: copy T(sNx+1,sNy ) into T(sNx+1,sNy+1) e.g. C | C x T(sNx,sNy+1)=====>x T(sNx+1,sNy+1) C | C ----------------|----------- C | C x T(sNx,sNy) | C | C o SE corner: copy T(sNx+1,1 ) into T(sNx+1,0 ) e.g. C | C x T(sNx,1) | C | C ----------------|----------- C | C x T(sNx,0) =====>x T(sNx+1, 0) IF ( southWestCorner ) THEN DO j=1,Oly DO i=1,Olx uFld( 1-i , 1-j ) = vFld( j , 1-i ) vFld( 1-i , 1-j ) = uFld( j , 1-i )*negOne ENDDO ENDDO ENDIF IF ( southEastCorner ) THEN DO j=1,Oly DO i=1,Olx uFld(sNx+i, 1-j ) = vFld(sNx+1-j, 1-i )*negOne vFld(sNx+i, 1-j ) = uFld(sNx+1-j, 1-i ) ENDDO ENDDO ENDIF IF ( northWestCorner ) THEN DO j=1,Oly DO i=1,Olx uFld( 1-i ,sNy+j) = vFld( j ,sNy+i)*negOne vFld( 1-i ,sNy+j) = uFld( j ,sNy+i) ENDDO ENDDO ENDIF IF ( northEastCorner ) THEN DO j=1,Oly DO i=1,Olx uFld(sNx+i,sNy+j) = vFld(sNx+1-j,sNy+i) vFld(sNx+i,sNy+j) = uFld(sNx+1-j,sNy+i)*negOne ENDDO ENDDO ENDIF C- End of Y direction case. ENDIF C-- End useCubedSphereExchange ENDIF RETURN END