pkg/cost/cost_final.F

C $Header: /u/gcmpack/MITgcm/pkg/cost/cost_final.F,v 1.8 2003/11/14 23:07:30 heimbach Exp $

#include "PACKAGES_CONFIG.h"
#include "COST_CPPOPTIONS.h"

      subroutine cost_final( mythid )

c     ==================================================================
c     SUBROUTINE cost_final
c     ==================================================================
c
c     o Sum of all cost function contributions.
c
c     started: Christian Eckert eckert@mit.edu 30-Jun-1999
c     changed: Christian Eckert eckert@mit.edu 25-Feb-2000
c     heimbach@mit.edu 05-Nov-2003 Modularize cost package
c
c     ==================================================================
c     SUBROUTINE cost_final
c     ==================================================================

      implicit none

c     == global variables ==

#include "EEPARAMS.h"
#include "SIZE.h"
#include "PARAMS.h"

#include "cost.h"
#include "ctrl.h"

c     == routine arguments ==

      integer mythid

#ifdef ALLOW_COST
c     == local variables ==

      integer i,j,k
      integer bi,bj
      integer itlo,ithi
      integer jtlo,jthi

c     == end of interface ==

      jtlo = mybylo(mythid)
      jthi = mybyhi(mythid)
      itlo = mybxlo(mythid)
      ithi = mybxhi(mythid)

#if (defined (ALLOW_ECCO))
      CALL ECCO_COST_FINAL (myThid)

#elif (defined (ALLOW_COST_VECTOR))
      CALL COST_VECTOR (myThid)

#elif (defined (ALLOW_COST_STATE_FINAL))
      CALL COST_STATE_FINAL (myThid)

#endif /* above stuff undef */

#ifndef ALLOW_ECCO

# ifdef ALLOW_COST_TEST
      CALL COST_TEST (myThid)
# endif
# ifdef ALLOW_COST_ATLANTIC_HEAT
      CALL COST_ATLANTIC_HEAT (myThid)
# endif

c--   Sum up all contributions.
      do bj = jtlo,jthi
        do bi = itlo,ithi

          print*,' --> objf_test(bi,bj)   =',objf_test(bi,bj)
          print*,' --> objf_tracer(bi,bj) =',objf_tracer(bi,bj)
          print*,' --> objf_atl(bi,bj)    =',objf_atl(bi,bj)

          fc = fc
     &            + mult_test   * objf_test(bi,bj)
     &            + mult_tracer * objf_tracer(bi,bj)
     &            + mult_atl    * objf_atl(bi,bj)
        enddo
      enddo

      print*,' local fc = ', fc

c--   Do global summation.
      _GLOBAL_SUM_R8( fc , myThid )

      print*,' global fc = ', fc

#endif /* ALLOW_ECCO */

c--   set averaging freq. to zero to avoid re-write of
c--   averaged fields in reverse checkpointing loops
      taveFreq   = 0.
      dumpFreq   = 0.
      pChkptFreq = 0.

#endif /* ALLOW_COST */

      return
      end

1	C $Header: /u/gcmpack/MITgcm/pkg/cost/cost_final.F,v 1.8 2003/11/14 23:07:30 heimbach Exp $
2
3	#include "PACKAGES_CONFIG.h"
4	#include "COST_CPPOPTIONS.h"
5
6	subroutine cost_final( mythid )
7
8	c ==================================================================
9	c SUBROUTINE cost_final
10	c ==================================================================
11	c
12	c o Sum of all cost function contributions.
13	c
14	c started: Christian Eckert eckert@mit.edu 30-Jun-1999
15	c changed: Christian Eckert eckert@mit.edu 25-Feb-2000
16	c heimbach@mit.edu 05-Nov-2003 Modularize cost package
17	c
18	c ==================================================================
19	c SUBROUTINE cost_final
20	c ==================================================================
21
22	implicit none
23
24	c == global variables ==
25
26	#include "EEPARAMS.h"
27	#include "SIZE.h"
28	#include "PARAMS.h"
29
30	#include "cost.h"
31	#include "ctrl.h"
32
33	c == routine arguments ==
34
35	integer mythid
36
37	#ifdef ALLOW_COST
38	c == local variables ==
39
40	integer i,j,k
41	integer bi,bj
42	integer itlo,ithi
43	integer jtlo,jthi
44
45	c == end of interface ==
46
47	jtlo = mybylo(mythid)
48	jthi = mybyhi(mythid)
49	itlo = mybxlo(mythid)
50	ithi = mybxhi(mythid)
51
52	#if (defined (ALLOW_ECCO))
53	CALL ECCO_COST_FINAL (myThid)
54
55	#elif (defined (ALLOW_COST_VECTOR))
56	CALL COST_VECTOR (myThid)
57
58	#elif (defined (ALLOW_COST_STATE_FINAL))
59	CALL COST_STATE_FINAL (myThid)
60
61	#endif /* above stuff undef */
62
63	#ifndef ALLOW_ECCO
64
65	# ifdef ALLOW_COST_TEST
66	CALL COST_TEST (myThid)
67	# endif
68	# ifdef ALLOW_COST_ATLANTIC_HEAT
69	CALL COST_ATLANTIC_HEAT (myThid)
70	# endif
71
72	c-- Sum up all contributions.
73	do bj = jtlo,jthi
74	do bi = itlo,ithi
75
76	print*,' --> objf_test(bi,bj) =',objf_test(bi,bj)
77	print*,' --> objf_tracer(bi,bj) =',objf_tracer(bi,bj)
78	print*,' --> objf_atl(bi,bj) =',objf_atl(bi,bj)
79
80	fc = fc
81	& + mult_test * objf_test(bi,bj)
82	& + mult_tracer * objf_tracer(bi,bj)
83	& + mult_atl * objf_atl(bi,bj)
84	enddo
85	enddo
86
87	print*,' local fc = ', fc
88
89	c-- Do global summation.
90	_GLOBAL_SUM_R8( fc , myThid )
91
92	print*,' global fc = ', fc
93
94	#endif /* ALLOW_ECCO */
95
96	c-- set averaging freq. to zero to avoid re-write of
97	c-- averaged fields in reverse checkpointing loops
98	taveFreq = 0.
99	dumpFreq = 0.
100	pChkptFreq = 0.
101
102	#endif /* ALLOW_COST */
103
104	return
105	end
106