igmk / pamtra

Passive and Active Microwave TRAnsfer model
GNU General Public License v3.0
19 stars 16 forks source link

Scattering function not normalized #17

Open mariomech opened 4 years ago

mariomech commented 4 years ago

In CHECK_NORM4 routine in file rad_scat4.f90 there is a check for the normalization of the scattering function. This fails for the passive calculations when using SSRGA and too large size parameters. The value usually checked for is 1.d-6. When setting it to 1.d-4 most calculations work fine.

This needs to be sorted out at some time.

DaveOri commented 4 years ago

I can take care of this. Can you provide a workable script to use for testing? Also if I remember correctly this is a problem also for mie-sphere, it might be related to the exansion in legendre coefficients for all rt3-style scattering routines