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

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

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

revision 1.2 by heimbach, Thu Jan 17 17:03:34 2002 UTC revision 1.8 by jmc, Mon Nov 5 18:56:37 2007 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "COST_CPPOPTIONS.h"  #include "COST_CPPOPTIONS.h"
5    
# Line 20  C     == Global variables === Line 21  C     == Global variables ===
21  #include "PARAMS.h"  #include "PARAMS.h"
22  #include "DYNVARS.h"  #include "DYNVARS.h"
23  #include "GRID.h"  #include "GRID.h"
24  #include "TR1.h"  #ifdef ALLOW_COST_TRACER
25    # ifdef ALLOW_PTRACERS
26    #  include "PTRACERS_SIZE.h"
27    c# include "PTRACERS_PARAMS.h"
28    #  include "PTRACERS_FIELDS.h"
29    # endif
30    #endif
31    
32  #include "cost.h"  #include "cost.h"
33    
# Line 44  ce    some reference temperature Line 51  ce    some reference temperature
51        k=1        k=1
52        DO j=1,sNy        DO j=1,sNy
53           DO i=1,sNx           DO i=1,sNx
54    #ifdef ALLOW_PTRACERS
55              locfc = locfc + maskC(i,j,k,bi,bj)*              locfc = locfc + maskC(i,j,k,bi,bj)*
56       &           lambdaTr1ClimRelax*Tr1(i,j,k,bi,bj)*       &           lambdaTr1ClimRelax*ptracer(i,j,k,bi,bj,1)*
57       &           rA(i,j,bi,bj)*drF(k)*deltaTtracer       &           rA(i,j,bi,bj)*drF(k)*dTtracerLev(k)
58    #endif
59           ENDDO           ENDDO
60        ENDDO        ENDDO
61    
 cph(  
       print *, 'COST TRACER vor ', objf_tracer(bi,bj)  
 cph)  
   
62        objf_tracer(bi,bj) = objf_tracer(bi,bj) + locfc        objf_tracer(bi,bj) = objf_tracer(bi,bj) + locfc
63    
64  cph(  #ifdef ALLOW_PTRACERS
65        print *, 'COST TRACER nach', objf_tracer(bi,bj)        print *, 'COST TRACER nach', objf_tracer(bi,bj),
66  cph)       &     ptracer(83,33,1,1,1,1), ptracer(83,33,2,1,1,1)
67    #endif
68    
69  #endif /* ALLOW_COST_TRACER */  #endif /* ALLOW_COST_TRACER */
70    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22