fixes bug in returned values

pyatsoleiloutils/readmatfile.py

@@ -34,7 +34,6 @@ def readmomapmat(
     s_indexes=f['momap']['s_indexes'][:]
     s_array=f['momap']['s_array'][:]
     index_first_element = f['momap']['index_first_element'][:]
-    f.close()
 
-    return momap_pn, s_indexes, s_array, first_element_index
+    return momap_pn, s_indexes, s_array, index_first_element