Changeset 112 for trunk

Timestamp:

06/10/08 03:01:01 (16 years ago)

Author:

ericg

Message:

Added corr_thres in plt_def to define correlation lower limit for plotting regression maps

Location:

trunk

Files:

• procs/com_eg.pro (modified) (1 diff)
• procs/init_post_it.pro (modified) (1 diff)
• procs/macros/make_linfit.pro (modified) (2 diffs)
• procs/macros/make_linfitdom.pro (modified) (2 diffs)
• usr/plt_def.pro (modified) (1 diff)

trunk/procs/com_eg.pro

@@ -13 +13 @@
 ;
 COMMON formats, look, free_1d_minmax, lat_axis, cont_fill, multi_win, atmos_msk, line_thick, line_style, line_color, title_type, fill_space, vector_sample, symbol_style, symbol_color, symbol_families, contour_options, mean_sc_only, nc_grids_list, cont_real, symbol_size, xchartxt, ychartxt, linfit_sep, linfit_map, bin_interval,  var3_ybinx
-COMMON switches, cmd1_back,  cmd2_back, file_name1, file_name2, ncdf_db1, ncdf_db2, grid1_full, grid2_full
+COMMON switches, cmd1_back,  cmd2_back, file_name1, file_name2, ncdf_db1, ncdf_db2, grid1_full, grid2_full, corr_thres
 
 ;

trunk/procs/init_post_it.pro

@@ -27 +27 @@
    bin_interval = 0
    var3_ybinx = ""
+   corr_thresh = 0.2

trunk/procs/macros/make_linfit.pro

@@ -126 +126 @@
   ; only take points where correlation larger than 0.2
 
-   idc = where(abs(pt_corr) LT 0.12)
+   idc = where(abs(pt_corr) LT corr_thres)
    pt_linfit(idc) = !values.f_nan
+
+   corr_txt = '[r>'+strtrim(string(corr_thres, format = '(f5.2)'), 2)+']'
 
    varname =  varname+' '+ntxt
@@ -150 +152 @@
    ENDCASE 
 
-   field.legend = lintxt+' for '+ntxt+' in ['+cmdm.date1+'-'+cmdm.spec+'] -'
+   field.legend = lintxt+' for '+ntxt+' in ['+cmdm.date1+'-'+cmdm.spec+'] '+corr_txt+' -'
 
 ; additional computations (pac_5 nino_3 nino_4 averages) *** requires whole domain ****

trunk/procs/macros/make_linfitdom.pro

@@ -110 +110 @@
   ; only take points where ABS(correlation) larger than a value
 
-   idc = where(abs(pt_corr) LT 0.1)
+   idc = where(abs(pt_corr) LT corr_thres)
    pt_linfit(idc) = !values.f_nan
+   corr_txt = '[r>'+strtrim(string(corr_thres, format = '(f5.2)'), 2)+']'
 
    varname =  varname+' '+ntxt
@@ -134 +135 @@
    ENDCASE 
 
-   field.legend = lintxt+' for '+ntxt+' in ['+cmdm.date1+'-'+cmdm.spec+'] -'
+   field.legend = lintxt+' for '+ntxt+' in ['+cmdm.date1+'-'+cmdm.spec+'] '+corr_txt+' -'
 
 ; additional computations (pac_5 nino_3 nino_4 averages) *** requires whole domain ****

trunk/usr/plt_def.pro

@@ -170 +170 @@
                       ; set to -99999 to desactivate
    linfit_map =  ''  ; ''/'p'/'m' to compute fit for all/before/after linfit_sep (maps only)
-
+   corr_thres =  0.1 ; correlation threshold to plot regression maps (between 0.1 and 0.2)
                       ; bining interval (one number) or 
                       ; binning array [b1,b2,b3,b4,...]