Changeset 35

Timestamp:

01/16/08 01:30:08 (16 years ago)

Author:

ericg

Message:

Bug correction in c_correlate of fld and fld_prev if different sizes

File:

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

trunk/procs/trends.pro

@@ -8 +8 @@
 ; type = 't', 'xt', 'yt', 'zt' or 'xyt'
 
-   print, 'ENTER trends...'
+   print, '    Enter trends.pro'
 
    trend_type = strtrim(string(trend_int), 2)
@@ -162 +162 @@
 ; compute time lag correlation
             IF ioverchk EQ 2 AND c_correl EQ 1 THEN BEGIN
-               lag_array = findgen(2*lag_correl+1)-lag_correl
-               lag_correlation = C_CORRELATE(fld, fld_prev_t, lag_array)
-               print, '     1D time series lag correlation = ', lag_correlation
-               indx = where (lag_correlation EQ max(lag_correlation))
-               IF indx[0] NE -1 THEN BEGIN 
-                  print, '        Lag of max correlation =', lag_array(indx),lag_correlation(indx) 
-                  print, '        Lag correlation limit [-N,N] =', lag_correl 
+               IF (size(fld))(1) EQ (size(fld_prev_t))(1) THEN BEGIN 
+                  lag_array = findgen(2*lag_correl+1)-lag_correl
+                  lag_correlation = C_CORRELATE(fld, fld_prev_t, lag_array)
+                  print, '     1D time series lag correlation = ', lag_correlation
+                  indx = where (lag_correlation EQ max(lag_correlation))
+                  IF indx[0] NE -1 THEN BEGIN 
+                     print, '        Lag of max correlation =', lag_array(indx),lag_correlation(indx) 
+                     print, '        Lag correlation limit [-N,N] =', lag_correl 
+                  ENDIF 
                ENDIF 
             ENDIF 
@@ -374 +376 @@
    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