Ionization/Recombination rates for Ni_26, _27

[ionizationcalc]

I'm trying to plot out ionization/recombination rates for Ni. There are some 'ValueErrors' once temperatures are lower than 10^4K, but it works for other charge states such as Ni_25 or Ni_28. Here is an example:

Z = 28
z1 = 26
Te = numpy.logspace(3.7, 9, 53)
datacache = {}
ion, rec = pyatomdb.atomdb.get_ionrec_rate(Te, Z=Z, z1=z1, datacache=datacache)

The outputs are:

...
~/opt/anaconda3/lib/python3.7/site-packages/pyatomdb/atomdb.py in _calc_ionrec_dr(cidat, Te, extrap)
2605 ((cidat['par_type'] >= const.INTERP_IONREC_RATE_INC_MAX) &
2606 (cidat['par_type'] <= const.INTERP_IONREC_RATE_INC_MAX + const.MAX_IONREC))):
-> 2607 dr[idr] = _interpolate_ionrec_rate(cidat,Te[idr])
2608
2609 else:
ValueError: shape mismatch: value array of shape (4,1) could not be broadcast to indexing result of shape (4,)

[jagophile]

Hi,

Thanks for finding this error. This was something I should have caught earlier, and affected to any DR rates below the minimum value listed in the files.

I've fixed it now - please try downloading the newest version (0.10.7) from Github or pip and let me know if it works for you.

Thanks,
Adam

[ionizationcalc]

Hi Adam,

Thanks! It works for ionization rates. But I get NAN values for some recombination rates running on v0.10.7.

Te = np.logspace(3.7, 9, 13)
ion, rec = pyatomdb.atomdb.get_ionrec_rate(Te, Z=28, z1=26, datacache=datacache)
print(rec)
[ nan 4.95871026e-10 2.63527910e-10 1.37312280e-10
7.00233057e-11 3.46450748e-11 1.64201687e-11 7.45279043e-12
3.23254481e-12 1.65366426e-12 7.10469611e-13 2.25968988e-13
6.22427916e-14]

Chengcai