Opened 2 years ago

#351 new defect - convergence

Database H2, PDR Convergence failure at low density, high ionization limit

Reported by: mchatzikos Owned by: nobody
Priority: major Milestone: C17_branch
Component: ionization convergence Version: trunk
Keywords: Cc:

Description

The sim:

#
# commands controlling continuum
#
cmb redshift 3
table ism
table Trapezium
intensity -1.8, range 0.44 to 1 Ryd
extinguish 21
#
# commands for density & abundances
#
hden 6
abundances ISM
database H2
grains PAH 3
cosmic ray background
#
# other commands for details
#
iterate
stop temperature off
stop Av 10

fails to converge as given, producing errors such as:

 PROBLEM  ConvFail 1,  ionization not converged iteration 1 zone 931 fnzone 954.05 reason change SO BadConvIoniz0:0.0016665 [1]=0.00323345
 PROBLEM  ConvFail 2,  ionization not converged iteration 1 zone 932 fnzone 943.53 reason change C3 BadConvIoniz0:0.000369811 [1]=0.000211395
 PROBLEM  ConvFail 3,  ionization not converged iteration 1 zone 933 fnzone 934.99 reason change C3 BadConvIoniz0:0.000468988 [1]=0.00020337
 PROBLEM  ConvFail 4,  ionization not converged iteration 1 zone 934 fnzone 936.49 reason change SO BadConvIoniz0:0.00232437 [1]=0.00472749
 PROBLEM  ConvFail 5,  ionization not converged iteration 1 zone 935 fnzone 936.99 reason edn mole-grn BadConvIoniz0:0.0431219 [1]=0.0440445
 PROBLEM  ConvFail 6,  ionization not converged iteration 1 zone 936 fnzone 937.49 reason change SO BadConvIoniz0:0.00204379 [1]=0.00332147
 PROBLEM  ConvFail 7,  ionization not converged iteration 1 zone 937 fnzone 937.99 reason change C3 BadConvIoniz0:0.000833294 [1]=0.000353221
 PROBLEM  ConvFail 8,  ionization not converged iteration 1 zone 938 fnzone 939.99 reason change C3 BadConvIoniz0:0.000877815 [1]=0.00038686
 PROBLEM  ConvFail 9,  ionization not converged iteration 1 zone 939 fnzone 940.99 reason change C3 BadConvIoniz0:0.000912609 [1]=0.000414018
 PROBLEM  ConvFail 10,  ionization not converged iteration 1 zone 940 fnzone 940.99 reason change C3 BadConvIoniz0:0.000891184 [1]=0.000408287
 PROBLEM  ConvFail 11,  ionization not converged iteration 1 zone 941 fnzone 942.99 reason change C3 BadConvIoniz0:0.00102985 [1]=0.000506269
 PROBLEM  ConvFail 12,  ionization not converged iteration 1 zone 942 fnzone 943.49 reason change C3 BadConvIoniz0:0.00105679 [1]=0.000534282
 PROBLEM  ConvFail 13,  ionization not converged iteration 1 zone 943 fnzone 943.99 reason change C3 BadConvIoniz0:0.00113576 [1]=0.000595867
 PROBLEM  ConvFail 14,  ionization not converged iteration 1 zone 944 fnzone 947.49 reason change C3 BadConvIoniz0:0.00119276 [1]=0.000642296
 PROBLEM  ConvFail 15,  ionization not converged iteration 1 zone 945 fnzone 946.49 reason change C3 BadConvIoniz0:0.00124662 [1]=0.000709808
 PROBLEM  ConvFail 16,  ionization not converged iteration 1 zone 946 fnzone 949.99 reason edn mole-grn BadConvIoniz0:0.0416553 [1]=0.0440757
 PROBLEM  ConvFail 17,  ionization not converged iteration 1 zone 947 fnzone 954.49 reason edn mole-grn BadConvIoniz0:0.0402359 [1]=0.0441054
 PROBLEM  ConvFail 18,  ionization not converged iteration 1 zone 948 fnzone 949.49 reason change C3 BadConvIoniz0:0.00174068 [1]=0.0011573
 PROBLEM  ConvFail 19,  ionization not converged iteration 1 zone 949 fnzone 949.99 reason change C3 BadConvIoniz0:0.00191937 [1]=0.00133476
 PROBLEM  ConvFail 20,  ionization not converged iteration 1 zone 950 fnzone 952.99 reason change C3 BadConvIoniz0:0.00193357 [1]=0.00139257
 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.

However, when the command 'database H2' is commented out, Cloudy exits normally.

Change History (0)

Note: See TracTickets for help on using tickets.