findif.h

matlab编写的三维FDTD程序 pml边界条件 实现电磁波传播数值模拟和仿真

/* findif.h last modified 7/26/94 */
/* 
   This is an include file needed to compile fdmenu.c and fdtd3d.c .
*/
struct fddat {
	char filout[60];	/* ROOT NAME OF OUTPUT FILES */
	float centfreq;		/* NOT USED */
	float amp;		/* NOT USED */
	float perm;		/* NOT USED */
	float cond;		/* NOT USED */
	float dist;		/* NOT USED */
	float size;		/* NOT USED */
	float dx;		/* GRID CELL DIMENSION IN X-DIRECTION (m) */
	float dy;		/* GRID CELL DIMENSION IN Y-DIRECTION (m) */
	float dz;		/* GRID CELL DIMENSION IN Z-DIRECTION (m) */
	float dt;		/* TIME STEP INCREMENT (ns) */
	float modsiz;		/* PEAK VOLTAGE OF INPUT GAUSSIAN IMPULSE (V) */
	char pol;		/* NOT USED */
	short int timnum;	/* NUMBER OF TIME STEPS TO CALCULATE */
	float xmax;		/* X-DIMENSION OF GRID */
	float ymax;		/* Y-DIMENSION OF GRID */
	float zmax;		/* Z-DIMENSION OF GRID */
	short int bounswi;	/* TYPE OF ABC IMPLEMENTED */
        float tlength;		/* TIME LENGTH OF GAUSSIAN IMPULSE */
        float tcval;		/* DELAY IN ns BEFORE PEAK VOLTAGE */
        short int moddim;	/* NOT USED */
        short int indat;	/* NOT ACTIVATED */
        char infile[60];	/* NOT ACTIVATED */
	short int obj;		/* TYPE OF TARGET ACTIVATED IN GRID */	
	float rad;		/* RADIUS OF TARGET (IF CYLINDER, PIPE, OR SPHERE) (m) */
	float xpos;		/* X-POSITION OF TARGET IN GRID (m) */
	float ypos;		/* Y-POSITION OF TARGET IN GRID (m) */
	float zpos;		/* DEPTH OF TARGET BELOW INTERFACE IN GRID (m) */
	float xcen;		/* X-POSITION OF CENTER OF TRANSMIT ANTENNA IN GRID (m) */
	float ycen;		/* Y-POSITION OF CENTER OF TRANSMIT ANTENNA IN GRID (m) */
	float xlen;		/* LENGTH OF TRANSMIT ANTENNA (m) */
        float xang;		/* FLARE ANGLE OF TRANSMIT ANTENNA (degrees) */
	float xrcen;		/* X-POSITION OF CENTER OF RECEIVE ANTENNA IN GRID (m) */
	float yrcen;		/* Y-POSITION OF CENTER OF RECEIVE ANTENNA IN GRID (m) */
	float xrlen;		/* LENGTH OF RECEIVE ANTENNA (m) */
        float xrang;		/* FLARE ANGLE OF RECEIVE ANTENNA (degrees) */
        float anthei;		/* HEIGHT OF TRANSMIT AND RECEIVE ANTENNAS ABOVE INTERFACE (m)*/
        float zinter;		/* DISTANCE OF INTERFACE FROM TOP OF GRID (m) */
        float zcond;		/* CONDUCTIVITY OF MEDIUM BELOW INTERFACE (S/m) */
        float zperm;		/* RELATIVE PERMITTIVITY OF MEDIUM BELOW INTERFACE (DIMENSIONLESS) */
        short int antpat;	/* NOT USED */
	short int fullEx;	/* FLAG TO OUTPUT ALL EX VALUES IN GRID TO FILE AFTER fd1.timnum TIME STEPS */
	short int fullEy;	/* FLAG TO OUTPUT ALL Ey VALUES IN GRID TO FILE AFTER fd1.timnum TIME STEPS */
	short int fullEz;	/* FLAG TO OUTPUT ALL Ez VALUES IN GRID TO FILE AFTER fd1.timnum TIME STEPS */
	short int fullHx;	/* NOT ACTIVATED */
	short int fullHy;	/* NOT ACTIVATED */
	short int fullHz;	/* NOT ACTIVATED */
	short int HzEz;		/* NOT USED */
	short int ExEy;		/* NOT USED */
	short int HxHy;		/* NOT USED */
	float xlin;		/* NOT USED */		
	float ylin;		/* NOT USED */
	short int timpro;	/* TIME STEP INCREMENT BETWEEN WRITING TO OUTPUT FILES (IN TIMESTEPS) */	
	float timx;		/* X-POSITION OF VERTICAL PROFILE PARALLEL TO THE Y-AXIS (m) */
	float timy;		/* Y-POSITION OF VERTICAL PROFILE PARALLEL TO THE X-AXIS (m) */
        short int dump;		/* NOT ACTIVATED */
        char dumpfile[60];	/* NOT ACTIVATED */
        short int dumptime;	/* NOT ACTIVATED */
   	char filauto[60];	/* NOT ACTIVATED */
   	float zinc;		/* MOVEMENT OF TARGET IN GRID IN Z-DIRECTION BETWEEN TRACES IF fd1.numsteps>1 */
   	float xinc;		/* MOVEMENT OF TARGET IN GRID IN X-DIRECTION BETWEEN TRACES IF fd1.numsteps>1 */
   	float yinc;		/* MOVEMENT OF TARGET IN GRID IN Y-DIRECTION BETWEEN TRACES IF fd1.numsteps>1 */
   	short int numsteps;     /* NUMBER OF TRACES TO COLLECT IN MODEL RUN (OUTTERMOST LOOP) */
   	short int anttype;	/* TYPE OF ANTENNA TO USE (TRIANGULAR OR OSU) */
   	short int sheet;	/* THIN SHEET OPTION FLAG */
   	float objlen;		/* LENGTH OF TARGET (IF TARGET IS CUBE OR CYLINDER) */
   	float shtydim;		/* LENGTH OF THIN SHEET IN Y-DIRECTION (m) */
   	float shtxdim;		/* LENGTH OF THIN SHEET IN X-DIRECTION (m) */	
   	float shtthk;		/* THICKNESS OF THIN SHEET (m) */
   	float shtcond;		/* CONDUCTIVITY OF THIN SHEET (S/m) */
   	float shtperm;		/* RELATIVE PERMITTIVITY OF THIN SHEET (dimensionless) */
   	float shthei;		/* HEIGHT OF THIN SHEET ABOVE INTERFACE (m) */
   	short int patact;	/* FIELD PATTERN OUTPUT ACTIVATION FLAG */
   	float patrad;		/* RADIUS OF FIELD PATTERN */
   	float patphi;		/* AZIMUTH OF FIELD PATTERN */
   	float cylphi;		/* AZIMUTH OF CYLINDRICAL TARGET */
   	float cyldip; 		/* DIP OF CYLINDRICAL TARGET */
   	short int perconflg;	/* FLAG, IF = 1 THEN TARGET IS A PERFECT CONDUCTOR */
   	float targperm;		/* RELATIVE PERMITTIVITY OF TARGET (IF fd1.perconflg=0) (dimensionless) */	
   	float targcond;		/* CONDUCTIVITY OF TARGET (IF fd1.perconflg=0) (S/m) */
   	float zur;		/* RELATIVE MAGNETIC PERMEABILITY OF MEDIUM BELOW INTERFACE (dimensionless) */
   	short int tmproflg;	/* TIME PROFILE OPTION USED TO OUTPUT FILES (0,1,2,OR 3) */
   	short int absfreq;	/* FREQ. OF ELECTRICAL PROPERTIES OF OSU ANTENNA (KEEP SET AT 125 FOR OSU ANTENNA OPTION) */
   	short int pdflag;	/* ACTIVATION FLAG FOR INFINITESIMAL DIPOLE */	
   	short int attcoax;	/* NOT USED */
   	float coaxurt;		/* RELATIVE PERMEABILITY OF INTERIOR OF COAXIAL CABLE ATTACHED TO TRANSMIT ANTENNA (dimensionless)*/
   	float coaxst;		/* CONDUCTIVITY OF INTERIOR OF COAXIAL CABLE ATTACHED TO TRANSMIT ANTENNA (S/m) */
   	float coaxert;		/* RELATIVE PERMITTIVITY OF INTERIOR OF COAXIAL CABLE ATTACHED TO TRANSMIT ANTENNA (dimensionless)*/
   	float coaxurr;		/* RELATIVE PERMEABILITY OF INTERIOR OF COAXIAL CABLE ATTACHED TO RECEIVE ANTENNA (dimensionless)*/
   	float coaxsr;		/* CONDUCTIVITY OF INTERIOR OF COAXIAL CABLE ATTACHED TO RECEIVE ANTENNA (S/m) */
   	float coaxerr;		/* RELATIVE PERMITTIVITY OF INTERIOR OF COAXIAL CABLE ATTACHED TO RECEIVE ANTENNA (dimensionless)*/
   	float prdepth;		/* DEPTH OF PROBE OR SPOT E-FIELD MEASUREMENT POINT BELOW INTERFACE (m) */
   	float pvala;		/* DEPTH OF HORIZONTAL PLANE WHERE E-FIELD VALUES ARE WRITTEN TO OUTPUT FILE EVERY fd1.timpro
					TIMESTEPS (m) */ 
   	float pvalb;		/* DEPTH OF SECOND HORIZONTAL PLANE WHERE E-FIELD VALUES ARE WRITTEN TO OUTPUT FILE EVERY fd1.timpro
					TIMESTEPS (m) */
   	float pvalc;		/* DEPTH OF THIRD HORIZONTAL PLANE WHERE E-FIELD VALUES ARE WRITTEN TO OUTPUT FILE EVERY fd1.timpro
					TIMESTEPS (m) */
   	float delaya;		/* TIME DELAY BEFORE E-FIELD VALUES ALONG FIRST HORIZONTAL PLANE START TO BE WRITTEN TO 
					AN OUTPUT FILE (ns) */
   	float delayb;		/* TIME DELAY BEFORE E-FIELD VALUES ALONG SECOND HORIZONTAL PLANE START TO BE WRITTEN TO 
					AN OUTPUT FILE (ns) */
   	float delayc;		/* TIME DELAY BEFORE E-FIELD VALUES ALONG THIRD HORIZONTAL PLANE START TO BE WRITTEN TO 
					AN OUTPUT FILE (ns) */ 		
} fd1;