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

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

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

revision 1.1 by adcroft, Tue Jan 30 21:03:00 2001 UTC revision 1.2 by adcroft, Fri Feb 2 21:36:29 2001 UTC
# Line 0  Line 1 
1    C $Header$
2    C $Name$
3    
4    #include "OBCS_OPTIONS.h"
5    
6          SUBROUTINE OBCS_APPLY_W( bi, bj, K,
7         U                         wFld,
8         I                         myThid )
9    C     /==========================================================\
10    C     | S/R OBCS_APPLY_W                                         |
11    C     \==========================================================/
12          IMPLICIT NONE
13    C     == Global variables ==
14    #include "SIZE.h"
15    #include "EEPARAMS.h"
16    #include "PARAMS.h"
17    #include "GRID.h"
18    #include "OBCS.h"
19    
20    C     == Routine Arguments ==
21          INTEGER bi,bj,K
22          _RL wFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
23          INTEGER myThid
24    
25    #ifdef ALLOW_OBCS
26    
27    C     == Local variables ==
28          INTEGER I,J,I_obc,J_obc
29          _RL obc_mask
30    
31    #ifdef ALLOW_NONHYDROSTATIC
32    
33    C     Set model variables to OB values on North/South Boundaries
34          DO I=1-Olx,sNx+Olx
35    C Northern boundary
36           J_obc=OB_Jn(I,bi,bj)
37           IF (J_obc.NE.0) THEN
38            obc_mask = _maskS(I,J_obc,K,bi,bj)
39            wFld(I,J_obc,K,bi,bj)=OBNw(I,K,bi,bj) * obc_mask
40           ENDIF
41    C Southern boundary
42           J_obc=OB_Js(I,bi,bj)
43           IF (J_obc.NE.0) THEN
44            obc_mask = _maskS(I,J_obc+1,K,bi,bj)
45            wFld(I,J_obc,K,bi,bj)=OBSw(I,K,bi,bj) * obc_mask
46           ENDIF
47          ENDDO
48    
49    C     Set model variables to OB values on East/West Boundaries
50          DO J=1-Oly,sNy+Oly
51    C Eastern boundary
52           I_obc=OB_Ie(J,bi,bj)
53           IF (I_obc.NE.0) THEN
54            obc_mask = _maskW(I_obc,J,K,bi,bj)
55            wFld(I_obc,J,K,bi,bj)=OBEw(J,K,bi,bj) * obc_mask
56           ENDIF
57    C Western boundary
58           I_obc=OB_Iw(J,bi,bj)
59           IF (I_obc.NE.0) THEN
60            obc_mask = _maskW(I_obc+1,J,K,bi,bj)
61            wFld(I_obc,J,K,bi,bj)=OBWw(J,K,bi,bj) * obc_mask
62           ENDIF
63          ENDDO
64    
65    #endif
66    
67    #endif
68          RETURN
69          END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22