pkg/cost/cost_accumulate_mean.F

C $Header: /u/gcmpack/MITgcm/pkg/cost/cost_accumulate_mean.F,v 1.1 2004/09/27 18:00:19 heimbach Exp $
C $Name: checkpoint56 $

#include "CPP_OPTIONS.h"

      subroutine cost_accumulate_mean( myThid )
C     *==========================================================*
C     | subroutine cost_accumulate_mean                          |
C     | o accumulate mean state for cost evalualtion             |
C     *==========================================================*
C     |                                                          |
C     *==========================================================*
      IMPLICIT NONE

C     == Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "DYNVARS.h"
#include "GRID.h"

#include "cost.h"

C     == Routine arguments ==
C     myThid - Thread number for this instance of the routine.
      integer bi, bj
      integer myThid

#ifdef ALLOW_COST_ATLANTIC_HEAT
C     == Local variables
      _RL thetaRef

      integer i, j, k
      integer ig, jg
      integer itlo,ithi
      integer jtlo,jthi

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

C--   Calculate cost function on tile of this instance
      do bj = jtlo,jthi
        do bi = itlo,ithi
          do k = 1, Nr
            do j=1,sNy
              do i=1,sNx
                cMeanTheta(i,j,k,bi,bj) = cMeanTheta(i,j,k,bi,bj)
     &                                    + theta(i,j,k,bi,bj)
     &                                    /lastinterval*deltaTClock
                cMeanUVel(i,j,k,bi,bj) = cMeanUVel(i,j,k,bi,bj)
     &                                    + uVel(i,j,k,bi,bj)
     &                                    /lastinterval*deltaTClock
                cMeanVVel(i,j,k,bi,bj) = cMeanVVel(i,j,k,bi,bj)
     &                                    + vVel(i,j,k,bi,bj)
     &                                    /lastinterval*deltaTClock
              end do
            end do
          end do
        end do
      end do

#endif /* ALLOW_COST_ATLANTIC_HEAT */

      END
1	C $Header: /u/gcmpack/MITgcm/pkg/cost/cost_accumulate_mean.F,v 1.1 2004/09/27 18:00:19 heimbach Exp $
2	C $Name: checkpoint56 $
3
4	#include "CPP_OPTIONS.h"
5
6	subroutine cost_accumulate_mean( myThid )
7	C ==========================================================
8	C \| subroutine cost_accumulate_mean \|
9	C \| o accumulate mean state for cost evalualtion \|
10	C ==========================================================
11	C \| \|
12	C ==========================================================
13	IMPLICIT NONE
14
15	C == Global variables ===
16	#include "SIZE.h"
17	#include "EEPARAMS.h"
18	#include "PARAMS.h"
19	#include "DYNVARS.h"
20	#include "GRID.h"
21
22	#include "cost.h"
23
24	C == Routine arguments ==
25	C myThid - Thread number for this instance of the routine.
26	integer bi, bj
27	integer myThid
28
29	#ifdef ALLOW_COST_ATLANTIC_HEAT
30	C == Local variables
31	_RL thetaRef
32
33	integer i, j, k
34	integer ig, jg
35	integer itlo,ithi
36	integer jtlo,jthi
37
38	jtlo = mybylo(mythid)
39	jthi = mybyhi(mythid)
40	itlo = mybxlo(mythid)
41	ithi = mybxhi(mythid)
42
43	C-- Calculate cost function on tile of this instance
44	do bj = jtlo,jthi
45	do bi = itlo,ithi
46	do k = 1, Nr
47	do j=1,sNy
48	do i=1,sNx
49	cMeanTheta(i,j,k,bi,bj) = cMeanTheta(i,j,k,bi,bj)
50	& + theta(i,j,k,bi,bj)
51	& /lastinterval*deltaTClock
52	cMeanUVel(i,j,k,bi,bj) = cMeanUVel(i,j,k,bi,bj)
53	& + uVel(i,j,k,bi,bj)
54	& /lastinterval*deltaTClock
55	cMeanVVel(i,j,k,bi,bj) = cMeanVVel(i,j,k,bi,bj)
56	& + vVel(i,j,k,bi,bj)
57	& /lastinterval*deltaTClock
58	end do
59	end do
60	end do
61	end do
62	end do
63
64	#endif /* ALLOW_COST_ATLANTIC_HEAT */
65
66	END