/[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.1 by heimbach, Fri Jul 13 13:37:45 2001 UTC revision 1.5 by jmc, Tue Jul 13 18:08:06 2004 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "COST_CPPOPTIONS.h"
5    
6        subroutine cost_tracer( bi, bj, myThid )        subroutine cost_tracer( bi, bj, myThid )
7  C     /==========================================================\  C     /==========================================================\
# 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    # if (defined (ALLOW_PASSIVE_TRACER))
26    #  include "TR1.h"
27    # elif (defined (ALLOW_PTRACERS))
28    #  include "PTRACERS_SIZE.h"
29    #  include "PTRACERS.h"
30    # endif
31    #endif
32    
33  #include "cost.h"  #include "cost.h"
34    
# Line 44  ce    some reference temperature Line 52  ce    some reference temperature
52        k=1        k=1
53        DO j=1,sNy        DO j=1,sNy
54           DO i=1,sNx           DO i=1,sNx
55    #if (defined (ALLOW_PASSIVE_TRACER))
56              locfc = locfc + maskC(i,j,k,bi,bj)*              locfc = locfc + maskC(i,j,k,bi,bj)*
57       &           lambdaTr1ClimRelax*Tr1(i,j,k,bi,bj)*       &           lambdaTr1ClimRelax*Tr1(i,j,k,bi,bj)*
58       &           rA(i,j,bi,bj)*drF(k)*deltaTtracer       &           rA(i,j,bi,bj)*drF(k)*deltaTtracer
59    #elif (defined (ALLOW_PTRACERS))
60                locfc = locfc + maskC(i,j,k,bi,bj)*
61         &           lambdaTr1ClimRelax*ptracer(i,j,k,bi,bj,1)*
62         &           rA(i,j,bi,bj)*drF(k)*deltaTtracer
63    #endif
64           ENDDO           ENDDO
65        ENDDO        ENDDO
66    
 cph(  
       print *, 'COST TRACER vor ', objf_tracer(bi,bj)  
 cph)  
   
67        objf_tracer(bi,bj) = objf_tracer(bi,bj) + locfc        objf_tracer(bi,bj) = objf_tracer(bi,bj) + locfc
68    
69  cph(  #if (defined (ALLOW_PASSIVE_TRACER))
70        print *, 'COST TRACER nach', objf_tracer(bi,bj)        print *, 'COST TRACER nach', objf_tracer(bi,bj),
71  cph)       &     Tr1(83,33,1,1,1), Tr1(83,33,2,1,1)
72    #elif (defined (ALLOW_PTRACERS))
73          print *, 'COST TRACER nach', objf_tracer(bi,bj),
74         &     ptracer(83,33,1,1,1,1), ptracer(83,33,2,1,1,1)
75    #endif
76    
77  #endif /* ALLOW_COST_TRACER */  #endif /* ALLOW_COST_TRACER */
78    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22