tsGetPOTAndRLargest: debug the case of 1 input percentiles

menta78 · Jun 16, 2016 · 3e7a704 · 3e7a704
1 parent 1288864
commit 3e7a704
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/tsGetPOTAndRlargest.m b/tsGetPOTAndRlargest.m
@@ -28,7 +28,7 @@
 if length(pcts) == 1
     % testing at least 2 percentages, to be able to compute the error on
     % the percentage.
-    pcts = [pcts(1), pcts(1) - 3];
+    pcts = [pcts(1) - 3, pcts(1)];
 end
 
 numperyear=nan(length(pcts),1);
@@ -51,7 +51,7 @@
     nperYear=tsGetNumberPerYear(ms,locs);
     minnumperyear(ipp)=nanmin(nperYear);
 
-    if length(pks)/nyears<desiredEventsPerYear & nanmin(nperYear)<desiredEventsPerYear
+    if (ipp > 1) && (length(pks)/nyears<desiredEventsPerYear) && (nanmin(nperYear)<desiredEventsPerYear)
 
         break