Skip to content

normalise Mueller matrix error in 3D mode #172

Description

@jiema-astro

I encountered the following error when running basic mcfost para command:
Computing dust properties ... ERROR: normalise_Mueller_matrix: exact normalisation is smaller than numerical integration

Minimal setup:

#Grid geometry and size
  1                          
  100 90 100 10            

#Grain properties
 1                              
 DHS  1 1 0.  0.8  0.9        
 Draine_Si_sUV.dat 1.0          
 1                                 
 0.3 1000.0 3.5 50          

The error appears only in 3D mode. The same setup works in 2D mode with 100 90 1 10 in geometry, or using a smaller distribution of dust sizes 0.3 30.0 3.5 50 for example.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions