/[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.1 by adcroft, Tue Sep 4 14:53:11 2001 UTC revision 1.2 by heimbach, Wed Sep 5 17:46:03 2001 UTC
# Line 53  C     == Local variables == Line 53  C     == Local variables ==
53          psiM=d0+d1*thetaM          psiM=d0+d1*thetaM
54          psiM=max(0., min(min(1.,psiM),(1.-cfl)/cfl*thetaM) )          psiM=max(0., min(min(1.,psiM),(1.-cfl)/cfl*thetaM) )
55          uT(i,j)=          uT(i,j)=
 c    &   0.5*(uTrans(i,j)+abs(uTrans(i,j)))  
 c    &      *( Tracer(i-1,j) + d0*Rj + d1*Rjm )  
 c    &  +0.5*(uTrans(i,j)-abs(uTrans(i,j)))  
 c    &      *( Tracer( i ,j) - d0*Rj + d1*Rjp )  
56       &   0.5*(uTrans(i,j)+abs(uTrans(i,j)))       &   0.5*(uTrans(i,j)+abs(uTrans(i,j)))
57       &      *( Tracer(i-1,j) + psiP*Rj )       &      *( Tracer(i-1,j) + psiP*Rj )
58       &  +0.5*(uTrans(i,j)-abs(uTrans(i,j)))       &  +0.5*(uTrans(i,j)-abs(uTrans(i,j)))
59       &      *( Tracer( i ,j) - psiM*Rj )       &      *( Tracer( i ,j) - psiM*Rj )
60    cph *note* put these comments after end of continued line
61    cph        to ensure TAMC compatibility
62    c    &   0.5*(uTrans(i,j)+abs(uTrans(i,j)))
63    c    &      *( Tracer(i-1,j) + d0*Rj + d1*Rjm )
64    c    &  +0.5*(uTrans(i,j)-abs(uTrans(i,j)))
65    c    &      *( Tracer( i ,j) - d0*Rj + d1*Rjp )
66    
67         ENDDO         ENDDO
68        ENDDO        ENDDO
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22