c ionconst.inc
       common /ionconst/s107,hallmin,pedmin,sigma_ratio,ped0,
     $     p_hardy,c_hardy,h_hardy,iter,itmax,errmax,
     $     ion_skip,looptries,ionmodel
	logical*1 ionmodel