User physical src term cannot be enrolled when restarting about athena

[tomidakn]

Hi Wendy,

Please make sure you have the latest version. EnrollUserSourceTermFunction is a
member of Mesh, not MeshBlock. So I believe this should work as intended, and
you should be able to use the InitUserMeshProperties function for restarting
initialization.

Cheers,
Kengo

On 2016/02/10 13:11, Wenhua Ju wrote:

With the new boundary condition interface, I see that the new
�$B!H�(BMesh::InitUserMeshProperties�$B!I�(B function is added to problem generators which is
called when we start a new run or restart. This is nice, since user boundary
conditions in this function can be enrolled when we restart.

However, in this case, user physical source term still cannot be enrolled when
we restart. The �$B!H�(BEnrollSrcTermFunction�$B!I�(B is under MeshBlock. It�$B!G�(Bs not a property
for the whole mesh. If we call this function in �$B!H�(BMesh::InitUserMeshProperties�$B!I�(B,
we have to iterate over every meshblock. But if we call this function in
�$B!H�(BProblemGenerator�$B!I�(B, it is not called when we restart a run.

Before the new BC interface was written, my old way to solve this problem was
adding an extra function of "RestartInitialization" to the problem generator
where I enroll user src term function (like athena 4.2). I let the main function
call this "RestartInitialization" function if res_flag==1. This way of hacking
is not good for future users.

�$B!=�(B
Reply to this email directly or view it on GitHub
#18.

Kengo TOMIDA [email protected]
Department of Astrophysical Sciences, Princeton University

[wenhuaju]

Hi Kengo,

Thanks! Sorry I missed your update a few weeks ago. Now everything works.

Best,
Wendy