readTsgDataTxt.m

function [error] = readTsgDataTxt( hMainFig, filename)
% Function to read the TSG data. Should be a NetCDF file
%
% Input
% -----
% hMainFig ............ Handle to the main user interface
% filename ........... Data filename
%
% Output
% ------
% error .............. 1: OK - -1 : an error occured
%
% The data are store using setappdata - Variable name : 'tsg_data'
%
%
% Caution : replace the fill-value with NaN
% $Id$

% Get the data from the application GUI
% -------------------------------------
tsg = getappdata( hMainFig, 'tsg_data');

% Display read file info on console
% ---------------------------------
fprintf('\nREAD_ASCII_FILE\n'); tic;

% Open the file
% -------------
fid = fopen( filename, 'r' );

% Check file
% -----------
if fid == -1
  msg_error = ['TSG_GOSUD file_lecture : Open file error : ' filename];
  warndlg( msg_error, 'ASCII error dialog');
  sprintf('...cannot locate %s\n', filename);
  error = -1;
  return;
end

% Display more info about read file on console
% --------------------------------------------
fprintf('...reading %s : ', filename);

% Read the file
% -------------
tsgData = fscanf(fid, '%d/%d/%d %d:%d:%d %f %f %f %f %d %f %f', [13 Inf])';

% Every variable are put in a structure
% -------------------------------------
tsg.DAYD = datenum(tsgData(:,3), tsgData(:,2),tsgData(:,1), ...
                   tsgData(:,4),tsgData(:,5),tsgData(:,6));

% save original date
% ------------------
tsg.DATE                = datestr(  tsg.DAYD, 'yyyymmddHHMMSS' );
tsg.LATX                = tsgData(:,7);
tsg.LONX                = tsgData(:,8);
tsg.SSJT                = tsgData(:,9);
tsg.SSPS                = tsgData(:,10);
tsg.SSPS_QC             = tsgData(:,11);
tsg.SSPS_ADJUSTED       = tsgData(:,12);
tsg.SSPS_ADJUSTED_ERROR = tsgData(:,13);

tsg.SSJT_QC             = tsg.SSPS_QC;

% populate tsg.file structure
% ---------------------------
[tsg.file.pathstr, tsg.file.name, tsg.file.ext, tsg.file.versn] = ...
                                                         fileparts(filename);
%tsg.file.name = filename;
tsg.file.type = 'ASCII';

% Keep somme information for the log file
% ---------------------------------------
tsg.report.tsgfile = filename;

% Save the data in the application GUI
% ------------------------------------
setappdata( hMainFig, 'tsg_data', tsg );

% Perform somme automatic tests
% -----------------------------
automaticQC( hMainFig );

% Close the file
% --------------
fclose( fid );

% Clear the Workspace
% -------------------
clear tsgdata

% Display time to read file on console
% ------------------------------------
t = toc; fprintf('...done (%6.2f sec).\n\n',t);

% Everything OK
% -------------
error = 1;