Good that you emailed. There is a old and new system for PSDs. pnd_fieldCalcFromscat_speciesFields uses the old system. We have now decided that the old system will be removed. Has not yet happen due to lack of time. The new system supports retrievals, in contrast to the old one.

There could be a lack of a demo case for the new system.
The new PSDs should all be found in The new system uses pnd_fieldCalcFromParticleBulkProps

This is all I have time to write now. I can try to explain more carefully on Monday. I can also send you example cfiles.



On 2019-09-27 17:25, Fox, Stuart wrote:
Hi all,

I would like to add some further PSD options to ARTS (specifically ones for rain and graupel that are consistent with the single-moment schemes used in the Met Office NWP model). I will be making use of them by defining hydrometeor mass density fields and then using pnd_fieldCalcFromscat_speciesFields.

However, I’m a little bit confused as to the correct way to implement these (and there appears to be incomplete implementations of some of the existing options e.g. Abel & Boutle 2012 for rain, which happens to be one of the ones I’d like to use). The guidance in the arts developer guide also doesn’t seem to follow what’s actually in the code for some cases.

So far I have:

-added logic to pnd_fieldCalcFromscat_speciesFields to call pnd_fieldXXX for each of the new parametrizations

-updated the documentation for pnd_fieldCalcFromscat_speciesFields in to include the new parametrizations

-added new pnd_fieldXXX function to (not as suggested by the developer guide) to calculate the pnd field according to the raw PSD function psd_XXX and the scattering meta-data (and added this to microphysics.h as well)

-added new “raw” psd calculation psd_XXX to

This seems to be all that is required to make my use-case work, but I can see that it is not quite complete. In particular, I believe that I should add a new workspace method dNdD_XXX to allow a direct calculation of the raw PSD. Should this go in (and be added to This seems to be where other ones are, but again the developer guide suggests it should be in

What is the purpose of the psd_XXX functions in Are these also required?

Thanks for your help,


Dr Stuart Fox  Radiation Research Manager

*Met Office*FitzRoy Road  Exeter  Devon  EX1 3PB  United Kingdom
Tel: +44 (0)330 135 2480  Fax: +44 (0)1392 885681
Email:  Website:
See our guide to climate change at

arts_dev.mi mailing list

arts_dev.mi mailing list

Reply via email to