[ET Trac] [Einstein Toolkit] #1841: GRHydro tracers broken
Einstein Toolkit
trac-noreply at einsteintoolkit.org
Mon Nov 23 10:20:55 CST 2015
#1841: GRHydro tracers broken
------------------------------------+---------------------------------------
Reporter: I.Hawke@… | Owner:
Type: defect | Status: new
Priority: unset | Milestone:
Component: EinsteinToolkit thorn | Version: development version
Resolution: | Keywords: GRHydro
------------------------------------+---------------------------------------
Comment (by I.Hawke@…):
This is not the only problem with tracers, as I've since found.
Currently all reconstruction of tracers is done on the primitive
variables, but the conserved are never updated. This can be fixed by
adding
if (evolve_tracer .ne. 0) then
call Prim2ConservativeTracer(CCTK_PASS_FTOF)
end if
at the end of GRHydro_Reconstruct.F90, but *only* if the use_cxx_code
parameter is set to no.
To truly fix tracers, the C++ version GRHydro_Reconstruct.cc will also
need a P2C_Tracer routine, which needs calling at the end.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/1841#comment:1>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
More information about the Trac
mailing list