From f2834d835beab1e8987c91f5f21ea8bdc5f2b5c2 Mon Sep 17 00:00:00 2001
From: Donovan Allum <dallum@uwaterloo.ca>
Date: Thu, 26 Jan 2023 14:10:03 -0500
Subject: [PATCH] Updating the baroclinic vorticity calculation in
derivatives.cpp to include a factor of 1/rho_0
---
src/cases/derivatives/derivatives.cpp | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/cases/derivatives/derivatives.cpp b/src/cases/derivatives/derivatives.cpp
index 2165711..ab383f0 100644
--- a/src/cases/derivatives/derivatives.cpp
+++ b/src/cases/derivatives/derivatives.cpp
@@ -221,19 +221,19 @@ class userControl : public BaseCase {
// Store Temperature in T, it is free
init_tracer_restart("t",u);
compute_baroclinic_vort(deriv_var, *temp4, u, gradient_op, grid_type, v_exist);
- deriv_var=deriv_var*g;
+ deriv_var=deriv_var*g/rho_0;
write_array(deriv_var,"bar",plotnum);
if (master()) fprintf(stdout,"Completed the write for bar.%d\n",plotnum);
if ( v_exist ) {
compute_baroclinic_vort_x(deriv_var, u, gradient_op, grid_type);
- deriv_var=deriv_var*g;
+ deriv_var=deriv_var*g/rho_0;
write_array(deriv_var,"barx",plotnum);
if (master()) fprintf(stdout,"Completed the write for barx.%d\n",plotnum);
}
compute_baroclinic_vort_y(deriv_var, u, gradient_op, grid_type);
- deriv_var=deriv_var*g;
+ deriv_var=deriv_var*g/rho_0;
write_array(deriv_var,"bary",plotnum);
if (master()) fprintf(stdout,"Completed the write for bary.%d\n",plotnum);
--
GitLab