/*CoolArgo compute argon cooling */
#include "cddefines.h"
#include "coolheavy.h"
#include "tepowers.h"
#include "phycon.h"
#include "taulines.h"
#include "fivel.h"
#include "tcool.h"
#include "ligbar.h"
#include "coladd.h"
#include "level2.h"
#include "putcs.h"
#include "level3.h"
#include "pop3.h"
#include "coolmetals.h"
#include "popexc.h"
#include "ionfracs.h"

void CoolArgo()
{
	float a12, 
	  a13, 
	  a14, 
	  a15, 
	  a23, 
	  a24, 
	  a25, 
	  a34, 
	  a35, 
	  a45, 
	  cs, 
	  cs12, 
	  cs13, 
	  cs14, 
	  cs15, 
	  cs23, 
	  cs24, 
	  cs25, 
	  cs34, 
	  cs35, 
	  cs45, 
	  pop2, 
	  popn3;
	double cs2s2p, 
	  cs2s3p, 
	  p[5];
	static double gAr4[5]={4.,4.,6.,2.,4.};
	static double exAr4[4]={21090.4,128.9,13636.2,177.1};

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

	/* Argon II 6.98 micron
	 * trans prob from 
	 * >>refer Nussbaumer, H., & Storey, P.J. 1988, A&A, 200, L25
	 * Y(ik) from
	 * >>refer  Pelan, J., & Berrington, K.A. 1995, A&A Suppl, 110, 209 */
	PutCS(3.1,&TauLines[ipTAr7]);
	level2(&TauLines[ipTAr7]);

	/* A III 7136, 7751, 3109, 5192 CS
	 * >>refer Galavis, M.E., Mendoza, C., & Zeippen, C.J. 1995, A&AS, 111, 347
	 * c7136 = popexc(4.825,9.,5.,0.403,1.955e4,xIonFracs(18,3))*2.7e-12
	 * call coladd( 'Ar 3' , 7136 , C7136  )
	 * >>chng 97 mar 19, break up into three level atom */
	cs12 = 4.825f;
	cs13 = 0.841f;
	cs23 = (float)MIN2(1.30,3.296/(TePowers.te10*TePowers.te01/TePowers.te001/
	  TePowers.te001));
	a12 = 0.3963f;
	a23 = 2.59f;
	a13 = 3.952f;
	/* POP3(G1,G2,G3,O12,O13,O23,A21,A31,A32,E12,E23,P2,ABUND,GAM2) */
	popn3 = (float)(pop3(9.,5.,1.,cs12,cs13,cs23,a12,a13,a23,1.955e4,2.770e4,
	  &pop2,xIonFracs[3][17],0.));
	CoolHeavy.c7136 = pop2*a12*2.7e-12;
	CoolHeavy.c5192 = popn3*a23*3.83e-12;
	CoolHeavy.c3109 = popn3*a13*6.40e-12;
	coladd("Ar 3",7136,CoolHeavy.c7136);
	tcool.dCooldT += CoolHeavy.c7136*1.955e4*tcool.tsq1;
	coladd("Ar 3",5192,CoolHeavy.c5192);
	coladd("Ar 3",3109,CoolHeavy.c3109);

	/* Ar III 21.8(J=0,1), 9.0 (J=1,2) mircon lines,
	 * CS 
	 * >>refer Galavis, M.E., Mendoza, C., & Zeippen, C.J. 1995, A&AS, 111, 347 */
	PutCS(3.1,&TauLines[ipTAr9]);
	cs = (float)MIN2(1.384,3.110/(TePowers.te10/TePowers.te001*TePowers.te001));
	PutCS(cs,&TauLines[ipTAr22]);
	if( phycon.te < 1e4 )
	{
		cs = 0.671f;
	}
	else
	{
		cs = (float)MIN2(0.906,0.150*TePowers.te20/TePowers.te02/TePowers.te02*
		  TePowers.te003);
	}
	PutCS(cs,&TauDummy);
	level3(&TauLines[ipTAr9],&TauLines[ipTAr22],&TauDummy);

	/* Argon IV 4711+4740, 7335 lines (O II like)
	 * CS from
	 * >>refer Zeippen, C.J., Le Bourlot, J., Butler, K. 1987, A&A, 188, 251
	 * expression is 1-3 by itself, convert to whole mult by *2.5
	 * cs12 = MIN(3., 93.094/(sqrte/te10*te02) ) * 2.5
	 * POP3(G1,G2,G3,O12,O13,O23,A21,A31,A32,E12,E23,P2,ABUND,GAM2)
	 * C7335 = POP3( 4.,10.,6., cs12 ,0.476,10.2 , 0.0103,1.7,1.16,
	 *  1 3.13e4 , 1.98e4 , pop2 ,xIonFracs(18,4),0.) * 1.16*2.75E-12
	 * c4740 = pop2 * 0.0103 * 4.2e-12
	 *
	 * >>chng 97 jan 31, increase to full 5 level atom
	 * Ar IV, cs data from 
	 * >>refer Ramsbottom, C.A., Bell., K.L., & Keenan, F.P., 1997,
	 * >>refer MNRAS 284, 754
	 * differs by 2-3x from older values
	 * temp dependence form 
	 *>>refer Ramsbottom, C.A., & Bell, K.L. 1997, At. Data Nucl. Data Tables, 66, 65 */
	cs12 = (float)MAX2(0.761,0.481*TePowers.te05);
	cs12 = (float)MIN2(0.853,cs12);
	a12 = 2.23e-2f;

	cs13 = (float)MAX2(1.14,0.719*TePowers.te05);
	cs13 = (float)MIN2(1.3,cs13);
	a13 = 1.77e-3f;

	cs14 = (float)MAX2(0.39,0.108*TePowers.te10*TePowers.te02*TePowers.te02);
	cs14 = (float)MIN2(0.5,cs14);
	a14 = 0.862f;

	cs15 = (float)MAX2(0.78,0.216*TePowers.te10*TePowers.te02*TePowers.te02);
	cs15 = (float)MIN2(1.0,cs15);
	a15 = 2.11f;

	cs23 = 7.06f;
	a23 = 2.30e-5f;

	cs24 = (float)MAX2(1.53,0.346*TePowers.te10*TePowers.te05);
	cs24 = (float)MIN2(1.96,cs24);
	a24 = 0.603f;

	cs25 = (float)MAX2(2.18,0.664*TePowers.te10*TePowers.te02);
	cs25 = (float)MIN2(2.65,cs25);
	a25 = 0.789f;

	cs34 = (float)MAX2(1.56,0.475*TePowers.te10*TePowers.te02);
	cs34 = (float)MIN2(1.89,cs34);
	a34 = 0.119f;

	cs35 = (float)MAX2(4.01,1.00*TePowers.te10*TePowers.te02*TePowers.te02);
	cs35 = (float)MIN2(5.03,cs35);
	a35 = 0.598f;

	cs45 = (float)(0.0359*TePowers.te20*TePowers.te20*TePowers.te02*TePowers.te02);
	a45 = 4.94e-5f;

	/* FIVEL( G(1-5) , ex(wn,1-5), cs12,cs13,14,15,23,24,25,34,35,45,
	 *  A21,31,41,51,32,42,52,43,53,54, pop(1-5), abund) */
	fivel(gAr4,exAr4,cs12,cs13,cs14,cs15,cs23,cs24,cs25,cs34,cs35,
	  cs45,a12,a13,a14,a15,a23,a24,a25,a34,a35,a45,p,xIonFracs[4][17]);
	CoolHeavy.Ar4740 = p[1]*a12*4.20e-12;
	CoolHeavy.Ar4711 = p[2]*a13*4.20e-12;
	CoolHeavy.Ar2868 = p[3]*a14*6.94e-12;
	CoolHeavy.Ar2854 = p[4]*a15*6.94e-12;
	CoolHeavy.Ar7263 = p[3]*a24*2.74e-12;
	CoolHeavy.Ar7171 = p[4]*a25*2.74e-12;
	CoolHeavy.Ar7331 = p[3]*a34*2.74e-12;
	CoolHeavy.Ar7237 = p[4]*a35*2.74e-12;
	tcool.dCooldT +=  (CoolHeavy.Ar4740 + CoolHeavy.Ar4711)*(3.04e4*
	  tcool.tsq1 - tcool.halfte) ;
	tcool.dCooldT += (CoolHeavy.Ar2868 + CoolHeavy.Ar2854 + CoolHeavy.Ar7263 + 
	  CoolHeavy.Ar7171 + CoolHeavy.Ar7331 + CoolHeavy.Ar7237)*(5.02e4*
	  tcool.tsq1 - tcool.halfte);
	coladd("Ar 4",7335,CoolHeavy.Ar7263+CoolHeavy.Ar7171+CoolHeavy.Ar7331+
	  CoolHeavy.Ar7237);
	coladd("Ar 4",4740,CoolHeavy.Ar4740+CoolHeavy.Ar4711);
	coladd("Ar 4",2868,CoolHeavy.Ar2868+CoolHeavy.Ar2854);

	/* Argon V 6435+7007, A from 
	 * >>refer Mendoza, C., & Zeippen, C.J. 1982, MNRAS, 199, 1025
	 * CS from 
	 * >>refer Galavis, M.E., Mendoza, C., & Zeippen, C.J. 1995, A&AS, 111, 347
	 * POPEXC( O12,g1,g2,A21,excit,abund) ; result already*a21, excit in Kelvin */
	if( phycon.te < 1e4 )
	{
		cs = 3.09f;
	}
	else
	{
		cs = (float)MIN2(4.454,0.634*TePowers.te20/TePowers.te03*TePowers.te001*
		  TePowers.te001);
	}
	CoolHeavy.c7007 = popexc(cs,9.,5.,0.68,2.055e4,xIonFracs[5][17])*
	  2.84e-12;
	coladd("Ar 5",7007,CoolHeavy.c7007);

	/* Ar V 3P fine structure lines , A from
	 * >>refer Mendoza, C. 1982, in Planetary Nebulae, IAU Symp No. 103,
	 * >>refer ed by D.R. Flower, (D. Reidel: Holland), 143
	 * cs=
	 * >>refer Galavis, M.E., Mendoza, C., & Zeippen, C.J. 1995, A&AS, 111, 347
	 * >>chng 96 july 16, bug upward change in collision strengths from 
	 * >>refer Mendoza, C. 1982, in Planetary Nebulae, IAU Symp No. 103,
	 * >>refer ed by D.R. Flower, (D. Reidel: Holland), 143 */
	cs = (float)MIN2(3.26,26.27/(TePowers.te20*TePowers.te03*TePowers.te01*
	  TePowers.te003));
	PutCS(cs,&TauLines[ipTAr13]);
	cs = (float)MIN2(8.47,44.31/(TePowers.te20/TePowers.te003/TePowers.te003));
	PutCS(cs,&TauLines[ipTAr8]);
	cs = (float)MIN2(1.95,7.280/(TePowers.te10*TePowers.te05*TePowers.te005));
	PutCS(cs,&TauDummy);
	level3(&TauLines[ipTAr13],&TauLines[ipTAr8],&TauDummy);

	/* [Ar VI] 4.53 micron, cs from 
	 * >>refer Saraph, H.E., & Storey, P.J. A&AS, 115, 151
	 * call PutCS( .8 , Ar06453 )
	 * >>chng 96 Dec 11, cs should have been 6.33, caught by Simon Casassus */
	PutCS(6.33,&TauLines[ipAr06453]);
	level2(&TauLines[ipAr06453]);
	/* fs4p5 = popexc( 0.8, 2.,4., 9.7e-2, 3174.3,xIonFracs(18,6))
	 *  1 * 4.39E-13
	 * call coladd( 'Ar 6' , 4 , FS4P5  )
	 *
	 * [ArX] 5533.4A, cs 
	 * >>refer Saraph, H.E. & Tully, J.A. 1994, A&AS, 107, 29 */
	cs = (float)MIN2(0.573,20.05/(TePowers.te30*TePowers.te03*TePowers.te005));
	PutCS(cs,&TauLines[ipAr1055]);
	level2(&TauLines[ipAr1055]);

	/* [Ar XI] 2.60, microns, 6917A */
	cs = (float)MIN2(0.207,2.685e-3*TePowers.te20*TePowers.te20*
	  TePowers.te001*TePowers.te001);
	cs = (float)MAX2(0.09,cs);
	PutCS(cs,&TauLines[ipAr1126]);

	cs = (float)MIN2(0.64,0.0127*TePowers.te30*TePowers.te05*TePowers.te001*
	  TePowers.te001);
	cs = (float)MAX2(0.25,cs);
	PutCS(cs,&TauLines[ipAr1178]);

	cs = (float)MIN2(0.17,2.789e-3*TePowers.te30*TePowers.te05*
	  TePowers.te02/TePowers.te003);
	cs = (float)MAX2(0.06,cs);
	PutCS(cs,&TauDummy);

	level3(&TauLines[ipAr1126],&TauLines[ipAr1178],&TauDummy);

	/* Ar X IV 4413, wavelength+a from 
	 * >>refer Froese Fischer, C. 1983, J.Phys. B, 16, 157
	 * cs is guess */
	CoolHeavy.fs4413 = popexc(0.2,2.,4.,105.,3.259e4,xIonFracs[14][17])*
	  4.51e-12;
	coladd("Ar14",4413,CoolHeavy.fs4413);

	/* Ar 15 409A, Be seq interpolated lambda, cs, A. */
	CoolHeavy.c409 = popexc(0.1,1.,9.,1.1e6,3.52e5,xIonFracs[15][17])*
	  4.86e-11;
	coladd("Ar15",409,CoolHeavy.c409);

	/* Ar 16 365, 25, Li Seq
	 * >>refer Cochrane, D.M., & McWhirter, R.W.P. 1983, PhyS, 28, 25 */
	ligbar(18,&TauLines[ipT354],&TauLines[ipT25],&cs2s2p,&cs2s3p);

	PutCS(cs2s2p,&TauLines[ipT354]);
	level2(&TauLines[ipT354]);

	/* funny factor (should have been 0.5) due to energy change */
	PutCS(cs2s2p*0.454,&TauLines[ipT389]);
	level2(&TauLines[ipT389]);

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


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