Merge remote-tracking branch 'origin/devel' into SOdevel

.github/CONTRIBUTING.md

@@ -21,6 +21,7 @@ Git
 2. Routines in [master](https://github.com/irfu/irfu-matlab/tree/master) branch should do not less and not more than written in their help!
 3. Our [devel](https://github.com/irfu/irfu-matlab/tree/devel) branch includes latest delivered development changes for the next release.
 4. Semantic release numbering https://semver.org/
+5. Write clear and useful git commit messages https://chris.beams.io/posts/git-commit/
 
 Common data types
 -----------------

.github/pom.xml

@@ -0,0 +1,17 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<project xmlns="http://maven.apache.org/POM/4.0.0"
+  xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"
+  xsi:schemaLocation="http://maven.apache.org/POM/4.0.0 http://maven.apache.org/xsd/maven-4.0.0.xsd">
+  <modelVersion>4.0.0</modelVersion>
+  <groupId>se.irfu</groupId>
+  <artifactId>irfu-matlab</artifactId>
+  <version>1.14.2</version>
+  <dependencies>
+    <!-- https://mvnrepository.com/artifact/org.xerial/sqlite-jdbc -->
+    <dependency>
+      <groupId>org.xerial</groupId>
+      <artifactId>sqlite-jdbc</artifactId>
+      <version>3.28.0</version>
+    </dependency>
+  </dependencies>
+</project>

@dataobj/dataobj.m

@@ -111,7 +111,11 @@
               end
             else
               data{iVar(i)}(data{iVar(i)}==0) = NaN; % fillvalue timeline
-              data{iVar(i)} = irf_time(data{iVar(i)},'cdfepoch>epoch');
+              dTmp = data{iVar(i)};
+              for iCol=1:size(dTmp,2)
+                dTmp(:,iCol) = irf_time(dTmp(:,iCol),'cdfepoch>epoch');
+              end
+              data{iVar(i)} = dTmp; clear dTmp
               % bug fix for spdfcdfread (time comes out as row vector)
               if size(data{iVar(i)},1)~=info.Variables{iVar(i),3}
                 data{iVar(i)} = data{iVar(i)}';

@dataobj/get_ts.m

@@ -32,7 +32,10 @@
 
 %userData
 ud = data; 
-ud = rmfield(ud,{'DEPEND_0','data','nrec','dim','name','variance','UNITS'});
+
+field_list={'DEPEND_0','data','nrec','dim','name','variance','UNITS'};
+ud = rmfield(ud,field_list(isfield(ud,field_list)));
+
 repres = [];
 if isfield(data,'TENSOR_ORDER') % CAA data has TENSOR_ORDER>=1
   tensorOrder = data.TENSOR_ORDER; ud = rmfield(ud,'TENSOR_ORDER');