/[MITgcm]/MITgcm/pkg/cost/cost_tracer.F
ViewVC logotype

Contents of /MITgcm/pkg/cost/cost_tracer.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.4 - (show annotations) (download)
Fri Jun 27 01:54:20 2003 UTC (20 years, 10 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint52l_pre, hrcube4, hrcube5, checkpoint52d_pre, checkpoint52j_pre, checkpoint51o_pre, checkpoint51l_post, checkpoint52l_post, checkpoint52k_post, checkpoint54, checkpoint53, checkpoint52, checkpoint52f_post, checkpoint51f_post, checkpoint51d_post, checkpoint51t_post, checkpoint51n_post, checkpoint52i_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint51s_post, checkpoint51j_post, checkpoint52e_pre, checkpoint52e_post, checkpoint51n_pre, checkpoint53d_post, checkpoint52b_pre, checkpoint54b_post, checkpoint51l_pre, checkpoint52m_post, checkpoint51q_post, checkpoint52b_post, checkpoint52c_post, checkpoint51h_pre, checkpoint52f_pre, checkpoint54a_pre, checkpoint53c_post, branchpoint-genmake2, checkpoint54a_post, checkpoint51r_post, checkpoint51i_post, checkpoint51b_post, checkpoint51c_post, checkpoint53a_post, checkpoint52d_post, checkpoint53g_post, checkpoint52a_pre, checkpoint52i_post, checkpoint51i_pre, checkpoint52h_pre, checkpoint53f_post, checkpoint52j_post, branch-netcdf, checkpoint52n_post, checkpoint53b_pre, checkpoint51e_post, checkpoint51o_post, checkpoint51f_pre, checkpoint53b_post, checkpoint52a_post, checkpoint51g_post, ecco_c52_e35, checkpoint51m_post, checkpoint53d_pre, checkpoint51p_post, checkpoint51u_post
Branch point for: branch-genmake2, branch-nonh, tg2-branch, netcdf-sm0, checkpoint51n_branch
Changes since 1.3: +24 -9 lines 
extend cost, ctr, to deal with ALLOW_PTRACERS
1 C $Header: /u/gcmpack/MITgcm/pkg/cost/cost_tracer.F,v 1.3 2002/07/13 03:01:39 heimbach Exp $
2
3 #include "COST_CPPOPTIONS.h"
4 #ifdef ALLOW_PTRACERS
5 # include "PTRACERS_OPTIONS.h"
6 #endif
7
8 subroutine cost_tracer( bi, bj, myThid )
9 C /==========================================================\
10 C | subroutine cost_tracer |
11 C | o this routine computes the cost function for the tiles |
12 C | of this processor |
13 C |==========================================================|
14 C | |
15 C | Notes |
16 C | ===== |
17 C \==========================================================/
18 IMPLICIT NONE
19
20 C == Global variables ===
21 #include "SIZE.h"
22 #include "EEPARAMS.h"
23 #include "PARAMS.h"
24 #include "DYNVARS.h"
25 #include "GRID.h"
26 #ifdef ALLOW_COST_TRACER
27 # if (defined (ALLOW_PASSIVE_TRACER))
28 # include "TR1.h"
29 # elif (defined (ALLOW_PTRACERS))
30 # include "PTRACERS.h"
31 # endif
32 #endif
33
34 #include "cost.h"
35
36 C == Routine arguments ==
37 C myThid - Thread number for this instance of the routine.
38 integer bi, bj
39 integer myThid
40
41 #ifdef ALLOW_COST_TRACER
42 C == Local variables
43 _RL thetaRef
44 _RL locfc
45
46 integer i, j, k
47 integer ig, jg
48 ce some reference temperature
49 thetaRef = 24.0D0
50
51 locfc = 0. _d 0
52
53 k=1
54 DO j=1,sNy
55 DO i=1,sNx
56 #if (defined (ALLOW_PASSIVE_TRACER))
57 locfc = locfc + maskC(i,j,k,bi,bj)*
58 & lambdaTr1ClimRelax*Tr1(i,j,k,bi,bj)*
59 & rA(i,j,bi,bj)*drF(k)*deltaTtracer
60 #elif (defined (ALLOW_PTRACERS))
61 locfc = locfc + maskC(i,j,k,bi,bj)*
62 & lambdaTr1ClimRelax*ptracer(i,j,k,bi,bj,1)*
63 & rA(i,j,bi,bj)*drF(k)*deltaTtracer
64 #endif
65 ENDDO
66 ENDDO
67
68 objf_tracer(bi,bj) = objf_tracer(bi,bj) + locfc
69
70 #if (defined (ALLOW_PASSIVE_TRACER))
71 print *, 'COST TRACER nach', objf_tracer(bi,bj),
72 & Tr1(83,33,1,1,1), Tr1(83,33,2,1,1)
73 #elif (defined (ALLOW_PTRACERS))
74 print *, 'COST TRACER nach', objf_tracer(bi,bj),
75 & ptracer(83,33,1,1,1,1), ptracer(83,33,2,1,1,1)
76 #endif
77
78 #endif /* ALLOW_COST_TRACER */
79
80 RETURN
81 END
  ViewVC Help
Powered by ViewVC 1.1.22