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

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

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

revision 1.2 by jmc, Tue Oct 9 00:02:50 2007 UTC revision 1.3 by jmc, Mon Aug 11 22:29:55 2008 UTC
# Line 234  cgg    All these points need to be wet. Line 234  cgg    All these points need to be wet.
234              enddo              enddo
235    
236              do k = 1,nr              do k = 1,nr
237                call find_rho(bi,bj,imin,imax,jmin,jmax,k,k,  
238       &                      tbar,sbar,rholoc,mythid)  C-- jmc: both calls below are wrong if more than 1 tile => stop here
239           IF ( bi.NE.1 .OR. bj.NE.1 )
240         &  STOP 'COST_OBCS_AGEOS wrong with more than 1 tile/proc'
241                  call find_rho_2d(
242         I                         iMin, iMax, jMin, jMax, k,
243         I                         tbar(1-OLx,1-OLy,k,bi,bj),
244         I                         sbar(1-OLx,1-OLy,k,bi,bj),
245         O                         rholoc,
246         I                         k, bi, bj, myThid )
247                _EXCH_XY_R8(rholoc , myThid)                _EXCH_XY_R8(rholoc , myThid)
248    
249  cgg   Compute centered difference horizontal gradient on bdy.  cgg   Compute centered difference horizontal gradient on bdy.
# Line 294  cgg    All these points need to be wet. Line 302  cgg    All these points need to be wet.
302    
303              do k = 1,nr              do k = 1,nr
304    
305                call find_rho(bi,bj,imin,imax,jmin,jmax,k,k,                call find_rho_2d(
306       &                      tbar,sbar,rholoc,mythid)       I                         iMin, iMax, jMin, jMax, k,
307         I                         tbar(1-OLx,1-OLy,k,bi,bj),
308         I                         sbar(1-OLx,1-OLy,k,bi,bj),
309         O                         rholoc,
310         I                         k, bi, bj, myThid )
311    
312                _EXCH_XY_R8(rholoc , myThid)                _EXCH_XY_R8(rholoc , myThid)
313    
314  cgg   Compute centered difference horizontal gradient on bdy.  cgg   Compute centered difference horizontal gradient on bdy.
# Line 353  cgg    All these points need to be wet. Line 366  cgg    All these points need to be wet.
366    
367              do k = 1,nr              do k = 1,nr
368    
369                call find_rho(bi,bj,imin,imax,jmin,jmax,k,k,                call find_rho_2d(
370       &                      tbar,sbar,rholoc,mythid)       I                         iMin, iMax, jMin, jMax, k,
371         I                         tbar(1-OLx,1-OLy,k,bi,bj),
372         I                         sbar(1-OLx,1-OLy,k,bi,bj),
373         O                         rholoc,
374         I                         k, bi, bj, myThid )
375                _EXCH_XY_R8(rholoc , myThid)                _EXCH_XY_R8(rholoc , myThid)
376    
377  cgg   Compute centered difference horizontal gradient on bdy.  cgg   Compute centered difference horizontal gradient on bdy.
# Line 412  cgg    All these points need to be wet. Line 429  cgg    All these points need to be wet.
429    
430              do k = 1,nr              do k = 1,nr
431    
432                call find_rho(bi,bj,imin,imax,jmin,jmax,k,k,                call find_rho_2d(
433       &                      tbar,sbar,rholoc,mythid)       I                         iMin, iMax, jMin, jMax, k,
434         I                         tbar(1-OLx,1-OLy,k,bi,bj),
435         I                         sbar(1-OLx,1-OLy,k,bi,bj),
436         O                         rholoc,
437         I                         k, bi, bj, myThid )
438                _EXCH_XY_R8(rholoc , myThid)                _EXCH_XY_R8(rholoc , myThid)
439    
440  cgg   Compute centered difference horizontal gradient on bdy.  cgg   Compute centered difference horizontal gradient on bdy.
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22