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 |