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

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

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

revision 1.5 by jmc, Fri Sep 24 16:53:45 2004 UTC revision 1.6 by heimbach, Fri Aug 19 22:19:35 2005 UTC
# Line 8  C !ROUTINE: GAD_FLUXLIMIT_ADV_X Line 8  C !ROUTINE: GAD_FLUXLIMIT_ADV_X
8    
9  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
10        SUBROUTINE GAD_FLUXLIMIT_ADV_X(        SUBROUTINE GAD_FLUXLIMIT_ADV_X(
11       I           bi,bj,k,deltaT,       I           bi,bj,k,deltaTloc,
12       I           uTrans, uVel,       I           uTrans, uVel,
13       I           maskLocW, tracer,       I           maskLocW, tracer,
14       O           uT,       O           uT,
# Line 40  C  uVel                 :: zonal flow Line 40  C  uVel                 :: zonal flow
40  C  tracer               :: tracer field  C  tracer               :: tracer field
41  C  myThid               :: thread number  C  myThid               :: thread number
42        INTEGER bi,bj,k        INTEGER bi,bj,k
43        _RL deltaT        _RL deltaTloc
44        _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
45        _RL uVel  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL uVel  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
46        _RS maskLocW(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RS maskLocW(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
# Line 94  c       uFld = uVel(i,j,k,bi,bj) Line 94  c       uFld = uVel(i,j,k,bi,bj)
94       &   uTrans(i,j)*(Tracer(i,j)+Tracer(i-1,j))*0.5 _d 0       &   uTrans(i,j)*(Tracer(i,j)+Tracer(i-1,j))*0.5 _d 0
95       &   -0.5*(       &   -0.5*(
96       &        (1-Cr)*ABS(uTrans(i,j))       &        (1-Cr)*ABS(uTrans(i,j))
97       &        +uTrans(i,j)*uFld*deltaT       &        +uTrans(i,j)*uFld*deltaTloc
98       &         *recip_dxC(i,j,bi,bj)*Cr       &         *recip_dxC(i,j,bi,bj)*Cr
99       &        )*Rj       &        )*Rj
100         ENDDO         ENDDO
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22