#!/usr/bin/perl
use Statistics::OLS;
use Math::FFT;
use Math::Stat;
use Math::Spline qw(spline linsearch binsearch);
use Math::Derivative qw(Derivative2);
use Math::Approx;
use Astro::FITS::CFITSIO qw( :longnames :constants );
use PDL;
use PDL::Fit::Polynomial; 
use PDL::Filter::Linear;
use PGPLOT;  # Load PGPLOT module
use PDL::Fit::Gaussian;
use PDL::Slatec;
use PDL::Image2D;
#use PDL::Matrix;


$DIR_DATA="/disk-a/sanchez/ppak/legacy/DATA/V500/reduced_v1.3c/";
#$DIR_DATA="/disk-b/sanchez/ppak/legacy/DATA/V500/reduced_v1.3c/";

$max_size=10;
$frac=0.9;


$n=0;
open(DIR,"ls HII.*.rss.fits.gz |");
while($file=<DIR>) {
    chop($file);
    $name=$file;
    $cut=".rss.fits.gz";
    $name =~ s/$cut//;
    $cut="HII.";
    $name =~ s/$cut//;
    $FILE[$n]=$file;
    $NAME[$n]=$name;
    $n++;    
}
close(DIR);

open(LOG,">ana_seg_gas.log");
for ($i=0;$i<$n;$i++) {    
    $call="./ana_single_gas.pl ".$NAME[$i];
    mycall($call);
    $do_ftp=1;
    if ($do_ftp==1) {
	open(FTP,">ftp.input");
	print FTP "open ftp://mycalifa:myafi07lAc\@ftp.caha.es\n";
	print FTP "cd /TMP/GAS_tests/FIT3D/SEG6_v1.3c/maps/\n";
	print FTP "mput -f map.*_*.$NAME[$i]_*.fits.gz\n";
	close(FTP);
	$call="ncftp < ftp.input &";
	mycall($call);
    }

}
close(LOG);
exit;

sub mycall {                                      
                                                                   my $call=$_[0];

														     print "$call | PROCESSING\n";                                                                                 
    system($call);                                                                                                
    print LOG "$call\n";                                                                                          
    print "DONE\n";                                                                                              
}