diff --git a/src/Science/compute_baroclinic_vort.cpp b/src/Science/compute_baroclinic_vort.cpp
index d0bf1d73d8ed0b3a4ddbabb87b3c773db59395c3..a97aa642241b1a2a5ca42c865953c23ba5e34a53 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);