pkg/cost/cost_initvaria.F

C $Header: /u/gcmpack/MITgcm/pkg/cost/cost_initvaria.F,v 1.3 2003/10/02 21:34:45 heimbach Exp $

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

      subroutine cost_initvaria( mythid )

c     ==================================================================
c     SUBROUTINE cost_initvaria
c     ==================================================================
c
c     o Initialise the variable cost function part.
c
c     started: Christian Eckert eckert@mit.edu 30-Jun-1999
c
c     changed: Christian Eckert eckert@mit.edu 18-Apr-2000
c
c              - Restructured the code in order to create a package
c                for the MITgcmUV.
c
c     ==================================================================
c     SUBROUTINE cost_initvaria
c     ==================================================================

      implicit none

c     == global variables ==

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

#include "cost.h"

c     == routine arguments ==

      integer mythid

c     == local variables ==

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

c     == external functions ==

c     == end of interface ==
      jtlo = mybylo(mythid)
      jthi = mybyhi(mythid)
      itlo = mybxlo(mythid)
      ithi = mybxhi(mythid)

c--   Initialize the tiled cost function contributions.
      do bj = jtlo,jthi
        do bi = itlo,ithi
          objf_atl(bi,bj)  = 0. _d 0
          objf_test(bi,bj) = 0. _d 0
          objf_tracer(bi,bj) = 0. _d 0
#ifdef ALLOW_COST_VECTOR
          do i=1,sNx
            objf_vector(i,bi,bj) = 0. _d 0
          end do
#endif
#ifdef ALLOW_COST_STATE_FINAL
c          do k=1,nr
            do j=1,sNy
              do i=1,sNx
                objf_state_final(i,j,bi,bj) = 0. _d 0
cph No init. of cost_state_final here,
cph because we need it in ADM*TLM
              end do
            end do
c          end do
#endif
       enddo
      enddo

#ifdef ALLOW_ECCO
      call ecco_cost_initvaria( mythid )
#endif

c--   Initialise the "global" parts of the cost function.
      _BEGIN_MASTER( mythid )
        fc         = 0. _d 0
      _END_MASTER( mythid )

      _BARRIER

      return
      end

1	C $Header: /u/gcmpack/MITgcm/pkg/cost/cost_initvaria.F,v 1.3 2003/10/02 21:34:45 heimbach Exp $
2
3	#include "PACKAGES_CONFIG.h"
4	#include "COST_CPPOPTIONS.h"
5
6	subroutine cost_initvaria( mythid )
7
8	c ==================================================================
9	c SUBROUTINE cost_initvaria
10	c ==================================================================
11	c
12	c o Initialise the variable cost function part.
13	c
14	c started: Christian Eckert eckert@mit.edu 30-Jun-1999
15	c
16	c changed: Christian Eckert eckert@mit.edu 18-Apr-2000
17	c
18	c - Restructured the code in order to create a package
19	c for the MITgcmUV.
20	c
21	c ==================================================================
22	c SUBROUTINE cost_initvaria
23	c ==================================================================
24
25	implicit none
26
27	c == global variables ==
28
29	#include "EEPARAMS.h"
30	#include "SIZE.h"
31
32	#include "cost.h"
33
34	c == routine arguments ==
35
36	integer mythid
37
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 == external functions ==
46
47	c == end of interface ==
48	jtlo = mybylo(mythid)
49	jthi = mybyhi(mythid)
50	itlo = mybxlo(mythid)
51	ithi = mybxhi(mythid)
52
53	c-- Initialize the tiled cost function contributions.
54	do bj = jtlo,jthi
55	do bi = itlo,ithi
56	objf_atl(bi,bj) = 0. _d 0
57	objf_test(bi,bj) = 0. _d 0
58	objf_tracer(bi,bj) = 0. _d 0
59	#ifdef ALLOW_COST_VECTOR
60	do i=1,sNx
61	objf_vector(i,bi,bj) = 0. _d 0
62	end do
63	#endif
64	#ifdef ALLOW_COST_STATE_FINAL
65	c do k=1,nr
66	do j=1,sNy
67	do i=1,sNx
68	objf_state_final(i,j,bi,bj) = 0. _d 0
69	cph No init. of cost_state_final here,
70	cph because we need it in ADM*TLM
71	end do
72	end do
73	c end do
74	#endif
75	enddo
76	enddo
77
78	#ifdef ALLOW_ECCO
79	call ecco_cost_initvaria( mythid )
80	#endif
81
82	c-- Initialise the "global" parts of the cost function.
83	_BEGIN_MASTER( mythid )
84	fc = 0. _d 0
85	_END_MASTER( mythid )
86
87	_BARRIER
88
89	return
90	end
91