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

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

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

revision 1.1 by adcroft, Mon Sep 10 13:09:04 2001 UTC revision 1.1.6.2 by heimbach, Fri Mar 7 03:55:23 2003 UTC
# Line 33  C     == Routine arguments == Line 33  C     == Routine arguments ==
33        INTEGER myThid        INTEGER myThid
34    
35  C     == Local variables ==  C     == Local variables ==
36    C     wFld     :: velocity, vertical component
37        INTEGER i,j,kp1,km1,km2,bi,bj        INTEGER i,j,kp1,km1,km2,bi,bj
38        _RL Rjm,Rj,Rjp,cfl,d0,d1        _RL Rjm,Rj,Rjp,cfl,d0,d1
39        _RL psiP,psiM,thetaP,thetaM        _RL psiP,psiM,thetaP,thetaM
40          _RL wFld
41    
42        IF (.NOT. multiDimAdvection) THEN        IF (.NOT. multiDimAdvection) THEN
43  C      If using the standard time-stepping/advection schemes (ie. AB-II)  C      If using the standard time-stepping/advection schemes (ie. AB-II)
# Line 57  C      for maskC(...) and wVel(...) Line 59  C      for maskC(...) and wVel(...)
59        DO j=1-Oly,sNy+Oly        DO j=1-Oly,sNy+Oly
60         DO i=1-Olx,sNx+Olx         DO i=1-Olx,sNx+Olx
61          Rjp=(tracer(i,j,k,bi,bj)-tracer(i,j,kp1,bi,bj))          Rjp=(tracer(i,j,k,bi,bj)-tracer(i,j,kp1,bi,bj))
62       &         *maskC(i,j,kp1,bi,bj)       &         *maskC(i,j,kp1,bi_arg,bj_arg)
63          Rj =(tracer(i,j,km1,bi,bj)-tracer(i,j,k,bi,bj))          Rj =(tracer(i,j,km1,bi,bj)-tracer(i,j,k,bi,bj))
64       &         *maskC(i,j,k,bi,bj)*maskC(i,j,km1,bi,bj)       &         *maskC(i,j,k,bi_arg,bj_arg)*maskC(i,j,km1,bi_arg,bj_arg)
65          Rjm=(tracer(i,j,km2,bi,bj)-tracer(i,j,km1,bi,bj))          Rjm=(tracer(i,j,km2,bi,bj)-tracer(i,j,km1,bi,bj))
66       &         *maskC(i,j,km1,bi,bj)       &         *maskC(i,j,km1,bi_arg,bj_arg)
67    
68          cfl=abs(wVel(i,j,k,bi,bj)*dTarg*recip_drc(k))  c       wFld = wVel(i,j,k,bi_arg,bj_arg)
69            wFld = rTrans(i,j)*recip_rA(i,j,bi_arg,bj_arg)
70            cfl=abs(wFld*dTarg*recip_drC(k))
71          d0=(2.-cfl)*(1.-cfl)*oneSixth          d0=(2.-cfl)*(1.-cfl)*oneSixth
72          d1=(1.-cfl*cfl)*oneSixth          d1=(1.-cfl*cfl)*oneSixth
73  c       thetaP=0.  c       thetaP=0.
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.1.6.2
  ViewVC Help
Powered by ViewVC 1.1.22