NaN values encountered when trying to setup user-defined auxiliary variables

First of all, thanks a lot for developing Castro ; the code is very useful to me!

I recently tried to use auxiliary variables (described in this page of the documentation, accessed with UFX, etc.). More specifically, I did a small test where I started from the 1D Sedov example and added an auxiliary variable that is supposed to mirror exactly eint_e(as a way to check that the variable behaves as expected): see changes in the commit below
https://github.com/RemiLehe/Castro/commit/ccfd30fb9d6f236c493bc771305b8edc1fe7caea

However, when running this example, I get an error message:

amrex::Abort::0::State has NaNs in the rho_myenerg component::check_for_nan() !!!

Is this expected? Is there a way to work around this?

How to easily reproduce this issue

Clone my fork from Castro and switch to the branch test_aux_variables (this branch is identical to the current main branch of Castro, but with the additional commit listed above). Then run the 1d Sedov cylindrical test.

git clone https://github.com/RemiLehe/Castro
cd Castro
git checkout test_aux_variables
git submodule update --init
cd Exec/hydro_tests/Sedov
make DIM=1 -j 4
./Castro1d.gnu.MPI.ex inputs.1d.cyl