From 5c66a1fff0b6c33addfdca5c2abdc9cc9ffb712c Mon Sep 17 00:00:00 2001 From: Orso Meneghini Date: Thu, 23 Jan 2025 23:54:51 -0800 Subject: [PATCH] init trims data before the first equilibrium --- src/actors/compound/dynamic_plasma_actor.jl | 4 ++-- src/ddinit/init.jl | 3 +++ 2 files changed, 5 insertions(+), 2 deletions(-) diff --git a/src/actors/compound/dynamic_plasma_actor.jl b/src/actors/compound/dynamic_plasma_actor.jl index c877c35c3..76cc021b0 100644 --- a/src/actors/compound/dynamic_plasma_actor.jl +++ b/src/actors/compound/dynamic_plasma_actor.jl @@ -131,8 +131,8 @@ function _step(actor::ActorDynamicPlasma; n_steps::Int=0) prog = ProgressMeter.Progress(Nt * step_calls_per_2loop; dt=0.0, showspeed=true, enabled=par.verbose) old_logging = actor_logging(dd, false) - # remove timeslices after global_time - IMAS.trim_time!(dd, t0; trim_pulse_schedule=false) + # remove time dependent data after global_time + IMAS.trim_time!(dd, (-Inf, t0); trim_pulse_schedule=false) try for (kk, tt) in enumerate(range(t0, t1, 2 * Nt + 1)[2:end]) diff --git a/src/ddinit/init.jl b/src/ddinit/init.jl index 49da2f6ce..90e1a0294 100644 --- a/src/ddinit/init.jl +++ b/src/ddinit/init.jl @@ -226,6 +226,9 @@ function init( end end + # trim any data before the first equilibrium since things are really not robust against that + IMAS.trim_time!(dd, (dd.equilibrium.time[1], Inf)); + return dd end end