15 |
I myThid ) |
I myThid ) |
16 |
|
|
17 |
C !DESCRIPTION: |
C !DESCRIPTION: |
18 |
C Calculates the vertical flux due to down-gradient diffusion of a tracer: |
C Calculates the area integrated vertical flux due to down-gradient |
19 |
|
C diffusion of a tracer: |
20 |
C \begin{equation*} |
C \begin{equation*} |
21 |
C F^r_{diff} = - \kappa_r \partial_r \theta |
C F^r_{diff} = - A^r \kappa_r \frac{1}{\Delta r_c} \delta_k \theta |
22 |
C \end{equation*} |
C \end{equation*} |
23 |
|
|
24 |
C !USES: =============================================================== |
C !USES: =============================================================== |
47 |
C i,j :: loop indices |
C i,j :: loop indices |
48 |
C km1 :: =k+1 for k<Nr, =Nr for k>-Nr |
C km1 :: =k+1 for k<Nr, =Nr for k>-Nr |
49 |
INTEGER i,j,km1 |
INTEGER i,j,km1 |
50 |
|
CEOP |
51 |
|
|
52 |
km1=max(1,k-1) |
km1=max(1,k-1) |
53 |
|
|