Tsteinmetz98 commented 9 months ago

Hello, I am running models using V4 using the following input file:

I. GEOMETRY %(available options: sphere\sphere_matrix\slab)

geometry = sphere

II. PHYSICAL PARAMETERS 1) Central source %(available options: on\off)

    central = on

1.1) Shape: %(available options: black_body\engelkd_marengo\power_law\file_lambda_f_lambda\file_f_lambda\file_f_nu)

        Spectral shape = black_body  
        Number of BB = 1
        Temperature = 3200 K

1.2) Scale: %(available options: flux\Lum_r1\energy_den\dilutn_fac\T1)

Scale = Lum_r1   % Td at the inner boundary
L = 8.2E3
d = 1.496e13 cm % 1 AU

2) External source %(available options: on\off)

external = off

3) Dust Properties

3.1 Chemical composition %(available options: common_grain_composite\common_and_addl_grain\tabulated)

        optical properties index = common_grain_composite
    Abundances for supported grain types:
        Sil-Ow  Sil-Oc  Sil-DL  grf-DL  amC-Hn  SiC-Pg 
   x =  0.50    0.50    0.00    0.00    0.00    0.00
    SIZE DISTRIBUTION = MRN
    Tsub = 3000.

4) Density Distribution %(available options: powd\expd\rdw\rdwa\usr_suppld) density type = RDW Y = 200

5) Optical Depth: %(available options: linear\logarithmic\user_supplied) grid type = logarithmic % log grid lambda0 = .55 micron % fiducial wavelength % minimal optical depth @ fiducial wavelength
tau(min) = 10.0
% maximal optical depth @ fiducial wavelength tau(max) = 100.0 number of models = 10

III. NUMERICS

accuracy for flux conservation = 0.10

The code runs without issue for the first 4 models, producing the desired output files. For the 5th model, I get the following error

model = 5, tau= 27.8256 going to Solve In Solve: 1 iteration over fbol Calculating with nY= 57 nP= 147 nCav= 38 Done with transmitted radiation. Done with initial dust temperature. ......................... Done with finding dust temperature after 25 iterations errT: 6.925E-04 errU: 8.345E-02 Done with finding energy density and diffuse flux. Done with radiative transfer. Achieved error in bolometric Flux: 15.178% After Flux_Cons nY= 86 In Solve: 2 iteration over fbol Calculating with nY= 86 nP= 238 nCav= 57 Done with transmitted radiation. Done with initial dust temperature. ............................. Done with finding dust temperature after 29 iterations errT: 3.943E-04 errU: 8.754E-02 Done with finding energy density and diffuse flux. Done with radiative transfer. Achieved error in bolometric Flux: 10.741% After Flux_Cons nY=129 In Solve: 3 iteration over fbol Calculating with nY= 129 nP= 321 nCav= 86 Done with transmitted radiation. Done with initial dust temperature. ......................................... Done with finding dust temperature after 41 iterations errT: 2.674E-04 errU: 8.450E-02 Done with finding energy density and diffuse flux. Done with radiative transfer. Achieved error in bolometric Flux: 12.794% After Flux_Cons nY=196 In Solve: 4 iteration over fbol Calculating with nY= 196 nP= 471 nCav= 100 Done with transmitted radiation. Done with initial dust temperature. ................................ Done with finding dust temperature after 32 iterations errT: 2.652E-04 errU: 7.574E-02 Done with finding energy density and diffuse flux. Done with radiative transfer. Achieved error in bolometric Flux: 10.570% After Flux_Cons nY=294 In Solve: 5 iteration over fbol Calculating with nY= 294 nP= 647 nCav= 100 Done with transmitted radiation. Done with initial dust temperature. ............................................... Done with finding dust temperature after 47 iterations errT: 1.510E-04 errU: 8.156E-02 Done with finding energy density and diffuse flux. Done with radiative transfer. Achieved error in bolometric Flux: 6.319% After Flux_Cons nY=294 1 . iteration over eta Doing Dynamics free(): invalid pointer

Program received signal SIGABRT: Process abort signal.

Backtrace for this error:

0 0x7f3bef621768 in ???

9 0x7f3bef37ae36 in free

10 0x55b4bec134ef in analysis_

11 0x55b4bec3a21d in solve_

12 0x55b4bec43c7c in kernel_

13 0x55b4bec56f89 in rundusty_

14 0x55b4bec57ca0 in MAIN__

15 0x55b4bebf839e in main

Abort (core dumped)

I couldn't find the exact files in the traceback (there doesn't seem to be a directory under https://github.com/ivezic/dusty/issues) so I cannot investigate myself what the issues are. Any help with this would be much appreciated!!!

ivezic commented 9 months ago

Hi, I am not sure when I will have time to take a detailed look (probably not before Dec 8). Meanwhile, can you please try to run version V2 (which is much more stable than V4) with the same inputs?

Tsteinmetz98 commented 8 months ago

hello! I ran the same properties, using the following parameters: I PHYSICAL PARAMETERS
1) External radiation:
Spectrum = 1
Number of BB = 1 Temperature = 3200 K

 2) Dust Properties 

   2.1 Chemical composition                       
       optical properties index = 1
       Abundances for supported grain types:
           Sil-Ow  Sil-Oc  Sil-DL  grf-DL  amC-Hn  SiC-Pg 
       x =  0.50    0.50    0.00    0.00    0.00    0.00

   2.2 Grain size distribution

    - size distribution = 1  % standard MRN                         

   2.3 Dust temperature on inner boundary:

    - temperature = 1200 K                         

 3) Density Distribution
    - density type = 3                    
    - Y = 1.e3 
    central star radius is 300 Rsun, and dust envelope is 10000 Rsun. assume inner radius is 5 R*, therefore Y = 7

 4) Optical Depth     
    - grid type = 1                    % linear grid
    - lambda0 = 10.0 micron            % optical depth specified 
    - tau(min) = 1.0; tau(max) = 50.0   % for the visual wavelength   
    - number of models = 50

II NUMERICS

 - accuracy for flux conservation = 0.05

This runs without issue. When I re-ran the previous input file but with lambda0=10.0 micron, the same error appeared at tau=3.

ivezic / dusty

Error in model runs #9