Dear Everybody,
Thank you very much for your reply, but I still have some problems.
For example, if I want to extract data for the particle with Dmax closest to 
1mm, and save it in xml format.Based on your suggestion and the code from the 
file of '/DataInterfaces/DataInterfaces/Matlab/demo_ssp4arts.m', I typed the 
following command in atmlab:

D=load('6-BulletRosette');
dmax=[D.M.diameter_max];
[~,i]=min(abs(dmax-1e-3));
S=D.S(i);
M=D.M(i);
Sc_data{1}{1}=S(1);
Sc_meta{1}{1}=M(1);
outfolder=pwd;
sdname='6-BulletRosette_dmax1_data.xml';
smname='6-BulletRosette_dmax1_meta.xml';
xmlStore( fullfile( outfolder, sdname ), Sc_data, ...
            'ArrayOfArrayOfSingleScatteringData', 'binary' );
xmlStore( fullfile( outfolder, smname ), Sc_meta, ...
            'ArrayOfArrayOfScatteringMetaData', 'binary' );

Then there is such a mistake:
Error using xmlWriteScatteringMetaData (line 22)
Illegal ScatteringMetaData version number

Error in xmlWriteArrayOf (line 36)
  feval (func, fid, fidb, data{i}, precision);

Error in xmlWriteArrayOfArrayOf (line 31)
  xmlWriteArrayOf (fid, fidb, data{i}, type, precision);

Error in xmlStoreImpl (line 53)
  xmlWriteArrayOfArrayOf (fid, fidb, data, type(15:length(type)), precision);

Error in xmlStore (line 30)
  xmlStoreImpl(filename, '', data, type, precision)

I don't know how to solve this, and then how to use atmlab to import the xml 
file needed by arts?
I used to use this simple commands to load the scattering data in atmlab:

C.SCAT_DATA{1} = fullfile( arts_scadata_path, 'hongdda_column_dmax2um.xml' );
or C.SCAT_DATA{1} = xmlLoad( 'MieSphere_R1.00000e+02um.xml' );


Looking forward to hearing from you, and sorry for disturbing you.
_______________________________________________
arts_users.mi mailing list
arts_users.mi@lists.uni-hamburg.de
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi

Reply via email to