55 |
psiM=d0+d1*thetaM |
psiM=d0+d1*thetaM |
56 |
psiM=max(0., min(min(1.,psiM),(1.-cfl)/cfl*thetaM) ) |
psiM=max(0., min(min(1.,psiM),(1.-cfl)/cfl*thetaM) ) |
57 |
vT(i,j)= |
vT(i,j)= |
|
c & 0.5*(vTrans(i,j)+abs(vTrans(i,j))) |
|
|
c & *( Tracer(i,j-1) + d0*Rj + d1*Rjm ) |
|
|
c & +0.5*(vTrans(i,j)-abs(vTrans(i,j))) |
|
|
c & *( Tracer(i, j ) - d0*Rj + d1*Rjp ) |
|
58 |
& 0.5*(vTrans(i,j)+abs(vTrans(i,j))) |
& 0.5*(vTrans(i,j)+abs(vTrans(i,j))) |
59 |
& *( Tracer(i,j-1) + psiP*Rj ) |
& *( Tracer(i,j-1) + psiP*Rj ) |
60 |
& +0.5*(vTrans(i,j)-abs(vTrans(i,j))) |
& +0.5*(vTrans(i,j)-abs(vTrans(i,j))) |
61 |
& *( Tracer(i, j ) - psiM*Rj ) |
& *( Tracer(i, j ) - psiM*Rj ) |
62 |
|
cph *note* put these comments after end of continued line |
63 |
|
cph to ensure TAMC compatibility |
64 |
|
c & 0.5*(vTrans(i,j)+abs(vTrans(i,j))) |
65 |
|
c & *( Tracer(i,j-1) + d0*Rj + d1*Rjm ) |
66 |
|
c & +0.5*(vTrans(i,j)-abs(vTrans(i,j))) |
67 |
|
c & *( Tracer(i, j ) - d0*Rj + d1*Rjp ) |
68 |
|
|
69 |
ENDDO |
ENDDO |
70 |
ENDDO |
ENDDO |