Changeset 21

Timestamp:

11/28/07 16:08:51 (16 years ago)

Author:

kolasinski

Message:

Bugfix in trends : corrections of the errors of the standard deviation of 1D time-series

File:

• trunk/procs/trends.pro (modified) (3 diffs)

trunk/procs/trends.pro

@@ -8 +8 @@
 ; type = 't', 'xt', 'yt', 'zt' or 'xyt'
 
-   print, '    Enter trends.pro'
+   print, 'ENTER trends...'
 
    trend_type = strtrim(string(trend_int), 2)
@@ -144 +144 @@
 
             ; Compute std error-bar using running bootstrap window (width in plt_def)
-            temperr = stat_error(fld, boot_win(0))
-            err_std_1 = sqrt( total((temperr-sqrt((moment(fld))[1]) )^2)/n_elements(temperr) )
-            temperr = stat_error(fld, boot_win(1))
-            err_std_2 = sqrt( total((temperr-sqrt((moment(fld))[1]) )^2)/n_elements(temperr) )
-            temperr = stat_error(fld, boot_win(2))
-            err_std_3 = sqrt( total((temperr-sqrt((moment(fld))[1]) )^2)/n_elements(temperr) )
+            IF jpt GT boot_win(0) THEN BEGIN 
+               temperr = stat_error(fld, boot_win(0))
+               err_std_1 = sqrt( total((temperr-sqrt((moment(fld))[1]) )^2)/n_elements(temperr) )
+            ENDIF ELSE err_std_1 = 0
+            IF jpt GT boot_win(1) THEN BEGIN 
+               temperr = stat_error(fld, boot_win(1))
+               err_std_2 = sqrt( total((temperr-sqrt((moment(fld))[1]) )^2)/n_elements(temperr) )
+            ENDIF  ELSE err_std_2 = 0
+            IF jpt GT boot_win(2) THEN BEGIN 
+               temperr = stat_error(fld, boot_win(2))
+               err_std_3 = sqrt( total((temperr-sqrt((moment(fld))[1]) )^2)/n_elements(temperr) )
+            ENDIF  ELSE err_std_3 = 0
 
             print, '     1D times-serie anomaly info : var/stddev +- error (24,36,48)', (moment(fld))[1], sqrt((moment(fld))[1]), err_std_1, err_std_2, err_std_3 
@@ -368 +374 @@
    ENDIF 
 
+   IF debug_w THEN print, 'EXIT trends...'
+
    return, fld
-   
-   IF debug_w THEN print, '   Vargrid : ',  vargrid
-   IF debug_w THEN print, '   Leaving trends.pro'
 
 END