double free or corruption (fasttop): 0x0000000002faf5c0 in atomic line transfer

[Fitz-Hu]

mcfost ref3.0_3D.para -phantom tde_02000 ref3.0_3D.para.gz tde_02000.gz

You are running MCFOST 3.0.42
 Git SHA = e387e13c5de244c46e7f291ded95bec248a53c48
 WARNING: polarization is turned off in ray-traced SEDs
 it can be turned back on with -rt2
 Input file read successfully
 Thermal equilibrium calculation

 Temperature calculation under LTE approximation
 Parallelized code on   4 processors

Thu Oct 20 11:43:46 AEDT 2022
 Creating directory ././data_th
 Forcing 3D mode
 Using ray-tracing method 1
 Reading /fred/oz015/mcfost/utils/Dust/Draine_Si_sUV.dat
 Number of regions detected: 1
 zone 1 --> region= 1 : R=1.00 to 300.00 AU
 Performing phantom2mcfost setup
 Reading phantom density file: tde_02000
 ---- Reading header file #           1
 ERROR: could not find ndusttypes in header
 ---- Done
 ---- Reading data file #           1
 FT:Phantom:2022.0.1:0a99d82 (hydro+grav): 10/10/2022 13:58:01.9
 ERROR: could not find ndusttypes in header
  npart =       100000  ntypes =            8  ndusttypes =            0
  nptmass =            0
 ERROR: could not find gmw in header
 ---- Done
 Using temperature from phantom temperature writeout
 Found           0 point masses in the phantom file

 *** Using gas particles for Voronoi mesh ***           1

 SPH gas mass is    0.7447277     Msun
 SPH dust mass is   0.0000000E+00 Msun

 # Sink particles:

 Updating the stellar properties:
 There are now           0 stars in the model
 Using    74473 particles from Phantom file
 Phantom dump file processed ok
 Done

 # Farthest particules :
 x =  -9133.89212133848        19389.1208883494     
 y =  -10111.8746346924        18174.7031468849     
 z =  -18241.0836476695        16649.5387915112     
 Found       74473  SPH particles with            0 dust grains
 Selecting spatial range which contains
   99.90000     % of particles in each dimension
 # Model limits :
 x =  -6299.88378906250        8234.03125000000     
 y =  -5732.88867187500        11012.8945312500     
 z =  -7345.58398437500        9002.25292968750     
 Finding            6 walls
 Performing Voronoi tesselation on        74264 SPH particles
 mcfost will require ~   3.116241     MB of temporary memory for the tesselation
 Trying to allocate   4.725174     MB for neighbours list
 sha1 has changed
 Tesselation needs to be recomputed, deallocating neighbours list
 Trying to allocate   8.498840     MB for temporary neighbours list
 100% |==================================================| 
 Trying to allocate   4.374683     MB for neighbours list
 Freeing temporary neighbours list
 Tesselation done
 Building the kd-trees for the model walls
 Tesselation Time =  1.96s
 Found       74264 cells
 Neighbours list size =     1146797
 Average number of neighbours =   15.44217    
 Voronoi volume =   3897143833209.67     
 Using n_cells =       74264
 Using gas-to-dust ratio in mcfost parameter file
 Total  gas mass in model :  0.7426378      Msun
 Total dust mass in model :  7.4263779E-03  Msun
 Using       74264 cells
 Using a Voronoi mesh
 Using scattering method 2
 Reading dust properties from previous run ... Done
 lambda =  0.871355626679430     
 Integ tau dans plan eq. =    9.150955    
  Column density (g/cm?)   =   2.0694783E-02
 Integ tau (i =75.0 deg)   =  0.79326E+01
  Column density (g/cm?)   =   1.7939512E-02
 No dark zone
 Initializing thermal properties ... Done
 Error: wl #           4  No energy
 Error: wl #           2  No energy
 Error: wl #           3  No energy
 Error: wl #           1  No energy
 Exiting
 Exiting
 Exiting
 Exiting
*** Error in `/fred/oz015/mcfost/mcfost': double free or corruption (fasttop): 0x0000000002faf5c0 ***
*** Error in `/fred/oz015/mcfost/mcfost': double free or corruption (!prev): 0x0000000002fc34a0 ***
======= Backtrace: =========
[0x17ebc6c]
[0x17bcb00]
[0x17bcb25]
[0x14edfdb]
[0x8e8915]
[0xc9b6b6]
[0x1710ba3]
[0x16c4700]
[0x16c5cf4]
[0x1695964]
[0xc9d113]
[0xd35116]
[0x4077ce]
[0x17b1e54]
[0x17b20d1]
[0x4076a6]
======= Memory map: ========
00400000-01cfb000 r-xp 00000000 7e2:ef5ac 144121794574419632             /fred/oz015/mcfost/mcfost
01efa000-01f58000 rw-p 018fa000 7e2:ef5ac 144121794574419632             /fred/oz015/mcfost/mcfost
01f58000-0219a000 rw-p 00000000 00:00 0 
02f9c000-03163000 rw-p 00000000 00:00 0                                  [heap]
14f6b0000000-14f6b0050000 rw-p 00000000 00:00 0 
14f6b0050000-14f6b4000000 ---p 00000000 00:00 0 
14f6b54fe000-14f6b8000000 rw-p 00000000 00:00 0 
14f6b8000000-14f6b8960000 rw-p 00000000 00:00 0 
14f6b8960000-14f6bc000000 ---p 00000000 00:00 0 
14f6bc67d000-14f6c0000000 rw-p 00000000 00:00 0 
14f6c0000000-14f6c09e6000 rw-p 00000000 00:00 0 
14f6c09e6000-14f6c4000000 ---p 00000000 00:00 0 
14f6c4855000-14f6c777d000 rw-p 00000000 00:00 0 
14f6c7914000-14f6c7ffe000 rw-p 00000000 00:00 0 
14f6c7ffe000-14f6c7fff000 ---p 00000000 00:00 0 
14f6c7fff000-14f6e8000000 rw-p 00000000 00:00 0 
14f6e8000000-14f6e89e6000 rw-p 00000000 00:00 0 
14f6e89e6000-14f6ec000000 ---p 00000000 00:00 0 
14f6ec0e0000-14f6ed8d8000 rw-p 00000000 00:00 0 
14f6edaa7000-14f6edaa8000 rw-p 00000000 00:00 0 
14f6edab2000-14f6edab3000 rw-p 00000000 00:00 0 
14f6edab4000-14f6edef6000 rw-p 00000000 00:00 0 
14f6edef6000-14f6edef7000 ---p 00000000 00:00 0 
14f6edef7000-14f70def8000 rw-p 00000000 00:00 0 
14f70def8000-14f70def9000 ---p 00000000 00:00 0 
14f70def9000-14f72e14a000 rw-p 00000000 00:00 0 
7ffcead13000-7ffceadaa000 rw-p 00000000 00:00 0                          [stack]
7ffceadc2000-7ffceadc6000 r--p 00000000 00:00 0                          [vvar]
7ffceadc6000-7ffceadc8000 r-xp 00000000 00:00 0                          [vdso]
ffffffffff600000-ffffffffff601000 r-xp 00000000 00:00 0                  [vsyscall]
*** longjmp causes uninitialized stack frame ***: /fred/oz015/mcfost/mcfost terminated
======= Backtrace: =========
[0x17ebc6c]
[0x15275bb]
[0x14ee1ca]
[0x14edfd3]
[0x8e8915]
[0xc9b6b6]
[0x1710ba3]
[0x16c4700]
[0x16c395a]
[0x1710f59]
[0x17a9395]

[cpinte]

Hi Fitz, There is no sink particle, nor potential in your dump, so mcfost does not see a star, and cannot run a dust continuum calculation.

[Fitz-Hu]

How should i deal with this then? Does this mean mcfost cannot run with these dumps?

[cpinte]

Not really a bug. Models with no sink need to be run with -atom