Add files via upload

bengida1989 · Dec 4, 2020 · ae2a268 · ae2a268
1 parent 4e8a7c4
commit ae2a268
Showing 1 changed file with 17 additions and 19 deletions.
diff --git a/vorticity_threshold.m b/vorticity_threshold.m
@@ -19,18 +19,17 @@
  vort = reshape(tmp, nRows, nColumns);
 
  elseif MASK(1) == 1 && threshold>0 % Regular masking
- S = swirl;
- for i=1:nRows
- for j=1:nColumns
- if abs(S(i,j))<=threshold*max(abs(S(:)))
- vorticity(i,j)=0;
- end;
- end;
- end;
+ P = swirl;
+ P_max = max(max(abs(P)));
+ P_thres = threshold*P_max;
+ vorticity(abs(P)<=P_thres) = 0;
  vort = vorticity;
+
  else
  vort = vorticity;
- end;
+
+ end
+
 
  case 'Vorticity'
  if MASK(1) == 0 && threshold>0% Gaussian masking
@@ -48,16 +47,15 @@
  vort = reshape(tmp, nRows, nColumns);
 
  elseif MASK(1) == 1 && threshold>0% Regular masking
- V = vorticity;
- for i=1:nRows
- for j=1:nColumns
- if abs(V(i,j))<=threshold*max(abs(V(:)))
- vorticity(i,j)=0;
- end;
- end;
- end;
+ P = vorticity;
+ P_max = max(max(abs(P)));
+ P_thres = threshold*P_max;
+ vorticity(abs(P)<=P_thres) = 0;
  vort = vorticity;
+
  else
  vort = vorticity;
- end;
-end;
+
+ end
+
+end