/[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.1 by heimbach, Wed Feb 6 15:48:08 2002 UTC
# Line 57  C      for maskC(...) and wVel(...) Line 57  C      for maskC(...) and wVel(...)
57        DO j=1-Oly,sNy+Oly        DO j=1-Oly,sNy+Oly
58         DO i=1-Olx,sNx+Olx         DO i=1-Olx,sNx+Olx
59          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))
60       &         *maskC(i,j,kp1,bi,bj)       &         *maskC(i,j,kp1,bi_arg,bj_arg)
61          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))
62       &         *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)
63          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))
64       &         *maskC(i,j,km1,bi,bj)       &         *maskC(i,j,km1,bi_arg,bj_arg)
65    
66          cfl=abs(wVel(i,j,k,bi,bj)*dTarg*recip_drc(k))          cfl=abs(wVel(i,j,k,bi_arg,bj_arg)*dTarg*recip_drc(k))
67          d0=(2.-cfl)*(1.-cfl)*oneSixth          d0=(2.-cfl)*(1.-cfl)*oneSixth
68          d1=(1.-cfl*cfl)*oneSixth          d1=(1.-cfl*cfl)*oneSixth
69  c       thetaP=0.  c       thetaP=0.
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.1.6.1
  ViewVC Help
Powered by ViewVC 1.1.22