For the LwSrc, in step3d_t.F near lines 1146 we need to change
cff=dt(ng)*pm(i,j)*pn(i,j)
to
cff=dt(ng)*pm(Isrc,Jsrc)*pn(Isrc,Jsrc)
and the same near lines 1343 from
cff=dt(ng)*pm(i,j)*pn(i,j)
to
cff=dt(ng)*pm(Isrc,Jsrc)*pn(Isrc,Jsrc)
-j
LwSrc i,j bug
- arango
- Site Admin
- Posts: 1361
- Joined: Wed Feb 26, 2003 4:41 pm
- Location: DMCS, Rutgers University
- Contact:
Re: LwSrc i,j bug
Thanks, John, great catch. Weird, my aggressive compiling/linking flags didn't pick up those out-of-bounce indices in parallel. I will update the code.