this file was moved to tsg_tools directory

tsg_util/netcdfTOtsgqc.m

-% Conversion d'une liste de fichier NetCDF TSG en fichier ASCCI
-%
-% La liste des fichiers netCDF doit tre plac dans un dichier ASCII.
-% un nom de fichier par ligne
-%
-clear all;
-clc;
-
-% Repertoire des fichiers NetCdf en lecture
-% -----------------------------------------
-dirIn = 'D:\temp\';
-
-% Repertoire des fichiers ASCII en criture
-% -----------------------------------------
-dirOut = dirIn;
-
-% Nom du fichier qui liste les fichiers netCdf
-% --------------------------------------------
-fichierListe = 'list_nc.txt';
-
-% Liste des paramtres  crire - Create empty structure array
-% ------------------------------------------------------------
-PARA = { 'CNDC'; 'CNDC_CAL'; ...
-  'SSPS'; 'SSPS_QC'; 'SSPS_CAL';...
-  'SSPS_ADJUSTED'; 'SSPS_ADJUSTED_QC'; 'SSPS_ADJUSTED_ERROR'; ...
-  'SSJT'; 'SSJT_QC'; 'SSJT_CAL';...
-  'SSJT_ADJUSTED'; 'SSJT_ADJUSTED_QC'; 'SSJT_ADJUSTED_ERROR'; ...
-  'SSTP'; 'SSTP_QC'; 'SSTP_CAL';...
-  'SSTP_ADJUSTED'; 'SSTP_ADJUSTED_QC'; 'SSTP_ADJUSTED_ERROR'};
-nPara = size( PARA, 1 );
-
-% Lecture du nom des fichiers  convertir. Ceux-ci se trouvent dans
-% le fichier liste_nc.txt
-% -----------------------------------------------------------------
-list_nc = textread([dirIn fichierListe],'%s');
-
-% Conversion du nom des fichiers en caracteres
-% --------------------------------------------
-list_nc = char(list_nc);
-
-% Boucle sur les fichiers  convertir
-% ----------------------------------
-nfiles  = size(list_nc,1);
-for i = 1 : nfiles
-
-  % Create structure array
-  % ----------------------
-  tsg = struct;
-
-  fileName = deblank(list_nc(i,:));
-  display( ['Traitement du fichier : ' fileName]);
-
-  % Ouverture d'un fichier NetCdf en lecture
-  % ----------------------------------------
-  nc = netcdf([dirIn fileName],'read');
-
-  if ~isempty(nc)
-
-    % Creation et ouverture du fichier en criture .tsgqc
-    % ---------------------------------------------------
-    ind         = find(fileName == '.');
-    fileNameOut = fileName(1 : ind-1);
-    fid         = fopen([dirOut fileNameOut '.tsgqc'],'wt');
-    display( ['          Ecriture du fichier : ' fileNameOut]);
-
-    % --------------------------------
-    % Get global attributes: meta data
-    % --------------------------------
-    global_att = att(nc);
-    for iatt = 1:length(global_att)
-
-      % extract name and value from netcdf globals attributes
-      % -----------------------------------------------------
-      attribute = name(global_att{iatt});
-      value     = global_att{iatt}(:);
-
-      % assign globals attributes in base workspace
-      % -------------------------------------------
-      assignin('base', attribute, value);
-
-      % populate tsg structure with globals attributes
-      % ----------------------------------------------
-      tsg.(attribute) = value;
-
-    end % end of attributes loop
-
-    % ----------------------------
-    % Get variables describing TSG
-    % ----------------------------
-    variables = var(nc);
-    for ivar = 1:length(variables)
-
-      % extract name and value from netcdf variables
-      % --------------------------------------------
-      variable = name(variables{ivar});
-      nv       = nc{ivar};
-
-      % use autonan mode, remplace fillValue with NaN
-      % ---------------------------------------------
-      nv = autonan(nv, true);
-
-      % populate tsg structure with netcdf variables
-      % --------------------------------------------
-      tsg.(variable) = nv(:);
-
-      % assign netcdf variables in base workspace
-      % -----------------------------------------
-      assignin('base', variable, nv(:));
-
-    end % end of variables loop
-
-    % Initialise le tableau des variables  crire
-    % --------------------------------------------
-    nbRecords = length( tsg.DAYD );
-    tsg_data = NaN * ones(nPara+8, nbRecords);
-
-    % Selectionne les paramtres et cre le tableau en criture
-    % --------------------------------------------------------
-    if nbRecords ~= 0
-
-      % Met en forme les dates et positions
-      % ------------------------------------
-       year  = str2num(tsg.DATE(:,1:4));
-       month = str2num(tsg.DATE(:,5:6));
-       day   = str2num(tsg.DATE(:,7:8));
-       hour  = str2num(tsg.DATE(:,10:11));
-       min   = str2num(tsg.DATE(:,12:13));
-       sec   = str2num(tsg.DATE(:,14:15));
-
-      %[year, month, day, hour, min, sec] = datevec( julianToDatenum(tsg.DAYD) );
-      tsg_data = [year month day hour min fix(sec) tsg.LATX mod(tsg.LONX+180,360)-180 ];
-      header = 'YEAR MNTH DAYX hh mi ss LATX LONX';
-      format = '%04d %02d %02d %02d %02d %02d %12.7f %12.7f';
-
-      % Boucle sur les paramtres  sauvegarder
-      % ---------------------------------------
-      for iPara = 1 : nPara
-
-        % Slectionne le nom du paramtre
-        % -------------------------------
-        para = PARA{iPara};
-
-        if isempty(tsg.(para))
-          tsg.(para) = NaN * ones( size( tsg.DAYD ) );
-        end
-
-        if findstr( '_QC', para)
-          
-          % Test pour contourner un bug de TSGQC. Les QC devaient tre casts
-          % en byte sinon la valeur 0 devenait gale  -127
-          ind = 0;
-          ind = find( tsg.(para) == -127 );
-          if ~isempty(ind)
-            tsg.(para)(ind) = 0;
-          end
-          tsg_data = [tsg_data double(tsg.(para))];
-          format   = [format ' %1d'];
-        else
-          tsg_data = [tsg_data tsg.(para)];
-          format   = [format ' %6.3f'];
-        end
-        header   = [header ' ' para];
-
-      end % Fin de boucle sur les paramtres  sauvegarder
-
-      % Write some information
-      % ----------------------
-      if ~isempty( tsg.DATE_CREATION )
-        fprintf( fid, '%%DATE_CREATION %s\n', tsg.DATE_CREATION );
-      end
-      if ~isempty( tsg.CYCLE_MESURE )
-        fprintf( fid, '%%CYCLE_MESURE %s\n', tsg.CYCLE_MESURE);
-      end
-      if ~isempty( tsg.PLATFORM_NAME )
-        fprintf( fid, '%%PLATFORM_NAME %s\n', tsg.PLATFORM_NAME);
-      end
-      if ~isempty( tsg.SHIP_CALL_SIGN )
-        fprintf( fid, '%%SHIP_CALL_SIGN %s\n', tsg.SHIP_CALL_SIGN);
-      end
-      if ~isempty( tsg.PI_NAME )
-        fprintf( fid, '%%PI_NAME %s\n', tsg.PI_NAME);
-      end
-      if ~isempty( tsg.DATA_ACQUISITION )
-        fprintf( fid, '%%DATA_ACQUISITION %s\n', tsg.DATA_ACQUISITION);
-      end
-
-      % Write the file
-      % --------------
-      fprintf( fid, '%%HEADER %s\n', header);
-      fprintf( fid, [format '\n'], tsg_data');
-
-      % Everything OK
-      % -------------
-      error = 1;
-    end   % Selectionne les paramtres et cre le tableau en criture
-
-    fclose( fid );
-
-    % Clear the Workspace
-    % -------------------
-    clear tsg tsg_data;
-
-  else
-
-    msg_error = 'ok';
-    if isempty(nc)
-      msg_error = ['Problem  opening file : ' fileName];
-    end
-    
-  end
-end   % Find Boucle sur les fichiers  convertir