/*CoolSili compute silicon cooling */
#include "cddefines.h"
#include "coolheavy.h"
#include "taulines.h"
#include "tepowers.h"
#include "phycon.h"
#include "poplevls.h"
#include "h.h"
#include "edsqte.h"
#include "ionfracs.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 "popexc.h"
#include "coolmetals.h"
#include "expion.h"
#include "embesq.h"

void CoolSili()
{
	double cs, 
	  cs2s2p, 
	  cs2s3p, 
	  cs01, 
	  cs02, 
	  cs12;
	float
	  p2;
	

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

	/* Si I 2518 */
	MakeCS(&TauLines[ipSii2518]);
	level2(&TauLines[ipSii2518]);

	/* Si I 2215 */
	MakeCS(&TauLines[ipSii2215]);
	level2(&TauLines[ipSii2215]);

	/* Silicon II 35 micron; elec cs from 
	 *>>refer Dufton, P.L., & Kingston, A.E. 1994, At. Data Nucl. Data Tables,
	 *>>refer 57, 273 */
	cs = 5.77;
	/* hydrogen collision strength from 
	 * >>refer Tielens, A.G.G., & Hollenbach, D. 1985, ApJ, 291, 722
	 * they give rate de-ex 6.5E-10 cm^3 s^-1, indep of temp */
	cs += 6.5e-10/edsqte.cdsqte*4.*h.hi;
	PutCS(cs,&TauLines[ipTSi35]);
	level2(&TauLines[ipTSi35]);

	/* Si II  2337, atomic data from 
	 * >>refer Dumont, A.M., & Mathez, G. 1981, A&A, 102, 1
	 * coll str 
	 * >>refer Dufton, P.L., & Kingston, A.E. 1991, MNRAS, 248, 827
	 * also 
	 * >>refer Mendoza, C. 1982, in Planetary Nebulae, IAU Symp No. 103,
	 * >>refer ed by D.R. Flower, (D. Reidel: Holland), 143 */
	cs = MIN2(6.0,15.836/(TePowers.te10*TePowers.te02/TePowers.te003));
	PutCS(5.5,&TauLines[ipT2335]);
	level2(&TauLines[ipT2335]);

	/* Si II 1808, permitted resonance line,
	 * osc str from morton et al 88 (apj sup); 
	 * all si ii collision data (following 4 lines) are from
	 * >>refer Dufton, P.L., & Kingston, A.E. 1991, MNRAS, 248, 827
	 * following assumes there is typo in table 1 of dufton and kingston
	 * and that they meant 2s 2p^2 ^2D instead of 3d */
	PutCS(13.01,&TauLines[ipT1808]);
	level2(&TauLines[ipT1808]);

	PutCS(3.61,&TauLines[ipT1527]);
	level2(&TauLines[ipT1527]);

	PutCS(2.89,&TauLines[ipT1305]);
	level2(&TauLines[ipT1305]);

	PutCS(12.25,&TauLines[ipT1260]);
	level2(&TauLines[ipT1260]);

	/* permitted Si III 1206.5, collision strength from 
	 * >>refer Callaway, J. 1994, At. Data Nucl. Data Tables, 57, 9 */
	cs = MIN2(7.0,1.442*TePowers.te10*TePowers.te03*TePowers.te03/
	  TePowers.te01);
	PutCS(cs,&TauLines[ipT1207]);
	level2(&TauLines[ipT1207]);

	/* Si III] 1895, CS=
	 * >>refer Dufton, P.L., & Kingston, A.E. 1989, MNRAS, 241, 209
	 * >>refer Callaway, J. 1994, At. Data Nucl. Data Tables, 57, 9
	 * a=
	 * >>refer Nussbaumer, H. 1986, A&A , 155, 205
	 * grnd 3s^2 ^1S, upper lev 3p ^3P^o j=0,1,2 */
	cs = 106./(TePowers.te10*TePowers.te10*TePowers.te10*TePowers.te02);
	PutCS(cs,&TauLines[ipT1895]);
	beseq(1.8,3.6,10.4,&TauLines[ipT1895],.013);
	embesq.em1895 = (float)(PopLevls.PopLevels[3]*0.013*1.05e-11);

	/* Si IV 1394, 1403, data from 
	 * >>refer Mendoza, C. 1982, in Planetary Nebulae, IAU Symp No. 103,
	 * >>refer ed by D.R. Flower, (D. Reidel: Holland), 143
	 * cs from 
	 * >>refer Dufton, P.L., & Kingston, A.E. 1987, J.Phys. B, 20, 3899 */
	cs = 6.37*TePowers.te10;
	PutCS(cs*0.667,&TauLines[ipT1394]);
	PutCS(cs*0.333,&TauLines[ipT1403]);
	PutCS(1.0,&TauDummy);
	level3(&TauLines[ipT1403],&TauDummy,&TauLines[ipT1394]);

	/* Si VI 1.96 micron
	 * >>refer Saraph, H.E. & Tully, J.A. 1994, A&AS, 107, 29
	 * >>chng 96 july 16 had been constant */
	cs = MIN2(0.43,0.0448*TePowers.te20/TePowers.te003/TePowers.te003);
	cs = MAX2(0.3,cs);
	PutCS(cs,&TauLines[ipSi619]);
	level2(&TauLines[ipSi619]);

	/* Si VII 2148- OIII like, 
	 * >>refer Kafatos, M., & Lynch, J.P. 1980, ApJS, 42, 611 */
	CoolHeavy.c2148 = 
		popexc(0.4,9.,5.,15.,6.7e4,xIonFracs[7][13])*9.26e-12;
	coladd("Si 7",2148,CoolHeavy.c2148);

	/* Si VII ground term, 2.48, 6.51 microns
	 * cs 
	 * >>refer Butler, K., & Zeippen, C.J. 1994, A&AS, 108, 1 */
	cs = MIN2(0.217,0.0904*TePowers.te05*TePowers.te03/TePowers.te003/
	  TePowers.te001);
	PutCS(cs,&TauLines[ipTSi65]);
	cs = MIN2(0.70,8.79e-2*TePowers.te10*TePowers.te10/TePowers.te02);
	PutCS(cs,&TauLines[ipTSi25]);
	cs = MIN2(0.20,9.751e-3*TePowers.te20*TePowers.te03*TePowers.te03/
	  TePowers.te003);
	PutCS(cs,&TauDummy);
	level3(&TauLines[ipTSi25],&TauLines[ipTSi65],&TauDummy);

	/* Si 8 1446, 3727-like, 
	 * >>refer Kafatos, M., & Lynch, J.P. 1980, ApJS, 42, 611 */
	CoolHeavy.c1446 = 
		popexc(0.4,4.,10.,1.,9.97e4,xIonFracs[8][13])*
	  1.39e-11;
	coladd("Si 8",1446,CoolHeavy.c1446);

	/* Si 9 1985, 2150
	 * cs, As from
	 * >>refer Aggarwal, K.M. 1983, J.Phys. B, 16, L59
	 * >>refer Baluja, K.L. 1985, J.Phys. B, 18, L413 */
	CoolHeavy.c949 = 
		pop3(9.,5.,1.,0.5913,0.0757,0.225,26.3,214.,5.16,
	  7.62e4,7.902e4,&p2,xIonFracs[9][13],0.)*214.*2.096e-11;
	CoolHeavy.c1815 = CoolHeavy.c949*1.912*0.0516;
	CoolHeavy.c1985 = p2*26.3*1.0e-11;
	coladd("Si 9",949,CoolHeavy.c949);
	coladd("Si 9",1815,CoolHeavy.c1815);
	coladd("Si 9",1985,CoolHeavy.c1985);

	/* Si 9 3P fine structure lines, A=
	 * >>refer Baluja, K.L. 1985, J.Phys. B, 18, L413
	 * 2.583, 3.9microns
	 * CS=
	 * >>refer Lennon, D.J. Burke, V.M. 1994, A&AS, 103, 273 */
	cs01 = MIN2(0.98,28.51/(TePowers.te10*TePowers.te10*TePowers.te10*
	  TePowers.te10/TePowers.te03*TePowers.te003*TePowers.te001*TePowers.te001));

	cs12 = MIN2(2.7,81.21/(TePowers.te10*TePowers.te10*TePowers.te10*
	  TePowers.te10/TePowers.te01/TePowers.te01/TePowers.te001/TePowers.te001));

	cs02 = MIN2(0.70,19.67/(TePowers.te10*TePowers.te10*TePowers.te10*
	  TePowers.te10/TePowers.te03*TePowers.te001));

	PutCS(cs01,&TauLines[ipTSi4]);
	PutCS(cs12,&TauLines[ipTSi3]);
	PutCS(cs02,&TauDummy);

	level3(&TauLines[ipTSi4],&TauLines[ipTSi3],&TauDummy);

	/* 5S0 - 3P, cs from, A=guess
	 * >>refer Aggarwal, K.M. 1984, ApJS, 54, 1 */
	CoolHeavy.c691 = popexc(40.6/TePowers.sqrte*TePowers.te10,9.,5.,
	  1e4,2.081e5,xIonFracs[9][13])*2.88e-11;
	coladd("Si 9",691,CoolHeavy.c691);

	/* Si 10 606, actually four lines clumped together.
	 * atomic data 
	 * >>refer Saha, H.P., & Trefftz, E. 1982, A&A, 116, 224 */
	CoolHeavy.c606 = 
		0.10*1.42e-16*expion(2.4e5,xIonFracs[10][13]);
	coladd("Si10",606,CoolHeavy.c606);

	/* Si X 1.43m, A from 
	 * >>refer Chandra, S. 1982, SoPh, 75, 133
	 * cs from 
	 * >>refer Zhang, H.L., Graziani, M., Pradhan, A.K. 1994, A&A, 283, 319 */
	if( phycon.te <= 40500. )
	{
		cs = 0.190*TePowers.te20/TePowers.te001;
	}
	else
	{
		cs = 24.93/(TePowers.te20*TePowers.te03*TePowers.te01*TePowers.te003*
		  TePowers.te003);
	}
	PutCS(cs,&TauLines[ipSi10143]);
	level2(&TauLines[ipSi10143]);

	/* SI 11 582.9, 1909-LIKE, CS=
	 * >>refer Berrington, K.A., Burke, P.G., Dufton, P.L., Kingston, A.E.
	 * >>refer 1985, At. Data Nucl. Data Tables, 33, 195
	 * A=
	 * >>refer Muhlethaler, H.P., & Nussbaumer, H. 1976, A&A 48, 109 */
	CoolHeavy.c583 = 
		popexc(0.10,1.,9.,1e5,2.47e5,xIonFracs[11][13])*
	  3.4e-11;
	coladd("Si11",583,CoolHeavy.c583);

	/* li seq 2s2p and 2s3p, Si XII 499, 521
	 * >>refer Cochrane, D.M., & McWhirter, R.W.P. 1983, PhyS, 28, 25 */
	ligbar(14,&TauLines[ipTSi499],&TauLines[ipTSi41],&cs2s2p,&cs2s3p);
	PutCS(cs2s2p,&TauLines[ipTSi499]);
	PutCS(cs2s2p*0.5,&TauLines[ipTSi521]);
	PutCS(1.0,&TauDummy);
	level3(&TauLines[ipTSi521],&TauDummy,&TauLines[ipTSi499]);

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


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