#! /usr/bin/perl

####please update the directories before use

# Directory to run the code
$homedir = "**";  

#ET Archive under the folder of /data/ESI, which contains files "FPPM_SMN365.bin"
$climinputfd = "../data/ESI";

#ET climatology; this directory should be consistent with GETD_composite_ESI.cfg
$climfd = "../data/static/ ETClim";

$filekey="FPPM_SMN365.bin";

$tmpfd = "$homedir/tmp";

system("rm -rf $climfd");
system("mkdir $climfd");

$climyearstart = 2013;
$climyearend = 2023;
$imgW = 2466;
$imgH = 2520;

#build climatology
for($doy=1;$doy<=365;$doy++){
print "Processing $doy\n";

if($doy<10){
$doystr="00$doy";}
elsif($doy<100){
$doystr="0$doy";
}else{
$doystr = "$doy";
}

system("rm -rf $tmpfd");
system("mkdir $tmpfd");
for($year=$climyearstart; $year<=$climyearend; $year++){
for($ddd=${doy}-6;$ddd<=$doy;$ddd++){
if($ddd<=0){
$ddn=$ddd+365;
$yyyyddd = ($year-1) * 1000 + $ddn;
}else
{
$yyyyddd = $year * 1000 + $ddd;
}
print "$yyyyddd\n";
$infile = "$climinputfd/$yyyyddd/${filekey}";
$linkfile = "$tmpfd/data_$yyyyddd.bin";
symlink($infile, $linkfile);
}
}
$climmeanfile="$climfd/AVG_SMN_${doystr}.dat";
$climstdfile="$climfd/SD_SMN_${doystr}.dat";
$climnfile="$climfd/N_${doystr}.dat";
system("$homedir/calclim $imgW $imgH $tmpfd $climmeanfile $climstdfile $climnfile");

}