pkg/cost/cost_final.F

C $Header: /home/ubuntu/mnt/e9_copy/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.

#endif /* ALLOW_COST */

      return
      end

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