summaryrefslogtreecommitdiff
path: root/797/CH5/EX5.9.s/5_09_solution.sce
diff options
context:
space:
mode:
Diffstat (limited to '797/CH5/EX5.9.s/5_09_solution.sce')
-rw-r--r--797/CH5/EX5.9.s/5_09_solution.sce17
1 files changed, 17 insertions, 0 deletions
diff --git a/797/CH5/EX5.9.s/5_09_solution.sce b/797/CH5/EX5.9.s/5_09_solution.sce
new file mode 100644
index 000000000..6432c6b74
--- /dev/null
+++ b/797/CH5/EX5.9.s/5_09_solution.sce
@@ -0,0 +1,17 @@
+//Solution 5-09
+WD=get_absolute_file_path('5_09_solution.sce');
+datafile=WD+filesep()+'5_09_example.sci';
+clc;
+exec(datafile)
+V_A = V_A * 1000 / 3600; //conversion from [kmph] to [m/s]
+P_atmair = P_atmair / 100; //conversion from [cm of Hg] to [m og Hg]
+P_air = P_air / 100; //conversion from [cm of Hg] to [m of Hg]
+//(a)
+h_1 = rho_Hg / rho_sw * (P_atmair - P_air); //from pressure realtion P=rho*g*h
+printf("Ocean swell at point 3 is %1.2f m", h_1);
+//(b)
+h_air = (V_A)^2 / (2 * g); //Bernoulli equation application between A and B
+rho_air = P_air / P_atmair * rho_atmair //from ideal gas equation
+h_dynamic = rho_air / rho_sw * h_air; //surge of point 2 from point 3
+h_2 = h_1 + h_dynamic; //total surge at point 2
+printf("\nOcean swell at point 2 is %1.2f m", h_2);