use lib "../lib"; 

use Physics::Water::SoundSpeed;

use Data::Dumper; 

my $obj = new Physics::Water::SoundSpeed();

my $ss = $obj->sound_speed_t(5);
print "Sound Speed : $ss\n"; 

my $ss = $obj->sound_speed( 10, .902 );
print "Sound Speed : $ss\n"; 

my $ss = $obj->sound_speed( [10,20,30] );
print Dumper $ss; 

my $ss = $obj->sound_speed( [10,20,30], .101 );
print Dumper $ss; 

my $ss = $obj->sound_speed( [10,20,30], [.101, .201, .302 ] );
print Dumper $ss; 


my $pr = $obj->d2p_sea( 0 );
print "Pressure at surface $pr\n"; 

my $dp = [0,50,100,150,200,250,300,350,400,450,500,600,700,800,900,1000,2000,3000,4000,5000,6000];
my $pr = $obj->d2p_sea( $dp );
print Dumper $pr; 

my $ss = $obj->sound_speed( 35, $pr );

for (my $i=0; $i<=$#$dp; $i++)
{
   print "$pr->[$i] $dp->[$i] $ss->[$i]\n";
}

my $ss = $obj->sound_speed_sea_tps( 5, 1,35);
print "Speed Sound in Sea at 5, 1, 35 -> $ss\n"; 

my $ss = $obj->sound_speed_sea_tps( 4, 7, 55);
print "Speed Sound in Sea at 4, 7, 55 -> $ss\n"; 

my $obj_us = new Physics::Water::SoundSpeed('units'=>'US');

my $ss_us = sprintf("%5.3f", $obj_us->sound_speed_t(41));
my $ss_si = sprintf("%5.3f", $obj->sound_speed_t(5) / 0.3048) ;
print "Sound Speed : $ss_us / $ss_si\n";