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

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

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

revision 1.2 by adcroft, Fri Feb 2 21:36:30 2001 UTC revision 1.3 by adcroft, Mon Feb 5 15:45:47 2001 UTC
# Line 34  C     == Routine arguments == Line 34  C     == Routine arguments ==
34  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
35    
36  C     == Local variables ==  C     == Local variables ==
37        INTEGER I, J ,K        INTEGER I, J , K, I_obc, J_obc
38    
39  C     Eastern OB  C     Eastern OB
40        IF (useOrlanskiEast) THEN        IF (useOrlanskiEast) THEN
# Line 45  C     Eastern OB Line 45  C     Eastern OB
45        ELSE        ELSE
46          DO K=1,Nr          DO K=1,Nr
47            DO J=1-Oly,sNy+Oly            DO J=1-Oly,sNy+Oly
48              OBEu(J,K,bi,bj)=0.              I_obc=OB_Ie(J,bi,bj)
49              OBEv(J,K,bi,bj)=0.              IF (I_obc.ne.0) THEN
50              OBEt(J,K,bi,bj)=tRef(K)                OBEu(J,K,bi,bj)=0.
51              OBEs(J,K,bi,bj)=sRef(K)                OBEv(J,K,bi,bj)=0.
52                  OBEt(J,K,bi,bj)=tRef(K)
53                  OBEs(J,K,bi,bj)=sRef(K)
54  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
55              OBEw(J,K,bi,bj)=0.                OBEw(J,K,bi,bj)=0.
56  #endif  #endif
57                ENDIF
58            ENDDO            ENDDO
59          ENDDO          ENDDO
60        ENDIF        ENDIF
# Line 65  C     Western OB Line 68  C     Western OB
68        ELSE        ELSE
69          DO K=1,Nr          DO K=1,Nr
70            DO J=1-Oly,sNy+Oly            DO J=1-Oly,sNy+Oly
71              OBWu(J,K,bi,bj)=0.              I_obc=OB_Iw(J,bi,bj)
72              OBWv(J,K,bi,bj)=0.              IF (I_obc.ne.0) THEN
73              OBWt(J,K,bi,bj)=tRef(K)                OBWu(J,K,bi,bj)=0.
74              OBWs(J,K,bi,bj)=sRef(K)                OBWv(J,K,bi,bj)=0.
75                  OBWt(J,K,bi,bj)=tRef(K)
76                  OBWs(J,K,bi,bj)=sRef(K)
77  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
78              OBWw(J,K,bi,bj)=0.                OBWw(J,K,bi,bj)=0.
79  #endif  #endif
80                ENDIF
81            ENDDO            ENDDO
82          ENDDO          ENDDO
83        ENDIF        ENDIF
84    
85  C         Northern OB, template for forcing  C         Northern OB
86        IF (useOrlanskiNorth) THEN        IF (useOrlanskiNorth) THEN
87          CALL ORLANSKI_NORTH(          CALL ORLANSKI_NORTH(
88       &          bi, bj, futureTime,       &          bi, bj, futureTime,
# Line 85  C         Northern OB, template for forc Line 91  C         Northern OB, template for forc
91        ELSE        ELSE
92          DO K=1,Nr          DO K=1,Nr
93            DO I=1-Olx,sNx+Olx            DO I=1-Olx,sNx+Olx
94              OBNv(I,K,bi,bj)=0.              J_obc=OB_Jn(I,bi,bj)
95              OBNu(I,K,bi,bj)=0.              IF (J_obc.ne.0) THEN
96              OBNt(I,K,bi,bj)=tRef(K)                OBNv(I,K,bi,bj)=0.
97              OBNs(I,K,bi,bj)=sRef(K)                OBNu(I,K,bi,bj)=0.
98                  OBNt(I,K,bi,bj)=tRef(K)
99                  OBNs(I,K,bi,bj)=sRef(K)
100  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
101              OBNw(I,K,bi,bj)=0.                OBNw(I,K,bi,bj)=0.
102  #endif  #endif
103                ENDIF
104            ENDDO            ENDDO
105          ENDDO          ENDDO
106        ENDIF        ENDIF
107    
108  C         Southern OB, template for forcing  C         Southern OB
109    #ifdef ALLOW_ORLANSKI
110        IF (useOrlanskiSouth) THEN          IF (useOrlanskiSouth) THEN  
111          CALL ORLANSKI_SOUTH(          CALL ORLANSKI_SOUTH(
112       &          bi, bj, futureTime,       &          bi, bj, futureTime,
113       &          uVel, vVel, wVel, theta, salt,       &          uVel, vVel, wVel, theta, salt,
114       &          myThid )       &          myThid )
115        ELSE        ELSE
116    #endif
117          DO K=1,Nr          DO K=1,Nr
118            DO I=1-Olx,sNx+Olx            DO I=1-Olx,sNx+Olx
119              OBSu(I,K,bi,bj)=0.              J_obc=OB_Js(I,bi,bj)
120              OBSv(I,K,bi,bj)=0.              IF (J_obc.ne.0) THEN
121              OBSt(I,K,bi,bj)=tRef(K)                OBSu(I,K,bi,bj)=0.
122              OBSs(I,K,bi,bj)=sRef(K)                OBSv(I,K,bi,bj)=0.
123                  OBSt(I,K,bi,bj)=tRef(K)
124                  OBSs(I,K,bi,bj)=sRef(K)
125  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
126              OBSw(I,K,bi,bj)=0.                OBSw(I,K,bi,bj)=0.
127  #endif  #endif
128                ENDIF
129            ENDDO            ENDDO
130          ENDDO          ENDDO
131        ENDIF        ENDIF
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22