NODES and PROCS variables in export headers

CODE/etc/rawformats.conf

@@ -22,6 +22,7 @@ winston|EARTHWORM|EarthWorm Winston Wave Server data request||host:port
 gipsy|GNSS|JPL GIPSY-OASIS .tdp file|ANTENNA,RECEIVER,XYZ|fullpath of root directory containing YYYY/FID/*.tdp files
 gipsyx|GNSS|JPL GipsyX .tdp file|ANTENNA,RECEIVER,XYZ|fullpath of root directory containing YYYY/FID/*.tdp files
 globkval|GNSS|MIT GAMIT/GLOBK VAL file||fullpath of directory containing the .VAL file
+spotgins-enu-v2|GNSS|SPOTGINS solutions - version 2 (<2025-08)||fullpath of file(s) with bash wildcard facilities, possible $FID, or URL to single file
 spotgins-ippp|GNSS|SPOTGINS IPPP time series||fullpath of file(s) with bash wildcard facilities, possible $FID, or URL to single file
 gamit-pos|GNSS|MIT GAMIT/GLOBL POS time series||fullpath of file(s) with bash wildcard facilities, possible $FID, or URL to single file
 pbogps-pos|GNSS|PBO GPS POS time series||fullpath of file(s) with bash wildcard facilities, possible $FID, or URL to single file

CODE/matlab/add_export_metadata.m

@@ -0,0 +1,31 @@
+function E = add_export_metadata(E, NP, export_header_keylist, np_type)
+%ADD_EXPORT_METADATA Add metadata to export structure based on node and process parameters
+%
+% INPUTS:
+%   E - Export structure to modify
+%   NP - Node or Proc configuration structure
+%   export_header_keylist - Cell array of node/proc field names to export
+%   np_type - String indicating the type ('NODE' or 'PROC')
+%
+% OUTPUT:
+%   E - Modified export structure with added metadata
+%
+%   Author: Pierre Sakic / WEBOBS, IPGP
+%   Created: 2025-08-22
+%   Updated: 2025-08-22
+
+    % Initialize meta field if it doesn't exist
+    if ~isfield(E, 'meta')
+        E.meta = struct();
+    end
+
+    % Add node metadata
+    if ~isempty(export_header_keylist)
+        for iexport = 1:length(export_header_keylist)
+            fieldname = export_header_keylist{iexport};
+            if isfield(NP, fieldname) 
+                E.meta.([np_type, '.', fieldname]) = any2str(NP.(fieldname));
+            end
+        end
+    end
+end

CODE/matlab/plotorbit.m

@@ -19,7 +19,9 @@
 if size(d,2) > 1
 	set(gca,'Ylim',get(gca,'YLim'))	% freezes Y axis (error bars can overflow)
 	h = plot(repmat(t,[1,2])',(repmat(d(:,1),[1,2])+d(:,2)*[-1,1])','-','LineWidth',.1,'Color',.6*[1,1,1]);
-    uistack(h,'bottom')
+	if exist('uistack','file') == 2 % checks if uistack function exists, (only MATLAB, not octave compatible)
+		uistack(h,'bottom')
+	end
 end
 
 % overwrites non-final orbits

CODE/matlab/readfmtdata.m

@@ -65,7 +65,7 @@
     case {'miniseed','seedlink','arclink','combined','fdsnws-dataselect'}
         D(n) = readfmtdata_miniseed(WO,P,N(n),F);
 
-    case {'globkval','gipsy','gipsyx','gipsy-tdp','usgs-rneu','ies-neu','ogc-neu','ingv-gps','sbe37-ascii','spotgins-ippp','gamit-pos','pbogps-pos'}
+    case {'globkval','gipsy','gipsyx','gipsy-tdp','usgs-rneu','ies-neu','ogc-neu','ingv-gps','sbe37-ascii','spotgins-enu-v2','spotgins-ippp','gamit-pos','pbogps-pos'}
         D(n) = readfmtdata_gnss(WO,P,N(n),F);
 
     case {'hyp71sum2k','fdsnws-event','scevtlog-xml'}
@@ -99,7 +99,8 @@
         D(n) = readfmtdata_mc3(WO,P,N(n),F);
 
     otherwise
-        D(n).t = [];
+        fprintf('%s: ** WARNING ** unknown format "%s", will try to read it anyway...\n',wofun,F.fmt);
+    	D(n).t = [];
         [D(n).d,D(n).CLB] = calib([],[],N(n).CLB);
         D(n).e = ones(size(D(n).d));
         fprintf('%s: ** WARNING ** unknown format "%s". Nothing to do!\n',wofun,F.fmt);

CODE/matlab/readfmtdata_gnss.m

@@ -34,7 +34,13 @@
 %		data format: extract from SmoothFinal.tdp output file (grep "Station.SSSS.State.pos.[XYZ]" lines)
 %		node calibration: no .CLB file or 4 components (East, North, Up) in meters and (Orbit)
 %
-%	format 'gins-ippp'
+%	format 'spotgins-enu-v2'
+%		type: SPOTGINS solutions - version 2 (<2025-08)
+%		filename/url: P.RAWDATA (use $FID to point the right file/url)
+%		data format: jjjjj.jj E N V dE dN dV yyyymmddhhmmss yyyy.yyyyyyyyy
+%		node calibration: no .CLB file or 4 components (East, North, Up) in meters and (Orbit)
+%
+%	format 'spotgins-ippp'
 %		type: GINS IPPP solutions
 %		filename/url: P.RAWDATA (use $FID to point the right file/url)
 %		data format: yyyymmdd hhmmss yyyy.yyyyyyyyy jjjjj.jj X Y Z dX dY dZ E N V dE dN dV
@@ -300,8 +306,56 @@
 	end
 	%D.ITRF_YEAR = 'ITRF08';
 
+
+% -----------------------------------------------------------------------------
+case 'spotgins-enu-v2'
+        % format exemple
+        %#jjjjj.jjjjjjjj         _____E         _____N         _____U         ____dE         ____dN         ____dU  yyyymmddHHMMSS  yyyy.yyyyyyy  Const  Dateofexe       GinsVersion
+        % 52670.83876160       0.055822       0.051638       0.005578       0.001263       0.001163       0.004899  20030131200749  2003.0844898  G      250404_185253   VALIDE_24_2
+        % 52671.50195600       0.057207       0.054240      -0.004722       0.000705       0.000619       0.002600  20030201120249  2003.0863067  G      250404_185253   VALIDE_24_2
+
+	fdat = sprintf('%s/%s.dat',F.ptmp,N.ID);
+	wosystem(sprintf('rm -f %s',fdat),P);
+	for a = 1:length(F.raw)
+		fraw = F.raw{a};
+                cmd0 = sprintf('awk ''/^[^#]/ {print}'' >> %s',fdat); % removes header lines
+		if strncmpi('http',fraw,4)
+			s  = wosystem(sprintf('curl -s -S "%s" | %s',fraw,cmd0),P);
+			if s ~= 0
+				break;
+			end
+		else
+			s = wosystem(sprintf('cat %s | %s',fraw,cmd0),P);
+		end
+		if s ~= 0
+			fprintf('%s: ** WARNING ** Raw data "%s" not found.\n',wofun,fraw);
+		end
+	end
+
+
+        % load the file
+	if exist(fdat,'file')
+		dd = dlmread(fdat);
+	else
+		dd = [];
+	end
+	if ~isempty(dd)
+                t = dd(:,1) + 678941.5007; % converts MJD to datenum
+		d = [dd(:,2:4),zeros(size(dd,1),1)]; % North(mm),East(mm),Up(mm) => E(m),N(m),U(m),Orbit
+		e = dd(:,5:7);
+		e(e<min_error) = min_error;
+		fprintf('%d data imported.\n',size(dd,1));
+	else
+		fprintf('no data found!\n')
+		t = [];
+		d = [];
+		e = [];
+	end
+
+
 % -----------------------------------------------------------------------------
 case 'spotgins-ippp'
+        % From J.S. - ITES Strasbourg
 	% format example
 	% !yyyymmdd hhmmss yyyy.yyyyyyyyy  jjjjj.jj        X_position        Y_position        Z_position            dX            dY            dZ             E             N             V            dE            dN            dV
     %  20160723  65619 2016.558521561  57592.29    4182067.152057     570976.439258    4765940.539811      0.000611      0.000218      0.000673     -0.006574     -0.008848     -0.014844      0.000205      0.000307      0.000859

CODE/matlab/superprocs/gnss.m

@@ -81,6 +81,9 @@
 cmpnames = split(field2str(P,'COMPONENT_NAMELIST','Relative Eastern,Relative Northern,Relative Vertical'),',');
 disp_yscale = field2num(P,'DISP_YSCALE_M',0);
 
+export_header_proc_keylist = split(field2str(P,'EXPORT_HEADER_PROC_KEYLIST',''),',');
+export_header_node_keylist = split(field2str(P,'EXPORT_HEADER_NODE_KEYLIST',''),',');	
+
 % Harmonic correction: period list (day), pairs of sine, cosine (mm) for each component
 harm_refdate = field2num(P,'HARMONIC_ORIGIN_DATE');
 harm_period = field2num(P,'HARMONIC_PERIOD_DAY',0);
@@ -306,7 +309,7 @@
 modeltime_markersize = pi*(field2num(P,'MODELTIME_MARKERSIZE',10,'notempty')/2)^2; % scatter needs marker size as a surface (πr²)
 
 
-geo = [cat(1,N.LAT_WGS84),cat(1,N.LON_WGS84),cat(1,N.ALTITUDE)];
+geo = [cat(1,N.LAT_WGS84),cat(1,N.LAT_WGS84),cat(1,N.ALTITUDE)];
 
 V.name = P.NAME;
 V.velref = itrf;
@@ -653,7 +656,13 @@
 				E.header = [E.header,{'East_treat(m)','North_treat(m)','Up_treat(m)'}];
 			end
 			E.title = sprintf('%s {%s}',P.GTABLE(r).GTITLE,upper(N(n).ID));
+
+			E.meta = {};
+			E = add_export_metadata(E,N(n),export_header_node_keylist,"NODE");
+			E = add_export_metadata(E,P,export_header_proc_keylist,"PROC");
+
 			mkexport(WO,sprintf('%s_%s',N(n).ID,P.GTABLE(r).TIMESCALE),E,P.GTABLE(r));
+			E.meta = {}; % meta are erased after the export
 		end
 	end
 
@@ -2328,6 +2337,7 @@
                     E.infos = cat(2,E.infos,sprintf('Time period #%d = %g days (%s)',m,modeltime_period(m),days2h(modeltime_period(m),'round')));
                 end
                 mkexport(WO,sprintf('%s_VECTORS_%s_%s',summary,lower(N(kn(s)).ID),P.GTABLE(r).TIMESCALE),E,P.GTABLE(r));
+				E.meta = {}; % meta are erased after the export
             end
 
             % modeltime results
@@ -2364,6 +2374,7 @@
 end
 end
 
+
 if P.REQUEST
     mkendreq(WO,P);
 end