Re: [arts-users] scat_meta_single to scat_meta

2023-02-15 Thread Galligani
Hello! Perfect. I adapted everything to use only the bulk approach. Thanks for the help. Cheers, Vito On 03.02.2023 04:08, Patrick Eriksson wrote: Vito, Sorry, but partly filling pnd_field manually and partly filling it by using the PSD methods will likely never work. At least this

Re: [arts-users] scat_meta_single to scat_meta

2023-02-02 Thread Patrick Eriksson
Vito, Sorry, but partly filling pnd_field manually and partly filling it by using the PSD methods will likely never work. At least this deviates from how ARTS is designed. There are two main options here: 1: To fill pnd_field "manually". But this is likely more easily done outside of ARTS.

Re: [arts-users] scat_meta_single to scat_meta

2023-02-02 Thread Galligani
Hello ! Yes, for example this: ReadXML( pnd_field, "hail_pnddata1.xml" ) Tensor4Create(pnd_field2) ReadXML( pnd_field2, "hail_pnddata2.xml" ) Append( pnd_field, pnd_field2) # RWC ScatSpeciesSizeMassInfo( species_index=5, x_unit="dveq" ) Copy( psd_size_grid, scat_species_x ) Copy(

Re: [arts-users] scat_meta_single to scat_meta

2023-02-01 Thread Patrick Eriksson
Vito, Can you explain how it failed? Bye, Patrick On 2023-02-01 16:34, Galligani wrote: Hi Patrick! Thanks that worked out! I'm now running a Frankeinstein of code with pnd_data saved outside ARTS for different hail species, and using the ARTS psdModifiedGammaMass. Basically I'm trying to

Re: [arts-users] scat_meta_single to scat_meta

2023-02-01 Thread Galligani
Hi Patrick! Thanks that worked out! I'm now running a Frankeinstein of code with pnd_data saved outside ARTS for different hail species, and using the ARTS psdModifiedGammaMass. Basically I'm trying to append the RWC pnd_data from psdModifiedGammaMass pndFromPsdBasic to the pnd_field from

Re: [arts-users] scat_meta_single to scat_meta

2023-01-25 Thread Patrick Eriksson
Vito, scat_data and scat_data_raw are of type ArrayOfArrayOfSingleScatteringData. The first Array dimension matches scat_species and the second one size. So you need multiple Append to scat_data_raw and scat_meta. If I get it right, you should do this for each habit: ReadXML(

Re: [arts-users] scat_meta_single to scat_meta

2023-01-24 Thread Galligani
Hello everyone, I re-open this question because I am adding species to my simulations now (before I was running with only one species) I've run into the following error: 'Number of scattering species specified by scat_species does not agree with number of scattering species in scat_data:

Re: [arts-users] scat_meta_single to scat_meta

2022-12-19 Thread vito . galligani
Hello! Yes! I can use my own pnd_field. No problem. Thank you! Vito On 15.12.2022 18:39, Patrick Eriksson wrote: Hi, Thanks Oliver, I looked in the wrong function! Hm, it seems that I forgot about fixing the full chain for the mono PSDs, and that nobody has tried to use them until now. My

Re: [arts-users] scat_meta_single to scat_meta

2022-12-15 Thread Patrick Eriksson
Hi, Thanks Oliver, I looked in the wrong function! Hm, it seems that I forgot about fixing the full chain for the mono PSDs, and that nobody has tried to use them until now. My mistake is a bit shameful, but I feel reluctant to fix it now, as the present PSD system hopefully will be replaced

Re: [arts-users] scat_meta_single to scat_meta

2022-12-14 Thread Lemke, Oliver
Hi Patrick, hi Vito, > On 14. Dec 2022, at 20:04, Patrick Eriksson > wrote: > >> pndFromPsdBasic >> } >> Append( scat_species, species_id_string ) >> Append( pnd_agenda_array_input_names, pnd_agenda_input_names ) >> But the error message is that: >> The method requires that length of

Re: [arts-users] scat_meta_single to scat_meta

2022-12-14 Thread Patrick Eriksson
Vito, When it comes to: ArrayOfAgendaAppend( pnd_agenda_array ){   VectorSet(pnd_size_grid, [1])   VectorSet(psd_size_grid, [1])   psdMonoDispersive(psd_data, dpsd_data_dx, pnd_agenda_input_t, pnd_agenda_input,    pnd_agenda_input_names, dpnd_data_dx_names, scat_meta,

Re: [arts-users] scat_meta_single to scat_meta

2022-12-14 Thread vito . galligani
Hi Patrick! Thanks for your reply. I think I made some progress, but I'm still running into error messages (now with PSD) The single scattering data I'm happy with. Basically I'm doing this: ArrayOfScatteringMetaDataCreate( scat_meta_1habit ) ArrayOfSingleScatteringDataCreate(

Re: [arts-users] scat_meta_single to scat_meta

2022-12-14 Thread Patrick Eriksson
Hi, Yes, this is not totally streamlined. And nothing we will change now as the handling of scattering data is being rewritten from scratch. You must do a bit more than you included, as scat_dataCalc requires scat_data_raw. Not clear from where you get that. I think you can use Append to