pkg/cost/cost_tracer.F

C $Header: /u/gcmpack/MITgcm/pkg/cost/cost_tracer.F,v 1.11 2011/08/06 19:08:58 heimbach Exp $
C $Name:  $

#include "COST_OPTIONS.h"

      subroutine cost_tracer( bi, bj, myThid )
C     /==========================================================\
C     | subroutine cost_tracer                                   |
C     | o this routine computes the cost function for the tiles  |
C     |   of this processor                                      |
C     |==========================================================|
C     |                                                          |
C     | Notes                                                    |
C     | =====                                                    |
C     \==========================================================/
      IMPLICIT NONE

C     == Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "DYNVARS.h"
#include "GRID.h"
#ifdef ALLOW_COST_TRACER
# ifdef ALLOW_PTRACERS
#  include "PTRACERS_SIZE.h"
#  include "PTRACERS_PARAMS.h"
#  include "PTRACERS_FIELDS.h"
# endif
#endif
#include "cost.h"

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

#ifdef ALLOW_COST_TRACER
C     == Local variables
      _RL thetaRef
      _RL locfc

      integer i, j, k
      integer ig, jg

ce    some reference temperature
      thetaRef = 24.0D0

      locfc = 0. _d 0

      k=1
      DO j=1,sNy
         DO i=1,sNx
#ifdef ALLOW_PTRACERS
            locfc = locfc + hFacC(i,j,k,bi,bj)*
     &           lambdaTr1ClimRelax*ptracer(i,j,k,bi,bj,1)*
     &           rA(i,j,bi,bj)*drF(k)*dTtracerLev(k)
#endif
         ENDDO
      ENDDO

      objf_tracer(bi,bj) = objf_tracer(bi,bj) + locfc

#endif /* ALLOW_COST_TRACER */

      RETURN
      END
1	jmc	1.12	C $Header: /u/gcmpack/MITgcm/pkg/cost/cost_tracer.F,v 1.11 2011/08/06 19:08:58 heimbach Exp $
2	jmc	1.5	C $Name: $
3	heimbach	1.1
4	jmc	1.12	#include "COST_OPTIONS.h"
5	heimbach	1.1
6			subroutine cost_tracer( bi, bj, myThid )
7			C /==========================================================\
8			C \| subroutine cost_tracer \|
9			C \| o this routine computes the cost function for the tiles \|
10			C \| of this processor \|
11			C \|==========================================================\|
12			C \| \|
13			C \| Notes \|
14			C \| ===== \|
15			C \==========================================================/
16			IMPLICIT NONE
17
18			C == Global variables ===
19			#include "SIZE.h"
20			#include "EEPARAMS.h"
21			#include "PARAMS.h"
22			#include "DYNVARS.h"
23			#include "GRID.h"
24	heimbach	1.4	#ifdef ALLOW_COST_TRACER
25	heimbach	1.6	# ifdef ALLOW_PTRACERS
26	jmc	1.5	# include "PTRACERS_SIZE.h"
27	jmc	1.9	# include "PTRACERS_PARAMS.h"
28	jmc	1.8	# include "PTRACERS_FIELDS.h"
29	heimbach	1.4	# endif
30			#endif
31	heimbach	1.1	#include "cost.h"
32
33			C == Routine arguments ==
34			C myThid - Thread number for this instance of the routine.
35			integer bi, bj
36			integer myThid
37
38			#ifdef ALLOW_COST_TRACER
39			C == Local variables
40			_RL thetaRef
41			_RL locfc
42
43			integer i, j, k
44			integer ig, jg
45	heimbach	1.10
46	heimbach	1.1	ce some reference temperature
47			thetaRef = 24.0D0
48
49			locfc = 0. _d 0
50
51			k=1
52			DO j=1,sNy
53			DO i=1,sNx
54	heimbach	1.6	#ifdef ALLOW_PTRACERS
55	heimbach	1.10	locfc = locfc + hFacC(i,j,k,bi,bj)*
56	heimbach	1.4	& lambdaTr1ClimRelaxptracer(i,j,k,bi,bj,1)
57	jmc	1.7	& rA(i,j,bi,bj)drF(k)dTtracerLev(k)
58	heimbach	1.4	#endif
59	heimbach	1.1	ENDDO
60			ENDDO
61
62			objf_tracer(bi,bj) = objf_tracer(bi,bj) + locfc
63
64			#endif /* ALLOW_COST_TRACER */
65
66			RETURN
67			END