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

Diff of /MITgcm/pkg/obcs/obcs_apply_ts.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_TS( bi, bj, K,
7         U                          tFld, sFld,
8         I                          myThid )
9    C     /==========================================================\
10    C     | S/R OBCS_APPLY_TS                                        |
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 tFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
23          _RL sFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
24          INTEGER myThid
25    
26    #ifdef ALLOW_OBCS
27    
28    C     == Local variables ==
29          INTEGER I,J,I_obc,J_obc
30          _RL obc_mask
31    
32    C     Set model variables to OB values on North/South Boundaries
33          DO I=1-Olx,sNx+Olx
34    C Northern boundary
35           J_obc = OB_Jn(I,bi,bj)
36           IF (J_obc.NE.0) THEN
37            obc_mask = _maskS(I,J_obc,K,bi,bj)
38            tFld(I,J_obc,K,bi,bj)=OBNt(I,K,bi,bj)*obc_mask
39            sFld(I,J_obc,K,bi,bj)=OBNs(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            tFld(I,J_obc,K,bi,bj)=OBSt(I,K,bi,bj)*obc_mask
46            sFld(I,J_obc,K,bi,bj)=OBSs(I,K,bi,bj)*obc_mask
47           ENDIF
48          ENDDO
49    
50    C     Set model variables to OB values on East/West Boundaries
51          DO J=1-Oly,sNy+Oly
52    C Eastern boundary
53           I_obc = OB_Ie(J,bi,bj)
54           IF (I_obc.NE.0) THEN
55            obc_mask = _maskW(I_obc,J,K,bi,bj)
56            tFld(I_obc,J,K,bi,bj)=OBEt(J,K,bi,bj)*obc_mask
57            sFld(I_obc,J,K,bi,bj)=OBEs(J,K,bi,bj)*obc_mask
58           ENDIF
59    C Western boundary
60           I_obc=OB_Iw(J,bi,bj)
61           IF (I_obc.NE.0) THEN
62            obc_mask = _maskW(I_obc+1,J,K,bi,bj)
63            tFld(I_obc,J,K,bi,bj)=OBWt(J,K,bi,bj) * obc_mask
64            sFld(I_obc,J,K,bi,bj)=OBWs(J,K,bi,bj) * obc_mask
65           ENDIF
66          ENDDO
67    
68    #endif
69          RETURN
70          END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22