use lib "../lib"; 

use Physics::Water::SoundSpeed;

use Data::Dumper; 

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

my $pr = [1,5,10,15,20,25];
my $tm = [1,5,10,15,20,25,30,35,40,45,50,55];

for (my $j=0; $j<=$#$pr; $j++)
{
   my $pri = $pr->[$j];
   
   my $ss = $obj->sound_speed( $tm, $pri );
   print "\npressure $pri\n\n";
   open(FID,">isob$pri.txt");
   for (my $i=0; $i<=$#$ss; $i++)
    {
        print FID "$tm->[$i] $ss->[$i]\n";
    }
    close FID;
}