diff --git a/Exec/ParticleFilterTest/DarkMatterParticleContainer.cpp b/Exec/ParticleFilterTest/DarkMatterParticleContainer.cpp index b610a5b1..4e655eb8 100644 --- a/Exec/ParticleFilterTest/DarkMatterParticleContainer.cpp +++ b/Exec/ParticleFilterTest/DarkMatterParticleContainer.cpp @@ -122,7 +122,7 @@ struct ShellFilter ylen = src.m_aos[i].rdata(1+1+3)+(jdir)*(m_phi[1]-m_plo[1]) - m_center[1]; zlen = src.m_aos[i].rdata(2+1+3)+(kdir)*(m_phi[2]-m_plo[2]) - m_center[2]; Real mag = sqrt(xlen*xlen+ylen*ylen+zlen*zlen); - result+=int(mag>m_radius_inner && mag -1000000.0 and zlen < 1000000.0); + result+=int(mag>m_radius_inner && magm_radius_inner && mag -1000000.0 and zlen < 1000000.0)) + if(int(mag>m_radius_inner && mag0&&radius_outer>radius_inner) { int write_hdf5=0; #ifdef AMREX_USE_HDF5