/[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.7 - (show annotations) (download)
Sat Dec 4 00:14:26 2004 UTC (19 years, 6 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint57m_post, checkpoint57g_pre, checkpoint57s_post, checkpoint58b_post, checkpoint57b_post, checkpoint57g_post, checkpoint57y_post, checkpoint57r_post, checkpoint57d_post, checkpoint57i_post, checkpoint59, checkpoint58, checkpoint57, checkpoint58f_post, checkpoint57n_post, checkpoint58d_post, checkpoint58a_post, checkpoint57z_post, checkpoint58y_post, checkpoint58t_post, checkpoint58m_post, checkpoint57l_post, checkpoint57t_post, checkpoint57v_post, checkpoint57f_post, checkpoint57a_post, checkpoint57h_pre, checkpoint58w_post, checkpoint57h_post, checkpoint57y_pre, checkpoint58o_post, checkpoint57c_post, checkpoint58p_post, checkpoint58q_post, checkpoint58e_post, checkpoint57c_pre, checkpoint58r_post, checkpoint58n_post, checkpoint57e_post, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59i, checkpoint59h, checkpoint57p_post, checkpint57u_post, checkpoint57q_post, eckpoint57e_pre, checkpoint58k_post, checkpoint58v_post, checkpoint58l_post, checkpoint57h_done, checkpoint57j_post, checkpoint57f_pre, checkpoint58g_post, checkpoint58x_post, checkpoint58h_post, checkpoint58j_post, checkpoint57a_pre, checkpoint57o_post, checkpoint57k_post, checkpoint57w_post, checkpoint58i_post, checkpoint57x_post, checkpoint58c_post, checkpoint58u_post, checkpoint58s_post
Changes since 1.6: +2 -2 lines 
depth convergence accelerator: replace deltaTtracer by dTtracerLev(k)
1 C $Header: /u/gcmpack/MITgcm/pkg/cost/cost_tracer.F,v 1.6 2004/09/17 22:57:12 heimbach 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 & lambdaTr1ClimRelax*ptracer(i,j,k,bi,bj,1)*
56 & rA(i,j,bi,bj)*drF(k)*dTtracerLev(k)
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
  ViewVC Help
Powered by ViewVC 1.1.22