Index: /issm/trunk-jpl/src/m/contrib/defleurian/netCDF/export_netCDF.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/defleurian/netCDF/export_netCDF.m	(revision 26967)
+++ /issm/trunk-jpl/src/m/contrib/defleurian/netCDF/export_netCDF.m	(revision 26968)
@@ -1,5 +1,5 @@
 function export_netCDF(md,filename)
 %verbosity of the code, 0 is no messages, 5 is chatty
-	verbose = 0;
+	verbose = 5;
 	if exist(filename),
 		delete(filename)
@@ -182,5 +182,5 @@
 					end
 					for k=1:length(subfields),
-						if ~ismember(subfields{k}, {'errlog', 'outlog'})
+						if ~ismember(subfields{k}, {'errlog', 'outlog', 'SolutionType'})
 							StackedVar=restable();
 							for l=1:Listsize,
@@ -201,4 +201,11 @@
 							if ~isempty(varid),
 								FillVar(StackedVar,subgroupID,varid);
+							end
+						elseif ismember(subfields{k}, {'SolutionType'})
+							%We just add solution type once as an attribute
+							Var = md.(groups{i}).(fields{j})(1).(subfields{k});
+							[DimSize,DimValue,varid]=CreateVar(ncid,Var,subgroupID,subfields{k},DimSize,DimValue);
+							if ~isempty(varid),
+								FillVar(Var,subgroupID,varid);
 							end
 
Index: /issm/trunk-jpl/src/m/contrib/defleurian/netCDF/export_netCDF.py
===================================================================
--- /issm/trunk-jpl/src/m/contrib/defleurian/netCDF/export_netCDF.py	(revision 26967)
+++ /issm/trunk-jpl/src/m/contrib/defleurian/netCDF/export_netCDF.py	(revision 26968)
@@ -36,8 +36,14 @@
                     maxsize.append(np.nanmax(datasize[:, dim]))
             except IndexError:
-                maxsize.append(np.nanmax(datasize[:]))
+                if datasize.ndim == 0:
+                    maxsize.append(datasize)
+                else:
+                    maxsize.append(np.nanmax(datasize[:]))
             findim = np.insert(maxsize, 0, rows)
             #first check if all steps are the same size
-            SameSize = np.sum(np.abs(datasize - datasize[0])) == 0
+            if datasize.ndim == 0:
+                SameSize = True
+            else:
+                SameSize = np.sum(np.abs(datasize - datasize[0])) == 0
             if SameSize:
                 #same size for all steps, just reshape
Index: /issm/trunk-jpl/src/m/io/loadvars.py
===================================================================
--- /issm/trunk-jpl/src/m/io/loadvars.py	(revision 26967)
+++ /issm/trunk-jpl/src/m/io/loadvars.py	(revision 26968)
@@ -199,5 +199,11 @@
                                 if timeindex:
                                     if NewFormat:
-                                        if vardim == 1:
+                                        if vardim == 0:
+                                            try:
+                                                Tree[0].__dict__[str(var)] = varval[timeindex].data
+                                            except IndexError:
+                                                print('WARNING: No data on index {} for {} reverting to last time.'.format(timeindex, str(var)))
+                                                Tree[0].__dict__[str(var)] = varval[-1].data
+                                        elif vardim == 1:
                                             try:
                                                 Tree[0].__dict__[str(var)] = varval[timeindex].data