2D calculation of the baroclinic vorticity in compute_baroclinic_vort.cpp was...

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

@@ -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);