/[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.3 by jmc, Mon Aug 11 22:29:55 2008 UTC revision 1.4 by jmc, Tue Apr 28 18:13:28 2009 UTC
# Line 206  c--     Read time averages and the month Line 206  c--     Read time averages and the month
206    
207  cgg     Minor problem : grad T,S needs overlap values.  cgg     Minor problem : grad T,S needs overlap values.
208          _BARRIER          _BARRIER
209          _EXCH_XYZ_R8(tbar,myThid)          _EXCH_XYZ_RL(tbar,myThid)
210          _EXCH_XYZ_R8(sbar,myThid)          _EXCH_XYZ_RL(sbar,myThid)
211    
212    
213          do bj = jtlo,jthi          do bj = jtlo,jthi
# Line 244  C-- jmc: both calls below are wrong if m Line 244  C-- jmc: both calls below are wrong if m
244       I                         sbar(1-OLx,1-OLy,k,bi,bj),       I                         sbar(1-OLx,1-OLy,k,bi,bj),
245       O                         rholoc,       O                         rholoc,
246       I                         k, bi, bj, myThid )       I                         k, bi, bj, myThid )
247                _EXCH_XY_R8(rholoc , myThid)                _EXCH_XY_RL(rholoc , myThid)
248    
249  cgg   Compute centered difference horizontal gradient on bdy.  cgg   Compute centered difference horizontal gradient on bdy.
250                do i = imin, imax                do i = imin, imax
# Line 309  cgg    All these points need to be wet. Line 309  cgg    All these points need to be wet.
309       O                         rholoc,       O                         rholoc,
310       I                         k, bi, bj, myThid )       I                         k, bi, bj, myThid )
311    
312                _EXCH_XY_R8(rholoc , myThid)                _EXCH_XY_RL(rholoc , myThid)
313    
314  cgg   Compute centered difference horizontal gradient on bdy.  cgg   Compute centered difference horizontal gradient on bdy.
315                 do i = imin, imax                 do i = imin, imax
# Line 372  cgg    All these points need to be wet. Line 372  cgg    All these points need to be wet.
372       I                         sbar(1-OLx,1-OLy,k,bi,bj),       I                         sbar(1-OLx,1-OLy,k,bi,bj),
373       O                         rholoc,       O                         rholoc,
374       I                         k, bi, bj, myThid )       I                         k, bi, bj, myThid )
375                _EXCH_XY_R8(rholoc , myThid)                _EXCH_XY_RL(rholoc , myThid)
376    
377  cgg   Compute centered difference horizontal gradient on bdy.  cgg   Compute centered difference horizontal gradient on bdy.
378                do j = jmin, jmax                do j = jmin, jmax
# Line 435  cgg    All these points need to be wet. Line 435  cgg    All these points need to be wet.
435       I                         sbar(1-OLx,1-OLy,k,bi,bj),       I                         sbar(1-OLx,1-OLy,k,bi,bj),
436       O                         rholoc,       O                         rholoc,
437       I                         k, bi, bj, myThid )       I                         k, bi, bj, myThid )
438                _EXCH_XY_R8(rholoc , myThid)                _EXCH_XY_RL(rholoc , myThid)
439    
440  cgg   Compute centered difference horizontal gradient on bdy.  cgg   Compute centered difference horizontal gradient on bdy.
441                do j = jmin, jmax                do j = jmin, jmax
# Line 498  c--         Print cost function for each Line 498  c--         Print cost function for each
498    
499  #ifdef ECCO_VERBOSE  #ifdef ECCO_VERBOSE
500  c--     Print cost function for all tiles.  c--     Print cost function for all tiles.
501          _GLOBAL_SUM_R8( fcthread , myThid )          _GLOBAL_SUM_RL( fcthread , myThid )
502          write(msgbuf,'(a)') ' '          write(msgbuf,'(a)') ' '
503          call print_message( msgbuf, standardmessageunit,          call print_message( msgbuf, standardmessageunit,
504       &                      SQUEEZE_RIGHT , mythid)       &                      SQUEEZE_RIGHT , mythid)
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22