32 |
_RL mask_Kp1,mask_Km1,wBarXm,wBarXp |
_RL mask_Kp1,mask_Km1,wBarXm,wBarXp |
33 |
_RL uZm,uZp |
_RL uZm,uZp |
34 |
LOGICAL rAdvAreaWeight |
LOGICAL rAdvAreaWeight |
|
PARAMETER( rAdvAreaWeight =.TRUE.) |
|
35 |
c _RL umask_Kp1,umask_K,umask_Km1 |
c _RL umask_Kp1,umask_K,umask_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 wBarXZ,uZbarZ |
c1 _RL wBarXZ,uZbarZ |
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. |