Commit cf482c33 by maaf

A bug regarding critical frequency is fixed.

parent 8f88f6df
......@@ -137,10 +137,10 @@ FirstOrder::FirstOrder(
// Potentials output
oti_ = 0;
otj_ = 0;
oti_ = 16;
otj_ = 17;
otk_ = 0;
ots_ = 0;
ots_ = 1;
// ---------------------------------------------
//
......@@ -182,10 +182,15 @@ FirstOrder::FirstOrder(
// NOTE : I wrote the following functions based on beta=0.
// The beta is irrelevant to the value of tau
double depth_k = find_depth_critical(0, h_, mwl_, 0, inputData.U, inputData.g, 1e10); // The limit wave number
tau_ = U_/g_ * ( sqrt(g_ * depth_k * tanh(depth_k * h_) ) - depth_k * U_);
tau_ = 0.;
if(U_ != 0)
{
double depth_k = find_depth_critical(0, h_, mwl_, 0, inputData.U, inputData.g, 1e10); // The limit wave number
tau_ = U_/g_ * ( sqrt(g_ * depth_k * tanh(depth_k * h_) ) - depth_k * U_);
}
} // End of the constructor
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or sign in to comment