28 |
# include "PTRACERS_FIELDS.h" |
# include "PTRACERS_FIELDS.h" |
29 |
# endif |
# endif |
30 |
#endif |
#endif |
|
|
|
31 |
#include "cost.h" |
#include "cost.h" |
32 |
|
|
33 |
C == Routine arguments == |
C == Routine arguments == |
42 |
|
|
43 |
integer i, j, k |
integer i, j, k |
44 |
integer ig, jg |
integer ig, jg |
45 |
|
|
46 |
ce some reference temperature |
ce some reference temperature |
47 |
thetaRef = 24.0D0 |
thetaRef = 24.0D0 |
48 |
|
|
52 |
DO j=1,sNy |
DO j=1,sNy |
53 |
DO i=1,sNx |
DO i=1,sNx |
54 |
#ifdef ALLOW_PTRACERS |
#ifdef ALLOW_PTRACERS |
55 |
locfc = locfc + maskC(i,j,k,bi,bj)* |
locfc = locfc + hFacC(i,j,k,bi,bj)* |
56 |
& lambdaTr1ClimRelax*ptracer(i,j,k,bi,bj,1)* |
& lambdaTr1ClimRelax*ptracer(i,j,k,bi,bj,1)* |
57 |
& rA(i,j,bi,bj)*drF(k)*dTtracerLev(k) |
& rA(i,j,bi,bj)*drF(k)*dTtracerLev(k) |
58 |
#endif |
#endif |