my $maxProc=4;
my $pm = new Parallel::ForkManager($maxProc);

for(my $idy=0;$idy<=$#ySurf;$idy++)
{
	my $currY = $ySurf[$idy];
	for(my $alfaIdx=0;$alfaIdx<=$#alfas;$alfaIdx++)
	{
		$pm->start and next;
		
		my $alfa = $alfas[$alfaIdx];
		
		my $norVectX = cos($alfa);
		my $norVectY = 0;
		my $norVectZ = sin($alfa);

		my $norTransf=TransfNorVec($input->{az}, $input->{el}, $norVectX, $norVectY, $norVectZ, $Side);
		
		my $intSurfx = $input->{MiddleX} - $input->{r}*$norVectX;
		my $intSurfy = $currY;
		my $intSurfz = ($input->{MiddleZ}+$input->{zOffset})-$input->{r}*$norVectZ;
		
		my $center = Divide($input->{w}, $input->{h}, $input->{M}, $input->{N}, $intSurfx, $intSurfy, $intSurfz, $alfa);
		
		my $surf=0;

		for(my $idx=0;$idx<$center->{count};$idx++)
		{
			my $xT = $center->{Px}[$idx];
			my $yT = $center->{Py}[$idx];
			my $zT = $center->{Pz}[$idx];
			
			my $xTr = coordSystTransf($input->{az}, $input->{el}, $input->{MiddleX}, $input->{MiddleY}, $input->{MiddleZ}, $input->{ZeroX}, $input->{ZeroY}, $input->{ZeroZ}, $xT, $yT, $zT, $Side);
			
			my $sqDiag = ($input->{w}/$input->{M})**2+($input->{h}/$input->{N})**2;
			
			my $P = new Points();
			$P->{count}=0;
			
			for(my $dataIdx=0;$dataIdx<$data->{count};$dataIdx++)
			{
				my $inSphere = Check($sqDiag, $xTr->{xe}, $xTr->{ye}, $xTr->{ze}, $data->{abs_x}[$dataIdx], $data->{abs_y}[$dataIdx], $data->{abs_z}[$dataIdx]);
				
				if($inSphere == 1)
				{
					$P->{Px}[$P->{count}]=($data->{magX}[$dataIdx])*cos($data->{phX}[$dataIdx]);
					$P->{Py}[$P->{count}]=($data->{magY}[$dataIdx])*cos($data->{phY}[$dataIdx]);
					$P->{Pz}[$P->{count}]=($data->{magZ}[$dataIdx])*cos($data->{phZ}[$dataIdx]);							
					$P->{count}++;
				}
			}
			
			my $Poynt = new Trans();
			if(defined($P->{Px}[0]) && defined($P->{Py}[0]) && defined($P->{Pz}[0]))
			{
				$Poynt->{xe} = summe($P->{Px})/$P->{count};
				$Poynt->{ye} = summe($P->{Py})/$P->{count};
				$Poynt->{ze} = summe($P->{Pz})/$P->{count};
				
				my $innerProd = $Poynt->{xe}*$norTransf->{xe}+$Poynt->{ye}*$norTransf->{ye}+$Poynt->{ze}*$norTransf->{ze};
				$surf += ((($input->{w})/$input->{M})*(($input->{h})/$input->{N}))*$innerProd;
			}
		}
		
		
		print "\ny= ".$currY."; alfa= ".$alfa .";\n";
		
		if($surf == 0)
		{
			print 1e-100 ."dB\n";
		}else{
			print sprintf("%.4f", 20*log(abs($surf))) ." dB\n";
		}
		
		$pm->finish;
	}
	$pm->wait_all_children;
}