I happened to look at the code in the t3dmix2_s.h and met something I couldnt understand.

Code: Select all

```
cff=0.25_r8*(diff3d_r(i,j,k)+diff3d_r(i-1,j,k))*pmon_u(i,j)
FX(i,j)=cff*(Hz(i,j,k)+Hz(i-1,j,k))*(t(i,j,k,nrhs,itrc)-t(i-1,j,k,nrhs,itrc))
similar for FE
cff=dt(ng)*pm(i,j)*pn(i,j)
cff1=cff*FX(i+1,j)-FX(i,j))
cff2=cff*FE(i,j+1)-FE(i,j))
cff3=cff1+cff2
t(i,j,k,nnew,itrc)=t(i,j,k,nnew,itrc)+cff3
```

Thanks a lot