Skip to content
Snippets Groups Projects
Commit dbee209b authored by leprob001's avatar leprob001
Browse files

Updated save functions and tests to the new rfb input file.

parent ffcaaf9b
No related branches found
No related tags found
No related merge requests found
......@@ -11,8 +11,8 @@ rfb = artoa.load.rfb(dataFileName);
data = [];
data(:, 1) = rfb.DATA(1:13, rfb.VARIABLE_LIST.time_of_arrival(1));
data(:, 2) = str2num(cell2mat(rfb.SAT_DATA(:, 1)));
data(:, 3) = str2num(cell2mat(rfb.SAT_DATA(:, 2)));
data(:, 2) = rfb.SAT_DATA(1:13, rfb.SAT_FORMAT.lat_sat);
data(:, 3) = rfb.SAT_DATA(1:13, rfb.SAT_FORMAT.lon_sat);
data(:, 4) = rfb.DATA(1:13, rfb.VARIABLE_LIST.temperature);
data(:, 5) = rfb.DATA(1:13, rfb.VARIABLE_LIST.pressure);
......
......
......@@ -13,13 +13,14 @@ header = "Created in 2019" + newline + "a multiline test header";
data = [];
data(:, 1) = rfb.DATA(1:13, rfb.VARIABLE_LIST.time_of_arrival(1));
data(:, 2) = str2num(cell2mat(rfb.SAT_DATA(:, 1)));
data(:, 3) = str2num(cell2mat(rfb.SAT_DATA(:, 2)));
data(:, 2) = rfb.SAT_DATA(1:13, rfb.SAT_FORMAT.lat_sat);
data(:, 3) = rfb.SAT_DATA(1:13, rfb.SAT_FORMAT.lon_sat);
data(:, 4) = rfb.DATA(1:13, rfb.VARIABLE_LIST.temperature);
data(:, 5) = rfb.DATA(1:13, rfb.VARIABLE_LIST.pressure);
data(:, 6) = rfb.DATA(1:13, rfb.VARIABLE_LIST.pressure);
data(:, 7) = rfb.DATA(1:13, rfb.VARIABLE_LIST.temperature);
saveToFilename = fullfile(currentDirectory, '..', 'testdata', 'output', 'testOutput.trj4');
success = artoa.save.trj4(saveToFilename, header, data);
......
......
......@@ -95,6 +95,12 @@ floatCycleEndTimeRafosDays = num2str( ...
);
floatCycleEndTime = num2str(floatCycleEndTime);
if size(pFloatDetails.offset, 1) > 1
floatOffsetEnd = pFloatDetails.offset(2,end);
else
floatOffsetEnd = NaN;
end
header = { ...
['** Float: ' floatname], ...
['** Variables: InterpFlag LineNum RafosDay Temp Pres Lat Lon U V W'], ...
......@@ -109,7 +115,7 @@ header = { ...
['** Launch time : ' launchTime ' (RAFOS day ' launchTimeRafosDays ')'], ...
['** Cycle End time : ' floatCycleEndTime ' (RAFOS day ' floatCycleEndTimeRafosDays ')'], ...
['** First surface Position time : ' recoverTime ' (RAFOS day ' recoverTimeRafosDays ')'], ...
['** Float clock offset (init/final) : ' num2str(pFloatDetails.offset(1,end)) ' / ' num2str(pFloatDetails.offset(2,end)),' seconds'], ...
['** Float clock offset (init/final) : ' num2str(pFloatDetails.offset(1,end)) ' / ' num2str(floatOffsetEnd),' seconds'], ...
['** -------'] ...
};
......
......
......@@ -49,7 +49,7 @@ defaultValue = -999;
for i = 1:size(pData, 2)
nanIndex = find(isnan(pData(:, i)));
if ~isempty(nanIndex)
pData(nanIndex, i) = repmat(defaultValue, 1, size(pData, 1));
pData(nanIndex, i) = defaultValue;
end
clear nanIndex;
end
......
......
......@@ -45,7 +45,7 @@ defaultValue = -999;
for i = 1:size(pData, 2)
nanIndex = find(isnan(pData(:, i)));
if ~isempty(nanIndex)
pData(nanIndex, i) = repmat(defaultValue, 1, size(pData, 1));
pData(nanIndex, i) = defaultValue;
end
clear nanIndex;
end
......
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please to comment