/[MITgcm]/MITgcm/pkg/generic_advdiff/gad_fluxlimit_adv_y.F
ViewVC logotype

Diff of /MITgcm/pkg/generic_advdiff/gad_fluxlimit_adv_y.F

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

revision 1.5 by edhill, Mon Mar 29 03:33:51 2004 UTC revision 1.6 by jmc, Fri Sep 24 16:53:46 2004 UTC
# Line 10  C !INTERFACE: ========================== Line 10  C !INTERFACE: ==========================
10        SUBROUTINE GAD_FLUXLIMIT_ADV_Y(        SUBROUTINE GAD_FLUXLIMIT_ADV_Y(
11       I           bi,bj,k,deltaT,       I           bi,bj,k,deltaT,
12       I           vTrans, vVel,       I           vTrans, vVel,
13       I           tracer,       I           maskLocS, tracer,
14       O           vT,       O           vT,
15       I           myThid )       I           myThid )
16    
# Line 43  C  myThid            :: thread number Line 43  C  myThid            :: thread number
43        _RL deltaT        _RL deltaT
44        _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
45        _RL vVel  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL vVel  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
46          _RS maskLocS(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
47        _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
48        INTEGER myThid        INTEGER myThid
49    
# Line 73  CEOP Line 74  CEOP
74  c       vFld = vVel(i,j,k,bi,bj)  c       vFld = vVel(i,j,k,bi,bj)
75          vFld = vTrans(i,j)*recip_dxG(i,j,bi,bj)          vFld = vTrans(i,j)*recip_dxG(i,j,bi,bj)
76       &       *recip_drF(k)*recip_hFacS(i,j,k,bi,bj)       &       *recip_drF(k)*recip_hFacS(i,j,k,bi,bj)
77          Rjp=(tracer(i,j+1)-tracer(i,j))*maskS(i,j+1,k,bi,bj)          Rjp=(tracer(i,j+1)-tracer(i, j ))*maskLocS(i,j+1)
78          Rj=(tracer(i,j)-tracer(i,j-1))*maskS(i,j,k,bi,bj)          Rj =(tracer(i, j )-tracer(i,j-1))*maskLocS(i, j )
79          Rjm=(tracer(i,j-1)-tracer(i,j-2))*maskS(i,j-1,k,bi,bj)          Rjm=(tracer(i,j-1)-tracer(i,j-2))*maskLocS(i,j-1)
80    
81          IF (Rj.NE.0.) THEN          IF (Rj.NE.0.) THEN
82           IF (vTrans(i,j).GT.0) THEN           IF (vTrans(i,j).GT.0) THEN
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22