#include <malloc.h>
#include "pdp1.h"

// Function void mccdiag_init(void) {{{
void mccdiag_init(void)
{
	int i, j;

	if(ecollisional==0 && icollisional==0) {
		ndiag = 0;
		return;
	}

	/* default gas = argon */

	if ( (gas < HELIUM) || (gas > OXYGEN) ) gas = ARGON;

	switch (gas) {

		case HELIUM: // HELIUM {{{
			ndiag = 4;

			rate_title = (char **)malloc(ndiag*sizeof(char *));
			rate_title[0]= "MCC Elastic";
			rate_title[1]= "MCC Excitation";
			rate_title[2]= "MCC Ionization";
			rate_title[3]= "MCC Charge X";
			break;
			// }}}

		case ARGON: // case ARGON {{{
			ndiag = 5;

			rate_title = (char **)malloc(ndiag*sizeof(char *));
			rate_title[0]= "MCC Elastic";
			rate_title[1]= "MCC Excitation";
			rate_title[2]= "MCC Ionization";
			rate_title[3]= "MCC Scattering";
			rate_title[4]= "MCC Charge X";
			break;
			// }}}

		case NEON: // case NEON {{{
			ndiag = 7;

			rate_title = (char **)malloc(ndiag*sizeof(char *));
			rate_title[0]= "MCC Elastic";
			rate_title[1]= "MCC Excitation1";
			rate_title[2]= "MCC Excitation2";
			rate_title[3]= "MCC Excitation3";
			rate_title[4]= "MCC Excitation4";
			rate_title[5]= "MCC Ionization";
			rate_title[6]= "MCC Charge X";
			break;
			// }}}

		case OXYGEN: // case OXYGEN {{{
			ndiag = 21;

			rate_title = (char **)malloc(ndiag*sizeof(char *));
			rate_title[0]=  "MCC Elastic";
			rate_title[1]=  "MCC Rotational";
			rate_title[2]=  "MCC Vibrational1";
			rate_title[3]=  "MCC Vibrational2";
			rate_title[4]=  "MCC Vibrational3";
			rate_title[5]=  "MCC Vibrational4";
			rate_title[6]=  "MCC O2 Sing Delta";
			rate_title[7]=  "MCC O2 Sing Sigma";
			rate_title[8]=  "MCC Diss Attach";
			rate_title[9]=  "MCC E Loss1";
			rate_title[10]= "MCC E Loss2";
			rate_title[11]= "MCC E Loss3";
			rate_title[12]= "MCC E Loss4";
			rate_title[13]= "MCC Ionization";
			rate_title[14]= "MCC E Loss5";
			rate_title[15]= "MCC Diss Recombine";
			rate_title[16]= "MCC E Impact Detach";
			rate_title[17]= "MCC O- Elastic";
			rate_title[18]= "MCC O- Detachment";
			rate_title[19]= "MCC Mutual Neutral";
			rate_title[20]= "MCC Charge X";
			break;
			// }}}
	}

	/*  Allocate arrays for diagnostic rates */
	mccrate  = (float **)malloc(ndiag*sizeof(float *));
	rate_show= (float **)malloc(ndiag*sizeof(float *));
	for(i=0; i<ndiag; i++) {
		mccrate[i]  = (float *)malloc(ng*sizeof(float));
		rate_show[i]= (float *)malloc(ng*sizeof(float));
		for (j=0; j<ng; j++) mccrate[i][j] = rate_show[i][j] = 0.0;
	}
} // }}}

// Modeline for ViM {{{1
// vim:set ts=4
// vim600:fdm=marker fdl=0 fdc=3
// 1}}}