Opened 6 years ago

Last modified 3 months ago

#221 new defect - code aborts

compile grains fails with illegal extrapolation slope

Reported by: Gary J. Ferland Owned by: peter
Priority: critical Milestone:
Component: grains Version: trunk
Keywords: Cc:

Description

tests with the grain size distribution extended to 1 micron show that the compile grains code can fail with an error at r5877 of trunk and gcc 4.2.1. To reproduce use the slightly modified szd file and the compile grains file. The error is

 Starting size bin 10, amin=0.70423 amax=1.00000 micron
    absorption cs extrapolated below 9.9980e-06 Ryd
    absorption cs extrapolated above 7.3311e+01 Ryd
    scattering cs extrapolated below 9.9980e-06 Ryd
    scattering cs extrapolated above 7.3311e+01 Ryd
 Illegal value for slope in extrapolation 6.438356e-02
 [Stop in mie_repair at ../grains_mie.cpp:2629, something went wrong]

Attachments (2)

compile_grains.in (94 bytes) - added by Gary J. Ferland 6 years ago.
compile grains
orion.szd (268 bytes) - added by Gary J. Ferland 6 years ago.
extend size distribution to 1 micron causes abort

Download all attachments as: .zip

Change History (4)

Changed 6 years ago by Gary J. Ferland

Attachment: compile_grains.in added

compile grains

Changed 6 years ago by Gary J. Ferland

Attachment: orion.szd added

extend size distribution to 1 micron causes abort

comment:1 Changed 3 years ago by peter

Ticket #233 is a duplicate report with a different test case.

comment:2 Changed 3 months ago by peter

The illegal value for the slope is the result of the fact that the refractive index data do not extend to low or high enough frequencies. So the only solution is to extend the refractive index data, preferably over the entire frequency range supported by Cloudy. For the time being I do not have the tools to do that. It is also not clear how much that would extend the size range for which we can derive valid opacities as at some point the Mie code itself will start to fail...

Note: See TracTickets for help on using tickets.