/[MITgcm]/MITgcm/pkg/cost/cost_final.F
ViewVC logotype

Diff of /MITgcm/pkg/cost/cost_final.F

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

revision 1.10 by heimbach, Sat Dec 4 23:22:28 2004 UTC revision 1.13 by heimbach, Thu Aug 25 16:10:42 2005 UTC
# Line 49  c     == end of interface == Line 49  c     == end of interface ==
49        itlo = mybxlo(mythid)        itlo = mybxlo(mythid)
50        ithi = mybxhi(mythid)        ithi = mybxhi(mythid)
51    
52    #ifdef ALLOW_SEAICE
53          CALL SEAICE_COST_FINAL (myThid)
54    #endif
55    
56  #if (defined (ALLOW_ECCO))  #if (defined (ALLOW_ECCO))
57        CALL ECCO_COST_FINAL (myThid)        CALL ECCO_COST_FINAL (myThid)
58    
# Line 73  c--   Sum up all contributions. Line 77  c--   Sum up all contributions.
77        do bj = jtlo,jthi        do bj = jtlo,jthi
78          do bi = itlo,ithi          do bi = itlo,ithi
79    
80            print*,' --> objf_test(bi,bj)   =',objf_test(bi,bj)            write(standardmessageunit,'(A,D22.15)')
81            print*,' --> objf_tracer(bi,bj) =',objf_tracer(bi,bj)       &          ' --> objf_test(bi,bj)   = ', objf_test(bi,bj)
82            print*,' --> objf_atl(bi,bj)    =',objf_atl(bi,bj)            write(standardmessageunit,'(A,D22.15)')
83         &         ' --> objf_tracer(bi,bj) = ', objf_tracer(bi,bj)
84              write(standardmessageunit,'(A,D22.15)')
85         &         ' --> objf_atl(bi,bj)    = ', objf_atl(bi,bj)
86    
87            fc = fc            fc = fc
88       &            + mult_test   * objf_test(bi,bj)       &            + mult_test   * objf_test(bi,bj)
# Line 84  c--   Sum up all contributions. Line 91  c--   Sum up all contributions.
91          enddo          enddo
92        enddo        enddo
93    
94        print*,' local fc = ', fc        write(standardmessageunit,'(A,D22.15)') '  local fc = ', fc
95    
96  c--   Do global summation.  c--   Do global summation.
97        _GLOBAL_SUM_R8( fc , myThid )        _GLOBAL_SUM_R8( fc , myThid )
98    
99        print*,' global fc = ', fc        write(standardmessageunit,'(A,D22.15)') ' global fc = ', fc
100    
101  #endif /* ALLOW_ECCO */  #endif /* ALLOW_ECCO */
102    
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22