polaris-mcrt / polaris Goto Github PK

Hello, I recently switched from Polaris Version 4.0.4 to 4.0.6. I am using an OctTree grid which I created from FLASH input data (but also a self-written test cube with the same result). With 4.0.4 the TEMP and RAT calculations work fine, but with 4.0.6 I get permanently the message "ERROR: Wrong cell border!" using the standard ISRF source. In the end about 80% of all the photons are "killed". This reduces to about 50% when I decrease the starting radius from "2" to "1". In contrast, with Version 4.0.4 the number of killed photons was marginal.
The read-in of the OctTree seems to work fine, all the data (length, mass, density etc) it reports in the POLARIS output are correct.
Do you have any idea what to look for or to check, as this seems to be relatively deeply embedded in the code? If needed I could provide a small example. Thanks.

Hi all,
recently me and independently also Michael Küffmeier encountered a problem with the latest Polaris version when trying to calculate the dust temperature. It seems to hang at the line
return phIDs[0];
in the routine getPhaseFunctionID in Parameters.h .

Michael gets the error shown in the image below (which I think is identical):

To my understanding this comes from the fact that phIDs is not an array but a scalar. We both downloaded the clean version and just tried to ran a simple TEMP calculation

.

Dear developers,
would it be possible to include the position of stars/point sources used for the TEMP/RAT calculation in the screen output of POLARIS for later cross-checking? If I am not mistaken, so far it is only contained in the command file and if that file is changed at a later point, this information is lost.
Best,
Daniel Seifried