part_6.br

The Finite Difference Time Domain Method for Electromagnetics With MATLAB Simulations Atef Elshe

	// Open Datafile For Observation Port
	pFile2 = fopen("Port.txt","wt");

	// Do the requested number of iterations
	for(i=0; i<N; i++){

        // We use to 2 sets of field streams to have 
	// a "New" and "Old" set
           if(i%2==0) {

		process_field_H(Hx, Chxh, Chxe, Hy, Chyh, Chye, 
		  Ez, Khx, Khy, Bhx, Bhy, Chx, Chy, psiHxy, psiHyx, 
		  dx, dy, it, o_Hx, o_Hy, o_psiHxy, o_psiHyx);

		process_field_Ez(Ez, Ceze, Cezh, o_Hx, o_Hy, Cs, 
		  Kex, Key, Bex, Bey, Cex, Cey, psiEzx, psiEzy, 
		  dx, dy, gauss[i], it, o_Ez, o_psiEzx, o_psiEzy);

		// Save a single Ez point for use later
		// Ez point is here is 10 cells from center
		Copy(o_Ez.domain( int2( (xsize/2)-10 , ysize/2 ), 
			  int2( (xsize/2)-10 , ysize/2 ) ), Obs);

		streamWrite(Obs,outputPort);
		fprintf (pFile2,"%e\n",outputPort[0]);

	   } else {

		process_field_H(o_Hx, Chxh, Chxe, o_Hy, Chyh, Chye,
		  o_Ez, Khx, Khy, Bhx, Bhy, Chx, Chy, o_psiHxy, 
		  o_psiHyx, dx, dy, it, Hx, Hy, psiHxy, psiHyx);

		process_field_Ez(o_Ez, Ceze, Cezh, Hx, Hy, Cs, Kex, 
		  Key, Bex, Bey, Cex, Cey, o_psiEzx, o_psiEzy, dx, 
		  dy, gauss[i], it, Ez, psiEzx, psiEzy);

		// Save a single Ez point for use later
		// Ez point is here is 10 cells from center

		Copy(Ez.domain( int2( (xsize/2)-10 , ysize/2 ), 
		     int2( (xsize/2)-10 , ysize/2 ) ), Obs);

		streamWrite(Obs,outputPort);
		fprintf (pFile2,"%e\n",outputPort[0]);

	   }

	   fprintf (pFile,"%e %e\n",gauss[i],t[i]);
	}