% This script can be used for transfroming arts-data from the old
% format (ArrayOfTensor3) to the new format (GriddedField3).
%
% It's only a "quick hack", you have to adapt the paths in this file.
%
% 2004-02-24 Claudia Emde


% Initialize atmlab.

addpath('/smiles_local/atmlab/atmlab')
atmlab_init;

path = '/rinax/storage/users/claudia/MC_comparison/2101data/'

basename = 'jpl2101'

% midlatitude-winter, mitlatitude summer, tropical
% $$$ species = {'C2H2' 'C2H6' 'CH3Cl' 'CH4' 'ClO' 'CO2' 'COF2' 'CO' 'H2CO' ...
% $$$           'H2O2' 'H2O' 'HBr' 'HCl' 'HCN' 'HF' 'HI' 'HNO3' 'HOCl' ...
% $$$           'N2O' 'N2' 'NH3' 'NO2' 'NO' 'O2' 'O3' 'OH' 'PH3' 'SF6' ...
% $$$           'SO2' 't' 'z'}

%subarctic-summer, subarctic winter
species = {'CH3Cl' 'CH4' 'ClO' 'CO2' 'CO' 'H2CO' ...
     'H2O' 'HBr' 'HCl' 'HCN' 'HF' 'HI' 'HNO3' 'HOCl' ...
     'N2O' 'N2' 'NH3' 'NO2' 'NO' 'O2' 'O3' 'OH'  ...
     'SO2' 't' 'z'}
	   
% jpl_data
species = {'H2O' 'N2O' 'O2' 'O3' 't' 'z'}

for i = 1:length(species)

disp(i);
  
filename = [basename, '.', species{i}, '.xml'];

field =  xmlLoad([path, filename]);

gridded_field.p_grid = field{1}(:,1,1);
gridded_field.lat_grid = [0];
gridded_field.lon_grid = [0];
gridded_field.data = field{4}(:,:,:);

xmlStore([path, '/new_jpl_data/', filename], gridded_field, 'GriddedField3')

end