#!/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";}
      }
 }