search

will-henney / proplyd-cloudy

Ionization/dynamic models of proplyds
http://deprecated.github.com/proplyd-cloudy/
2 stars 0 forks source link

Convergence failure with Adal Models #13

Closed nahieflores closed 12 years ago

nahieflores commented 12 years ago

El modelo usando:

Teff = 39000 atm = WMBasic r0 = 1.91 e 15 cm D = 1.5824 e 17 cm Q0 = 1 e 49 ph s-1 Phi(H) = 3.178 e13 ph s-1 cm-2

que es lo apropiado para p177-341, quebro pero ya no por un problema de programacion, sino al encontrar los limites para buscar la convergencia. Esto fue en el angulo de 67 grados, supongo que teniendo los otros 3 puedo tener un aproximado de los cocientes de lineas. Abajo el mensaje del modelo

Current brackets: None 5.90638155099 I-front radius = 0.9787, target = 0.9786 Running cloudy with x0 = -2.31, dR = -0.00839 Calculating proplyd-auto-models/WM039000-phi13.50-r15.28/th67/it15n1.211e+03x0-2.31 MainProcess finished with it15n1.211e+03x0-2.31 Current brackets: -2.31326285635 5.90638155099

Traceback (most recent call last): File "auto-models.py", line 377, in args=(cmdargs, extras, x0monitor)) File "/opt/epd/7.1/lib/python2.7/site-packages/scipy/optimize/zeros.py", line 364, in brentq r = _zeros._brentq(f,a,b,xtol,maxiter,args,full_output,disp) ValueError: xtol must be >= 0 28036.697u 52.481s 7:48:27.69 99.9% 0+0k 5336+169640io 1pf+0w

will-henney commented 12 years ago

This is clearly a bug in auto-models.py. Before it crashes, you have negative values for x0 and dR, which makes no sense. Perhaps the initial guess was too inaccurate. I will look into it.

Please update this issue with full details of the command line arguments that you used, and which triggered the error.

On Wed, Jan 25, 2012 at 8:09 PM, nahieflores reply@reply.github.com wrote:

El modelo usando:

 Teff = 39000  atm = WMBasic  r0 = 1.91 e 15 cm  D = 1.5824 e 17 cm Q0 = 1 e 49 ph s-1  Phi(H) = 3.178 e13 ph s-1 cm-2

que es lo apropiado para p177-341, quebro pero ya no por un problema de programacion, sino al encontrar los limites para buscar la convergencia. Esto fue en el angulo de 67 grados, supongo que teniendo los otros 3 puedo tener un aproximado de los cocientes de lineas. Abajo el mensaje del modelo

Current brackets:  None 5.90638155099 I-front radius = 0.9787, target = 0.9786 Running cloudy with x0 = -2.31, dR = -0.00839 Calculating proplyd-auto-models/WM039000-phi13.50-r15.28/th67/it15n1.211e+03x0-2.31 MainProcess  finished with  it15n1.211e+03x0-2.31 Current brackets:  -2.31326285635 5.90638155099

Traceback (most recent call last):  File "auto-models.py", line 377, in    args=(cmdargs, extras, x0monitor))  File "/opt/epd/7.1/lib/python2.7/site-packages/scipy/optimize/zeros.py", line 364, in brentq    r = _zeros._brentq(f,a,b,xtol,maxiter,args,full_output,disp) ValueError: xtol must be >= 0 28036.697u 52.481s 7:48:27.69 99.9%     0+0k 5336+169640io 1pf+0w

Reply to this email directly or view it on GitHub: https://github.com/deprecated/proplyd-cloudy/issues/13

  Dr William Henney, Centro de Radioastronomía y Astrofísica,   Universidad Nacional Autónoma de México, Campus Morelia

nahieflores commented 12 years ago

time python auto-models.py --cloudyexec /fs/astro04/other0/nahiely/Cloudy/c10.00/source/cloudy.exe --ntheta 5 --Tstar 39000 --logPhi 13.502 --r0 1.91e15

This is exactly the command line that I used

will-henney commented 12 years ago

I ran that exact model with current head of emission branch and it worked fine, although it took nearly 5 hours.

The th67 angle converged after 15 iterations, with the last few looking like this:

Current brackets:  None 7.01865865074
I-front radius = 0.9758, target = 0.9747
Running cloudy with x0 = 6.30, dR = 0.02267
Calculating proplyd-auto-models/WM039000-phi13.50-r15.28/th67/it13n1.219e+03x06.30
MainProcess  finished with  it13n1.219e+03x06.30
Current brackets:  None 6.29746557473
I-front radius = 0.9776, target = 0.9773
Running cloudy with x0 = 6.09, dR = 0.02191
Calculating proplyd-auto-models/WM039000-phi13.50-r15.28/th67/it14n1.219e+03x06.09
MainProcess  finished with  it14n1.219e+03x06.09
Current brackets:  None 6.0859128842
I-front radius = 0.9781, target = 0.9781
Running cloudy with x0 = 6.07, dR = 0.02186
Calculating proplyd-auto-models/WM039000-phi13.50-r15.28/th67/it15n1.219e+03x06.07
MainProcess  finished with  it15n1.219e+03x06.07
Current brackets:  6.07096184169 6.0859128842

It is strange that the density in my version 1.219e+03 is different from yours (1.211e+03). Another thing that worries me is that one of the brackets is always None. I will carry on looking into it.

will-henney commented 12 years ago

Since this works for me, I am closing the issue for now