Use correct size for BGField that is allocated with size 1:BGDataSize

piclas-framework · Mar 7, 2022 · 3c911f0 · 3c911f0
1 parent 865c5cc
commit 3c911f0
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/particles/pic/interpolation/pic_interpolation_tools.f90 b/src/particles/pic/interpolation/pic_interpolation_tools.f90
@@ -281,9 +281,9 @@ PPURE FUNCTION GetEMFieldDW(ElemID, PartPos_loc)
 ! Evaluate the electro-(magnetic) field using the reference position and return the field
 !===================================================================================================================================
 ! MODULES
-USE MOD_Mesh_Vars     ,ONLY: Elem_xGP
+USE MOD_Mesh_Vars             ,ONLY: Elem_xGP
 USE MOD_PICInterpolation_Vars ,ONLY: useBGField
-USE MOD_Interpolation_Vars    ,ONLY: BGField,BGType
+USE MOD_Interpolation_Vars    ,ONLY: BGField,BGType,BGDataSize
 USE MOD_Globals
 USE MOD_PreProc
 #if ! (USE_HDG)
@@ -384,7 +384,7 @@ PPURE FUNCTION GetEMFieldDW(ElemID, PartPos_loc)
   END SELECT
   ! Add contribution of the magnetic field
   DO k = 0, PP_N; DO l=0, PP_N; DO m=0, PP_N
-    GetEMFieldDW(ind1:ind2) = GetEMFieldDW(ind1:ind2) + PartDistDepo(k,l,m)/DistSum*BGField(ind1:ind2,k,l,m,ElemID)
+    GetEMFieldDW(ind1:ind2) = GetEMFieldDW(ind1:ind2) + PartDistDepo(k,l,m)/DistSum*BGField(1:BGDataSize,k,l,m,ElemID)
   END DO; END DO; END DO
 END IF ! useBGField