diff --git a/src/pe_coupling/momentum_exchange_method/restoration/Reconstructor.h b/src/pe_coupling/momentum_exchange_method/restoration/Reconstructor.h
index 7f9ce1017f6cc13c703ee18711c792b1d07f8c85..b8fb7affbfc2bb720b7eb0930419d221ad5343c2 100644
--- a/src/pe_coupling/momentum_exchange_method/restoration/Reconstructor.h
+++ b/src/pe_coupling/momentum_exchange_method/restoration/Reconstructor.h
@@ -120,6 +120,8 @@ void EquilibriumReconstructor< LatticeModel_T, BoundaryHandling_T >
 
    real_t cx, cy, cz;
    blockStorage_->getBlockLocalCellCenter( *block, Cell(x,y,z), cx, cy, cz );
+
+   WALBERLA_ASSERT_NOT_NULLPTR( (*bodyField)(x,y,z) );
    const auto velocity = (*bodyField)(x,y,z)->velFromWF(cx,cy,cz);
 
    pdfField->setToEquilibrium( x, y, z, Vector3< real_t >( velocity[0], velocity[1], velocity[2] ), averageDensity );