Changeset 120

Timestamp:

09/11/08 15:59:55 (16 years ago)

Author:

ericg

Message:

Added +/- 1 stddev curves in bining plots

File:

• trunk/procs/ybinx.pro (modified) (8 diffs)

trunk/procs/ybinx.pro

@@ -5 +5 @@
          ; ----------------------
          ;
+               IF debug_w THEN print, "    "
                IF debug_w THEN print, "    Enter ybinx..."
+               IF debug_w THEN print, "    "
                
                ; organise density bins: either read array or build form single number and min/max in fld_glo_mmx.def
@@ -19 +21 @@
 
                IF debug_w THEN print, '    bin_plt= ', bin_plt 
+               IF debug_w THEN print, '      nbins= ', nbins
 
                ; 3rd field to bin regression ?
@@ -27 +30 @@
 
                IF sw3 THEN BEGIN 
-                  idxmsk = where(fld GT valmask/20.)
+                  idxmsk = where(fld GT valmask/10.)
                   IF idxmsk(0) NE -1 THEN fld(idxmsk) = valmask
                   IF idxmsk(0) NE -1 THEN fld3(idxmsk) = valmask
                ENDIF 
+
+               ; print min/max of field for debug
+               idxmskp = where(fld LE valmask/10.)
+               IF debug_w THEN print, '    Min/max fld= ', min(fld(idxmskp)), max(fld(idxmskp))
+               IF debug_w THEN print, '    Min/max fld2= ', min(fld2(idxmskp)), max(fld2(idxmskp))
 
                ; remove mean seasonal cycle if required
@@ -102 +110 @@
                max_idx = max(binpop)
 
-               IF debug_w THEN print, '  max bin size ', max_idx
-
+               IF debug_w THEN print, '    Max bin size ', max_idx
                fldy = fltarr(nbins+1, max_idx)
                fldy[*, *] = !values.f_nan 
@@ -152 +159 @@
                   WHILE ib LE nbins DO BEGIN 
                      
+                     binsz = binpop(ib)
                      
-                     binsz = binpop(ib)
-                     sfc_tot = total(surfb(ib, 0:binsz-1))
-                     IF sfc_tot GT 0 THEN BEGIN  
+                     IF binsz GT 1 THEN BEGIN  
+
+                        sfc_tot = total(surfb(ib, 0:binsz-1))
                         yplt(ib) = total(fldys(ib, 0:binsz-1))/sfc_tot
+                        yerr(ib) = sqrt((moment(fldy(ib, 0:binsz-1)))[1])
+
                         mean_fld = mean_fld + yplt(ib)*float(binpop(ib)) 
+
+                        ; print bin info
+                        IF ib GT 0 AND ib LT nbins THEN print, '    Bin size, occurence, average: ', bin_interval(ib-1), bin_interval(ib), binpop(ib), (binpop(ib)/total(binpop))*100., yplt(ib)
+                        IF ib EQ 0 THEN print, '    Bin size, occurence, average:      min' , bin_interval(ib), binpop(ib), (binpop(ib)/total(binpop))*100. , yplt(ib)
+                        IF ib EQ nbins THEN print, '    Bin size, occurence, average: ', bin_interval(ib-1),'     max ', binpop(ib), (binpop(ib)/total(binpop))*100. , yplt(ib)
                      ENDIF ELSE yplt(ib) = !values.f_nan
                   
@@ -184 +199 @@
 
                      IF binsz GT 1 THEN BEGIN 
+
                         idx1 = where(fldy(ib, 0:binsz-1) NE valmask)
                         idx2 = where(fldy2(ib, 0:binsz-1) NE valmask)
@@ -193 +209 @@
 
                         mean_fld = mean_fld + yplt(ib)*float(binpop(ib)) 
+
+                        ; print bin info
                         IF ib GT 0 AND ib LT nbins THEN print, '    Bin size, occurence, regress.: ', bin_interval(ib-1), bin_interval(ib), binpop(ib), (binpop(ib)/total(binpop))*100., yplt(ib)
                         IF ib EQ 0 THEN print, '    Bin size, occurence, regress.:      min' , bin_interval(ib), binpop(ib), (binpop(ib)/total(binpop))*100. , yplt(ib)
@@ -228 +246 @@
                res = execute(pltcmd)      
 
-               ; plot error bars for regression
-
-               IF sw3 THEN BEGIN 
-                  pltcmd = 'pltsc,yplt-yerr,bin_plt,minc,maxc,minc2,maxc2,varlegend'+com_strplt+',ov1d=1,COLOR=1, thick=1, STY1D=-1,subtitle=""'
-                  res = execute(pltcmd)      
-                  pltcmd = 'pltsc,yplt+yerr,bin_plt,minc,maxc,minc2,maxc2,varlegend'+com_strplt+',ov1d=1,COLOR=1, thick=1, STY1D=-1,subtitle=""'
-                  res = execute(pltcmd)      
-               ENDIF 
+               ; plot +/- 1 stdedv for field binnig
+
+               pltcmd = 'pltsc,yplt-yerr,bin_plt,minc,maxc,minc2,maxc2,varlegend'+com_strplt+',ov1d=1,COLOR=1, thick=1, STY1D=-1,subtitle=""'
+               res = execute(pltcmd)      
+               pltcmd = 'pltsc,yplt+yerr,bin_plt,minc,maxc,minc2,maxc2,varlegend'+com_strplt+',ov1d=1,COLOR=1, thick=1, STY1D=-1,subtitle=""'
+               res = execute(pltcmd)      
 
                IF debug_w THEN print, "    ... Exit ybinx"