32 |
_RL mask_Kp1,mask_Km1,wBarYm,wBarYp |
_RL mask_Kp1,mask_Km1,wBarYm,wBarYp |
33 |
_RL vZm,vZp |
_RL vZm,vZp |
34 |
LOGICAL rAdvAreaWeight |
LOGICAL rAdvAreaWeight |
|
PARAMETER( rAdvAreaWeight =.TRUE.) |
|
35 |
c _RL vmask_Kp1,vmask_K,vmask_Km1 |
c _RL vmask_Kp1,vmask_K,vmask_Km1 |
36 |
c LOGICAL freeslipK,noslipK |
c LOGICAL freeslipK,noslipK |
37 |
c PARAMETER(freeslipK=.TRUE.) |
c PARAMETER(freeslipK=.TRUE.) |
41 |
c PARAMETER(noslip1=.NOT.freeslip1) |
c PARAMETER(noslip1=.NOT.freeslip1) |
42 |
c1 _RL wBarYZ,vZbarZ |
c1 _RL wBarYZ,vZbarZ |
43 |
|
|
44 |
|
rAdvAreaWeight =.TRUE. |
45 |
|
C- Area-weighted average either in KE or in vert. advection: |
46 |
|
IF ( selectKEscheme.EQ.1 .OR. selectKEscheme.EQ.3 ) |
47 |
|
& rAdvAreaWeight =.FALSE. |
48 |
|
|
49 |
Kp1=min(K+1,Nr) |
Kp1=min(K+1,Nr) |
50 |
mask_Kp1=1. |
mask_Kp1=1. |
51 |
IF (K.EQ.Nr) mask_Kp1=0. |
IF (K.EQ.Nr) mask_Kp1=0. |