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

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

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

revision 1.9 by adcroft, Thu Jul 8 16:54:43 2004 UTC revision 1.10 by heimbach, Mon Sep 20 23:22:57 2004 UTC
# Line 36  C     == Local variables == Line 36  C     == Local variables ==
36    
37          DO K=1,Nr          DO K=1,Nr
38           DO I=1-Olx,sNx+Olx           DO I=1-Olx,sNx+Olx
39            OBNu(I,K,bi,bj)=0.  #ifdef ALLOW_OBCS_NORTH
40            OBNv(I,K,bi,bj)=0.            OBNu(I,K,bi,bj)=0. _d 0
41            OBNt(I,K,bi,bj)=0.            OBNv(I,K,bi,bj)=0. _d 0
42            OBNs(I,K,bi,bj)=0.            OBNt(I,K,bi,bj)=0. _d 0
43            OBSu(I,K,bi,bj)=0.            OBNs(I,K,bi,bj)=0. _d 0
44            OBSv(I,K,bi,bj)=0.  # ifdef ALLOW_NONHYDROSTATIC
45            OBSt(I,K,bi,bj)=0.            OBNw(I,K,bi,bj)=0. _d 0
46            OBSs(I,K,bi,bj)=0.  # endif
47  #ifdef ALLOW_NONHYDROSTATIC  # ifdef ALLOW_OBCS_PRESCRIBE
48            OBNw(I,K,bi,bj)=0.            OBNu0(I,K,bi,bj)=0. _d 0
49            OBSw(I,K,bi,bj)=0.            OBNv0(I,K,bi,bj)=0. _d 0
50  #endif            OBNt0(I,K,bi,bj)=0. _d 0
51              OBNs0(I,K,bi,bj)=0. _d 0
52              OBNu1(I,K,bi,bj)=0. _d 0
53              OBNv1(I,K,bi,bj)=0. _d 0
54              OBNt1(I,K,bi,bj)=0. _d 0
55              OBNs1(I,K,bi,bj)=0. _d 0
56    # endif
57    #endif /* ALLOW_OBCS_NORTH */
58    
59    #ifdef ALLOW_OBCS_SOUTH
60              OBSu(I,K,bi,bj)=0. _d 0
61              OBSv(I,K,bi,bj)=0. _d 0
62              OBSt(I,K,bi,bj)=0. _d 0
63              OBSs(I,K,bi,bj)=0. _d 0
64    # ifdef ALLOW_NONHYDROSTATIC
65              OBSw(I,K,bi,bj)=0. _d 0
66    # endif
67    # ifdef ALLOW_OBCS_PRESCRIBE
68              OBSu0(I,K,bi,bj)=0. _d 0
69              OBSv0(I,K,bi,bj)=0. _d 0
70              OBSt0(I,K,bi,bj)=0. _d 0
71              OBSs0(I,K,bi,bj)=0. _d 0
72              OBSu1(I,K,bi,bj)=0. _d 0
73              OBSv1(I,K,bi,bj)=0. _d 0
74              OBSt1(I,K,bi,bj)=0. _d 0
75              OBSs1(I,K,bi,bj)=0. _d 0
76    # endif
77    #endif /* ALLOW_OBCS_SOUTH */
78           ENDDO           ENDDO
79    
80           DO J=1-Oly,sNy+Oly           DO J=1-Oly,sNy+Oly
81            OBEu(J,K,bi,bj)=0.  #ifdef ALLOW_OBCS_EAST
82            OBEv(J,K,bi,bj)=0.            OBEu(J,K,bi,bj)=0. _d 0
83            OBEt(J,K,bi,bj)=0.            OBEv(J,K,bi,bj)=0. _d 0
84            OBEs(J,K,bi,bj)=0.            OBEt(J,K,bi,bj)=0. _d 0
85            OBWu(J,K,bi,bj)=0.            OBEs(J,K,bi,bj)=0. _d 0
86            OBWv(J,K,bi,bj)=0.  # ifdef ALLOW_NONHYDROSTATIC
87            OBWt(J,K,bi,bj)=0.            OBEw(J,K,bi,bj)=0. _d 0
88            OBWs(J,K,bi,bj)=0.  # endif
89  #ifdef ALLOW_NONHYDROSTATIC  # ifdef ALLOW_OBCS_PRESCRIBE
90            OBEw(J,K,bi,bj)=0.            OBEu0(J,K,bi,bj)=0. _d 0
91            OBWw(J,K,bi,bj)=0.            OBEv0(J,K,bi,bj)=0. _d 0
92  #endif            OBEt0(J,K,bi,bj)=0. _d 0
93              OBEs0(J,K,bi,bj)=0. _d 0
94              OBEu1(J,K,bi,bj)=0. _d 0
95              OBEv1(J,K,bi,bj)=0. _d 0
96              OBEt1(J,K,bi,bj)=0. _d 0
97              OBEs1(J,K,bi,bj)=0. _d 0
98    # endif
99    #endif /* ALLOW_OBCS_EAST */
100    
101    #ifdef ALLOW_OBCS_WEST
102              OBWu(J,K,bi,bj)=0. _d 0
103              OBWv(J,K,bi,bj)=0. _d 0
104              OBWt(J,K,bi,bj)=0. _d 0
105              OBWs(J,K,bi,bj)=0. _d 0
106    # ifdef ALLOW_NONHYDROSTATIC
107              OBWw(J,K,bi,bj)=0. _d 0
108    # endif
109    # ifdef ALLOW_OBCS_PRESCRIBE
110              OBWu0(J,K,bi,bj)=0. _d 0
111              OBWv0(J,K,bi,bj)=0. _d 0
112              OBWt0(J,K,bi,bj)=0. _d 0
113              OBWs0(J,K,bi,bj)=0. _d 0
114              OBWu1(J,K,bi,bj)=0. _d 0
115              OBWv1(J,K,bi,bj)=0. _d 0
116              OBWt1(J,K,bi,bj)=0. _d 0
117              OBWs1(J,K,bi,bj)=0. _d 0
118    # endif
119    #endif /* ALLOW_OBCS_WEST */
120           ENDDO           ENDDO
121          ENDDO          ENDDO
122    
# Line 106  C--   Apply OBCS values to initial condi Line 161  C--   Apply OBCS values to initial condi
161         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
162          CALL OBCS_CALC( bi, bj, startTime, nIter0,          CALL OBCS_CALC( bi, bj, startTime, nIter0,
163       &        uVel, vVel, wVel, theta, salt, myThid )       &        uVel, vVel, wVel, theta, salt, myThid )
164           ENDDO
165          ENDDO
166    
167          DO bj = myByLo(myThid), myByHi(myThid)
168           DO bi = myBxLo(myThid), myBxHi(myThid)
169          DO K=1,Nr          DO K=1,Nr
170           CALL OBCS_APPLY_UV( bi, bj, k, uVel, vVel, myThid )           CALL OBCS_APPLY_UV( bi, bj, k, uVel, vVel, myThid )
171           CALL OBCS_APPLY_TS( bi, bj, k, theta, salt, myThid )           CALL OBCS_APPLY_TS( bi, bj, k, theta, salt, myThid )
# Line 114  C--   Apply OBCS values to initial condi Line 174  C--   Apply OBCS values to initial condi
174        ENDDO        ENDDO
175    
176  #endif /* ALLOW_OBCS */  #endif /* ALLOW_OBCS */
177    
178        RETURN        RETURN
179        END        END
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22