20 |
#include "GRID.h" |
#include "GRID.h" |
21 |
#include "DYNVARS.h" |
#include "DYNVARS.h" |
22 |
#include "CG2D.h" |
#include "CG2D.h" |
23 |
|
#ifdef ALLOW_OBCS |
24 |
#include "OBCS.h" |
#include "OBCS.h" |
25 |
|
#endif |
26 |
|
|
27 |
C === Routine arguments === |
C === Routine arguments === |
28 |
C myThid - Thread no. that called this routine. |
C myThid - Thread no. that called this routine. |
70 |
ENDDO |
ENDDO |
71 |
ENDDO |
ENDDO |
72 |
ENDDO |
ENDDO |
73 |
|
#ifdef ALLOW_OBCS |
74 |
IF (openBoundaries) THEN |
IF (openBoundaries) THEN |
75 |
DO I=1,sNx |
DO I=1,sNx |
76 |
IF (OB_Jn(I,bi,bj).NE.0) aS2d(I,OB_Jn(I,bi,bj),bi,bj)=0. |
IF (OB_Jn(I,bi,bj).NE.0) aS2d(I,OB_Jn(I,bi,bj),bi,bj)=0. |
85 |
IF (OB_Iw(J,bi,bj).NE.0) aW2d(OB_Iw(J,bi,bj),J,bi,bj)=0. |
IF (OB_Iw(J,bi,bj).NE.0) aW2d(OB_Iw(J,bi,bj),J,bi,bj)=0. |
86 |
ENDDO |
ENDDO |
87 |
ENDIF |
ENDIF |
88 |
|
#endif |
89 |
DO J=1,sNy |
DO J=1,sNy |
90 |
DO I=1,sNx |
DO I=1,sNx |
91 |
myNorm = MAX(ABS(aW2d(I,J,bi,bj)),myNorm) |
myNorm = MAX(ABS(aW2d(I,J,bi,bj)),myNorm) |