/[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.4 by jmc, Tue Apr 28 18:13:28 2009 UTC revision 1.5 by jmc, Tue May 24 20:45:48 2011 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "COST_CPPOPTIONS.h"  #include "COST_CPPOPTIONS.h"
 #ifdef ALLOW_OBCS  
 # include "OBCS_OPTIONS.h"  
 #endif  
5    
6        subroutine cost_obcs_ageos( myiter, mytime, mythid )        subroutine cost_obcs_ageos( myiter, mytime, mythid )
7    
# Line 33  c     == global variables == Line 30  c     == global variables ==
30  #include "DYNVARS.h"  #include "DYNVARS.h"
31  #include "PARAMS.h"  #include "PARAMS.h"
32  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
33  # include "OBCS.h"  # include "OBCS_GRID.h"
34  #endif  #endif
35    
36  #include "cal.h"  #include "cal.h"
# Line 236  cgg    All these points need to be wet. Line 233  cgg    All these points need to be wet.
233              do k = 1,nr              do k = 1,nr
234    
235  C-- jmc: both calls below are wrong if more than 1 tile => stop here  C-- jmc: both calls below are wrong if more than 1 tile => stop here
236         IF ( bi.NE.1 .OR. bj.NE.1 )         IF ( bi.NE.1 .OR. bj.NE.1 )
237       &  STOP 'COST_OBCS_AGEOS wrong with more than 1 tile/proc'       &  STOP 'COST_OBCS_AGEOS wrong with more than 1 tile/proc'
238                call find_rho_2d(                call find_rho_2d(
239       I                         iMin, iMax, jMin, jMax, k,       I                         iMin, iMax, jMin, jMax, k,
240       I                         tbar(1-OLx,1-OLy,k,bi,bj),       I                         tbar(1-OLx,1-OLy,k,bi,bj),
241       I                         sbar(1-OLx,1-OLy,k,bi,bj),       I                         sbar(1-OLx,1-OLy,k,bi,bj),
242       O                         rholoc,       O                         rholoc,
243       I                         k, bi, bj, myThid )       I                         k, bi, bj, myThid )
# Line 302  cgg    All these points need to be wet. Line 299  cgg    All these points need to be wet.
299    
300              do k = 1,nr              do k = 1,nr
301    
302    C-- jmc: both calls below are wrong if more than 1 tile => stop here
303           IF ( bi.NE.1 .OR. bj.NE.1 )
304         &  STOP 'COST_OBCS_AGEOS wrong with more than 1 tile/proc'
305                call find_rho_2d(                call find_rho_2d(
306       I                         iMin, iMax, jMin, jMax, k,       I                         iMin, iMax, jMin, jMax, k,
307       I                         tbar(1-OLx,1-OLy,k,bi,bj),       I                         tbar(1-OLx,1-OLy,k,bi,bj),
308       I                         sbar(1-OLx,1-OLy,k,bi,bj),       I                         sbar(1-OLx,1-OLy,k,bi,bj),
309       O                         rholoc,       O                         rholoc,
310       I                         k, bi, bj, myThid )       I                         k, bi, bj, myThid )
# Line 366  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           IF ( bi.NE.1 .OR. bj.NE.1 )
370         &  STOP 'COST_OBCS_AGEOS wrong with more than 1 tile/proc'
371                call find_rho_2d(                call find_rho_2d(
372       I                         iMin, iMax, jMin, jMax, k,       I                         iMin, iMax, jMin, jMax, k,
373       I                         tbar(1-OLx,1-OLy,k,bi,bj),       I                         tbar(1-OLx,1-OLy,k,bi,bj),
374       I                         sbar(1-OLx,1-OLy,k,bi,bj),       I                         sbar(1-OLx,1-OLy,k,bi,bj),
375       O                         rholoc,       O                         rholoc,
376       I                         k, bi, bj, myThid )       I                         k, bi, bj, myThid )
# Line 429  cgg    All these points need to be wet. Line 431  cgg    All these points need to be wet.
431    
432              do k = 1,nr              do k = 1,nr
433    
434           IF ( bi.NE.1 .OR. bj.NE.1 )
435         &  STOP 'COST_OBCS_AGEOS wrong with more than 1 tile/proc'
436                call find_rho_2d(                call find_rho_2d(
437       I                         iMin, iMax, jMin, jMax, k,       I                         iMin, iMax, jMin, jMax, k,
438       I                         tbar(1-OLx,1-OLy,k,bi,bj),       I                         tbar(1-OLx,1-OLy,k,bi,bj),
439       I                         sbar(1-OLx,1-OLy,k,bi,bj),       I                         sbar(1-OLx,1-OLy,k,bi,bj),
440       O                         rholoc,       O                         rholoc,
441       I                         k, bi, bj, myThid )       I                         k, bi, bj, myThid )
# Line 523  c--   End of loop over records. Line 527  c--   End of loop over records.
527    
528        return        return
529        end        end
   
   
   
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22