#!/usr/bin/perl -w
##############################################################################################
#############      DELOVANJE PROGRAMA  !!!!!!!!!!!!!!!!!!           #################
#Program je namenjen temu, da prebere datoteko mbmb1.dat, v kateri je za vsak termin poleg casa podana tudi temperatura in
# vlaga. Program zgenerira datoteko razredi.txt. V nadaljevanju spravim stolpec temperatur in vlage v vektor temperature in
# vlage. Temperature razdeli v razrede in pregleda, v kater ustrezni razred vlage pade. Tako dobim porazdelitev temperatur
# pri doloceni vlagi. 
#############      DELOVANJE PROGRAMA  !!!!!!!!!!!!!!!!!!           #################
##############################################################################################
print "\nProgram bo prebral datoteko in podatke razdelil v podrazrede\n";


open (VHODKOOR,"<datoteka.dat") || die "ne morem odpreti datoteke mbmb1.dat"; #odprem datoteko s koordinatami
open (PISIV,">razredi.txt") || die "ne morem narediti datoteke razredi.txt";

####Preberem vrstice koordinat
$ij=1;
while ( defined($koordinati=<VHODKOOR>)){
        chomp($koordinati);      
        $xy[$ij]=$koordinati;
	$ij++;}
$vrsticakoord=$ij;

print"$vrsticakoord\n";    
close VHODKOOR;

##################################
##Definiram vektor temperature in vlage
##################################

for($i=1;$i<$vrsticakoord;$i++){ ###gledam vse vrstice v radarskih koordinatah

#Branje temperatur
	$T=1*substr($xy[$i],20,7);
	$temp[$i]=$T;
        $H=1*substr($xy[$i],28,4);
        $vlaga[$i]=$H;

	}
##################################
##Definiram vektor temperature in vlage ########KONEC###############
##################################




#definiram razrede temperatur, po katerih gledam
for($t=-16.5;$t<=39.5;$t++){
	$k=0;
	$tmin=$t;
	$tmax=$t+1;
	for($x=1;$x<=84;$x++){$vrstica[$x]=0;} #vektor pogostosti vlage v posameznem intervalu temperature
	
#Berem vektor temperatur in si zapomnim mesta, kjer se za ustrezni interval nahajajo temperature
	for($i=1;$i<$vrsticakoord;$i++){
               	if($temp[$i]>$tmin && $temp[$i]<=$tmax){$k++;$pravi[$k]=$i;}
		} #konec for($i=1;$i<$vrsticakoord;$i++){  #zapomnim mesta temperatur iz intervala

#V vektorju vlag bom preveril ali se nahajajo v pravem intervalu za vlago
	for($v=1;$v<=84;$v++){
		$vmin=16.5+$v-1;
		$vmax=16.5+$v;
		for($n=1;$n<=$k;$n++){#preverjam z vektorjem temperatur, ki so v pravem intervalu		
							
			if($vlaga[$pravi[$n]]>$vmin && $vlaga[$pravi[$n]]<=$vmax){$vrstica[$v]=$vrstica[$v]+1;}
				}#konec for($n=1;$n<=$k;$k++){

			} #konec for($v=16.5;$v<=100.5;$v++){
	print PISIV "$tmin|$tmax|";
	for($j=1;$j<=84;$j++){	
	print PISIV "$vrstica[$j]|";	
	
	} #konec for($j=1;$j<=84;$j++){
	print PISIV "\n";
	}# konec for($t=-16.5;$t<39.5;$t++){


close PISIV;