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

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

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

revision 1.3 by heimbach, Mon Mar 28 23:49:49 2005 UTC revision 1.4 by heimbach, Thu Oct 26 00:34:05 2006 UTC
# Line 215  c--       Print cost function for each t Line 215  c--       Print cost function for each t
215       &                          SQUEEZE_RIGHT , mythid)       &                          SQUEEZE_RIGHT , mythid)
216            write(msgbuf,'(a,d22.15)')            write(msgbuf,'(a,d22.15)')
217       &        '     cost function (wvel) = ',       &        '     cost function (wvel) = ',
218       &        fctilet+fctiles       &        fcthread_wdrift
219            call print_message( msgbuf, standardmessageunit,            call print_message( msgbuf, standardmessageunit,
220       &                      SQUEEZE_RIGHT , mythid)       &                      SQUEEZE_RIGHT , mythid)
221            write(msgbuf,'(a)') ' '            write(msgbuf,'(a)') ' '
# Line 228  c--       Print cost function for each t Line 228  c--       Print cost function for each t
228    
229  #ifdef ECCO_VERBOSE  #ifdef ECCO_VERBOSE
230  c--     Print cost function for all tiles.  c--     Print cost function for all tiles.
231          _GLOBAL_SUM_R8( fcthread_tdrift , myThid )          _GLOBAL_SUM_R8( fcthread_wdrift , myThid )
         _GLOBAL_SUM_R8( fcthread_tdrifs , myThid )  
232          write(msgbuf,'(a)') ' '          write(msgbuf,'(a)') ' '
233          call print_message( msgbuf, standardmessageunit,          call print_message( msgbuf, standardmessageunit,
234       &                      SQUEEZE_RIGHT , mythid)       &                      SQUEEZE_RIGHT , mythid)
235          write(msgbuf,'(a,i8.8)')          write(msgbuf,'(a,i8.8)')
236       &    ' cost_Drift: irec = ',irec       &    ' cost_Driftw: irec = ',irec
237          call print_message( msgbuf, standardmessageunit,          call print_message( msgbuf, standardmessageunit,
238       &                      SQUEEZE_RIGHT , mythid)       &                      SQUEEZE_RIGHT , mythid)
239          write(msgbuf,'(a,a,d22.15)')          write(msgbuf,'(a,a,d22.15)')
240       &    '  cost function value',       &    '  cost function value',
241       &    ' (wvel) = ',fcthread_tdrifw       &    ' (wvel) = ',fcthread_wdrift
242          call print_message( msgbuf, standardmessageunit,          call print_message( msgbuf, standardmessageunit,
243       &                      SQUEEZE_RIGHT , mythid)       &                      SQUEEZE_RIGHT , mythid)
244          write(msgbuf,'(a)') ' '          write(msgbuf,'(a)') ' '
# Line 247  c--     Print cost function for all tile Line 246  c--     Print cost function for all tile
246       &                      SQUEEZE_RIGHT , mythid)       &                      SQUEEZE_RIGHT , mythid)
247  #endif  #endif
248    
 #else  
 c--   Do not enter the calculation of the drift contribution to  
 c--   the final cost function.  
   
       _BEGIN_MASTER( mythid )  
         write(msgbuf,'(a)') ' '  
         call print_message( msgbuf, standardmessageunit,  
      &                      SQUEEZE_RIGHT , mythid)  
         write(msgbuf,'(a,a)')  
      &    ' cost_Driftw: no contribution of drift between the first ',  
      &                 'and the last year to cost function.'  
         call print_message( msgbuf, standardmessageunit,  
      &                      SQUEEZE_RIGHT , mythid)  
         write(msgbuf,'(a,a,i9.8)')  
      &    ' cost_Driftw: number of records that would have',  
      &                ' been processed: ',nmonsrec  
         call print_message( msgbuf, standardmessageunit,  
      &                      SQUEEZE_RIGHT , mythid)  
         write(msgbuf,'(a)') ' '  
         call print_message( msgbuf, standardmessageunit,  
      &                      SQUEEZE_RIGHT , mythid)  
       _END_MASTER( mythid )  
249  #endif  #endif
250    
251        return        return
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22