search

cloudy-astrophysics / bug-tracker-migration-test

Trial run for importing the nublado.org Trac tickets as GitHub issues
0 stars 0 forks source link

grain convergence failures (trac #330) #332

Closed cloudy-bot closed 4 years ago

cloudy-bot commented 9 years ago

reported by: @CloudyLex

In both c13_branch and the trunk the following sim aborts before the first zone:

########### shape of considered galaxy
table AGN T=1.5e5 k, a(ox)=-1.4, a(uv)=-0.5, a(x)=-1
Q(h) 48 % vary
#
########### gas
hden 15
metallicity -1
abundance ism, no qheat
radius 22 cm        % 2kpc from qso
#
########### grid
iterate
#
stop zone 1

Typical output follows:

 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.68
 PROBLEM  ConvFail 18, a grain failure occurred iteration 1 zone 0 fnzone  13.68 NONE!!!!! 0 0
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.69
 PROBLEM  ConvFail 19, a grain failure occurred iteration 1 zone 0 fnzone  13.69 NONE!!!!! 0 0
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  ConvFail 20, a grain failure occurred iteration 1 zone 0 fnzone  13.70 NONE!!!!! 0 0
 Stop due to excessive convergence failures - there have been 20 so far. 
 This limit can be reset with the FAILURES command.
 ConvFail sets lgAbort since nTotalFailures=20 is >= LimFail=20
 This limit can be reset with the FAILURES command.
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70

Migrated from https://www.nublado.org/ticket/330

{
    "status": "closed",
    "changetime": "2019-02-04T12:08:14Z",
    "_ts": "1549282094687722",
    "description": "In both c13_branch and the trunk the following sim aborts before the first zone:\n{{{\n########### shape of considered galaxy\ntable AGN T=1.5e5 k, a(ox)=-1.4, a(uv)=-0.5, a(x)=-1\nQ(h) 48 % vary\n#\n########### gas\nhden 15\nmetallicity -1\nabundance ism, no qheat\nradius 22 cm        % 2kpc from qso\n#\n########### grid\niterate\n#\nstop zone 1\n}}}\n\nTypical output follows:\n{{{\n PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.68\n PROBLEM  ConvFail 18, a grain failure occurred iteration 1 zone 0 fnzone  13.68 NONE!!!!! 0 0\n PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.69\n PROBLEM  ConvFail 19, a grain failure occurred iteration 1 zone 0 fnzone  13.69 NONE!!!!! 0 0\n PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70\n PROBLEM  ConvFail 20, a grain failure occurred iteration 1 zone 0 fnzone  13.70 NONE!!!!! 0 0\n Stop due to excessive convergence failures - there have been 20 so far. \n This limit can be reset with the FAILURES command.\n ConvFail sets lgAbort since nTotalFailures=20 is >= LimFail=20\n This limit can be reset with the FAILURES command.\n PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70\n PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70\n PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70\n PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70\n PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70\n PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70\n}}}\n",
    "reporter": "gary",
    "cc": "",
    "resolution": "wontfix",
    "time": "2015-06-30T20:26:55Z",
    "component": "grains",
    "summary": "grain convergence failures",
    "priority": "blocker",
    "keywords": "",
    "version": "trunk",
    "milestone": "c17 branch",
    "owner": "peter",
    "type": "defect - convergence"
}
cloudy-bot commented 9 years ago

Attachment: projectV1.in

cloudy-bot commented 9 years ago

@CloudyLex changed milestone from "" to "C15_branch"

cloudy-bot commented 9 years ago

@peter-van-hoof-noaccount commented:

This sim is very pathological. The spectrum peaks at ~10 micron. It has an extremely low ionization parameter:

U(1.0----):2.654E-23 U(4.0----):3.042E-24

When I disable grains, it fails on lgNucleiConserved().

When I subsequently disable molecules, it fails on Te dropping below 2.8 K at the illuminated face.

I wonder if it is worthwhile fixing such an extreme model.

cloudy-bot commented 9 years ago

In both c13_branch and the trunk the following sim aborts before the first zone:

########### shape of considered galaxy
table AGN T=1.5e5 k, a(ox)=-1.4, a(uv)=-0.5, a(x)=-1
Q(h) 48 % vary
#
########### gas
hden 15
metallicity -1
abundance ism, no qheat
radius 22 cm        % 2kpc from qso
#
########### grid
iterate
#
stop zone 1

Typical output follows:

 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.68
 PROBLEM  ConvFail 18, a grain failure occurred iteration 1 zone 0 fnzone  13.68 NONE!!!!! 0 0
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.69
 PROBLEM  ConvFail 19, a grain failure occurred iteration 1 zone 0 fnzone  13.69 NONE!!!!! 0 0
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  ConvFail 20, a grain failure occurred iteration 1 zone 0 fnzone  13.70 NONE!!!!! 0 0
 Stop due to excessive convergence failures - there have been 20 so far. 
 This limit can be reset with the FAILURES command.
 ConvFail sets lgAbort since nTotalFailures=20 is >= LimFail=20
 This limit can be reset with the FAILURES command.
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70

to:

In both c13_branch and the trunk the following sim aborts before the first zone:

########### shape of considered galaxy
table AGN T=1.5e5 k, a(ox)=-1.4, a(uv)=-0.5, a(x)=-1
Q(h) 48 % vary
#
########### gas
hden 15
metallicity -1
abundance ism, no qheat
radius 22 cm        % 2kpc from qso
#
########### grid
iterate
#
stop zone 1

Typical output follows:

 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.68
 PROBLEM  ConvFail 18, a grain failure occurred iteration 1 zone 0 fnzone  13.68 NONE!!!!! 0 0
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.69
 PROBLEM  ConvFail 19, a grain failure occurred iteration 1 zone 0 fnzone  13.69 NONE!!!!! 0 0
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  ConvFail 20, a grain failure occurred iteration 1 zone 0 fnzone  13.70 NONE!!!!! 0 0
 Stop due to excessive convergence failures - there have been 20 so far. 
 This limit can be reset with the FAILURES command.
 ConvFail sets lgAbort since nTotalFailures=20 is >= LimFail=20
 This limit can be reset with the FAILURES command.
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
 PROBLEM  charge/temperature not converged for  sil-ism04 zone 13.70
cloudy-bot commented 9 years ago

Milestone renamed