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

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

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

revision 1.5 by jmc, Wed Mar 6 01:29:36 2002 UTC revision 1.6 by jmc, Fri Sep 24 16:53:46 2004 UTC
# Line 6  C $Name$ Line 6  C $Name$
6        SUBROUTINE GAD_DST3FL_ADV_Y(        SUBROUTINE GAD_DST3FL_ADV_Y(
7       I           bi,bj,k,deltaT,       I           bi,bj,k,deltaT,
8       I           vTrans, vVel,       I           vTrans, vVel,
9       I           tracer,       I           maskLocS, tracer,
10       O           vT,       O           vT,
11       I           myThid )       I           myThid )
12  C     /==========================================================\  C     /==========================================================\
# Line 26  C     == Routine arguments == Line 26  C     == Routine arguments ==
26        _RL deltaT        _RL deltaT
27        _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
28        _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)
29          _RS maskLocS(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
30        _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
31        _RL vT    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL vT    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
32        INTEGER myThid        INTEGER myThid
# Line 43  C     vFld    :: velocity [m/s], meridio Line 44  C     vFld    :: velocity [m/s], meridio
44        ENDDO        ENDDO
45        DO j=1-Oly+2,sNy+Oly-1        DO j=1-Oly+2,sNy+Oly-1
46         DO i=1-Olx,sNx+Olx         DO i=1-Olx,sNx+Olx
47          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)
48          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 )
49          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)
50    
51  c       vFld = vVel(i,j,k,bi,bj)  c       vFld = vVel(i,j,k,bi,bj)
52          vFld = vTrans(i,j)*recip_dxG(i,j,bi,bj)          vFld = vTrans(i,j)*recip_dxG(i,j,bi,bj)
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22