/*CoolAlum compute aluminum cooling */
#include "cddefines.h"
#include "poplevls.h"
#include "embesq.h"
#include "tepowers.h"
#include "taulines.h"
#include "logte.h"
#include "coolheavy.h"
#include "ligbar.h"
#include "makecs.h"
#include "coladd.h"
#include "level2.h"
#include "putcs.h"
#include "beseq.h"
#include "level3.h"
#include "pop3.h"
#include "ionfracs.h"
#include "coolmetals.h"

void CoolAlum()
{
	double cs ,
	  cs2s2p, 
	 cs2s3p ;
	float p2;

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

	/* Al I 3957 */
	MakeCS( &TauLines[ipAlI3957] );
	level2( &TauLines[ipAlI3957]);

	/* Al I 3090 */
	MakeCS(&TauLines[ipAlI3090]);
	level2(&TauLines[ipAlI3090]);

	/* AlII  1670.787
	 * >>chng 96 may 09 put back into level 2 for array processor speed up
	 * cs from 
	 * >>refer Tayal, S.S., Burke, P.G., Kingston, A.E. 1985, J.Phys. B, 18, 4321
	 * cs = MIN( 5.0 , 0.0125 * sqrte*te10*te003 )
	 * call PutCS( cs , al1671 )
	 * call level2( al1671 )
	 *
	 * Al 2, all data summarized in 
	 * >>refer Keenan, F.P., Harra, L.K., Aggarwal, K.M., Feibelman, W.A. 1992,
	 * >>refer ApJ, 385, 375
	 * doublet at 2660, 2669 */
	PutCS(3.56,&TauLines[ipT2670]);

	/* C2670 = beseq( 1.67,2.00,6.54, T2670 , 3.67E-3 )  * 7.45E-12 */
	beseq(1.67,2.00,6.54,&TauLines[ipT2670],3.67e-3);
	embesq.em2669 = (float)(PopLevls.PopLevels[3]*3.67e-3*7.45e-12);

	/* Aluminum al 3, 1854, 1862 doublet, 3s ^2 S gnd, ^2P^o 1/2 3/5 exc
	 * f=0.854 from 
	 * >>refer Dufton, P.L., Brown, P.J.F., Lennon, D.J., Lynas-Gray, A.E. 1986,
	 * >>refer MNRAS, 222, 713
	 * cs from 
	 * >>refer Dufton, P.L., & Kingston, A.E. 1987, J.Phys. B, 20, 3899 */
	cs = 4.407*TePowers.te10*TePowers.te03*TePowers.te01;
	cs = MIN2(25.0,cs);
	PutCS(cs*0.667,&TauLines[ipT1855]);
	PutCS(cs*0.333,&TauLines[ipT1863]);
	PutCS(1.0,&TauDummy);
	level3(&TauLines[ipT1863],&TauDummy,&TauLines[ipT1855]);

	/* Al V 2.905, cs 
	 * >>refer Saraph, H.E. & Tully, J.A. 1994, A&AS, 107, 29 */
	cs = MIN2(0.524,1.113/(TePowers.te10/TePowers.te02/TePowers.te003));
	PutCS(cs,&TauLines[ipAl529]);

	level2(&TauLines[ipAl529]);

	/* Al VI 3.66, 9.12 microns */
	cs = 639.1/(TePowers.sqrte*pow(TePowers.te03,TePowers.te003)*TePowers.te001);
	cs = MIN2(5.5 , cs);
	PutCS(cs,&TauLines[ipAl6366]);

	cs = MIN2(1.10,49.37/(TePowers.sqrte/TePowers.te10*TePowers.te02/
	  TePowers.te001));
	PutCS(cs,&TauLines[ipAl6912]);

	cs = MIN2(2.0,319.11/(TePowers.sqrte*TePowers.te10/TePowers.te02/
	  TePowers.te001));
	PutCS(cs,&TauDummy);

	level3(&TauLines[ipAl6366],&TauLines[ipAl6912],&TauDummy);

	/* [Al VI] 2428.4, 2601, 1169.8, 2124.9
	 * POP3(G1,G2,G3,O12,O13,O23,A21,A31,A32,E12,E23,P2,ABUND,GAM2)
	 * cs from 
	 * >>refer Butler, K., & Zeippen, C.J. 1994, A&AS, 108, 1
	 * >>refer As from NIST */
	CoolHeavy.c1170 = pop3(9.,5.,1.,1.044,0.145,0.463,6.63,72.9,7.79,
	  5.92e4,6.767e4,&p2,xIonFracs[6][12],0.)*72.9*1.70e-11;

	CoolHeavy.c2428 = p2*6.63*8.19e-12;
	CoolHeavy.c2125 = CoolHeavy.c1170*(7.79/72.9)*(1169.5/2124.9);
	coladd("Al 6",1170,CoolHeavy.c1170);
	coladd("Al 6",2428,CoolHeavy.c2428);
	coladd("Al 6",2125,CoolHeavy.c2125);

	/* Al VIII 5.85, 3.72 microns
	 * collision strength 
	 * >>refer Lennon, D.J. Burke, V.M. 1994, A&AS, 103, 273 */
	cs = MIN2(0.39,0.0459*TePowers.te20/TePowers.te003/TePowers.te003);
	PutCS(cs,&TauLines[ipAl8575]);
	cs = MIN2(1.062,0.0407*TePowers.te30/TePowers.te003/TePowers.te003);
	PutCS(cs,&TauLines[ipAl8370]);
	cs = MIN2(0.27,2.694e-3*TePowers.te20*TePowers.te20*TePowers.te01*
	  TePowers.te003);
	PutCS(cs,&TauDummy);
	level3(&TauLines[ipAl8575],&TauLines[ipAl8370],&TauDummy);

	/* [Al IX] 2.04 micron, no collision strength, A NIST */
	PutCS(1.,&TauLines[ipAl09204]);
	level2(&TauLines[ipAl09204]);

	/* Al 10 639, CS 
	 * >>refer Keenan, F.P. Berrington, K.A., Burke, P.G., Dufton, P.L.,
	 * >>refer Kingston, A.E. 1986, PhyS 34, 216
	 * A is extrapolation along iso seq */
	cs = 0.73492 - 0.16964*logte.alogte + 0.0096631*POW2(logte.alogte);
	cs = MAX2(0.01,cs);
	PutCS(cs,&TauLines[ipT639]);
	level2(&TauLines[ipT639]);

	/* Al 11 Li seq 2s2p 556
	 * >>refer Cochrane, D.M., & McWhirter, R.W.P. 1983, PhyS, 28, 25 */
	ligbar(13,&TauLines[ipTAl550],&TauLines[ipTAl48],&cs2s2p,&cs2s3p);
	PutCS(cs2s2p,&TauLines[ipTAl550]);
	PutCS(cs2s2p*0.5,&TauLines[ipTAl568]);
	PutCS(1.0,&TauDummy);
	level3(&TauLines[ipTAl568],&TauDummy,&TauLines[ipTAl550]);

	PutCS(cs2s3p,&TauLines[ipTAl48]);
	level2(&TauLines[ipTAl48]);

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