brittonsmith / cloudy_cooling_tools

Create an N-dimensional grid of models with the Cloudy photoionization code.
MIT License
5 stars 4 forks source link

New cloudy #1

Closed loikki closed 5 years ago

loikki commented 5 years ago

Hi,

I just want to report a bug with the last version of Cloudy (a previous one works fine).

When running the example sample_cooling.par, I obtain the error given at the end of the message. I tried with the version c13.05 and it seems to be fine.

I am not a Cloudy user, therefore I have no idea if the bug is due to this repository or Cloudy.

Best, Loic

                                                          Cloudy 17.01
                                                         www.nublado.org

                       **************************************17Jun01**************************************
                       *                                                                                 *
                       * stop zone 1                                                                     *
                       * iterate to convergence                                                          *
                       * hden -10                                                                        *
                       * coronal equilibrium T = 10.000000 linear                                        *
                       * punch last cooling file = "./cooling_run_run1.cooling.temp"                     *
                       * punch last heating file = "./cooling_run_run1.heating.temp"                     *
                       * punch last abundance file = "./cooling_run_run1.abundance.temp"                 *
                       * punch last ionization means file = "./cooling_run_run1.ionization.temp"         *
                       * punch last physical conditions file = "./cooling_run_run1.physical.temp"        *
                       *                                                                                 *
                       ***********************************************************************************

 NOTE
 NOTE The simulation is going into neutral gas but cosmic rays are not included.
 NOTE Ion-molecule chemistry will not occur without a source of ionization.
 NOTE The chemistry network may collapse deep in molecular regions.
 NOTE Consider adding galactic background cosmic rays with the COSMIC RAYS BACKGROUND command.
 NOTE You may be making a BIG mistake.
 NOTE

 NOTE Is the entered value of the hydrogen density (1.00e-10) reasonable?
 NOTE It seems pretty low to me.

 NOTE Simulation may crash because of extreme density.  The value was 1.00e-10

 NOTE Setcon: continuum has zero intensity starting at  3.0447e-09 Ryd.
  8228 cells in the incident continuum have zero intensity.  Problems???

 NOTE There is no hydrogen-ionizing radiation.
 Was this intended?

 NOTE There is no Balmer continuum radiation.<<<<
 Was this intended?

           8228CellPeak3.04E-09   Lo 3.04e-09=  29.98m   Hi-Con:7.13E+06 Ryd   E(hi):7.35E+06Ryd     E(hi):     100.00 MeV
           I(nu>1ryd):   0.0000   Average nu:0.000E+00   I( X-ray):   0.0000   I(BalC): -30.0000     Phi(BalmrC): -30.0000
           phi(1.0-1.8): 0.0000   phi(1.8-4.0):  0.000   phi(4.0-20):  0.000   phi(20--):  0.000     Ion pht flx:0.000E+00
           I(gam ray):   0.0000   phi(gam r):   0.0000   I(Infred):   0.0000   Alf(ox):   0.0000     Total inten: -35.9298
           U(1.0----):0.000E+00   U(4.0----):0.000E+00   T(En-Den):0.000E+00   T(Comp):0.000E+00     nuJnu(912A):0.000E+00
           Occ(FarIR):0.000E+00   Occ(H n=6):0.000E+00   Occ(1Ryd):0.000E+00   Occ(4R):0.000E+00     Occ (Nu-hi):0.000E+00
           Tbr(FarIR):0.000E+00   Tbr(H n=6):0.000E+00   Tbr(1Ryd):0.000E+00   Tbr(4R):0.000E+00     Tbr (Nu-hi):0.000E+00

 PROBLEM Network conservation error Li H 9.3113e-55 6.24476e-09 2.04e-19 0

 PROBLEM DISASTER
 An assert has been thrown, this is bad.
 Failed: lgNucleiConserved(c)
 It happened in the file mole_eval_balance.cpp at line number 141
 This is iteration 1, nzone 0, fzone 0.26, lgSearch=T.

           vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv 
          > PROBLEM DISASTER PROBLEM DISASTER.      <
          > Sorry, something bad has happened.      <
          > Please post this on the Cloudy web site <
          > discussion board at www.nublado.org     <
          > Please send all following information:  <
           ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 

 Cloudy version number is 17.01
 Cloudy compiled on Sep 17 2019 in OS Linux (AMD64) using the g++ 90100 compiler. Mode I32LP64, denormalized float: F double: F.

    0 warnings,  0 cautions,  0 temperature failures.  Messages follow.
 Input commands follow:
c ======================
stop zone 1
iterate to convergence
hden -10
coronal equilibrium T = 10.000000 linear
punch last cooling file = "./cooling_run_run1.cooling.temp"
punch last heating file = "./cooling_run_run1.heating.temp"
punch last abundance file = "./cooling_run_run1.abundance.temp"
punch last ionization means file = "./cooling_run_run1.ionization.temp"
punch last physical conditions file = "./cooling_run_run1.physical.temp"
c ======================
brittonsmith commented 5 years ago

Hi @loikki, this is an issue with Cloudy I'm afraid. Basically, the code is crashing because the conditions are too extreme, i.e., there is no radiation strong enough to affect the gas and the density is so low that nothing is really happening. It's likely that later versions of Cloudy have simply tightened the tolerances for detecting these sorts of issues. One way to get around this would be to increase the minimum density a few orders of magnitude. This would probably be fine since it's quite low.

If you want, you can also post these conditions on the Cloudy mailing list and ask if there is a way to keep the code running, but it's quite likely they're going to tell you that the conditions are unrealistic.

loikki commented 5 years ago

Thanks, I will reach the cloudy community