Small question about the using UpdateArg1 function

[ruohai0925]

Small question about the using UpdateArg1 function

[ajnonaka]

Maybe there is a factor of dt hidden somewhere in the code I haven't spotted yet, but I think I agree with you. This is the initial sync projection (composite over all levels) to compute the initial pressure. The quantity we are projecting is Vproj=(unew-uold)/dt. Thus the pressure increment should be pnew = pold + phi, where phi is the solution to the elliptic problem. @asalmgren do you agree? I am going offline for a few hours but can do some numerical tests to confirm later. Basically what you want to do it set "ns.init_iter" to a very large number and see if the pressure converges to something. Try it with "1/dt" and "1" in the UpdateArg1 weighting and see if anything interesting happens.

[WeiqunZhang]

dt = 1

[ruohai0925]

Aaha, I really apologize for asking this stupid question. But indeed I have not found it before. The dt is an input parameter, which equals to 1. I should be more careful next time.

Thanks, everyone.

[ajnonaka]

Doh! The use of "dt" as the input parameter name is a bit misleading. I'll change it to "factor".