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

Diff of /MITgcm/pkg/generic_advdiff/gad_dst3_adv_x.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 53  C     == Local variables == Line 53  C     == Local variables ==
53  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
54  C  i,j               :: loop indices  C  i,j               :: loop indices
55  C  uLoc              :: velocity [m/s], zonal component  C  uLoc              :: velocity [m/s], zonal component
56  C  cfl               :: Courant-Friedrich-Levy number  C  uCFL              :: Courant-Friedrich-Levy number
57        INTEGER i,j        INTEGER i,j
58        _RL uLoc        _RL uLoc
59        _RL Rjm,Rj,Rjp,cfl,d0,d1        _RL Rjm,Rj,Rjp,uCFL,d0,d1
60  #ifdef OLD_DST3_FORMULATION  #ifdef OLD_DST3_FORMULATION
61        _RL psiP,psiM,thetaP,thetaM        _RL psiP,psiM,thetaP,thetaM
62        _RL smallNo        _RL smallNo
# Line 79  c     _RL Rjjm,Rjjp Line 79  c     _RL Rjjm,Rjjp
79          Rjm=(tracer(i-1,j)-tracer(i-2,j))*maskLocW(i-1,j)          Rjm=(tracer(i-1,j)-tracer(i-2,j))*maskLocW(i-1,j)
80    
81          uLoc = uFld(i,j)          uLoc = uFld(i,j)
82  c       uLoc = uTrans(i,j)*recip_dyG(i,j,bi,bj)          uCFL = ABS( uLoc*deltaTloc
83  c    &       *recip_drF(k)*_recip_hFacW(i,j,k,bi,bj)       &                  *recip_dxC(i,j,bi,bj)*recip_deepFacC(k) )
84          cfl=ABS(uLoc*deltaTloc*recip_dxC(i,j,bi,bj))          d0=(2.-uCFL)*(1.-uCFL)*oneSixth
85          d0=(2.-cfl)*(1.-cfl)*oneSixth          d1=(1.-uCFL*uCFL)*oneSixth
         d1=(1.-cfl*cfl)*oneSixth  
86  #ifdef OLD_DST3_FORMULATION  #ifdef OLD_DST3_FORMULATION
87          IF ( ABS(Rj).LT.smallNo .OR.          IF ( ABS(Rj).LT.smallNo .OR.
88       &       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