bianjie.h

3D-FDTD垂直双阵子程序

 int bianjie()
  {int loop1,loop2,loop3;
   for(loop1=0;loop1<Nx;loop1++)
	   for(loop3=0;loop3<Nz+1;loop3++)       //Ex y界面
   { Exy[loop1][0][loop3]=0;
     Exz[loop1][0][loop3]=0;
	 Exy[loop1][Ny][loop3]=0;
     Exz[loop1][Ny][loop3]=0;
   }

  for(loop1=0;loop1<Nx;loop1++)
	   for(loop2=0;loop2<Ny+1;loop2++)       //Ex z界面
   { Exy[loop1][loop2][0]=0;
     Exz[loop1][loop2][0]=0;
	 Exy[loop1][loop2][Nz]=0;
     Exz[loop1][loop2][Nz]=0;
   }


  for(loop2=0;loop2<Ny;loop2++)
	   for(loop3=0;loop3<Nz+1;loop3++)       //Ey x界面
   { Eyx[0][loop2][loop3]=0;
     Eyz[0][loop2][loop3]=0;
	 Eyx[Nx][loop2][loop3]=0;
     Eyz[Nx][loop2][loop3]=0;
   }
  
 for(loop1=0;loop1<Nx+1;loop1++)
	   for(loop2=0;loop2<Ny;loop2++)       //Ey z界面
   { Eyx[loop1][loop2][0]=0;
     Eyz[loop1][loop2][0]=0;
	 Eyx[loop1][loop2][Nz]=0;
     Eyz[loop1][loop2][Nz]=0;
   }
 
 for(loop2=0;loop2<Ny+1;loop2++)
	   for(loop3=0;loop3<Nz;loop3++)       //Ez x界面
   { Ezx[0][loop2][loop3]=0;
     Ezy[0][loop2][loop3]=0;
	 Ezx[Nx][loop2][loop3]=0;
     Ezy[Nx][loop2][loop3]=0;
   }
  
 for(loop1=0;loop1<Nx+1;loop1++)
	   for(loop3=0;loop3<Nz;loop3++)       //Ez y界面
   { Ezx[loop1][0][loop3]=0;
     Ezy[loop1][0][loop3]=0;
	 Ezx[loop1][Ny][loop3]=0;
     Ezy[loop1][Ny][loop3]=0;
   }

 
        return(1);
  }