Dear WCS,

These scattering data require that you are using v2.3.

Could it be that you are using v2.2? If yes:

I have not considered this situation, but it could be quite easy to convert the data to a version working with v2.2. However, I don't have time right now to check the details. Anyone that remember the details here?

I hesitate to say that you should upgrade to v2.3. This as Atmlab v2.3 is not poorly synced with the development in ARTS around scattering solvers.

Bye,

Patrick


On 03/06/18 15:29, WCS wrote:
*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

_______________________________________________
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