diff --git a/src/include/ranges/mrange.cc.h b/src/include/ranges/mrange.cc.h
index 396d30b..5d62a51 100644
--- a/src/include/ranges/mrange.cc.h
+++ b/src/include/ranges/mrange.cc.h
@@ -6,6 +6,7 @@
 #include "operation/op_types.h"
 #include "operation/op_types.cc.h" // -> do .cc files
 #include "operation/op_utility.cc.h" // -> do .cc files
+#include "array/marray.h"
 
 namespace CNORXZ
 {
@@ -564,6 +565,13 @@ namespace CNORXZ
 	return mA[num];
     }
 
+    template <class... Ranges>
+    MArray<RangePtr> MRange<Ranges...>::sub() const
+    {
+	CXZ_ERROR("not implemented (SRange missing!!)"); //!!!
+	return MArray<RangePtr>();
+    }
+
     template <class... Ranges>
     SizeT MRange<Ranges...>::size() const
     {
diff --git a/src/include/ranges/mrange.h b/src/include/ranges/mrange.h
index a4b1f2e..18775f4 100644
--- a/src/include/ranges/mrange.h
+++ b/src/include/ranges/mrange.h
@@ -178,6 +178,7 @@ namespace CNORXZ
 	friend MRangeFactory<Ranges...>;
 
 	virtual RangePtr sub(SizeT num) const override final;
+	virtual MArray<RangePtr> sub() const override final;
 	virtual SizeT size() const override final;
 	virtual SizeT dim() const override final;
 	virtual String stringMeta(SizeT pos) const override final;
diff --git a/src/include/ranges/range_base.h b/src/include/ranges/range_base.h
index cc1a5f5..69a1803 100644
--- a/src/include/ranges/range_base.h
+++ b/src/include/ranges/range_base.h
@@ -52,6 +52,7 @@ namespace CNORXZ
 	
 	// virtual RangePtr sub() const; // Sptr<URange<RangePtr>> ; range of subranges (TODO!!!)
 	virtual RangePtr sub(SizeT num) const;
+	virtual MArray<RangePtr> sub() const;
 	virtual SizeT size() const = 0;
 	virtual SizeT dim() const = 0;
 	virtual const TypeInfo& type() const = 0;
@@ -120,6 +121,9 @@ namespace CNORXZ
     RangePack operator*(const RangePtr& a, const RangePack& b);
     RangePack operator*(const RangePack& a, const RangePtr& b);
     RangePack operator*(const RangePack& a, const RangePack& b);
+
+    RangePtr getSub(const RangePtr& r, SizeT num);
+    MArray<RangePtr> getSub(const RangePtr& r);
 }
 
 #endif
diff --git a/src/include/ranges/yrange.h b/src/include/ranges/yrange.h
index 2470b84..f674da4 100644
--- a/src/include/ranges/yrange.h
+++ b/src/include/ranges/yrange.h
@@ -120,6 +120,7 @@ namespace CNORXZ
 	friend YRangeFactory;
 
 	virtual RangePtr sub(SizeT i) const override final;
+	virtual MArray<RangePtr> sub() const override final;
 	virtual SizeT size() const override final;
 	virtual SizeT dim() const override final;
 	virtual String stringMeta(SizeT pos) const override final;
diff --git a/src/lib/ranges/range_base.cc b/src/lib/ranges/range_base.cc
index 52de19f..3b5d6fb 100644
--- a/src/lib/ranges/range_base.cc
+++ b/src/lib/ranges/range_base.cc
@@ -1,5 +1,6 @@
 
 #include "ranges/ranges.h"
+#include "array/array.h"
 
 namespace CNORXZ
 {
@@ -59,6 +60,11 @@ namespace CNORXZ
 	return nullptr;
     }
 
+    MArray<RangePtr> RangeBase::sub() const
+    {
+	return MArray<RangePtr>();
+    }
+
     bool RangeBase::operator==(const RangeBase& in) const
     {
 	return this == &in;
@@ -136,4 +142,18 @@ namespace CNORXZ
 	return o;
     }
 
+    RangePtr getSub(const RangePtr& r, SizeT num)
+    {
+	RangePtr o = r->sub(num);
+	CXZ_ASSERT(o != nullptr, "try to access empty subrange at position " << num);
+	return o;
+    }
+
+    MArray<RangePtr> getSub(const RangePtr& r)
+    {
+	MArray<RangePtr> o = r->sub();
+	CXZ_ASSERT(o.range() != nullptr, "try to access subrange of single range");
+	return o;
+    }
+
 } // end namespace CNORXZ
diff --git a/src/lib/ranges/yrange.cc b/src/lib/ranges/yrange.cc
index 9b041b8..e7b0e29 100644
--- a/src/lib/ranges/yrange.cc
+++ b/src/lib/ranges/yrange.cc
@@ -1,5 +1,6 @@
 
 #include "ranges/ranges.h"
+#include "array/array.h"
 
 namespace CNORXZ
 {
@@ -445,6 +446,14 @@ namespace CNORXZ
 	return mRVec[i];
     }
 
+    MArray<RangePtr> YRange::sub() const
+    {
+	if(mRVec.size() == 0){
+	    return MArray<RangePtr>();
+	}
+	return MArray<RangePtr>( CRangeFactory(mRVec.size()).create(), mRVec );
+    }
+
     SizeT YRange::size() const
     {
 	SizeT out = 1;