Skip to content

Commit

Permalink
update the nse_interp interface (#2650)
Browse files Browse the repository at this point in the history
this syncs up with the new Microphysics changes
  • Loading branch information
zingale authored Nov 10, 2023
1 parent a73a6ce commit edad6c8
Show file tree
Hide file tree
Showing 2 changed files with 8 additions and 2 deletions.
5 changes: 4 additions & 1 deletion Exec/reacting_tests/nse_test/problem_initialize_state_data.H
Original file line number Diff line number Diff line change
Expand Up @@ -67,8 +67,11 @@ void problem_initialize_state_data (int i, int j, int k,
Real abar;
Real dq;
Real dyedt;
Real dabardt;
Real dbeadt;
Real e_nu;

nse_interp(T, problem::rho0, ye, abar, dq, dyedt, burn_state.xn);
nse_interp(T, problem::rho0, ye, abar, dq, dyedt, dabardt, dbeadt, e_nu, burn_state.xn);
#elif defined(NSE_NET)
Real eps = 1.e-10_rt;
bool input_ye_is_valid = true;
Expand Down
5 changes: 4 additions & 1 deletion Exec/science/massive_star/problem_initialize_state_data.H
Original file line number Diff line number Diff line change
Expand Up @@ -82,10 +82,13 @@ void problem_initialize_state_data (int i, int j, int k,
Real abar;
Real dq;
Real dyedt;
Real dabardt;
Real dbeadt;
Real e_nu;
Real xn[NumSpec];

nse_interp(state(i,j,k,UTEMP), state(i,j,k,URHO),
state(i,j,k,UFX+AuxZero::iye), abar, dq, dyedt, xn);
state(i,j,k,UFX+AuxZero::iye), abar, dq, dyedt, dabardt, dbeadt, e_nu, xn);

state(i,j,k,UFX+AuxZero::iabar) = abar;
state(i,j,k,UFX+AuxZero::ibea) = dq;
Expand Down

0 comments on commit edad6c8

Please sign in to comment.