/[MITgcm]/MITgcm/pkg/obcs/obcs_apply_uv.F
ViewVC logotype

Diff of /MITgcm/pkg/obcs/obcs_apply_uv.F

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

revision 1.7 by jmc, Tue May 24 14:31:14 2011 UTC revision 1.8 by jmc, Tue Sep 18 20:09:17 2012 UTC
# Line 83  C--   Set Tangential component first: Line 83  C--   Set Tangential component first:
83  #ifdef ALLOW_OBCS_NORTH  #ifdef ALLOW_OBCS_NORTH
84          IF ( tileHasOBN(bi,bj) ) THEN          IF ( tileHasOBN(bi,bj) ) THEN
85  C Northern boundary  C Northern boundary
86           DO i=1-Olx,sNx+Olx           DO i=1-OLx,sNx+OLx
87            Jobc = OB_Jn(i,bi,bj)            Jobc = OB_Jn(i,bi,bj)
88            IF ( Jobc.NE.0 ) THEN            IF ( Jobc.NE.OB_indexNone ) THEN
89             DO k = kLo,kHi             DO k = kLo,kHi
90               uFld(i,Jobc,k,bi,bj) = OBNu(i,k,bi,bj)               uFld(i,Jobc,k,bi,bj) = OBNu(i,k,bi,bj)
91       &                            *_maskW(i,Jobc,k,bi,bj)       &                            *_maskW(i,Jobc,k,bi,bj)
# Line 97  C Northern boundary Line 97  C Northern boundary
97  #ifdef ALLOW_OBCS_SOUTH  #ifdef ALLOW_OBCS_SOUTH
98          IF ( tileHasOBS(bi,bj) ) THEN          IF ( tileHasOBS(bi,bj) ) THEN
99  C Southern boundary  C Southern boundary
100           DO i=1-Olx,sNx+Olx           DO i=1-OLx,sNx+OLx
101            Jobc = OB_Js(i,bi,bj)            Jobc = OB_Js(i,bi,bj)
102            IF ( Jobc.NE.0 ) THEN            IF ( Jobc.NE.OB_indexNone ) THEN
103             DO k = kLo,kHi             DO k = kLo,kHi
104               uFld(i,Jobc,k,bi,bj) = OBSu(i,k,bi,bj)               uFld(i,Jobc,k,bi,bj) = OBSu(i,k,bi,bj)
105       &                            *_maskW(i,Jobc,k,bi,bj)       &                            *_maskW(i,Jobc,k,bi,bj)
# Line 113  C     Set model variables to OB values o Line 113  C     Set model variables to OB values o
113  #ifdef ALLOW_OBCS_EAST  #ifdef ALLOW_OBCS_EAST
114          IF ( tileHasOBE(bi,bj) ) THEN          IF ( tileHasOBE(bi,bj) ) THEN
115  C Eastern boundary  C Eastern boundary
116           DO j=1-Oly,sNy+Oly           DO j=1-OLy,sNy+OLy
117            Iobc = OB_Ie(j,bi,bj)            Iobc = OB_Ie(j,bi,bj)
118            IF ( Iobc.NE.0 ) THEN            IF ( Iobc.NE.OB_indexNone ) THEN
119             DO k = kLo,kHi             DO k = kLo,kHi
120               vFld(Iobc,j,k,bi,bj) = OBEv(j,k,bi,bj)               vFld(Iobc,j,k,bi,bj) = OBEv(j,k,bi,bj)
121       &                            *_maskS(Iobc,j,k,bi,bj)       &                            *_maskS(Iobc,j,k,bi,bj)
# Line 127  C Eastern boundary Line 127  C Eastern boundary
127  #ifdef ALLOW_OBCS_WEST  #ifdef ALLOW_OBCS_WEST
128          IF ( tileHasOBW(bi,bj) ) THEN          IF ( tileHasOBW(bi,bj) ) THEN
129  C Western boundary  C Western boundary
130           DO j=1-Oly,sNy+Oly           DO j=1-OLy,sNy+OLy
131            Iobc = OB_Iw(j,bi,bj)            Iobc = OB_Iw(j,bi,bj)
132            IF ( Iobc.NE.0 ) THEN            IF ( Iobc.NE.OB_indexNone ) THEN
133             DO k = kLo,kHi             DO k = kLo,kHi
134               vFld(Iobc,j,k,bi,bj) = OBWv(j,k,bi,bj)               vFld(Iobc,j,k,bi,bj) = OBWv(j,k,bi,bj)
135       &                            *_maskS(Iobc,j,k,bi,bj)       &                            *_maskS(Iobc,j,k,bi,bj)
# Line 144  C--   Then set Normal component: Line 144  C--   Then set Normal component:
144  #ifdef ALLOW_OBCS_NORTH  #ifdef ALLOW_OBCS_NORTH
145          IF ( tileHasOBN(bi,bj) ) THEN          IF ( tileHasOBN(bi,bj) ) THEN
146  C Northern boundary  C Northern boundary
147           DO i=1-Olx,sNx+Olx           DO i=1-OLx,sNx+OLx
148            Jobc = OB_Jn(i,bi,bj)            Jobc = OB_Jn(i,bi,bj)
149            IF ( Jobc.NE.0 ) THEN            IF ( Jobc.NE.OB_indexNone ) THEN
150             DO k = kLo,kHi             DO k = kLo,kHi
151               vFld(i,Jobc,k,bi,bj)   = OBNv(i,k,bi,bj)               vFld(i,Jobc,k,bi,bj)   = OBNv(i,k,bi,bj)
152       &                              *_maskS(i,Jobc,k,bi,bj)       &                              *_maskS(i,Jobc,k,bi,bj)
# Line 161  C Northern boundary Line 161  C Northern boundary
161  #ifdef ALLOW_OBCS_SOUTH  #ifdef ALLOW_OBCS_SOUTH
162          IF ( tileHasOBS(bi,bj) ) THEN          IF ( tileHasOBS(bi,bj) ) THEN
163  C Southern boundary  C Southern boundary
164           DO i=1-Olx,sNx+Olx           DO i=1-OLx,sNx+OLx
165            Jobc = OB_Js(i,bi,bj)            Jobc = OB_Js(i,bi,bj)
166            IF ( Jobc.NE.0 ) THEN            IF ( Jobc.NE.OB_indexNone ) THEN
167             DO k = kLo,kHi             DO k = kLo,kHi
168               vFld(i,Jobc+1,k,bi,bj) = OBSv(i,k,bi,bj)               vFld(i,Jobc+1,k,bi,bj) = OBSv(i,k,bi,bj)
169       &                              *_maskS(i,Jobc+1,k,bi,bj)       &                              *_maskS(i,Jobc+1,k,bi,bj)
# Line 180  C     Set model variables to OB values o Line 180  C     Set model variables to OB values o
180  #ifdef ALLOW_OBCS_EAST  #ifdef ALLOW_OBCS_EAST
181          IF ( tileHasOBE(bi,bj) ) THEN          IF ( tileHasOBE(bi,bj) ) THEN
182  C Eastern boundary  C Eastern boundary
183           DO j=1-Oly,sNy+Oly           DO j=1-OLy,sNy+OLy
184            Iobc = OB_Ie(j,bi,bj)            Iobc = OB_Ie(j,bi,bj)
185            IF ( Iobc.NE.0 ) THEN            IF ( Iobc.NE.OB_indexNone ) THEN
186             DO k = kLo,kHi             DO k = kLo,kHi
187               uFld(Iobc,j,k,bi,bj)   = OBEu(j,k,bi,bj)               uFld(Iobc,j,k,bi,bj)   = OBEu(j,k,bi,bj)
188       &                              *_maskW(Iobc,j,k,bi,bj)       &                              *_maskW(Iobc,j,k,bi,bj)
# Line 197  C Eastern boundary Line 197  C Eastern boundary
197  #ifdef ALLOW_OBCS_WEST  #ifdef ALLOW_OBCS_WEST
198          IF ( tileHasOBW(bi,bj) ) THEN          IF ( tileHasOBW(bi,bj) ) THEN
199  C Western boundary  C Western boundary
200           DO j=1-Oly,sNy+Oly           DO j=1-OLy,sNy+OLy
201            Iobc = OB_Iw(j,bi,bj)            Iobc = OB_Iw(j,bi,bj)
202            IF ( Iobc.NE.0 ) THEN            IF ( Iobc.NE.OB_indexNone ) THEN
203             DO k = kLo,kHi             DO k = kLo,kHi
204               uFld(Iobc+1,j,k,bi,bj) = OBWu(j,k,bi,bj)               uFld(Iobc+1,j,k,bi,bj) = OBWu(j,k,bi,bj)
205       &                              *_maskW(Iobc+1,j,k,bi,bj)       &                              *_maskW(Iobc+1,j,k,bi,bj)
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22