/[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.3 by heimbach, Sat Jul 13 03:01:39 2002 UTC revision 1.4 by heimbach, Fri Jun 27 01:54:20 2003 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    
3  #include "COST_CPPOPTIONS.h"  #include "COST_CPPOPTIONS.h"
4    #ifdef ALLOW_PTRACERS
5    # include "PTRACERS_OPTIONS.h"
6    #endif
7    
8        subroutine cost_tracer( bi, bj, myThid )        subroutine cost_tracer( bi, bj, myThid )
9  C     /==========================================================\  C     /==========================================================\
# Line 20  C     == Global variables === Line 23  C     == Global variables ===
23  #include "PARAMS.h"  #include "PARAMS.h"
24  #include "DYNVARS.h"  #include "DYNVARS.h"
25  #include "GRID.h"  #include "GRID.h"
26  #include "TR1.h"  #ifdef ALLOW_COST_TRACER
27    # if (defined (ALLOW_PASSIVE_TRACER))
28    #  include "TR1.h"
29    # elif (defined (ALLOW_PTRACERS))
30    #  include "PTRACERS.h"
31    # endif
32    #endif
33    
34  #include "cost.h"  #include "cost.h"
35    
# Line 44  ce    some reference temperature Line 53  ce    some reference temperature
53        k=1        k=1
54        DO j=1,sNy        DO j=1,sNy
55           DO i=1,sNx           DO i=1,sNx
56    #if (defined (ALLOW_PASSIVE_TRACER))
57              locfc = locfc + maskC(i,j,k,bi,bj)*              locfc = locfc + maskC(i,j,k,bi,bj)*
58       &           lambdaTr1ClimRelax*Tr1(i,j,k,bi,bj)*       &           lambdaTr1ClimRelax*Tr1(i,j,k,bi,bj)*
59       &           rA(i,j,bi,bj)*drF(k)*deltaTtracer       &           rA(i,j,bi,bj)*drF(k)*deltaTtracer
60    #elif (defined (ALLOW_PTRACERS))
61                locfc = locfc + maskC(i,j,k,bi,bj)*
62         &           lambdaTr1ClimRelax*ptracer(i,j,k,bi,bj,1)*
63         &           rA(i,j,bi,bj)*drF(k)*deltaTtracer
64    #endif
65           ENDDO           ENDDO
66        ENDDO        ENDDO
67    
 cph(  
       print *, 'COST TRACER vor ', objf_tracer(bi,bj)  
 cph)  
   
68        objf_tracer(bi,bj) = objf_tracer(bi,bj) + locfc        objf_tracer(bi,bj) = objf_tracer(bi,bj) + locfc
69    
70  cph(  #if (defined (ALLOW_PASSIVE_TRACER))
71        print *, 'COST TRACER nach', objf_tracer(bi,bj)        print *, 'COST TRACER nach', objf_tracer(bi,bj),
72  cph)       &     Tr1(83,33,1,1,1), Tr1(83,33,2,1,1)
73    #elif (defined (ALLOW_PTRACERS))
74          print *, 'COST TRACER nach', objf_tracer(bi,bj),
75         &     ptracer(83,33,1,1,1,1), ptracer(83,33,2,1,1,1)
76    #endif
77    
78  #endif /* ALLOW_COST_TRACER */  #endif /* ALLOW_COST_TRACER */
79    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22