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 |
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)*dTtracerLev(k) |
57 |
#endif |
#endif |
58 |
ENDDO |
ENDDO |
59 |
ENDDO |
ENDDO |
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 |