/[MITgcm]/MITgcm/pkg/mom_vecinv/mom_vi_v_vertshear.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_vecinv/mom_vi_v_vertshear.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.8 by jmc, Tue Sep 27 23:42:06 2005 UTC revision 1.9 by jmc, Wed Sep 28 15:53:20 2005 UTC
# Line 32  C     == Local variables == Line 32  C     == Local variables ==
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.)
# Line 42  c     PARAMETER(freeslip1=.TRUE.) Line 41  c     PARAMETER(freeslip1=.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.
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22