/[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.16 by jmc, Mon Oct 8 23:59:21 2007 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "PACKAGES_CONFIG.h"  #include "PACKAGES_CONFIG.h"
5  #include "COST_CPPOPTIONS.h"  #include "COST_CPPOPTIONS.h"
# Line 28  c     == global variables == Line 29  c     == global variables ==
29  #include "PARAMS.h"  #include "PARAMS.h"
30    
31  #include "cost.h"  #include "cost.h"
32  #include "ctrl.h"  #ifdef ALLOW_CTRL
33    # include "ctrl.h"
34    #endif
35    
36  c     == routine arguments ==  c     == routine arguments ==
37    
# Line 49  c     == end of interface == Line 52  c     == end of interface ==
52        itlo = mybxlo(mythid)        itlo = mybxlo(mythid)
53        ithi = mybxhi(mythid)        ithi = mybxhi(mythid)
54    
55    #ifdef ALLOW_SEAICE
56          CALL SEAICE_COST_FINAL (myThid)
57    #endif
58    
59  #if (defined (ALLOW_ECCO))  #if (defined (ALLOW_ECCO))
60        CALL ECCO_COST_FINAL (myThid)        CALL ECCO_COST_FINAL (myThid)
61    
# Line 73  c--   Sum up all contributions. Line 80  c--   Sum up all contributions.
80        do bj = jtlo,jthi        do bj = jtlo,jthi
81          do bi = itlo,ithi          do bi = itlo,ithi
82    
83            print*,' --> objf_test(bi,bj)   =',objf_test(bi,bj)            write(standardmessageunit,'(A,D22.15)')
84            print*,' --> objf_tracer(bi,bj) =',objf_tracer(bi,bj)       &          ' --> objf_test(bi,bj)   = ', objf_test(bi,bj)
85            print*,' --> objf_atl(bi,bj)    =',objf_atl(bi,bj)            write(standardmessageunit,'(A,D22.15)')
86         &         ' --> objf_tracer(bi,bj) = ', objf_tracer(bi,bj)
87              write(standardmessageunit,'(A,D22.15)')
88         &         ' --> objf_atl(bi,bj)    = ', objf_atl(bi,bj)
89    #ifdef ALLOW_COST_TRANSPORT
90              write(standardmessageunit,'(A,D22.15)')
91         &         ' --> objf_transport(bi,bj) = ', objf_transport(bi,bj)
92    #endif
93    
94            fc = fc            fc = fc
95       &            + mult_test   * objf_test(bi,bj)       &            + mult_test   * objf_test(bi,bj)
96       &            + mult_tracer * objf_tracer(bi,bj)       &            + mult_tracer * objf_tracer(bi,bj)
97       &            + mult_atl    * objf_atl(bi,bj)       &            + mult_atl    * objf_atl(bi,bj)
98    #ifdef ALLOW_COST_TRANSPORT
99         &            + mult_transport * objf_transport(bi,bj)
100    #endif
101          enddo          enddo
102        enddo        enddo
103    
104        print*,' local fc = ', fc        write(standardmessageunit,'(A,D22.15)') '  local fc = ', fc
105    
106  c--   Do global summation.  c--   Do global summation.
107        _GLOBAL_SUM_R8( fc , myThid )        _GLOBAL_SUM_R8( fc , myThid )
108    
109        print*,' global fc = ', fc        write(standardmessageunit,'(A,D22.15)') ' global fc = ', fc
110    
111  #endif /* ALLOW_ECCO */  #endif /* ALLOW_ECCO */
112    
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22