Discussion:
[arts-users] problems of importing the ARTS Microwave Single Scattering Properties Database
WCS
2018-03-06 03:47:00 UTC
Permalink
Dear All,
Hello, I am a first year graduate student from China. I have encountered problems of importing data when using atmlab and arts to solve with scatterring.
Because I used the Hong scattering database before, each .xml file corresponds to one particle shape and size. ARTS Microwave Single Scattering Properties Database has a set of ARTS-ready and easy to use scattering data in the ARTS xml-format. A total of 16 habits are provided in StandardHabit.star.gz. After decompression, in the FullSet folder, there are four files for each particle shape, .xml, .meta, .xml.bin and matlab files, and I feel that the .xml file contains different particle sizes of each shape. So now I want to import a scatterring data of one kind of particle shape and one size in atmlab, how to enter the command in atmlab?

Look forward to your reply as soon as possible, thank you for helping me to solve the problem.
Sincerely
WCS
Patrick Eriksson
2018-03-06 08:55:36 UTC
Permalink
Hi,

To really just load one size, you need to work against the actual database.

However, if you are working in Matlab, it is probably still quicker to
use the .mat data, as these data can be imported quickly. Below is an
example how to use .mat, and extract data for the particle with Dmax
closest to 1mm.

As a general comment: Notice that the special habits are stored as
struct, in order to allow things like

dmax=[M.diameter_max];

However, to store the data to ARTS xml, you need to convert the data to
a cell array: Sc{i} = S(i) etc.
The details differ slightly if you want to store
ArrayOfArrayOfSingleScatteringData or ArrayOfSingleScatteringData.

The same applies to the meta data.

Bye,

Patrick
Post by WCS
D=load('8-ColumnAggregate');
dmax=[D.M.diameter_max];
[~,i] = min(abs(dmax-1e-3));
S=D.S(i);
M=D.M(i);
S
S =

struct with fields:

version: 3
ptype: 'totally_random'
description: 'Hexagonal column aggregate, using Hong
parametrization (Hong 2007).'
f_grid: [1x34 double]
T_grid: [190 210 230 250 270]
za_grid: [721x1 double]
aa_grid: []
pha_mat_data: [7-D double]
ext_mat_data: [34x5 double]
abs_vec_data: [34x5 double]
Post by WCS
M
M =

struct with fields:

version: 3
description: 'Meta data for Hexagonal column
aggregate, using Hong parametrization (Hong 2007).'
source: 'location: Chalmers, software:
Snowflake-toolkit (Rathsman 2016), parameterization: Hong 2007.'
refr_index: 'Matzler 2006'
mass: 6.4001e-08
diameter_max: 9.9258e-04
diameter_volume_equ: 5.1088e-04
diameter_area_equ_aerodynamical: 7.4820e-04
Post by WCS
Dear All,
        Hello, I am a first year graduate student from China. I have
encountered problems of importing data when using  atmlab and arts to
solve with scatterring.
        Because I used the Hong scattering database before,  each .xml
file corresponds to one particle shape and size. ARTS Microwave Single
Scattering Properties Database has a set of ARTS-ready and easy to use
scattering data in the ARTS xml-format. A total of 16 habits are
provided in StandardHabit.star.gz. After decompression, in the FullSet
folder, there are four files for each particle shape,  .xml, .meta,
.xml.bin and matlab files, and I feel that the .xml file contains
different particle sizes of each shape.  So now I want to import a
scatterring data of one kind of particle shape and one size in atmlab,
how to enter the command in atmlab?
        Look forward to your reply as soon as possible, thank you for
helping me to solve the problem.
                                                                                                                                                                                                                                 Sincerely
          WCS
_______________________________________________
arts_users.mi mailing list
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi
Loading...