22 |
C In the unstable case, A and B are non-zero and are chosen so as to |
C In the unstable case, A and B are non-zero and are chosen so as to |
23 |
C conserve total volume of T. |
C conserve total volume of T. |
24 |
|
|
25 |
C !USES: |
C !USES: |
26 |
IMPLICIT NONE |
IMPLICIT NONE |
|
C == Global data == |
|
27 |
#include "SIZE.h" |
#include "SIZE.h" |
28 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
29 |
#include "PARAMS.h" |
#include "PARAMS.h" |
30 |
#include "GRID.h" |
#include "GRID.h" |
31 |
|
|
32 |
C !INPUT/OUTPUT PARAMETERS: |
C !INPUT/OUTPUT PARAMETERS: |
|
C == Routine arguments == |
|
33 |
C bi,bj,k - indices |
C bi,bj,k - indices |
34 |
C rhoKm1 - rho in layer above |
C rhoKm1 - rho in layer above |
35 |
C rhoK - rho in layer below |
C rhoK - rho in layer below |
47 |
|
|
48 |
#ifdef INCLUDE_CONVECT_CALL |
#ifdef INCLUDE_CONVECT_CALL |
49 |
|
|
50 |
C !LOCAL VARIABLES: |
C !LOCAL VARIABLES: |
|
C == Local variables == |
|
51 |
C i,j :: Loop counter |
C i,j :: Loop counter |
52 |
C dS,d1,d2 :: Layer thicknesses |
C dS,d1,d2 :: Layer thicknesses |
53 |
INTEGER i,j |
INTEGER i,j |