Skip to content

Commit

Permalink
Fix source indexing when loading advection velocity for post-processing
Browse files Browse the repository at this point in the history
  • Loading branch information
@johnomotani
johnomotani committed Sep 20, 2024
1 parent 6c99965 commit ece2a94
Showing 1 changed file with 12 additions and 12 deletions.
24 changes: 12 additions & 12 deletions moment_kinetics/src/load_data.jl
View file
Original file line number Diff line number Diff line change
Expand Up @@ -4406,10 +4406,10 @@ function get_variable(run_info, variable_name; normalize_advection_speed_shape=t
dvth_dz=dvth_dz[:,:,:,it],
dqpar_dz=dqpar_dz[:,:,:,it],
vth=vth[:,:,:,it],
external_source_amplitude=external_source_amplitude[:,:,n_sources,it],
external_source_density_amplitude=external_source_density_amplitude[:,:,n_sources,it],
external_source_momentum_amplitude=external_source_momentum_amplitude[:,:,n_sources,it],
external_source_pressure_amplitude=external_source_pressure_amplitude[:,:,n_sources,it]),
external_source_amplitude=external_source_amplitude[:,:,:,it],
external_source_density_amplitude=external_source_density_amplitude[:,:,:,it],
external_source_momentum_amplitude=external_source_momentum_amplitude[:,:,:,it],
external_source_pressure_amplitude=external_source_pressure_amplitude[:,:,:,it]),
evolve_density=run_info.evolve_density,
evolve_upar=run_info.evolve_upar,
evolve_ppar=run_info.evolve_ppar)
Expand Down Expand Up @@ -4527,10 +4527,10 @@ function get_variable(run_info, variable_name; normalize_advection_speed_shape=t
dppar_dz=dppar_dz[:,:,it],
dqpar_dz=dqpar_dz[:,:,it],
dvth_dz=dvth_dz[:,:,it],
external_source_amplitude=external_source_amplitude[:,:,n_sources,it],
external_source_density_amplitude=external_source_density_amplitude[:,:,n_sources,it],
external_source_momentum_amplitude=external_source_momentum_amplitude[:,:,n_sources,it],
external_source_pressure_amplitude=external_source_pressure_amplitude[:,:,n_sources,it]),)
external_source_amplitude=external_source_amplitude[:,:,:,it],
external_source_density_amplitude=external_source_density_amplitude[:,:,:,it],
external_source_momentum_amplitude=external_source_momentum_amplitude[:,:,:,it],
external_source_pressure_amplitude=external_source_pressure_amplitude[:,:,:,it]),)
@views update_electron_speed_vpa!(advect, density[:,:,it], upar[:,:,it],
ppar[:,:,it], moments, run_info.vpa.grid,
run_info.external_source_settings.electron)
Expand Down Expand Up @@ -4665,10 +4665,10 @@ function get_variable(run_info, variable_name; normalize_advection_speed_shape=t
dvth_dz=dvth_dz[:,:,:,it],
dqz_dz=dqz_dz[:,:,:,it],
vth=vth[:,:,:,it],
external_source_amplitude=external_source_amplitude[:,:,n_sources,it],
external_source_density_amplitude=external_source_density_amplitude[:,:,n_sources,it],
external_source_momentum_amplitude=external_source_momentum_amplitude[:,:,n_sources,it],
external_source_pressure_amplitude=external_source_pressure_amplitude[:,:,n_sources,it]),
external_source_amplitude=external_source_amplitude[:,:,:,it],
external_source_density_amplitude=external_source_density_amplitude[:,:,:,it],
external_source_momentum_amplitude=external_source_momentum_amplitude[:,:,:,it],
external_source_pressure_amplitude=external_source_pressure_amplitude[:,:,:,it]),
evolve_density=run_info.evolve_density,
evolve_upar=run_info.evolve_upar,
evolve_ppar=run_info.evolve_ppar)
Expand Down

0 comments on commit ece2a94

Please sign in to comment.