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

Diff of /MITgcm/pkg/generic_advdiff/gad_dst3fl_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.2 by heimbach, Thu Nov 8 23:39:34 2001 UTC
# Line 64  C      for maskC(...) and wVel(...) Line 64  C      for maskC(...) and wVel(...)
64       &         *maskC(i,j,km1,bi,bj)       &         *maskC(i,j,km1,bi,bj)
65    
66          cfl=abs(wVel(i,j,k,bi,bj)*dTarg*recip_drc(k))          cfl=abs(wVel(i,j,k,bi,bj)*dTarg*recip_drc(k))
67          d0=(2.-cfl)*(1.-cfl)*oneSixth          d0=(2.D0-cfl)*(1.-cfl)*oneSixth
68          d1=(1.-cfl*cfl)*oneSixth          d1=(1.D0-cfl*cfl)*oneSixth
69  c       thetaP=0.  c       thetaP=0.D0
70  c       IF (Rj.NE.0.) thetaP=Rjm/Rj  c       IF (Rj.NE.0.D0) thetaP=Rjm/Rj
71          thetaP=Rjm/(1.D-20+Rj)          thetaP=Rjm/(1.D-20+Rj)
72          psiP=d0+d1*thetaP          psiP=d0+d1*thetaP
73          psiP=max(0.,min(min(1.,psiP),(1.-cfl)/(1.D-20+cfl)*thetaP))          psiP=max(0.D0,min(min(1.D0,psiP),
74         &       (1.D0-cfl)/(1.D-20+cfl)*thetaP))
75          thetaM=Rjp/(1.D-20+Rj)          thetaM=Rjp/(1.D-20+Rj)
76  c       thetaM=0.  c       thetaM=0.D0
77  c       IF (Rj.NE.0.) thetaM=Rjp/Rj  c       IF (Rj.NE.0.D0) thetaM=Rjp/Rj
78          psiM=d0+d1*thetaM          psiM=d0+d1*thetaM
79          psiM=max(0.,min(min(1.,psiM),(1.-cfl)/(1.D-20+cfl)*thetaM))          psiM=max(0.D0,min(min(1.D0,psiM),
80         &       (1.D0-cfl)/(1.D-20+cfl)*thetaM))
81          wT(i,j)=          wT(i,j)=
82       &   0.5*(rTrans(i,j)+abs(rTrans(i,j)))       &   0.5*(rTrans(i,j)+abs(rTrans(i,j)))
83       &      *( Tracer(i,j, k ,bi,bj) + psiM*Rj )       &      *( Tracer(i,j, k ,bi,bj) + psiM*Rj )
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22