#!/usr/bin/ruby
# File: lab14_v6t.rb
# Resources: Input14.txt
# => Opens file, prints variables to screen
# => Separates values from error
# => assigns 2D values to poi[i][j], pci[i][j], pof[i][j], pcj[i][j], ok[i][j], ol[i][j]
# => allows for more than 1 gas

lines = File.readlines("Input14.txt") #now an array of lines
name = lines[0]
experiment = lines[1]
partner = lines[2]
date = lines[3]
barometric_pressure = lines[4].split(',')
bap = barometric_pressure[0].to_f
ebap = barometric_pressure[1].to_f

dibutyl = lines[5].split(',')
#separate ddbp from eddbp
ddbp = dibutyl[0].to_f
puts "ddbp = #{ddbp}"
eddbp = dibutyl[1].to_f

mercury = lines[6].split(',')
#separate dhg from edhg
dhg = mercury[0].to_f
puts "dhg = #{dhg}"
edhg = mercury[1].to_f


ep = lines[7].to_f  #Uncertainty in U-Tube

ng = lines[8].to_i #number of gasses studied
linecount = 9
=begin
Q = lines[9]
m = lines[10].to_i
=end

=begin

puts "Name: #{name}"
puts "Experiment: #{experiment}"
puts "Partner: #{partner}"
puts "Date: #{date}"
puts "Barometric Pressure: #{bap}\tError: #{ebap}"

=end


$Q = Array.new(ng) #name of gasses
$sum1 = Array.new(ng){0.0}
$esum1 = Array.new(ng){0.0}
$MR = Array.new(ng){0.0}
$emR = Array.new(ng){0.0}


i=0

$Q[i] = lines[linecount]
linecount = linecount + 1
m = lines[linecount].to_i
linecount = linecount+1

$poi = Array.new(ng) {Array.new(m) {0.0}}
$pci = Array.new(ng) {Array.new(m) {0.0}}
$pof = Array.new(ng) {Array.new(m) {0.0}}
$pcf = Array.new(ng) {Array.new(m) {0.0}}
$ok = Array.new(ng) {Array.new(m) {0.0}}
$ol = Array.new(ng) {Array.new(m) {0.0}}

$Ci = Array.new(ng) {Array.new(m) {0.0}}
$ai = Array.new(ng) {Array.new(m) {0.0}}
$eai = Array.new(ng) {Array.new(m) {0.0}}
$p1 = Array.new(ng) {Array.new(m) {0.0}}
$ep1 = Array.new(ng) {Array.new(m) {0.0}}
$lnp1 = Array.new(ng) {Array.new(m) {0.0}}
$elnp1 = Array.new(ng) {Array.new(m) {0.0}}
$Cf = Array.new(ng) {Array.new(m) {0.0}}
$af = Array.new(ng) {Array.new(m) {0.0}}
$eaf = Array.new(ng) {Array.new(m) {0.0}}
$p3 = Array.new(ng) {Array.new(m) {0.0}}
$ep3 = Array.new(ng) {Array.new(m) {0.0}}
$lnp3 = Array.new(ng) {Array.new(m) {0.0}}
$elnp3 = Array.new(ng) {Array.new(m) {0.0}}
$n = Array.new(ng) {Array.new(m) {0.0}}
$en = Array.new(ng) {Array.new(m) {0.0}}
$d = Array.new(ng) {Array.new(m) {0.0}}
$ed = Array.new(ng) {Array.new(m) {0.0}}
$R = Array.new(ng) {Array.new(m) {0.0}}
$eR = Array.new(ng) {Array.new(m) {0.0}}


while i<ng
  j=0
  while j<m
    value = lines[linecount].split(',')
    $poi[i][j] = value[0].to_f 
    $pci[i][j] = value[1].to_f
    $pof[i][j] = value[2].to_f
    $pcf[i][j] = value[3].to_f

    $ok[i][j] = $poi[i][j] - $pci[i][j]
    $ol[i][j] = $pof[i][j] - $pcf[i][j]

    print "trial:#{j+1}\t#{$poi[i][j]} #{$pcf[i][j]} = "
    print "#{$ok[i][j]} mm DBP\t#{$pof[i][j]} #{$pcf[i][j]}=#{$ol[i][j]} mm DBP"
    puts

    j = j+1
    linecount = linecount + 1
  end
  puts
  i = i+1
end

lnp2 = Math.log(bap)
puts "lnp2 = #{lnp2}"

elnp2 = ebap/bap
puts "elnp2 = #{elnp2}"

puts "ddbp = #{ddbp}"
b = ddbp/dhg
puts "b = #{b}"
puts "dhg = #{dhg}"
eb = b*((eddbp/ddbp)**2 + (edhg/dhg)**2)**0.5
puts "eb = #{eb}"

ec = ep*2**0.5
puts "ec = #{ec}"

i = 0
while i<ng
  $sum1[i] = 0.0
  $esum1[i] = 0.0
  j=0
  while j<m
    $Ci[i][j] = $poi[i][j] - $pci[i][j]
    $ai[i][j] = b*$Ci[i][j]

    $eai[i][j] = $ai[i][j] * ((eb/b) **2 + (ec/$Ci[i][j]) **2) **0.5
    $p1[i][j] = bap + $ai[i][j]
    $ep1[i][j] = (ebap**2 + $eai[i][j]**2)**0.5
    $lnp1[i][j] = Math.log($p1[i][j])
    $elnp1[i][j] = $ep1[i][j]/$p1[i][j]
    $Cf[i][j] = $pof[i][j] - $pcf[i][j]
    $af[i][j] = b*$Cf[i][j]
    
    $eaf[i][j] = $af[i][j]* ((eb/b)**2 + (ec/$Cf[i][j])**2)**0.5
    $p3[i][j] = bap + $af[i][j]
    $ep3[i][j] = (ebap**2 + $eaf[i][j]**2)**0.5
  
    $lnp3[i][j] = Math.log($p3[i][j])
    $elnp3[i][j] = $ep3[i][j]/$p3[i][j]
    $n[i][j] = $lnp1[i][j] - lnp2
    $en[i][j] = ($elnp1[i][j]**2 + elnp2**2)**0.5
    $d[i][j] = $lnp1[i][j] - $lnp3[i][j]
    $ed[i][j] = ($elnp1[i][j]**2 + $elnp3[i][j]**2)**0.5
    $R[i][j] = $n[i][j]/$d[i][j]
    $eR[i][j] = $R[i][j]*(($en[i][j]/$n[i][j])**2 + ($ed[i][j]/$d[i][j])**2)**0.5
    $sum1[i] = $sum1[i] + $R[i][j]
    $esum1[i] = $esum1[i] + $eR[i][j]**2
    j = j+1
  end
  $MR[i] = $sum1[i]/m
  $emR[i] = $esum1[i]**0.5/m
  i = i+1
end
i=0
while i<ng
  puts "Gas: #{$Q[i]}"
  puts "Pressure, mm DBP,  Pressure, mmHG"
  puts "================,  =============="
  puts
  j=0
  while j<m
    puts "Trial = #{j+1}"
    print "Initial Equilibrium ="
    puts "#{$ok[i][j]}  = #{$p1[i][j]}"
    print "Final Equilibrium ="
    puts "#{$ol[i][j]}  = #{$p3[i][j]}"
    puts
    puts "Cp/Cv = #{$R[i][j]} +/- #{$eR[i][j]}"
    puts
    j = j+1
  end
  puts "Average Cp/Cv = #{$MR[i]} +/- #{$emR[i]}"
  puts
  puts
  i = i+1
end
=begin
#Check to see that they are there:

i= 0
puts "Gas: #{$Q[i]}"
puts "number of trials: #{m}"
while i<ng
  j=0
  while j<m
    puts "#{$poi[i][j]}\t#{$pci[i][j]}\t#{$pof[i][j]}\t#{$pcf[i][j]}"
    puts "#{$ok[i][j]}"
    j=j+1
  end
  i=i+1
end
=end

exit