Hi Patrick, Thanks for the suggestions.
1.) => I tried to check if Dmean was set to 0 value anywhere in the Bulkprop_fields when using "psdModifiedGammaMassXmean()" . The only place it is '0' is at the border of the cloudbox and outside the cloudbox. => This is how "pnd_fieldCalcFromParticleBulkProps()" expects to be set as given in the documentation. => When I try making the border values a positive non-zero value. It gives me this error. Exception: Call to ARTS WSM pnd_fieldCalcFromParticleBulkProps failed with error: *particle_bulkprop_field* allowed to contain non-zero values only inside the cloudbox. 2.) => So, this looks to me like a bug that if I make borders zero valued then "psdModifiedGammaMassXmean()" would have a problem, else if I make it non zero then "pnd_fieldCalcFromParticleBulkProps()" will give me problems. => Am I thinking this right? 3.) => Would you recommend me to change the source code in "PSD.CC" and change line 321 "if (ext_pars[1] <= 0) " and remove the "=" sign and build ARTS. Would this be a fix to the problem? Thanks, Renish Thomas ________________________________ From: Patrick Eriksson <patrick.eriks...@chalmers.se> Sent: Sunday, December 6, 2020 10:59 AM To: Thomas,Renish <renish.tho...@colostate.edu> Cc: arts users mi <arts_users.mi@lists.uni-hamburg.de> Subject: Re: Fwd: Re: [arts-users] ARTS internally generated PND fields Ranish, > 1.) So, if I must set the scattering particle temperature to an arbitrary > value after > calling "pnd_fieldCalcFromParticleBulkProps", do I just override and assign > 'pnd_agenda_input_t". You got me wrong. You don't need to do anything. pnd_agenda_input_t is set automatically where needed. > 2.) Can the particle temperature be set different than the ambient > temperature? Hardly. You can in principle try to "hack" pnd_agenda to set it to some other value. But why do you want to do it? And please note that the psdModifiedGamma PSD-s do not even have a temperature dependence. So in this case I really don't see any reason. The temperature also comes into play when interpolating scat_data, but that is handled by another part of ARTS. And over there, there is no option beside using the local temperature. > 3.) Do you recommend me to use internally generated PSD for MC or set it > externally? What you have started should work. So continue on this track. (The alternative is to set pnd_field already in Python. Then you can ignore particle_bulkprop_field ad pnd_agenda etc. I can also mention that the PSD system in ARTS v2.5 soon will change. Simon is making a complete revision of the handling of particles.) Bye, Patrick
_______________________________________________ arts_users.mi mailing list arts_users.mi@lists.uni-hamburg.de https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi