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

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

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

revision 1.10 by jmc, Sun Oct 22 01:08:04 2006 UTC revision 1.11 by jmc, Tue Dec 5 22:21:50 2006 UTC
# Line 52  C     == Local variables == Line 52  C     == Local variables ==
52  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
53  C  i,j               :: loop indices  C  i,j               :: loop indices
54  C  vLoc              :: velocity [m/s], meridional component  C  vLoc              :: velocity [m/s], meridional component
55  C  cfl               :: Courant-Friedrich-Levy number  C  vCFL              :: Courant-Friedrich-Levy number
56        INTEGER i,j        INTEGER i,j
57        _RL Rjm,Rj,Rjp,cfl,d0,d1        _RL Rjm,Rj,Rjp,vCFL,d0,d1
58        _RL vLoc        _RL vLoc
59  #ifdef OLD_DST3_FORMULATION  #ifdef OLD_DST3_FORMULATION
60        _RL psiP,psiM,thetaP,thetaM        _RL psiP,psiM,thetaP,thetaM
# Line 80  c     _RL Rjjm,Rjjp Line 80  c     _RL Rjjm,Rjjp
80          Rjm=(tracer(i,j-1)-tracer(i,j-2))*maskLocS(i,j-1)          Rjm=(tracer(i,j-1)-tracer(i,j-2))*maskLocS(i,j-1)
81    
82          vLoc = vFld(i,j)          vLoc = vFld(i,j)
83  c       vLoc = vTrans(i,j)*recip_dxG(i,j,bi,bj)          vCFL = ABS( vLoc*deltaTloc
84  c    &       *recip_drF(k)*_recip_hFacS(i,j,k,bi,bj)       &                  *recip_dyC(i,j,bi,bj)*recip_deepFacC(k) )
85          cfl=ABS(vLoc*deltaTloc*recip_dyC(i,j,bi,bj))          d0=(2.-vCFL)*(1.-vCFL)*oneSixth
86          d0=(2.-cfl)*(1.-cfl)*oneSixth          d1=(1.-vCFL*vCFL)*oneSixth
         d1=(1.-cfl*cfl)*oneSixth  
87  #ifdef OLD_DST3_FORMULATION  #ifdef OLD_DST3_FORMULATION
88          IF ( ABS(Rj).LT.smallNo .OR.          IF ( ABS(Rj).LT.smallNo .OR.
89       &       ABS(Rjm).LT.smallNo ) THEN       &       ABS(Rjm).LT.smallNo ) THEN
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22