From 0cfdb087e9225c5e26d4fe064b0a932d8af7b29b Mon Sep 17 00:00:00 2001 From: mlietzow Date: Thu, 30 Jan 2025 16:09:48 +0100 Subject: [PATCH] added missing bracket and fixed wrong stokes parameter --- src/DustComponent.cpp | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/DustComponent.cpp b/src/DustComponent.cpp index 9909b4a..8ce06cd 100644 --- a/src/DustComponent.cpp +++ b/src/DustComponent.cpp @@ -4932,9 +4932,6 @@ StokesVector CDustComponent::calcEmissivityEmi(CGridBasic * grid, { if(sizeIndexUsed(a, a_min, a_max)) { - // Get index of theta scattering - uint thID = phID == PH_MIE ? getScatThetaID(scattering_theta,a,w) : 0; - // Get cross sections and relative weight of the current dust grain size calcCrossSections(grid, pp, i_density, a, mag_field_theta, cs); @@ -4983,7 +4980,7 @@ StokesVector CDustComponent::calcEmissivityEmi(CGridBasic * grid, #if BENCHMARK == CAMPS // To perform Camps et. al (2015) benchmark. - tmp_stokes[a].addQ(cs.Cabs * pl; + tmp_stokes[a].addI(cs.Cabs * pl); #else // Add relative emissivity from this temperature tmp_stokes[a].addI(cs.Cabs * pl); @@ -5005,6 +5002,9 @@ StokesVector CDustComponent::calcEmissivityEmi(CGridBasic * grid, if(phID == PH_MIE) { + // Get index of theta scattering + uint thID = getScatThetaID(scattering_theta, a, w); + // Get scattering matrix const Matrix2D & mat_sca = getScatteringMatrix(a, w, 0, 0, thID);