pkg/cost/cost_tracer.F

C $Header: /u/gcmpack/MITgcm/pkg/cost/cost_tracer.F,v 1.5 2004/07/13 18:08:06 jmc Exp $
C $Name:  $

#include "COST_CPPOPTIONS.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.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 + maskC(i,j,k,bi,bj)*
     &           lambdaTr1ClimRelax*ptracer(i,j,k,bi,bj,1)*
     &           rA(i,j,bi,bj)*drF(k)*deltaTtracer
#endif
         ENDDO
      ENDDO

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

#ifdef ALLOW_PTRACERS
      print *, 'COST TRACER nach', objf_tracer(bi,bj),
     &     ptracer(83,33,1,1,1,1), ptracer(83,33,2,1,1,1)
#endif

#endif /* ALLOW_COST_TRACER */

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/cost/cost_tracer.F,v 1.5 2004/07/13 18:08:06 jmc Exp $
2	C $Name: $
3
4	#include "COST_CPPOPTIONS.h"
5
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	#ifdef ALLOW_COST_TRACER
25	# ifdef ALLOW_PTRACERS
26	# include "PTRACERS_SIZE.h"
27	# include "PTRACERS.h"
28	# endif
29	#endif
30
31	#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	ce some reference temperature
46	thetaRef = 24.0D0
47
48	locfc = 0. _d 0
49
50	k=1
51	DO j=1,sNy
52	DO i=1,sNx
53	#ifdef ALLOW_PTRACERS
54	locfc = locfc + maskC(i,j,k,bi,bj)*
55	& lambdaTr1ClimRelaxptracer(i,j,k,bi,bj,1)
56	& rA(i,j,bi,bj)drF(k)deltaTtracer
57	#endif
58	ENDDO
59	ENDDO
60
61	objf_tracer(bi,bj) = objf_tracer(bi,bj) + locfc
62
63	#ifdef ALLOW_PTRACERS
64	print *, 'COST TRACER nach', objf_tracer(bi,bj),
65	& ptracer(83,33,1,1,1,1), ptracer(83,33,2,1,1,1)
66	#endif
67
68	#endif /* ALLOW_COST_TRACER */
69
70	RETURN
71	END