diff --git a/src/solver_coriolis.cpp b/src/solver_coriolis.cpp index 0ab125f7..4a17897e 100644 --- a/src/solver_coriolis.cpp +++ b/src/solver_coriolis.cpp @@ -14,10 +14,10 @@ std::vector coriolis(std::vector velocity, precision_t rotation_rate, arma_cube lat_scgc) { std::vector coriolis_vec(3); - coriolis_vec[0] = - 2 * rotation_rate * velocity[1] % sin(lat_scgc) - - 2 * rotation_rate * velocity[2] % cos(lat_scgc); + coriolis_vec[0] = ( + 2 * rotation_rate * velocity[1] % sin(lat_scgc)); // - + //2 * rotation_rate * velocity[2] % cos(lat_scgc) ); coriolis_vec[1] = -2 * rotation_rate * velocity[0] % sin(lat_scgc); - coriolis_vec[2] = 2 * rotation_rate * cos(lat_scgc) % velocity[1]; + coriolis_vec[2] = 2 * rotation_rate * cos(lat_scgc) % velocity[1] * 0.0; return coriolis_vec; }