/[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.4 by heimbach, Fri Jun 27 01:54:20 2003 UTC revision 1.6 by heimbach, Fri Sep 17 22:57:12 2004 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "COST_CPPOPTIONS.h"  #include "COST_CPPOPTIONS.h"
 #ifdef ALLOW_PTRACERS  
 # include "PTRACERS_OPTIONS.h"  
 #endif  
5    
6        subroutine cost_tracer( bi, bj, myThid )        subroutine cost_tracer( bi, bj, myThid )
7  C     /==========================================================\  C     /==========================================================\
# Line 24  C     == Global variables === Line 22  C     == Global variables ===
22  #include "DYNVARS.h"  #include "DYNVARS.h"
23  #include "GRID.h"  #include "GRID.h"
24  #ifdef ALLOW_COST_TRACER  #ifdef ALLOW_COST_TRACER
25  # if (defined (ALLOW_PASSIVE_TRACER))  # ifdef ALLOW_PTRACERS
26  #  include "TR1.h"  #  include "PTRACERS_SIZE.h"
 # elif (defined (ALLOW_PTRACERS))  
27  #  include "PTRACERS.h"  #  include "PTRACERS.h"
28  # endif  # endif
29  #endif  #endif
# Line 53  ce    some reference temperature Line 50  ce    some reference temperature
50        k=1        k=1
51        DO j=1,sNy        DO j=1,sNy
52           DO i=1,sNx           DO i=1,sNx
53  #if (defined (ALLOW_PASSIVE_TRACER))  #ifdef ALLOW_PTRACERS
             locfc = locfc + maskC(i,j,k,bi,bj)*  
      &           lambdaTr1ClimRelax*Tr1(i,j,k,bi,bj)*  
      &           rA(i,j,bi,bj)*drF(k)*deltaTtracer  
 #elif (defined (ALLOW_PTRACERS))  
54              locfc = locfc + maskC(i,j,k,bi,bj)*              locfc = locfc + maskC(i,j,k,bi,bj)*
55       &           lambdaTr1ClimRelax*ptracer(i,j,k,bi,bj,1)*       &           lambdaTr1ClimRelax*ptracer(i,j,k,bi,bj,1)*
56       &           rA(i,j,bi,bj)*drF(k)*deltaTtracer       &           rA(i,j,bi,bj)*drF(k)*deltaTtracer
# Line 67  ce    some reference temperature Line 60  ce    some reference temperature
60    
61        objf_tracer(bi,bj) = objf_tracer(bi,bj) + locfc        objf_tracer(bi,bj) = objf_tracer(bi,bj) + locfc
62    
63  #if (defined (ALLOW_PASSIVE_TRACER))  #ifdef ALLOW_PTRACERS
       print *, 'COST TRACER nach', objf_tracer(bi,bj),  
      &     Tr1(83,33,1,1,1), Tr1(83,33,2,1,1)  
 #elif (defined (ALLOW_PTRACERS))  
64        print *, 'COST TRACER nach', objf_tracer(bi,bj),        print *, 'COST TRACER nach', objf_tracer(bi,bj),
65       &     ptracer(83,33,1,1,1,1), ptracer(83,33,2,1,1,1)       &     ptracer(83,33,1,1,1,1), ptracer(83,33,2,1,1,1)
66  #endif  #endif
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22