/*CoolMagn compute magnesium cooling */
#include "cddefines.h"
#include "coolheavy.h"
#include "poplevls.h"
#include "taulines.h"
#include "tepowers.h"
#include "phycon.h"
#include "mgexc.h"
#include "logte.h"
#include "ligbar.h"
#include "makecs.h"
#include "coladd.h"
#include "level2.h"
#include "putcs.h"
#include "level3.h"
#include "popexc.h"
#include "coolmetals.h"
#include "ionfracs.h"

void CoolMagn()
{
	float cs, 
	  csdum, 
	  csoi;
	double cs2s2p, 
	  cs2s3p;

#	ifdef DEBUG_FUN
	fputs( "<+>CoolMagn()\n", debug_fp );
#	endif

	/* Mg I 2853 
	 * fit to Dima's  integration of
	 * >>refer Leep, D., & Gallagher, A. 1976, Phys Rev A, 13, 148 */
	cs = (float)(5.21e-4*phycon.te/TePowers.te10);
	PutCS(cs,&TauLines[ipMgI2853]);
	level2(&TauLines[ipMgI2853]);

	/* Mg I 2026,  */
	MakeCS(&TauLines[ipMgI2026]);
	level2(&TauLines[ipMgI2026]);

	/* Mg 1 4571, data from Mendoza
	 * >>refer Mendoza, C. 1982, in Planetary Nebulae, IAU Symp No. 103,
	 * >>refer ed by D.R. Flower, (D. Reidel: Holland), 143
	 * cs set to O I */
	csoi = (float)(2.68e-5*phycon.te*(1. + 1.67e-6*phycon.te - 2.95e-10*phycon.te*
	  phycon.te));
	csoi = (float)MAX2(0.1,csoi);
	PutCS(csoi/9.,&TauLines[ipT4561]);
	level2(&TauLines[ipT4561]);

	/* Mg II 2798
	 * cs from 
	 * >>refer Sigut, A., & Pradhan, A.K., 1994, J Phys B sub
	 * refer not publ as of '97, cs agrees fairly well with Harrington et al.
	 * previous reference (~'82) */
	cs2s2p = 4.50*TePowers.te10;
	PutCS(cs2s2p,&TauLines[ipT2796]);
	PutCS(cs2s2p*0.5,&TauLines[ipT2804]);
	PutCS(1.0,&TauDummy);
	level3(&TauLines[ipT2804],&TauDummy,&TauLines[ipT2796]);
	/* call PutCS( cs , t2800 )
	 * call level2( t2800 )
	 *
	 * following used in MAGNES for photo destruction rate */
	if( PopLevls.PopLevels[0] > 0. )
	{
		mgexc.popmg2 = (PopLevls.PopLevels[2] + PopLevls.PopLevels[1])/
		  PopLevls.PopLevels[0];
	}
	else
	{
		mgexc.popmg2 = 0.;
	}

	/* MG IV 4.487 MIC 
	 * cs 
	 * >>refer Saraph, H.E. & Tully, J.A. 1994, A&AS, 107, 29 */
	cs = (float)MIN2(0.425,0.180*TePowers.te05*TePowers.te02);
	cs = (float)MAX2(0.356,cs);
	PutCS(cs,&TauLines[ipTMg4]);

	level2(&TauLines[ipTMg4]);

	/* MG V 5.61, 13.54 micron, 
	 * >>refer Butler, K., & Zeippen, C.J. 1994, A&AS, 108, 1
	 * >>chng 96 july 16 had been constant 0.3 */
	cs = (float)MIN2(0.311,0.11*TePowers.te10);
	PutCS(cs,&TauLines[ipTMg14]);
	cs = (float)MIN2(1.06,0.339*TePowers.te10);
	PutCS(cs,&TauLines[ipTMg6]);
	/* >>chng 96 july 16 had been constant 0.3 */
	cs = (float)MIN2(0.297,0.0745*TePowers.te10*TePowers.te02);
	PutCS(cs,&TauDummy);
	level3(&TauLines[ipTMg6],&TauLines[ipTMg14],&TauDummy);

	/* [Mg V] 2751+2893- Ne III-like, cs 
	 * >>refer Mendoza, C., & Zeippen, C.J. 1987, MNRAS, 224, 7p
	 * >>chng 96 aug 5 to three level atom
	 * c2751 = popexc(1.33,9.,5.,2.4,5.14e4,xmg(5))*7.11e-12
	 * dCooldT = dCooldT + c2751*5.14e4*tsq1
	 * call coladd( 'Mg 5' , 2751 , C2751 )
	 *
	 * following is 2-1 transition, both 2928 and 2783, 
	 * >>refer Butler, K., & Zeippen, C.J. 1994, A&AS, 108, 1 */
	PutCS(1.187,&TauLines[ipxMg52855]);

	cs = (float)MIN2(0.278,0.0171*TePowers.te20*TePowers.te05/
	  TePowers.te005/TePowers.te003);
	cs = (float)MAX2(0.182,cs);
	/* 3-2 transition, 2417.5 */
	PutCS(cs,&TauLines[ipxMg52417]);

	/* 3-1 transition, 1324.58 */
	PutCS(0.153,&TauLines[ipxMg51325]);

	level3(&TauLines[ipxMg52855],&TauLines[ipxMg52417],&TauLines[ipxMg51325]);

	/* Mg VI, 1806- OII like, data 
	 * >>refer Kafatos, M., & Lynch, J.P. 1980, ApJS, 42, 611 */
	CoolHeavy.c1806 = popexc(0.6,4.,10.,0.1,7.974e4,xIonFracs[6][11])*
	  1.11e-11;
	coladd("Mg 6",1806,CoolHeavy.c1806);

	/* [Mg VII] IR lines at 5.517 and 9.03 microns, 
	 * carbon-like, 
	 * >>refer Lennon, D.J. Burke, V.M. 1994, A&AS, 103, 273 */
	if( logte.alogte < 4.4 )
	{
		cs = (float)(0.027*TePowers.te30/TePowers.te03*TePowers.te003*TePowers.te001);
	}
	else
	{
		cs = 0.44f;
	}
	PutCS(cs,&TauLines[ipfsMg790]);
	if( logte.alogte < 4.6 )
	{
		cs = (float)(MIN2(1.456,0.0577*TePowers.te30*TePowers.te02/TePowers.te001/
		  TePowers.te001));
		csdum = (float)(8.275e-3*TePowers.sqrte/TePowers.te10/TePowers.te001);
	}
	else
	{
		cs = (float)(3.257/(TePowers.te05*TePowers.te02*TePowers.te003*TePowers.te003));
		csdum = (float)(1.456/(TePowers.te10*TePowers.te01*TePowers.te005));
	}
	PutCS(cs,&TauLines[ipfsMg755]);
	PutCS(csdum,&TauDummy);
	/* level3(  t10,t21,t20) */
	level3(&TauLines[ipfsMg790],&TauLines[ipfsMg755],&TauDummy);

	/* [mg vii] 2510, 2629
	 * c2596 = popexc( 1.7,9.,5.,10.,5.76e4,xIonFracs(12,7))*7.96e-12
	 * dCooldT = dCooldT + c2596 * 5.76e4*tsq1
	 * call coladd( 'Mg 7' , 2596 , C2596 )
	 *
	 * >>chng 96 aug 5, converted to 3 level atom */
	cs = (float)MIN2(0.22,0.3622/(TePowers.te05*TePowers.te02*TePowers.te003));
	PutCS(cs,&TauLines[ipxMg71190]);

	/* 2-1 transitions, 2509.2A+2629.1A together */
	cs = (float)MIN2(1.067,0.247*TePowers.te10*TePowers.te03*TePowers.te005);
	PutCS(cs,&TauLines[ipxMg72569]);

	/* 3-2 transition, 2261.5A */
	cs = (float)MIN2(0.542,3.863/(TePowers.te20*TePowers.te03*
	  TePowers.te01/TePowers.te003));
	cs = (float)MAX2(0.3735,cs);
	PutCS(cs,&TauLines[ipxMg72261]);

	level3(&TauLines[ipxMg72569],&TauLines[ipxMg72261],&TauLines[ipxMg71190]);
	/* level3(  t10,t21,t20)
	 *
	 * Mg VIII 3.03 micron, data from 
	 * >>refer Chandra, S. 1982, SoPh, 75, 133
	 * cs from 
	 * >>refer Zhang, H.L., Graziani, M., Pradhan, A.K. 1994, A&A, 283, 319 */
	PutCS(1.0,&TauLines[ipxMg08303]);
	level2(&TauLines[ipxMg08303]);
	/* fs303 = popexc(0.26,2.,4.,0.324,4752.,xmg(8))*6.58E-13
	 * call coladd( 'Mg 8' , 303 , FS303 )
	 *
	 * Mg IX 704.5, 1909-like, A from
	 * >>refer Muhlethaler, H.P., & Nussbaumer, H. 1976, A&A 48, 109
	 * beseq line, cs data from 
	 * >>refer Keenan, F.P. Berrington, K.A., Burke, P.G., Dufton, P.L.,
	 * >>refer Kingston, A.E. 1986, PhyS 34, 216
	 * A's 
	 * >>refer Fleming, J., Bell, K.L, Hibbert, A., Vaeck, N., Godefroid, M.R.
	 * >>refer 1996, MNRAS, 279, 1289
	 * 666 error! use beseq here */
	cs = (float)(0.98288 - 0.23766*logte.alogte + 0.014334*POW2(logte.alogte));
	cs = (float)MAX2(0.01,cs);
	PutCS(cs,&TauLines[ipT705]);
	level2(&TauLines[ipT705]);

	/* Mg X 610
	 * >>refer Cochrane, D.M., & McWhirter, R.W.P. 1983, PhyS, 28, 25 */
	ligbar(12,&TauLines[ipTMg610],&TauLines[ipT58],&cs2s2p,&cs2s3p);
	PutCS(cs2s2p,&TauLines[ipTMg610]);
	PutCS(cs2s2p*0.5,&TauLines[ipTMg625]);
	PutCS(1.0,&TauDummy);
	level3(&TauLines[ipTMg625],&TauDummy,&TauLines[ipTMg610]);

	/* Mg X 58A */
	PutCS(cs2s3p,&TauLines[ipT58]);
	level2(&TauLines[ipT58]);


#	ifdef DEBUG_FUN
	fputs( " <->CoolMagn()\n", debug_fp );
#	endif
	return;
}