#!/usr/bin/perl -w
##############################################################################################
#############      DELOVANJE PROGRAMA  !!!!!!!!!!!!!!!!!!           #################
#Program je namenjen temu, da prebere datoteko datoteka.dat, v kateri je za vsak termin poleg casa podana tudi temperatura in
# vlaga. Program zgenerira datoteko razredi6_14.txt. V nadaljevanju spravim stolpec temperatur in vlage v vektor temperature 
#in vlage, ki so bile izmerjene med 6 in 14h. Temperature razdeli v razrede in pregleda, v kater ustrezni razred vlage pade.
# Tako dobim porazdelitev temperatur pri doloceni vlagi med 6 in 14h.
#############      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,">razredi6_14h.txt") || die "ne morem narediti datoteke razredi.txt";

####Preberem vrstice koordinat
$ij=1;
while ( defined($koordinati=<VHODKOOR>)){
        chomp($koordinati);
#Vzamem le termine od 6 do 14h
	$ura=1*substr($koordinati,11,3);
	if(6<=$ura && $ura<14) {      
        	$xy[$ij]=$koordinati;
		$ij++;} #konec if(6<=$ura && $ura<14) {
	} #konec while ( defined($koordinati=<VHODKOOR>)){
$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;