19 |
#include "GRID.h" |
#include "GRID.h" |
20 |
#include "DYNVARS.h" |
#include "DYNVARS.h" |
21 |
#include "CG2D.h" |
#include "CG2D.h" |
22 |
|
#include "OBCS.h" |
23 |
|
|
24 |
C === Routine arguments === |
C === Routine arguments === |
25 |
C myThid - Thread no. that called this routine. |
C myThid - Thread no. that called this routine. |
67 |
ENDDO |
ENDDO |
68 |
ENDDO |
ENDDO |
69 |
ENDDO |
ENDDO |
70 |
|
IF (openBoundaries) THEN |
71 |
|
DO I=1,sNx |
72 |
|
IF (OB_Jn(I,bi,bj).NE.0) aS2d(I,OB_Jn(I,bi,bj),bi,bj)=0. |
73 |
|
IF (OB_Jn(I,bi,bj).NE.0) aS2d(I,OB_Jn(I,bi,bj)+1,bi,bj)=0. |
74 |
|
IF (OB_Js(I,bi,bj).NE.0) aS2d(I,OB_Js(I,bi,bj)+1,bi,bj)=0. |
75 |
|
IF (OB_Js(I,bi,bj).NE.0) aS2d(I,OB_Js(I,bi,bj),bi,bj)=0. |
76 |
|
ENDDO |
77 |
|
DO J=1,sNy |
78 |
|
IF (OB_Ie(J,bi,bj).NE.0) aW2d(OB_Ie(J,bi,bj),J,bi,bj)=0. |
79 |
|
IF (OB_Ie(J,bi,bj).NE.0) aW2d(OB_Ie(J,bi,bj)+1,J,bi,bj)=0. |
80 |
|
IF (OB_Iw(J,bi,bj).NE.0) aW2d(OB_Iw(J,bi,bj)+1,J,bi,bj)=0. |
81 |
|
IF (OB_Iw(J,bi,bj).NE.0) aW2d(OB_Iw(J,bi,bj),J,bi,bj)=0. |
82 |
|
ENDDO |
83 |
|
ENDIF |
84 |
DO J=1,sNy |
DO J=1,sNy |
85 |
DO I=1,sNx |
DO I=1,sNx |
86 |
myNorm = MAX(ABS(aW2d(I,J,bi,bj)),myNorm) |
myNorm = MAX(ABS(aW2d(I,J,bi,bj)),myNorm) |