Hi Frank, To initialize the SSD from a dictionary, you can use the from_data classmethod. Also note that using integer numbers for ptype is obsolete, they have been replaced by strings ("general", "totally_random", "azimuthally_random" for version 3 of the SSD).
Here's an example with data dimensions adjusted for the azimuthally random case: ------------------------------------ import typhon import numpy as np params = {'f_grid': [174000000000.0, 210000000000.0], 'T_grid': [190, 200, 210, 220.0, 230], 'za_grid': [0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 110, 120, 130, 140, 150, 160, 170, 180], 'aa_grid': [0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 110, 120, 130, 140, 150, 160, 170, 180], 'aspect_ratio': 1.5, 'NP': -1, 'precision': .001, 'phase': 'ice', 'ptype': 'azimuthally_random', 'equiv_radius': 100} scatt_data = typhon.arts.scattering.SingleScatteringData.from_data(params) scatt_data.pha_mat_data = np.zeros( (len(params['f_grid']), len(params['T_grid']), len(params['za_grid']), len(params['aa_grid']), len(params['za_grid']), 1, 16)) scatt_data.ext_mat_data = np.zeros( (len(params['f_grid']), len(params['T_grid']), len(params['za_grid']), 1, 3)) scatt_data.abs_vec_data = np.zeros( (len(params['f_grid']), len(params['T_grid']), len(params['za_grid']), 1, 2)) typhon.arts.xml.save(scatt_data, 'scattering.xml') ------------------------------------ Cheers, Oliver > On 15. Nov 2019, at 19:31, Werner, Frank (329D) <frank.wer...@jpl.nasa.gov> > wrote: > > Hi all, > > For some ARTS runs including scattering by cirrus I am trying to create and > write single scattering properties with the module ‘SingleScatteringData’ > from the typhon package. Unfortunately, I cannot seem to get the syntax right > (just calculating the data, not even writing it into an xml file). A very > simple example would be highly appreciated. Here is what I tried to do: > > import typhon > params = {'f_grid':[174000000000.0, 210000000000.0],'T_grid':[190,200,210, > 220.0,230], > 'za_grid':[0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 110, 120, > 130, 140, 150, 160, 170, 180], > 'aa_grid':[0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 110, 120, > 130, 140, 150, 160, 170, 180], > 'aspect_ratio':1.5,'NP':-1,'precision':.001,'phase':'ice','ptype':30,'equiv_radius':100} > > scatt_data = typhon.arts.scattering.SingleScatteringData (params) > > This yields the error: > > ‘TypeError: __init__() takes 1 positional argument but 2 were given’. This is > probably a trivial issue; the possible parameters are described in the user > manual. I just cannot seem to invoke the function correctly. > > Thanks for your help and best wishes, > Frank > > -- > Frank Werner > Mail Stop 183-701, Jet Propulsion Laboratory > 4800 Oak Grove Drive, Pasadena, California 91109, United States > Phone: +1 818 354-1918 > Fax: +1 818 393 5065 > _______________________________________________ > arts_users.mi mailing list > arts_users.mi@lists.uni-hamburg.de <mailto:arts_users.mi@lists.uni-hamburg.de> > https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi > <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