/[MITgcm]/MITgcm/pkg/ecco/cost_obcsvol.F
ViewVC logotype

Diff of /MITgcm/pkg/ecco/cost_obcsvol.F

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

revision 1.10 by jmc, Fri Aug 10 19:45:26 2012 UTC revision 1.11 by jmc, Tue Sep 18 20:16:34 2012 UTC
# Line 29  c     == global variables == Line 29  c     == global variables ==
29  #include "SIZE.h"  #include "SIZE.h"
30  #include "PARAMS.h"  #include "PARAMS.h"
31  #include "GRID.h"  #include "GRID.h"
 CML#include "DYNVARS.h"  
32  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
 c# include "OBCS_PARAMS.h"  
33  # include "OBCS_GRID.h"  # include "OBCS_GRID.h"
34  #endif  #endif
35    
# Line 206  c--         Determine the weights to be Line 204  c--         Determine the weights to be
204                do k = 1, Nr                do k = 1, Nr
205                  do i = imin,imax                  do i = imin,imax
206                    j = OB_Jn(I,bi,bj)                    j = OB_Jn(I,bi,bj)
207                      IF ( j.EQ.OB_indexNone ) j = 1
208  cgg    Barotropic velocity is stored in level 1.  cgg    Barotropic velocity is stored in level 1.
209                    tmpx = tmpfldxz(i,1,bi,bj)                    tmpx = tmpfldxz(i,1,bi,bj)
210                    if (maskS(i,j,k,bi,bj) .ne. 0.) then                    if (maskS(i,j,k,bi,bj) .ne. 0.) then
# Line 247  c--         Determine the weights to be Line 246  c--         Determine the weights to be
246                do k = 1, Nr                do k = 1, Nr
247                  do i = imin,imax                  do i = imin,imax
248                    j = OB_Js(I,bi,bj)                    j = OB_Js(I,bi,bj)
249                      IF ( j.EQ.OB_indexNone ) j = 1
250  cgg    Barotropic velocity is stored in level 1.  cgg    Barotropic velocity is stored in level 1.
251                    tmpx = tmpfldxz(i,1,bi,bj)                    tmpx = tmpfldxz(i,1,bi,bj)
252                    if (maskS(i,j+1,k,bi,bj) .ne. 0.) then                    if (maskS(i,j+1,k,bi,bj) .ne. 0.) then
# Line 289  c--         Determine the weights to be Line 289  c--         Determine the weights to be
289                do k = 1, Nr                do k = 1, Nr
290                  do j = jmin,jmax                  do j = jmin,jmax
291                    i = OB_Iw(j,bi,bj)                    i = OB_Iw(j,bi,bj)
292                      IF ( i.EQ.OB_indexNone ) i = 1
293  cgg    Barotropic velocity is stored in the level 1.  cgg    Barotropic velocity is stored in the level 1.
294                    tmpy = tmpfldyz(j,1,bi,bj)                    tmpy = tmpfldyz(j,1,bi,bj)
295                    if (maskW(i+1,j,k,bi,bj) .ne. 0.) then                    if (maskW(i+1,j,k,bi,bj) .ne. 0.) then
# Line 331  c--         Determine the weights to be Line 332  c--         Determine the weights to be
332                do k = 1, Nr                do k = 1, Nr
333                  do j = jmin,jmax                  do j = jmin,jmax
334                    i = OB_Ie(j,bi,bj)                    i = OB_Ie(j,bi,bj)
335                      IF ( i.EQ.OB_indexNone ) i = 1
336  cgg    Barotropic velocity stored in level 1.  cgg    Barotropic velocity stored in level 1.
337                    tmpy = tmpfldyz(j,1,bi,bj)                    tmpy = tmpfldyz(j,1,bi,bj)
338                    if (maskW(i,j,k,bi,bj) .ne. 0.) then                    if (maskW(i,j,k,bi,bj) .ne. 0.) then
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22