#!/usr/bin/perl

my $sac = 1.0;

print "<table border=1>\n";
print "<tr>\n";
print "<th>Depth</th>\n";
print "<th>Rock Bottom</th>\n";
print "<th>Al6</th>\n";
print "<th>Al9</th>\n";
print "<th>Al13</th>\n";
print "<th>Al19</th>\n";
print "<th>Al30</th>\n";
print "<th>Al40</th>\n";
print "</tr>\n";

foreach $d (1..13) {
        my $depth = $d * 10;
        my $ata = $depth / 33 + 1.0;

        my $probgas = $sac * $ata * 1.0;

        my $asctime = $depth / 30;
        my $ascgas  = $sac * ( $ata + 1 ) / 2 * $asctime;

        my $stopgas = $sac * ( 20 / 33 + 1.0 ) * 3.0;

        my $gas = $probgas + $ascgas + $stopgas;

        my $al6 = $gas / 6.1 * 3000.0;
        my $al9 = $gas / 8.2 * 3000.0;
        my $al13 = $gas / 13.2 * 3000.0;
        my $al19 = $gas / 19.9 * 3000.0;
        my $al30 = $gas / 30.0 * 3000.0;
        my $al40 = $gas / 40.0 * 3000.0;

        $al6 = 500 if $al6 < 500;
        $al9 = 500 if $al9 < 500;
        $al13 = 500 if $al13 < 500;
        $al19 = 500 if $al19 < 500;
        $al30 = 500 if $al30 < 500;
        $al40 = 500 if $al40 < 500;

        $al6 = "NA" if $al6 > 3000;
        $al9 = "NA" if $al9 > 3000;
        $al13 = "NA" if $al13 > 3000;
        $al19 = "NA" if $al19 > 3000;
        $al30 = "NA" if $al30 > 3000;
        $al40 = "NA" if $al40 > 3000;

        print "<tr>";
        print "<td>$depth fsw</td>";
        printf "<td>%5.2f cu ft</td>", $gas;
        ( $al6 != "NA" ) ? printf "<td>%d psi</td>", $al6 : print "<td>N/A</td>";
        ( $al9 != "NA" ) ? printf "<td>%d psi</td>", $al9 : print "<td>N/A</td>";
        ( $al13 != "NA" ) ? printf "<td>%d psi</td>", $al13 : print "<td>N/A</td>";
        ( $al19 != "NA" ) ? printf "<td>%d psi</td>", $al19 : print "<td>N/A</td>";
        ( $al30 != "NA" ) ? printf "<td>%d psi</td>", $al30 : print "<td>N/A</td>";
        ( $al40 != "NA" ) ? printf "<td>%d psi</td>", $al40 : print "<td>N/A</td>";
        print "</tr>\n";
}

print "</table>\n";