antenna.for

fortran code of dipole antenna

			exiyb1(i,ys,k)=exi(i,ys,k)
			exiyb1(i,ys-1,k)=exi(i,ys-1,k)
		end do
	end do

	do i=0,xs-1
		do j=0,ys
			exiza2(i,j,0)=exiza1(i,j,0)
			exiza2(i,j,1)=exiza1(i,j,1)
			exizb2(i,j,zs)=exizb1(i,j,zs)
			exizb2(i,j,zs-1)=exizb1(i,j,zs-1)
			exiza1(i,j,0)=exi(i,j,0)
			exiza1(i,j,1)=exi(i,j,1)
			exizb1(i,j,zs)=exi(i,j,zs)
			exizb1(i,j,zs-1)=exi(i,j,zs-1)

		end do
	end do

	do i=0,xs
		do j=0,ys-1
			eyiza2(i,j,0)=eyiza1(i,j,0)
			eyiza2(i,j,1)=eyiza1(i,j,1)
			eyizb2(i,j,zs)=eyizb1(i,j,zs)
			eyizb2(i,j,zs-1)=eyizb1(i,j,zs-1)
			eyiza1(i,j,0)=eyi(i,j,0)
			eyiza1(i,j,1)=eyi(i,j,1)
			eyizb1(i,j,zs)=eyi(i,j,zs)
			eyizb1(i,j,zs-1)=eyi(i,j,zs-1)
		end do
	end do

	do j=0,ys-1
		do k=0,zs
			eyixa2(0,j,k)=eyixa1(0,j,k)
			eyixa2(1,j,k)=eyixa1(1,j,k)
			eyixb2(xs,j,k)=eyixa1(xs,j,k)
			eyixb2(xs-1,j,k)=eyixa1(xs-1,j,k)
			eyixa1(0,j,k)=eyi(0,j,k)
			eyixa1(1,j,k)=eyi(1,j,k)
			eyixb1(xs,j,k)=eyi(xs,j,k)
			eyixb1(xs-1,j,k)=eyi(xs-1,j,k)
		end do
	end do

	do j=0,ys
		do k=0,zs-1
			ezixa2(0,j,k)=ezixa1(0,j,k)
			ezixa2(1,j,k)=ezixa1(1,j,k)
			ezixb2(xs,j,k)=ezixb1(xs,j,k)
			ezixb2(xs-1,j,k)=ezixb1(xs-1,j,k)
			ezixa1(0,j,k)=ezi(0,j,k)
			ezixa1(1,j,k)=ezi(1,j,k)
			ezixb1(xs,j,k)=ezi(xs,j,k)
			ezixb1(xs-1,j,k)=ezi(xs-1,j,k)
		end do
	end do

	do i=0,xs
		do k=0,zs-1
			eziya2(i,0,k)=eziya1(i,0,k)
			eziya2(i,1,k)=eziya1(i,1,k)
			eziyb2(i,ys,k)=eziyb1(i,ys,k)
			eziyb2(i,ys-1,k)=eziyb1(i,ys-1,k)
			eziya1(i,0,k)=ezi(i,0,k)
			eziya1(i,1,k)=ezi(i,1,k)
			eziyb1(i,ys,k)=ezi(i,ys,k)
			eziyb1(i,ys-1,k)=ezi(i,ys-1,k)
		end do
	end do

ccccccccccccccccccc E_FIELD  cccccccccccccccccccccccccccccccccccccccccccc
	do i=0,xs-1
		do j=1,ys-1
			do k=1,zs-1
				bbx=((hzi(i,j,k)-hzi(i,j-1,k))/dy)
     &			   -((hyi(i,j,k)-hyi(i,j,k-1))/dz)
				exi(i,j,k)=exi(i,j,k)+bbx*rh2
			end do
		end do
	end do

	do i=1,xs-1
		do j=0,ys-1
			do k=1,zs-1
				bby=((hxi(i,j,k)-hxi(i,j,k-1))/dz)
     &			   -((hzi(i,j,k)-hzi(i-1,j,k))/dx)
				eyi(i,j,k)=eyi(i,j,k)+bby*rh2
			end do
		end do
	end do

	do i=1,xs-1
		do j=1,ys-1
			do k=0,zs-1
				bbz=((hyi(i,j,k)-hyi(i-1,j,k))/dx)
     &			   -((hxi(i,j,k)-hxi(i,j-1,k))/dy)
				ezi(i,j,k)=ezi(i,j,k)+bbz*rh2
			end do
		end do
	end do

cccccccccccccccccc           PML        ccccccccccccccccccccccccccccccccccc

	do	i=1,xs-2
		do k=1,zs-1
			r1=exi(i,1,k)+exiya2(i,0,k)
			r2=exiya1(i,0,k)+exiya1(i,1,k)
			r3=exiya1(i,0,k+1)-2*exiya1(i,0,k)+exiya1(i,0,k-1)
     &		  +exiya1(i,1,k+1)-2*exiya1(i,0,k)+exiya1(i,1,k-1)
			r4=exiya1(i+1,0,k)-2*exiya1(i,1,k)+exiya1(i-1,1,k)
     &		  +exiya1(i+1,1,k)-2*exiya1(i,1,k)+exiya1(i-1,1,k)
			exi(i,0,k)=-exiya2(i,1,k)+uy1*r1+uy2*r2+uy3*r3+uy4*r4

			r1=exi(i,ys-1,k)+exiyb2(i,ys,k)
			r2=exiyb1(i,ys,k)+exiyb1(i,ys-1,k)
			r3=exiyb1(i,ys,k+1)-2*exiyb1(i,ys,k)+exiyb1(i,ys,k-1)
     &		 +exiyb1(i,ys-1,k+1)-2*exiyb1(i,ys-1,k)+exiyb1(i,ys-1,k-1)
			r4=exiyb1(i+1,ys,k)-2*exiyb1(i,ys,k)+exiyb1(i-1,ys,k)
     &		 +exiyb1(i+1,ys-1,k)-2*exiyb1(i,ys-1,k)+exiyb1(i-1,ys-1,k)
			exi(i,ys,k)=-exiyb2(i,ys-1,k)+uy1*r1+uy2*r2+uy3*r3+uy4*r4
	    end do
	end do

	do k=1,zs-1
	   r1=exi(0,1,k)-exiya1(0,0,k)
	   exi(0,0,k)=exiya1(0,1,k)+uy1*r1
         r1=exi(xs-1,1,k)-exiya1(xs-1,0,k)
	   exi(xs-1,0,k)=exiya1(xs-1,1,k)+uy1*r1
	   r1=exi(0,ys-1,k)-exiyb1(0,ys,k)
	   exi(0,ys,k)=exiyb1(0,ys-1,k)+uy1*r1
	   r1=exi(xs-1,ys-1,k)-exiyb1(xs-1,ys,k)
	   exi(xs-1,ys,k)=exiyb1(xs-1,ys-1,k)+uy1*r1
	end do

	do i=1,xs-2
		do j=1,ys-1
			r1=exi(i,j,k)+exiza2(i,j,0)
			r2=exiza1(i,j,0)+exiza1(i,j,1)
			r3=exiza1(i+1,j,0)-2*exiza1(i,j,0)+exiza1(i-1,j,0)
     &		  +exiza1(i+1,j,1)-2*exiza1(i,j,1)+exiza1(i-1,j,1)
			r4=exiza1(i,j+1,0)-2*exiza1(i,j,0)+exiza1(i,j-1,0)
     &		  +exiza1(i,j+1,1)-2*exiza1(i,j,1)+exiza1(i,j-1,1)
		    exi(i,j,0)=-exiza2(i,j,1)+uz1*r1+uz2*r2+uz3*r3+uz4*r4

	        r1=exi(i,j,zs-1)+exizb2(i,j,zs)

			r2=exizb1(i,j,zs)+exizb1(i,j,zs-1)
	        r3=exizb1(i+1,j,zs)-2*exizb1(i,j,zs)+exizb1(i-1,j,zs)
     &		  +exizb1(i+1,j,zs)-2*exizb1(i,j,zs-1)+exizb1(i-1,j,zs-1)
			r4=exizb1(i,j+1,zs)-2*exizb1(i,j,zs)+exizb1(i,j-1,zs)
     &	     +exizb1(i,j+1,zs-1)-2*exizb1(i,j,zs-1)+exizb1(i,j-1,zs-1)
			exi(i,j,zs)=-exizb2(i,j,zs-1)+uz1*r1+uz2*r2+uz3*r3+uz4*r4
	      end do
	end do

	do j=1,ys-1
	    r1=eyi(0,j,1)-eyixa1(0,j,0)
	    exi(0,j,0)=exiza1(0,j,1)+uz1*r1
		r1=exi(xs-1,j,1)-exiza1(xs-1,j,0)
		exi(xs-1,j,0)=exiza1(xs-1,j,1)+uz1*r1
		r1=exi(xs-1,j,zs-1)-exizb1(xs-1,j,zs)
		exi(xs-1,j,zs)=exizb1(xs-1,j,zs-1)+uz1*r1
	end do

	do i=0,xs-1
	   r1=exi(i,0,1)-exiza1(i,0,0)
	   exi(i,0,0)=exiza1(i,0,1)+uz1*r1
	   r1=exi(i,ys,1)-exiza1(i,ys,0)
	   exi(i,ys,0)=exiza1(i,ys,1)+uz1*r1
	   r1=exi(i,0,zs-1)-exizb1(i,0,zs)
	   exi(i,0,zs)=exizb1(i,0,zs-1)+uz1*r1
	   r1=exi(i,ys,zs-1)-exizb1(i,ys,zs)
	   exi(i,ys,zs)=exizb1(i,ys,zs-1)+uz1*r1
	end do
	
	do j=1,ys-2
	   do i=1,xs-1
	      r1=eyi(i,j,1)+eyiza2(i,j,0)
		  r2=eyiza1(i,j,0)+eyiza1(i,j,1)
		  r3=eyiza1(i+1,j,0)-2*eyiza1(i,j,0)+eyiza1(i-1,j,0)
     &	    +eyiza1(i+1,j,0)-2*eyiza1(i,j,1)+eyiza1(i-1,j,1)
            r4=eyiza1(i,j+1,0)-2*eyiza1(i,j,0)+eyiza1(i,j-1,0)
     &	    +eyiza1(i,j+1,1)-2*eyiza1(i,j,1)+eyiza1(i,j-1,1)
     		  eyi(i,j,0)=-eyiza2(i,j,1)+uz1*r1+uz2*r2+uz3*r3+uz4*r4
		  
		  r1=eyi(i,j,zs-1)+eyizb2(i,j,zs)
		  r2=eyizb1(i,j,zs)+eyizb1(i,j,zs-1)
		  r3=eyizb1(i+1,j,zs)-2*eyizb1(i,j,zs)+eyizb1(i-1,j,zs)
     &	    +eyizb1(i+1,j,zs-1)-2*eyizb1(i,j,zs-1)+eyizb1(i-1,j,zs-1)
     		  r4=eyizb1(i,j+1,zs)-2*eyizb1(i,j,zs)+eyizb1(i,j-1,zs)
     &	    +eyizb1(i,j,zs-1)-2*eyizb1(i,j,zs-1)+eyizb1(i,j-1,zs-1)
            eyi(i,j,zs)=-eyizb2(i,j,zs-1)+uz1*r1+uz2*r2+uz3*r3+uz4*r4
		end do
	end do  								

	do i=1,xs-1
	   r1=eyi(i,0,1)-eyiza1(i,0,0)
	   eyi(i,0,0)=eyiza1(i,0,1)+uz1*r1
	   r1=eyi(i,ys-1,1)-eyiza1(i,ys-1,0)
	   eyi(i,ys-1,0)=eyiza1(i,ys-1,1)+uz1*r1
	   r1=eyi(i,0,zs-1)-eyizb1(i,0,zs)
	   eyi(i,0,zs)=eyizb1(i,0,zs-1)+uz1*r1
	   r1=eyi(i,ys-1,zs-1)-eyizb1(i,ys-1,zs)
	   eyi(i,ys-1,zs)=eyizb1(i,ys-1,zs-1)+uz1*r1
	end do

      do j=1,ys-2
	   do i=1,zs-1
	      r1=eyi(1,j,k)+eyixa2(0,j,k)
		  r2=eyixa1(0,j,k)+eyixa1(1,j,k)
		  r3=eyixa1(0,j+1,k)-2*eyixa1(0,j,k)+eyixa1(0,j-1,k)
     &	    +eyixa1(1,j+1,k)-2*eyixa1(1,j,k)+eyixa1(1,j-1,k)
            r4=eyixa1(0,j,K+1)-2*eyixa1(0,j,k)+eyixa1(0,j,k-1)
     &	    +eyixa1(1,j,k+1)-2*eyixa1(1,j,k)+eyixa1(1,j,k-1)
     		  eyi(0,j,k)=-eyixa2(1,j,k)+ux1*r1+ux2*r2+ux3*r3+ux4*r4
		  
		  r1=eyi(xs-1,j,k)+eyizb2(xs,j,k)
		  r2=eyixb1(xs,j,k)+eyixb1(xs-1,j,k)
		  r3=eyixb1(xs,j+1,k)-2*eyixb1(xs,j,k)+eyixb1(xs,j-1,k)
     &	    +eyixb1(xs-1,j+1,k)-2*eyixb1(xs-1,j,k)+eyixb1(xs-1,j-1,k)
     		  r4=eyixb1(xs,j,k+1)-2*eyixb1(xs,j,k)+eyixb1(xs,j,k-1)
     &	    +eyixb1(xs-1,j,k+1)-2*eyixb1(xs-1,j,k)+eyixb1(xs-1,j,k-1)
            eyi(i,j,zs)=-eyizb2(i,j,zs-1)+uz1*r1+uz2*r2+uz3*r3+uz4*r4
		end do
	end do  			

      do k=1,zs-1
	   r1=eyi(1,0,k)-eyixa1(0,0,k)
	   eyi(0,0,k)=eyixa1(1,0,k)+ux1*r1
	   r1=eyi(1,ys-1,k)-eyixa1(0,ys-1,k)
	   eyi(0,ys-1,k)=eyixa1(1,ys-1,k)+ux1*r1
	   r1=eyi(xs-1,0,k)-eyixb1(xs,0,k)
	   eyi(xs,0,k)=eyixb1(xs-1,0,k)+ux1*r1
	   r1=eyi(xs-1,ys-1,k)-eyixb1(xs,ys-1,k)
	   eyi(xs,ys-1,k)=eyixb1(xs-1,ys-1,k)+ux1*r1
	end do

	do i=0,ys-1
	   r1=eyi(0,j,1)-eyiza1(0,j,0)
	   eyi(0,j,0)=eyiza1(0,j,1)+uz1*r1
	   r1=eyi(xs,j,1)-eyiza1(xs,j,0)
	   eyi(xs,j,0)=eyiza1(xs,j,1)+uz1*r1
	   r1=eyi(0,j,zs-1)-eyizb1(0,j,zs)
	   eyi(0,j,zs)=eyizb1(0,j,zs-1)+uz1*r1
	   r1=eyi(xs,j,zs-1)-eyizb1(xs,j,zs)
	   eyi(xs,j,zs)=eyizb1(xs,j,zs-1)+uz1*r1
	end do

	do k=1,zs-2
	   do j=1,ys-1
	      r1=ezi(1,j,k)+ezixa2(0,j,k)
		  r2=ezixa1(0,j,k)+ezixa1(1,j,k)
		  r3=ezixa1(0,j+1,k)-2*ezixa1(0,j,k)+eyixa1(0,j-1,k)
     &	    +ezixa1(1,j+1,k)-2*ezixa1(1,j,k)+eyixa1(1,j-1,k)
            r4=ezixa1(0,j,K+1)-2*ezixa1(0,j,k)+eyixa1(0,j,k-1)
     &	    +ezixa1(1,j,k+1)-2*ezixa1(1,j,k)+eyixa1(1,j,k-1)
     		  ezi(0,j,k)=-ezixa2(1,j,k)+ux1*r1+ux2*r2+ux3*r3+ux4*r4
		  
		  r1=ezi(xs-1,j,k)+ezixb2(xs,j,k)
		  r2=ezixb1(xs,j,k)+ezixb1(xs-1,j,k)
		  r3=ezixb1(xs,j+1,k)-2*ezixb1(xs,j,k)+ezixb1(xs,j-1,k)
     &	    +ezixb1(xs-1,j+1,k)-2*ezixb1(xs-1,j,k)+ezixb1(xs-1,j-1,k)
     		  r4=ezixb1(xs,j,k+1)-2*ezixb1(xs,j,k)+ezixb1(xs,j,k-1)
     &	    +ezixb1(xs-1,j,k+1)-2*ezixb1(xs-1,j,k)+ezixb1(xs-1,j,k-1)
            ezi(i,j,zs)=-ezixb2(i,j,zs-1)+uz1*r1+uz2*r2+uz3*r3+uz4*r4
		end do
	end do  	
	
	do j=1,ys-1
	   r1=ezi(1,j,0)-ezixa1(0,j,0)
	   ezi(0,j,0)=ezixa1(1,j,0)+ux1*r1
	   r1=eyi(1,j,zs-1)-ezixa1(0,j,zs-1)
	   ezi(0,j,zs-1)=ezixa1(1,j,zs-1)+ux1*r1
	   r1=ezi(xs-1,j,0)-ezixb1(xs,j,0)
	   ezi(xs,j,0)=ezixb1(xs-1,j,0)+ux1*r1
	   r1=ezi(xs-1,j,zs-1)-ezixb1(xs,j,zs-1)
	   ezi(xs,j,zs-1)=ezixb1(xs-1,j,zs-1)+ux1*r1
	end do

	do	k=1,zs-2
		do i=1,xs-1
			r1=ezi(i,1,k)+eziya2(i,0,k)
			r2=eziya1(i,0,k)+eziya1(i,1,k)
			r3=eziya1(i,0,k+1)-2*eziya1(i,0,k)+eziya1(i,0,k-1)
     &		  +eziya1(i,1,k+1)-2*eziya1(i,0,k)+eziya1(i,1,k-1)
			r4=eziya1(i+1,0,k)-2*eziya1(i,1,k)+eziya1(i-1,1,k)
     &		  +eziya1(i+1,1,k)-2*eziya1(i,1,k)+eziya1(i-1,1,k)
			ezi(i,0,k)=-eziya2(i,1,k)+uy1*r1+uy2*r2+uy3*r3+uy4*r4

			r1=ezi(i,ys-1,k)+eziyb2(i,ys,k)
			r2=eziyb1(i,ys,k)+eziyb1(i,ys-1,k)
			r3=eziyb1(i,ys,k+1)-2*eziyb1(i,ys,k)+eziyb1(i,ys,k-1)
     &		 +eziyb1(i,ys-1,k+1)-2*eziyb1(i,ys-1,k)+eziyb1(i,ys-1,k-1)
			r4=eziyb1(i+1,ys,k)-2*eziyb1(i,ys,k)+eziyb1(i-1,ys,k)
     &		 +eziyb1(i+1,ys-1,k)-2*eziyb1(i,ys-1,k)+eziyb1(i-1,ys-1,k)
			ezi(i,ys,k)=-eziyb2(i,ys-1,k)+uy1*r1+uy2*r2+uy3*r3+uy4*r4
	    end do
	end do
	
	do i=1,xs-1
	   r1=ezi(i,1,0)-eziya1(i,0,0)
	   ezi(i,0,0)=eziya1(i,1,0)+uy1*r1
	   r1=ezi(i,1,zs-1)-eziya1(i,0,zs-1)
	   ezi(i,0,zs-1)=eziya1(i,1,zs-1)+uy1*r1
	   r1=ezi(i,ys-1,0)-eziyb1(i,ys,0)
	   ezi(i,ys,0)=eziyb1(i,ys-1,0)+uy1*r1
	   r1=ezi(i,ys-1,zs-1)-eziyb1(i,ys,zs-1)
	   ezi(i,ys,zs-1)=eziyb1(i,ys-1,zs-1)+uy1*r1
	end do

	do k=0,zs-1
	   r1=ezi(0,1,k)-eziya1(0,0,k)
	   ezi(0,0,k)=eziya1(0,1,k)+uy1*r1
	   r1=ezi(xs,1,k)-eziya1(xs,0,k)
	   ezi(xs,0,k)=eziya1(xs,1,k)+uy1*r1
	   r1=ezi(0,ys-1,k)-eziyb1(0,ys,k)
	   ezi(0,ys,k)=eziyb1(0,ys-1,k)+uy1*r1
	   r1=ezi(xs,ys-1,k)-eziyb1(xs,ys,k)
	   ezi(xs,ys,k)=eziyb1(xs,ys-1,k)+uy1*r1
	end do
	return
	end

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC