/*CoolTita compute titanium cooling */
#include "cddefines.h"
#include "taulines.h"
#include "ionfracs.h"
#include "coolheavy.h"
#include "coladd.h"
#include "level2.h"
#include "putcs.h"
#include "pop3.h"
#include "coolmetals.h"

void CoolTita()
{
	float a21, 
	  a31, 
	  a32, 
	  p2, 
	  p3;

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

	/* Ti Titanium cooling
	 *
	 * these are 3 lines estimated by Jim Kingdon
	 * a's are bad, collision strengths just one */
	a21 = 0.015f;
	a31 = 0.032f;
	a32 = 0.002f;

	/* POP3(G1,G2,G3,O12,O13,O23,A21,A31,A32,E12,E23,P2,ABUND,GAM2)
	 * energies are in kelvin */
	p3 = (float)pop3(21.,5.,9.,21.,5.0,9.0,a21,a31,a32,11844.,3148.4,&p2,
	  xIonFracs[3][21],0.);

	/* multiplet at roughly 9594 */
	CoolHeavy.Ti3l31 = p3*a31*2.07e-12;
	/* multiplet at roughly 4.57 microns */

	CoolHeavy.Ti3l32 = p3*a32*4.35e-13;
	/* multiplet at roughly 1.21 microns */

	CoolHeavy.Ti3l21 = p2*a21*1.64e-12;
	coladd("Ti 3",9594,CoolHeavy.Ti3l31);
	coladd("Ti 3",4,CoolHeavy.Ti3l32);
	coladd("Ti 3",1,CoolHeavy.Ti3l21);

	/* [Ti VI] 1.7150 mic
	 * Y(ik) from 
	 * >>refer Pelan, J., & Berrington, K.A. 1995, A&A Suppl, 110, 209 */
	PutCS(3.48,&TauLines[ipTi06172]);
	level2(&TauLines[ipTi06172]);

	/* [Ti XIV] 2117.79, cs from 
	 * >>refer Saraph, H.E. & Tully, J.A. 1994, A&AS, 107, 29 */
	PutCS(0.23,&TauLines[ipTi14212]);
	level2(&TauLines[ipTi14212]);

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