From b14a7b52adbfdb53112a028f56c9f6fc68850307 Mon Sep 17 00:00:00 2001
From: Donovan Allum <dallum@uwaterloo.ca>
Date: Thu, 26 Jan 2023 14:02:07 -0500
Subject: [PATCH] 2D calculation of the baroclinic vorticity in
compute_baroclinic_vort.cpp was incorrectly clearly variables from a previous
timestep, this has been fixed
---
src/Science/compute_baroclinic_vort.cpp | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/Science/compute_baroclinic_vort.cpp b/src/Science/compute_baroclinic_vort.cpp
index d0bf1d7..a97aa64 100644
--- a/src/Science/compute_baroclinic_vort.cpp
+++ b/src/Science/compute_baroclinic_vort.cpp
@@ -18,11 +18,11 @@ using namespace Transformer;
void compute_baroclinic_vort(TArrayn::DTArray & barovort, TArrayn::DTArray & temp,
TArrayn::DTArray & T, TArrayn::Grad * gradient_op, const string * grid_type, bool v_exist ) {
if ( v_exist ) {
- compute_baroclinic_vort_x(temp, T, gradient_op, grid_type); }
+ compute_baroclinic_vort_x(temp, T, gradient_op, grid_type);
+ barovort = temp*temp; }
else {
barovort = 0;
}
- barovort = temp*temp;
compute_baroclinic_vort_y(temp, T, gradient_op, grid_type);
barovort = barovort + temp*temp;
barovort = sqrt(barovort);
--
GitLab