pdefs.h

美国berkeley大学开发的有界等离子体粒子1d3v计算机模拟程序,很实用

/* 
**	PDEFS.H
**	
**	Header file for all PDS1 modules.
*/ 

#include <stdlib.h>
#include <math.h>

#define TWOPI			6.28318530717959
#define FOURPI			12.5663706
#define PI				3.141592653589793
#define EPS0			8.854e-12			/* (F/m)	*/
#define NperTORR		8.3221e20		
#define MAXLEN 		5000
#define NSMAX 			4
#define NCMAX 			400
#define NEMAX 			500
#define HISTMAX 		512  					/* upper bound on histories */
#define DBL_MIN 		1E-200
#define frand()    	(rand()/32768.0)

/*******************************************************************/
#if defined(GLOBALORIGIN)
	#define GLOBAL
#else
	#define GLOBAL extern
#endif

GLOBAL float nc2p, r0, r1, area, rhoback, backj, dde,	epsilon,extr,
         extl, extc, w0, dcbias, acbias, extq, extq_1, extq_2,	extq_3,
 			exti, sigma, oldsigma, jp, dr, dt, vrscale, vscale, dttx,
			pressure, gtemp, gden, selsmax, elsengy0, elsengy1, elsengy2,
			sextmax, extengy0, extengy1, extengy2, sionmax, ionengy0,
			ionengy1, ionengy2, achrgx, bchrgx, ascat, bscat, seec[NSMAX],
			df, ramp, theta0, b0, jwall[NSMAX],	qm[NSMAX], q[NSMAX], m[NSMAX],
			jj0[NSMAX][2], v0[NSMAX][2],	vt[NSMAX][2], vc[NSMAX][2],
			vtpt[NSMAX], emin[NSMAX], de[NSMAX], dtheta[NSMAX],
			Escale[NSMAX],	risetime, enter[NSMAX][2];

GLOBAL int nsp, nc, ng, ntri, secondary, ecolsp, icolsp, secsp, ionsp, 
       ecollisional, icollisional, hist_hi,	thist_hi, freq_hi, interval,
		 np[NSMAX], index[1024], nfft, src, dcramped, nbin[NSMAX], reflux,
		 inject[NSMAX];

GLOBAL char freqbuf[30];

GLOBAL double t, Trf;

GLOBAL float far **r, far **vr, far **vth, far **vph, far **srho,	far *rho, 
         far *e, far *phi, far *r_array,	far *t_array, far *f_array,
			far *sels, far *sext, far *sion, far **np_hist,	far **jwall_hist,	
			far **phi_hist, far *sigma_hist, far *ese_hist,	far *ke_hist,
			far *te_hist, far **kes_hist, far *cur_hist,	far *com_cur_hist,
			far **com_phi_hist, far *cur_fft, far *phi_fft,	far *z_fft,
			far *mphi_fft, far *Local_t_array, far *ionrate, far *iontemp,
			far *chrgxrate, far *rg, far *rg3, far *avene, far *avenemp,
			far **fe, far **ftheta, far **e_array, far **th_array,
			far *Power_hist;

GLOBAL void circuit1(),	circuit2(),	circuit3(), circuit4(), (*circuitptr)();
GLOBAL float bc1(), bc2(), bc3(), bc4(), (*bcptr)();   

#undef GLOBAL