diff options
| -rw-r--r-- | src/multipole.cc | 17 |
1 files changed, 10 insertions, 7 deletions
diff --git a/src/multipole.cc b/src/multipole.cc index c79653a..1a2881a 100644 --- a/src/multipole.cc +++ b/src/multipole.cc @@ -133,17 +133,20 @@ static void output_1D(CCTK_ARGUMENTS, const variable_desc *v, CCTK_REAL rad, if (out_1d_every != 0 && (cctk_iteration) % out_1d_every == 0) { const char *real_name = CCTK_VarName(v->index); - const char *imag_name = CCTK_VarName(v->imag_index); - sprintf(name_tmp, "mp_%s_r%1.2f.th.asc", real_name, rad); Multipole_Output1D(CCTK_PASS_CTOC, name_tmp, array_size, th, ph, mp_theta, real); - sprintf(name_tmp, "mp_%s_r%1.2f.th.asc", imag_name, rad); - Multipole_Output1D(CCTK_PASS_CTOC, name_tmp, array_size, th, ph, mp_theta, imag); - sprintf(name_tmp, "mp_%s_r%1.2f.ph.asc", real_name, rad); Multipole_Output1D(CCTK_PASS_CTOC, name_tmp, array_size, th, ph, mp_phi, real); - sprintf(name_tmp, "mp_%s_r%1.2f.ph.asc", imag_name, rad); - Multipole_Output1D(CCTK_PASS_CTOC, name_tmp, array_size, th, ph, mp_phi, imag); + + if (v->imag_index != -1) + { + const char *imag_name = CCTK_VarName(v->imag_index); + sprintf(name_tmp, "mp_%s_r%1.2f.th.asc", imag_name, rad); + Multipole_Output1D(CCTK_PASS_CTOC, name_tmp, array_size, th, ph, mp_theta, imag); + + sprintf(name_tmp, "mp_%s_r%1.2f.ph.asc", imag_name, rad); + Multipole_Output1D(CCTK_PASS_CTOC, name_tmp, array_size, th, ph, mp_phi, imag); + } } } } |