#! /usr/bin/perl -w

#load some functions
use File::Copy qw(copy);
use File::Copy qw(move);
use Cwd;

print(">>>>>>>>>>>>Make sure you have netcdf libraries>>>>>>>>>>>>>\n");
print(">>>>>>>>>>>>Make sure the following excutables work on your machine:>>>>>>>>>>>>>\n");
print(">>>>>>>>>>>>    ./Nudge_gr3/gen_nudge2>>>>>>>>>>>>>\n");
print(">>>>>>>>>>>>    ./gen_nudge_from_hycom>>>>>>>>>>>>>\n");
print(">>>>>>>>>>>>Recompile them if necessary;>>>>>>>>>>>>>\n");
print(">>>>>>>>>>>>see an example compiling cmd in each source code.>>>>>>>>>>>>>\n");
print("\n");
#{
#  local( $| ) = ( 1 );
#  print "Press <Enter> or <Return> to continue: \n";
#  my $resp = <STDIN>;
#}

#dirs
$script_dir="../Grid_manipulation/";
$hycom_dir="../HYCOM_2015/";

$thisdir=cwd();
chdir("..");
$rundir = cwd();
chdir($thisdir);
system("ln -sf $hycom_dir/*.nc .");

#UTM grid
system("ln -sf ../hgrid.gr3 .");
system("ln -sf ../hgrid.ll .");
system("ln -sf ../vgrid.in .");

chdir("Nudge_gr3");
system("ln -sf $rundir/hgrid.ll hgrid.gr3");
system("./gen_nudge2 < gen_nudge2.in");
unlink("$rundir/SAL_nudge.gr3");
copy("nudge.gr3","$rundir/SAL_nudge.gr3");
unlink("$rundir/TEM_nudge.gr3");
copy("nudge.gr3","$rundir/TEM_nudge.gr3");



chdir($thisdir);
copy("./Nudge_gr3/nudge.gr3","TEM_nudge.gr3");

#set include zone, which is a bit larger than the actual nudging zone
#no longer needed after the update on gen_nudge_from_hycom, which directly reads TEM_nudge.gr3 and expand the region based on neighboring nodes within an element
#system("$script_dir/auto_edit_region 0 include.reg hgrid.gr3 1 0");
#move("out.gr3","include.gr3");

#generate nudging files
system("./gen_nudge_from_hycom.exe");

unlink("../SAL_nu.nc");
unlink("../TEM_nu.nc");

copy("SAL_nu.nc","../SAL_nu.nc");
copy("TEM_nu.nc","../TEM_nu.nc");

print(">>>>>>>>>>>>Done.>>>>>>>>>>>>>\n")