/[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.5 by jmc, Tue Jul 13 18:08:06 2004 UTC revision 1.6 by heimbach, Fri Sep 17 22:57:12 2004 UTC
# Line 22  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
 #  include "TR1.h"  
 # elif (defined (ALLOW_PTRACERS))  
26  #  include "PTRACERS_SIZE.h"  #  include "PTRACERS_SIZE.h"
27  #  include "PTRACERS.h"  #  include "PTRACERS.h"
28  # endif  # endif
# Line 52  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 66  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.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22