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

Diff of /MITgcm/pkg/generic_advdiff/gad_dst3fl_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 Jun 18 23:31:35 2006 UTC revision 1.11 by jmc, Mon Jun 19 14:40:43 2006 UTC
# Line 52  C       with no need to compute thetaM ( Line 52  C       with no need to compute thetaM (
52          Rj =(tracer( i ,j)-tracer(i-1,j))*maskLocW( i ,j)          Rj =(tracer( i ,j)-tracer(i-1,j))*maskLocW( i ,j)
53          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)
54    
55  c       uLoc = uFld(i,j)          uLoc = uFld(i,j)
56          uLoc = uTrans(i,j)*recip_dyG(i,j,bi,bj)  c       uLoc = uTrans(i,j)*recip_dyG(i,j,bi,bj)
57       &       *recip_drF(k)*_recip_hFacW(i,j,k,bi,bj)  c    &       *recip_drF(k)*_recip_hFacW(i,j,k,bi,bj)
58          cfl=abs(uLoc*deltaTloc*recip_dxC(i,j,bi,bj))          cfl=abs(uLoc*deltaTloc*recip_dxC(i,j,bi,bj))
59          d0=(2. _d 0 -cfl)*(1. _d 0 -cfl)*oneSixth          d0=(2. _d 0 -cfl)*(1. _d 0 -cfl)*oneSixth
60          d1=(1. _d 0 -cfl*cfl)*oneSixth          d1=(1. _d 0 -cfl*cfl)*oneSixth
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22