From 25e57e5cd17a178cb487cabb3a07ca4779c4cf68 Mon Sep 17 00:00:00 2001 From: curt Date: Tue, 27 May 2003 18:57:40 +0000 Subject: [PATCH] This script should be moved because it really isn't dafift specific. It builds the ILS database from several sources and uses an external compiled utility to crunch the math to straighten out the approaches whenever we can match them to specific runway headings and threshold locations. This should be *always* but our runway data is horribly outdated so that needs to be updated soon to get into sync. --- scripts/perl/dafif/build_ils.pl | 714 ++++++++++++++++++++++++++++++++ 1 file changed, 714 insertions(+) create mode 100755 scripts/perl/dafif/build_ils.pl diff --git a/scripts/perl/dafif/build_ils.pl b/scripts/perl/dafif/build_ils.pl new file mode 100755 index 000000000..b30db0075 --- /dev/null +++ b/scripts/perl/dafif/build_ils.pl @@ -0,0 +1,714 @@ +#!/usr/bin/perl + +######################################################################## +# Convert DAFIFT ARPT/ILS.TXT to FlightGear format. +######################################################################## + +use strict; + +my($faa_ils_file) = shift(@ARGV); +my($dafift_arpt_file) = shift(@ARGV); +my($dafift_ils_file) = shift(@ARGV); +my($fgfs_ils_file) = shift(@ARGV); +my($fgfs_apt_file) = shift(@ARGV); +my($output_file) = shift(@ARGV); + +my($calc_loc) + = "/home/curt/projects/FlightGear-0.9/source/src/Airports/calc_loc"; + +die "Usage: $0 " . + " \n" + if !defined($faa_ils_file) || !defined($dafift_arpt_file) + || !defined($dafift_ils_file) || !defined($fgfs_ils_file) + || !defined($fgfs_apt_file) || !defined($output_file); + +my( %CODES ); +my( %CodesByICAO ); +my( %ILS ); +my( %AIRPORTS ); +my( %RUNWAYS ); + + +&load_dafift( $dafift_arpt_file, $dafift_ils_file ); +&load_faa( $faa_ils_file ); +&load_fgfs( $fgfs_ils_file ); +&load_fgfs_airports( $fgfs_apt_file ); +&fix_localizer(); + +&write_result( $output_file ); + +exit; + + +######################################################################## +# Process DAFIFT data +######################################################################## + +sub load_dafift() { + my( $arpt_file ) = shift; + my( $ils_file ) = shift; + + my( $record ); + + my( $id, $rwy, $type ); + my( $has_dme, $has_gs, $has_loc, $has_im, $has_mm, $has_om ); + my( $dme_lon, $dme_lat, $dme_elev, $dme_bias ); + my( $gs_lon, $gs_lat, $gs_elev, $gs_angle ); + my( $loc_type, $loc_lon, $loc_lat, $loc_elev, $loc_freq, $loc_hdg, + $loc_width, $loc_id ); + my( $im_lon, $im_lat, $mm_lon, $mm_lat, $om_lon, $om_lat ); + + # load airport file so we can lookup ICAO from internal ID + + open( ARPT, "<$arpt_file" ) || die "Cannot open DAFIFT: $arpt_file\n"; + + ; # skip header line + + while ( ) { + chomp; + my(@F) = split(/\t/); + my($icao) = $F[3]; + if ( length($icao) < 3 ) { + if ( length( $F[4] ) >= 3 ) { + $icao = $F[4]; + } else { + $icao = "[none]"; + } + } + $CODES{$F[0]} = $icao; + $CodesByICAO{$icao} = 1; + # print "$F[0] - $icao\n"; + } + + # Load the DAFIFT ils file + + my( $last_id, $last_rwy ) = ("", ""); + + open( DAFIFT_ILS, "<$ils_file" ) || die "Cannot open DAFIFT: $ils_file\n"; + + ; # skip header line + + while ( ) { + chomp; + my @F = split(/\t/); + $id = $F[0]; + $rwy = $F[1]; + + if ( $last_id ne "" && ($last_id ne $id || $last_rwy ne $rwy) ) { + # just hist the start of the next record, dump the current data + + if ( ! $has_gs ) { + ( $gs_elev, $gs_angle, $gs_lat, $gs_lon ) = ( 0, 0, 0, 0 ); + } + if ( ! $has_dme ) { + ( $dme_lat, $dme_lon ) = ( 0, 0 ); + } + if ( ! $has_om ) { + ( $om_lat, $om_lon ) = ( 0, 0 ); + } + if ( ! $has_mm ) { + ( $mm_lat, $mm_lon ) = ( 0, 0 ); + } + if ( ! $has_im ) { + ( $im_lat, $im_lon ) = ( 0, 0 ); + } + if ( $ILS{$CODES{$last_id} . $last_rwy} eq "" ) { + print "DAFIFT adding: $CODES{$last_id} - $last_rwy\n"; + &safe_add_record( $CODES{$last_id}, $last_rwy, "ILS", + $loc_freq, $loc_id, $loc_hdg, $loc_lat, + $loc_lon, $gs_elev, $gs_angle, $gs_lat, + $gs_lon, $dme_lat, $dme_lon, $om_lat, + $om_lon, $mm_lat, $mm_lon, $im_lat, + $im_lon ); + } + + $has_dme = 0; + $has_gs = 0; + $has_loc = 0; + $has_im = 0; + $has_mm = 0; + $has_om = 0; + } + + $type = $F[2]; + if ( $type eq "D" ) { + # DME entry + $has_dme = 1; + $dme_lon = make_dcoord( $F[16] ); + $dme_lat = make_dcoord( $F[14] ); + $dme_elev = $F[10]; + if ( $dme_elev !~ m/\d/ ) { + $dme_elev = ""; + } else { + $dme_elev += 0; + } + $dme_bias = $F[27]; + # print "$id DME $dme_lon $dme_lat $dme_elev $dme_bias\n"; + } elsif ( $type eq "G" ) { + # GlideSlope entry + $has_gs = 1; + $gs_lon = make_dcoord( $F[16] ); + $gs_lat = make_dcoord( $F[14] ); + $gs_elev = $F[10]; + if ( $gs_elev !~ m/\d/ ) { + $gs_elev = ""; + } else { + $gs_elev += 0; + } + $gs_angle = $F[7]; + # print "$id GS $gs_lon $gs_lat $gs_elev $gs_angle\n"; + } elsif ( $type eq "Z" ) { + # Localizer entry + $has_loc = 1; + $loc_lon = make_dcoord( $F[16] ); + $loc_lat = make_dcoord( $F[14] ); + $loc_elev = $F[10]; + if ( $loc_elev !~ m/\d/ ) { + $loc_elev = ""; + } else { + $loc_elev += 0; + } + ($loc_freq) = $F[5] =~ m/(\d\d\d\d\d\d)/; + $loc_freq /= 1000.0; + my( $magvar ) = make_dmagvar( $F[22] ); + # print "mag var = $F[22] (" . $magvar . ")\n"; + $loc_hdg = $F[24] + make_dmagvar( $F[22] ); + $loc_width = $F[25]; + $loc_id = $F[18]; + if ( length( $loc_id ) >= 4 ) { + $loc_id =~ s/^I//; + } + # print "$id LOC $loc_lon $loc_lat $loc_elev $loc_freq $loc_hdg $loc_width\n"; + } elsif ( $type eq "I" ) { + # Inner marker entry + $has_im = 1; + $im_lon = make_dcoord( $F[16] ); + $im_lat = make_dcoord( $F[14] ); + # print "$id IM $im_lon $im_lat\n"; + } elsif ( $type eq "M" ) { + # Middle marker entry + $has_mm = 1; + $mm_lon = make_dcoord( $F[16] ); + $mm_lat = make_dcoord( $F[14] ); + # print "$id MM $mm_lon $mm_lat\n"; + } elsif ( $type eq "O" ) { + # Outer marker entry + $has_om = 1; + $om_lon = make_dcoord( $F[16] ); + $om_lat = make_dcoord( $F[14] ); + # print "$id OM $om_lon $om_lat\n"; + } + + $last_id = $id; + $last_rwy = $rwy; + # printf("%-5s %10.6f %11.6f\n", $F[0], $F[14], $F[16]); + } + + if ( ! $has_gs ) { + ( $gs_elev, $gs_angle, $gs_lat, $gs_lon ) = ( 0, 0, 0, 0 ); + } + if ( ! $has_dme ) { + ( $dme_lat, $dme_lon ) = ( 0, 0 ); + } + if ( ! $has_om ) { + ( $om_lat, $om_lon ) = ( 0, 0 ); + } + if ( ! $has_mm ) { + ( $mm_lat, $mm_lon ) = ( 0, 0 ); + } + if ( ! $has_im ) { + ( $im_lat, $im_lon ) = ( 0, 0 ); + } + if ( $ILS{$CODES{$last_id} . $last_rwy} eq "" ) { + print "DAFIFT adding (last): $CODES{$last_id} - $last_rwy\n"; + &safe_add_record( $CODES{$last_id}, $last_rwy, "ILS", $loc_freq, + $loc_id, $loc_hdg, $loc_lat, $loc_lon, + $gs_elev, $gs_angle, $gs_lat, $gs_lon, + $dme_lat, $dme_lon, $om_lat, $om_lon, $mm_lat, + $mm_lon, $im_lat, $im_lon ); + } +} + + +######################################################################## +# Process FAA data +######################################################################## + +sub load_faa() { + my( $file ) = shift; + + open( FAA_ILS, "<$file" ) || die "Cannot open FAA data: $file\n"; + + ; # skip header line + + while ( ) { + chomp; + + my ( $rec_type, $faa_id, $rwy, $type, $faa_date, + $faa_apt_name, $faa_city, $faa_st, $faa_state, + $faa_region, $id, $faa_len, $faa_wid, $faa_cat, + $faa_owner, $faa_operator, $faa_bearing, $faa_magvar, + $loc_type, $loc_id, $loc_freq, $faa_loc_latd, + $faa_loc_lats, $faa_loc_lond, $faa_loc_lons, $loc_width, + $faa_stop_dist, $faa_app_dist, $faa_gs_type, $gs_angle, + $faa_gs_freq, $faa_gs_latd, $faa_gs_lats, $faa_gs_lond, + $faa_gs_lons, $faa_gs_dist, $gs_elev, $faa_im_type, + $faa_im_latd, $faa_im_lats, $faa_im_lond, $faa_im_lons, + $faa_im_dist, $faa_mm_type, $faa_mm_id, $faa_mm_name, + $faa_mm_freq, $faa_mm_latd, $faa_mm_lats, $faa_mm_lond, + $faa_mm_lons, $faa_mm_dist, $faa_om_type, $faa_om_id, + $faa_om_name, $faa_om_freq, $faa_om_latd, $faa_om_lats, + $faa_om_lond, $faa_om_lons, $faa_om_dist, + $faa_om_backcourse, $faa_dme_channel, $faa_dme_latd, + $faa_dme_lats, $faa_dme_lond, $faa_dme_lons, $faa_dme_app_dist, + $faa_dme_stop_dist, $blank) + = $_ =~ + m/^(.{4})(.{11})(.{3})(.{10})(.{10})(.{42})(.{26})(.{2})(.{20})(.{3})(.{4})(.{5})(.{4})(.{9})(.{50})(.{50})(.{3})(.{3})(.{15})(.{5})(.{6})(.{14})(.{11})(.{14})(.{11})(.{5})(.{5})(.{6})(.{15})(.{4})(.{6})(.{14})(.{11})(.{14})(.{11})(.{6})(.{7})(.{15})(.{14})(.{11})(.{14})(.{11})(.{6})(.{15})(.{2})(.{5})(.{3})(.{14})(.{11})(.{14})(.{11})(.{6})(.{15})(.{2})(.{5})(.{3})(.{14})(.{11})(.{14})(.{11})(.{6})(.{9})(.{4})(.{14})(.{11})(.{14})(.{11})(.{6})(.{5})(.{34})/; + + $id = &strip_ws( $id ); + $rwy = &strip_ws( $rwy ); + $rwy =~ s/\/$//; + $rwy =~ s/\/$//; + $loc_id =~ s/^I-//; + my( $loc_hdg ) = $faa_bearing + make_dmagvar($faa_magvar); + my( $loc_lat ) = make_dcoord($faa_loc_lats) / 3600.0; + my( $loc_lon ) = make_dcoord($faa_loc_lons) / 3600.0; + # print "$loc_lon $loc_lat $faa_loc_lons $faa_loc_lats\n"; + my( $gs_lat ) = make_dcoord($faa_gs_lats) / 3600.0; + my( $gs_lon ) = make_dcoord($faa_gs_lons) / 3600.0; + my( $im_lat ) = make_dcoord($faa_im_lats) / 3600.0; + my( $im_lon ) = make_dcoord($faa_im_lons) / 3600.0; + my( $mm_lat ) = make_dcoord($faa_mm_lats) / 3600.0; + my( $mm_lon ) = make_dcoord($faa_mm_lons) / 3600.0; + my( $om_lat ) = make_dcoord($faa_om_lats) / 3600.0; + my( $om_lon ) = make_dcoord($faa_om_lons) / 3600.0; + my( $dme_lat ) = make_dcoord($faa_dme_lats) / 3600.0; + my( $dme_lon ) = make_dcoord($faa_dme_lons) / 3600.0; + + # my( $key ); + # print "$id - $rwy\n"; + # $key = $id . $rwy; + # print "-> $key -> $ILS{$key}\n"; + # $key = "K" . $id . $rwy; + # print "-> $key -> $ILS{$key}\n"; + + if ( $rec_type eq "ILS1" ) { + if ( length( $id ) < 4 ) { + if ( $CodesByICAO{"K" . $id} ) { + $id = "K" . $id; + } + } + if ( $ILS{$id . $rwy} ne "" ) { + print "FAA updating: $id - $rwy $type\n"; + &update_type( $id, $rwy, $type ); + } else { + print "FAA adding: $id - $rwy\n"; + &safe_add_record( $id, $rwy, $type, $loc_freq, $loc_id, + $loc_hdg, $loc_lat, $loc_lon, $gs_elev, + $gs_angle, $gs_lat, $gs_lon, $dme_lat, + $dme_lon, $om_lat, $om_lon, $mm_lat, + $mm_lon, $im_lat, $im_lon ); + } + } + } +} + + +######################################################################## +# Process FlightGear ILS data +######################################################################## + +sub load_fgfs() { + my( $ils_file ) = shift; + + open( FGILS, "zcat $ils_file|" ) || die "Cannot open FGFS: $ils_file\n"; + + ; # skip header line + + while ( ) { + chomp; + if ( ! m/\[End\]/ && length($_) > 1 ) { + # print "$_\n"; + my( $type_code, $type_name, $icao, $rwy, $loc_freq, $loc_id, + $loc_hdg, $loc_lat, $loc_lon, $gs_elev, $gs_angle, $gs_lat, + $gs_lon, $dme_lat, $dme_lon, $om_lat, $om_lon, $mm_lat, $mm_lon, + $im_lat, $im_lon ) = split(/\s+/); + my( $code ) = $icao; + $code =~ s/^K//; + if ( $ILS{$icao . $rwy} ne "" ) { + print "FGFS: Skipping $icao - $rwy - already exists\n"; + # skip approaches already in FAA or DAFIFT data + } elsif ( length( $icao ) < 4 || $icao =~ m/^K/ ) { + print "FGFS: Skipping $icao - $rwy - USA\n"; + # skip USA approaches not found in FAA or DAFIFT data + } else { + print "FGFS adding: $icao $rwy\n"; + &safe_add_record( $icao, $rwy, $type_name, $loc_freq, $loc_id, + $loc_hdg, $loc_lat, $loc_lon, $gs_elev, + $gs_angle, $gs_lat, $gs_lon, $dme_lat, + $dme_lon, $om_lat, $om_lon, $mm_lat, + $mm_lon, $im_lat, $im_lon ); + } + } else { + print "FGFS discarding: $_\n"; + } + } +} + + +######################################################################## +# Load the x-plane/flightgear airport/runway information +######################################################################## + +sub load_fgfs_airports() { + my( $infile ) = shift; + + my( $id ); + + open( IN, "zcat $infile|" ) || die "Cannot open: $infile\n"; + + ; # skip header line + + while ( ) { + chomp; + my(@F) = split(/\s+/); + if ( $F[0] eq "A" ) { + $id = $F[1]; + } elsif ( $F[0] eq "R" ) { + my($recip) = &rwy_recip( $F[1] ); + $RUNWAYS{ $id . $F[1] } = "FOR $_"; + $RUNWAYS{ $id . $recip } = "REV $_"; + } + } +} + + +######################################################################## +# Run through the final ils list and adjust the heading to match the runway +######################################################################## + +sub fix_localizer() { + my( $key ); + foreach $key ( sort (keys %ILS) ) { + my(@F) = split( /\s+/, $ILS{$key} ); + print "FIXING: $key $F[2] $F[3] $F[4]\n"; + if ( $RUNWAYS{$key} ) { + my(@G) = split( /\s+/, $RUNWAYS{$key} ); + print " LocPos = $F[7] $F[8]\n"; + print " RwyInfo = $G[0] $G[3] $G[4] $G[5] $G[6] $G[7]\n"; + my($cmd) + = "$calc_loc $F[7] $F[8] $G[0] $G[3] $G[4] $G[5] $G[6] $G[7]"; + open( CALC, "$cmd |" ) || die "Cannot open $calc_loc\n"; + my( $j, $lat, $lon ); + while ( ) { + if ( m/New localizer/ ) { + ($j, $j, $j, $lat, $lon) = split( /\s+/ ); + } + } + if ( $G[0] eq "REV" ) { + $G[5] += 180.0; + if ( $G[5] >= 360.0 ) { + $G[5] -= 360.0; + } + } + &update_loc( $F[2], $F[3], $lat, $lon, $G[5]); + } + } +} + + +######################################################################## +# Write out the accumulated combined result +######################################################################## + +sub write_result() { + my( $outfile ) = shift; + + open( OUT, ">$outfile" ) || die "Cannot write to: $outfile\n"; + + # dump out the final results + print OUT "// FlightGear ILS data, generated from DAFIFT ARPT/ILS.TXT and FAA data\n"; + + my( $key ); + foreach $key ( sort (keys %ILS) ) { + print OUT "$ILS{$key}\n"; + } + print OUT "[End]\n"; +} + + +######################################################################## +# Utility functions +######################################################################## + + +# add a record to the master list if it doesn't already exist + +sub safe_add_record() { + my( $apt_id ) = shift; + my( $rwy ) = shift; + my( $type ) = shift; + my( $loc_freq ) = shift; + my( $loc_id ) = shift; + my( $loc_hdg ) = shift; + my( $loc_lat ) = shift; + my( $loc_lon ) = shift; + my( $gs_elev ) = shift; + my( $gs_angle ) = shift; + my( $gs_lat ) = shift; + my( $gs_lon ) = shift; + my( $dme_lat ) = shift; + my( $dme_lon ) = shift; + my( $om_lat ) = shift; + my( $om_lon ) = shift; + my( $mm_lat ) = shift; + my( $mm_lon ) = shift; + my( $im_lat ) = shift; + my( $im_lon ) = shift; + + if ( $ILS{$apt_id . $rwy} eq "" ) { + # print "Safe adding (common): $apt_id - $rwy\n"; + &update_record( $apt_id, $rwy, $type, $loc_freq, $loc_id, + $loc_hdg, $loc_lat, $loc_lon, $gs_elev, + $gs_angle, $gs_lat, $gs_lon, $dme_lat, + $dme_lon, $om_lat, $om_lon, $mm_lat, + $mm_lon, $im_lat, $im_lon ); + } +} + + +# replace a record in the master list (or add it if it doesn't exist) + +sub update_record() { + my( $apt_id ) = shift; + my( $rwy ) = shift; + my( $type ) = shift; + my( $loc_freq ) = shift; + my( $loc_id ) = shift; + my( $loc_hdg ) = shift; + my( $loc_lat ) = shift; + my( $loc_lon ) = shift; + my( $gs_elev ) = shift; + my( $gs_angle ) = shift; + my( $gs_lat ) = shift; + my( $gs_lon ) = shift; + my( $dme_lat ) = shift; + my( $dme_lon ) = shift; + my( $om_lat ) = shift; + my( $om_lon ) = shift; + my( $mm_lat ) = shift; + my( $mm_lon ) = shift; + my( $im_lat ) = shift; + my( $im_lon ) = shift; + + my( $record ); + + # remap $type as needed + $type = &strip_ws( $type ); + if ( $type eq "LOCALIZER" ) { + $type = "LOC"; + } elsif ( $type eq "ILS/DME" ) { + $type = "ILS"; + } elsif ( $type eq "SDF/DME" ) { + $type = "SDF"; + } elsif ( $type eq "LOC/DME" ) { + $type = "ILS"; + } elsif ( $type eq "LOC/GS" ) { + $type = "LOC"; + } elsif ( $type eq "LDA/DME" ) { + $type = "LDA"; + } + + $record = sprintf( "%1s %-5s %-4s %-3s %06.2f %-4s %06.2f %10.6f %11.6f ", + substr( $type, 0, 1 ), $type, $apt_id, $rwy, + $loc_freq, $loc_id, $loc_hdg, $loc_lat, $loc_lon ); + $record .= sprintf( "%5d %5.2f %10.6f %11.6f ", + $gs_elev, $gs_angle, $gs_lat, $gs_lon ); + $record .= sprintf( "%10.6f %11.6f ", $dme_lat, $dme_lon ); + $record .= sprintf( "%10.6f %11.6f ", $om_lat, $om_lon ); + $record .= sprintf( "%10.6f %11.6f ", $mm_lat, $mm_lon ); + $record .= sprintf( "%10.6f %11.6f ", $im_lat, $im_lon ); + + # print "Updating (common): $apt_id - $rwy\n"; + $ILS{$apt_id . $rwy} = $record; + $AIRPORTS{$apt_id} = 1; +} + + +# update the $type of the record +sub update_type() { + my( $apt_id ) = shift; + my( $rwy ) = shift; + my( $new_type ) = shift; + + my( $record ); + + if ( $ILS{$apt_id . $rwy} ne "" ) { + my( $type_code, $type_name, $apt_id, $rwy, $loc_freq, $loc_id, + $loc_hdg, $loc_lat, $loc_lon, $gs_elev, $gs_angle, $gs_lat, + $gs_lon, $dme_lat, $dme_lon, $om_lat, $om_lon, $mm_lat, $mm_lon, + $im_lat, $im_lon ) = split( /\s+/, $ILS{$apt_id . $rwy} ); + # print "Updating type: $apt_id $rwy: $type_name -> $new_type\n"; + $type_name = $new_type; + &update_record( $apt_id, $rwy, $type_name, $loc_freq, $loc_id, + $loc_hdg, $loc_lat, $loc_lon, $gs_elev, + $gs_angle, $gs_lat, $gs_lon, $dme_lat, + $dme_lon, $om_lat, $om_lon, $mm_lat, + $mm_lon, $im_lat, $im_lon ); + } else { + die "Error, trying to update $apt_id - $rwy which doesn't exist\n"; + } +} + + +# update the localizer position of the record +sub update_loc() { + my( $apt_id ) = shift; + my( $rwy ) = shift; + my( $nloc_lat ) = shift; + my( $nloc_lon ) = shift; + my( $nloc_hdg ) = shift; + + my( $record ); + + if ( $ILS{$apt_id . $rwy} ne "" ) { + my( $type_code, $type_name, $apt_id, $rwy, $loc_freq, $loc_id, + $loc_hdg, $loc_lat, $loc_lon, $gs_elev, $gs_angle, $gs_lat, + $gs_lon, $dme_lat, $dme_lon, $om_lat, $om_lon, $mm_lat, $mm_lon, + $im_lat, $im_lon ) = split( /\s+/, $ILS{$apt_id . $rwy} ); + print " Old pos = " . $loc_lat . "," . $loc_lon . " " . $loc_hdg . "\n"; + print " New pos = " . $nloc_lat . "," . $nloc_lon . " " . $nloc_hdg . "\n"; + $loc_lat = $nloc_lat; + $loc_lon = $nloc_lon; + $loc_hdg = $nloc_hdg; + &update_record( $apt_id, $rwy, $type_name, $loc_freq, $loc_id, + $loc_hdg, $loc_lat, $loc_lon, $gs_elev, + $gs_angle, $gs_lat, $gs_lon, $dme_lat, + $dme_lon, $om_lat, $om_lon, $mm_lat, + $mm_lon, $im_lat, $im_lon ); + } else { + die "Error, trying to update $apt_id - $rwy which doesn't exist\n"; + } +} + + +# convert a lon/lat coordinate in various formats to signed decimal + +sub make_dcoord() { + my($coord) = shift; + my( $dir, $deg, $min, $sec ); + my( $value ) = 0.0; + + $coord = &strip_ws( $coord ); + + if ( $coord =~ m/^[WE]/ ) { + ( $dir, $deg, $min, $sec ) + = $coord =~ m/^([EW])(\d\d\d)(\d\d)(\d\d\d\d)/; + $value = $deg + $min/60.0 + ($sec/100)/3600.0; + if ( $dir eq "W" ) { + $value = -$value; + } + } elsif ( $coord =~ m/^[NS]/ ) { + ( $dir, $deg, $min, $sec ) + = $coord =~ m/^([NS])(\d\d)(\d\d)(\d\d\d\d)/; + $value = $deg + $min/60.0 + ($sec/100)/3600.0; + if ( $dir eq "S" ) { + $value = -$value; + } + } elsif ( $coord =~ m/[EW]$/ ) { + ($value, $dir) = $coord =~ m/([\d\s\.]+)([EW])/; + if ( $dir eq "W" ) { + $value = -$value; + } + } elsif ( $coord =~ m/[NS]$/ ) { + ($value, $dir) = $coord =~ m/([\d\s\.]+)([NS])/; + if ( $dir eq "S" ) { + $value = -$value; + } + } + # print "$dir $deg:$min:$sec = $value\n"; + return $value; +} + +# convert a magnetic variation in various formats to signed decimal + +sub make_dmagvar() { + my( $coord ) = shift; + my( $value ); + + if ( $coord =~ m/^[EW]/ ) { + my( $dir, $deg, $min, $date ) + = $coord =~ m/^([EW])(\d\d\d)(\d\d\d) (\d\d\d\d)/; + $value = $deg + ($min/10)/60.0; + if ( $dir eq "W" ) { + $value = -$value; + } + } elsif ( $coord =~ m/[EW]$/ ) { + my( $deg, $dir ) + = $coord =~ m/^(\d\d)([EW])/; + $value = $deg; + if ( $dir eq "W" ) { + $value = -$value; + } + } + # print "$dir $deg:$min = $value\n"; + + return $value; +} + +# strip white space off front and back of string + +sub strip_ws() { + my( $string ) = shift; + $string =~ s/^\s+//; + $string =~ s/\s+$//; + return $string; +} + + +# return the reciprical runway number +sub rwy_recip() { + my( $input ) = shift; + + my($num, $letter); + + if ( length($input) == 3 ) { + ($num, $letter) = $input =~ m/(\d\d)(.)/; + } elsif ( length($input) == 2 ) { + $num = $input; + $letter = ""; + } else { + $num = ""; + $letter = $input; + } + # print "RWY: $num - $letter <==> "; + + if ( $num ne "" ) { + $num += 18; + if ( $num > 35 ) { + $num -= 36; + } + if ( $num < 10 ) { + $num = "0$num"; + } + } + + if ( $letter eq "R" ) { + $letter = "L"; + } elsif ( $letter eq "L" ) { + $letter = "R"; + } elsif ( $letter eq "C" ) { + $letter = "C"; + } elsif ( $letter eq "N" ) { + $letter = "S"; + } elsif ( $letter eq "S" ) { + $letter = "N"; + } elsif ( $letter eq "E" ) { + $letter = "W"; + } elsif ( $letter eq "W" ) { + $letter = "E"; + } + # print "$num - $letter\n"; + + return "$num$letter"; +} -- 2.39.5