rendered paste body#!/usr/bin/perl
use strict;
use Geo::Coordinates::UTM;
open (TOWNSHIPS, "<townships");
open (OFILE, ">>combined");
my @qsectiondata;
my @townshipdata;
my @sectiondata;
my $townshipline;
my $sectionline;
my $qsectionline;
foreach $townshipline ( <TOWNSHIPS> )
{
# print "came back here\n";
@townshipdata = split(" ",$townshipline);
my $townshipPPID = $townshipdata[0];
my $meridian = $townshipdata[4];
my $range = $townshipdata[3];
my $township = $townshipdata[2];
open (SECTIONS, "<sections");
open (QSECTIONS, "<quartersections");
foreach $sectionline ( <SECTIONS> )
{
@sectiondata = split(" ",$sectionline);
my $sectionPPID = $sectiondata[0];
my $sectionPPID2 = $sectiondata[2];
my $sectionnum = $sectiondata[1]-1;
# print "TownshipPPID $townshipPPID SectionPPID $sectionPPID\n";
if ( $townshipPPID == $sectionPPID )
{
# print "SECTION PRINT $townshipPPID = $sectionPPID\n";
# print "sectionPPID check $sectiondata[7]";
foreach $qsectionline ( <QSECTIONS> )
{
@qsectiondata = split(" ",$qsectionline);
my $qsectionPPID2 = $qsectiondata[2];
# print "Qsectoin $qsectionPPID2 **********";
if ( $sectionPPID2 == $qsectionPPID2 )
{
# print "QSECTION PRINT $townshipPPID = $sectionPPID also qsections $qsectionPPID2 = sections $sectionPPID2\n";
my ($latitude,$longitude)=utm_to_latlon(5,'13V',$qsectiondata[3],$qsectiondata[4]);
printf OFILE " $meridian $range $township $sectionnum $qsectiondata[1] $latitude $longitude\n";
}
else
{
# print " no match";
}
}
}
# else{print "no match";}
}
}