48 |
gTracer(i,j,k,bi,bj)=tracer(i,j,k,bi,bj) |
gTracer(i,j,k,bi,bj)=tracer(i,j,k,bi,bj) |
49 |
& +deltaTloc*gTracer(i,j,k,bi,bj) |
& +deltaTloc*gTracer(i,j,k,bi,bj) |
50 |
#ifdef ALLOW_OBCS |
#ifdef ALLOW_OBCS |
51 |
& *maskInC(i,j,bi,bj) |
CML For the Stevens open boundary conditions, we need to step forward |
52 |
|
CML the tracer even on the boundary in order to include surface |
53 |
|
CML forcing and (vertical) diffusion contributions. For other cases |
54 |
|
CML this mask was meant to avoid unrealistic values on the open boundaries. |
55 |
|
CML Lets hope that we can live without this safety net. |
56 |
|
C & *maskInC(i,j,bi,bj) |
57 |
#endif |
#endif |
58 |
ENDDO |
ENDDO |
59 |
ENDDO |
ENDDO |