35 |
#include "PARAMS.h" |
#include "PARAMS.h" |
36 |
|
|
37 |
C !INPUT PARAMETERS: =================================================== |
C !INPUT PARAMETERS: =================================================== |
38 |
C bi_arg,bj_arg :: tile indices |
C bi_arg,bj_arg :: tile indices |
39 |
C k :: vertical level |
C k :: vertical level |
40 |
C rTrans :: vertical volume transport |
C rTrans :: vertical volume transport |
41 |
C wVel :: vertical flow |
C wVel :: vertical flow |
42 |
C tracer :: tracer field |
C tracer :: tracer field |
43 |
C myThid :: thread number |
C myThid :: thread number |
44 |
INTEGER bi_arg,bj_arg,k |
INTEGER bi_arg,bj_arg,k |
45 |
_RL dTarg |
_RL dTarg |
46 |
_RL rTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL rTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
49 |
INTEGER myThid |
INTEGER myThid |
50 |
|
|
51 |
C !OUTPUT PARAMETERS: ================================================== |
C !OUTPUT PARAMETERS: ================================================== |
52 |
C wT :: vertical advective flux |
C wT :: vertical advective flux |
53 |
_RL wT (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL wT (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
54 |
|
|
55 |
C !LOCAL VARIABLES: ==================================================== |
C !LOCAL VARIABLES: ==================================================== |
56 |
C i,j :: loop indices |
C i,j :: loop indices |
57 |
C kp1 :: =min( k+1 , Nr ) |
C kp1 :: =min( k+1 , Nr ) |
58 |
C km1 :: =max( k-1 , 1 ) |
C km1 :: =max( k-1 , 1 ) |
59 |
C km2 :: =max( k-2 , 1 ) |
C km2 :: =max( k-2 , 1 ) |
60 |
C bi,bj :: tile indices or (1,1) depending on use |
C bi,bj :: tile indices or (1,1) depending on use |
61 |
C Cr :: slope ratio |
C Cr :: slope ratio |
62 |
C Rjm,Rj,Rjp :: differences at i-1,i,i+1 |
C Rjm,Rj,Rjp :: differences at i-1,i,i+1 |
63 |
C wFld :: velocity, vertical component |
C wFld :: velocity, vertical component |
64 |
INTEGER i,j,kp1,km1,km2,bi,bj |
INTEGER i,j,kp1,km1,km2,bi,bj |
65 |
_RL Cr,Rjm,Rj,Rjp |
_RL Cr,Rjm,Rj,Rjp |
66 |
_RL wFld |
_RL wFld |