mccdiaginit.c

来自「一维等离子体静电粒子模拟程序」· C语言 代码 · 共 64 行

#include "pdc1.h"

void mccdiag_init(void)
{
  int i, j;
  
  if(ecollisional==0 && icollisional==0) 
  {
    ncolls = 0;
    return;
  }
  
/* default gas = argon */

  if ( (gas < ARGON_LK) || (gas > ARGON) ) gas = ARGON;

   switch (gas) 
   {
	case ARGON:
	   ncolls= 7;
	   rate_title = (char **)malloc(ncolls*sizeof(char *));
	   rate_title[0]= "MCC Elastic";
	   rate_title[1]= "MCC Excitation";
	   rate_title[2]= "MCC Ionization";
	   rate_title[3]= "MCC Charge X";
	   rate_title[4]= "MCC Scattering";
	   rate_title[5]= "MCC Ionization Ex";
	   rate_title[6]= "MCC Metastable";
	   break;

	case ARGON_LK:
	   ncolls= 7;
	   rate_title = (char **)malloc(ncolls*sizeof(char *));
	   rate_title[0]= "MCC Elastic";
	   rate_title[1]= "MCC Excitation";
	   rate_title[2]= "MCC Ionization";
	   rate_title[3]= "MCC Charge X";
	   rate_title[4]= "MCC Scattering";
	   rate_title[5]= "MCC Ionization Ex";
	   rate_title[6]= "MCC Metastable";
	   break;

  }
  
  /*  Allocate arrays for diagnostic rates */
  Pel       = (float *)malloc(ng*sizeof(float));
  Pel_show  = (float *)malloc(ng*sizeof(float));
  Pcx       = (float *)malloc(ng*sizeof(float));
  Pcx_show  = (float *)malloc(ng*sizeof(float));
  Pie       = (float *)malloc(ng*sizeof(float));
  Pie_show  = (float *)malloc(ng*sizeof(float));
  for (j=0; j<ng; j++) Pel[j] = Pel_show[j] = Pcx[j] = Pcx_show[j] = 
		       Pie[j] = Pie_show[j]=0.0;

  mccrate  = (float **)malloc(ncolls*sizeof(float *));
  rate_show= (float **)malloc(ncolls*sizeof(float *));
  for(i=0; i<ncolls; 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;
  }
}