Changeset 1007 for trunk/UTIL/PYTHON/planetoplot_v2

Timestamp:

Jul 19, 2013, 10:35:53 AM (11 years ago)

Author:

aslmd

Message:

UTIL PYTHON. various updates on example + mcd + plot scripts. nothing major.

Location:

trunk/UTIL/PYTHON/planetoplot_v2

Files:

• ppcompute.py (modified) (1 diff)
• ppplot.py (modified) (2 diffs)

trunk/UTIL/PYTHON/planetoplot_v2/ppcompute.py

@@ -95 +95 @@
            else:
               return np.array(field).sum(axis=axis)
+
+## compute mean over bins
+## author AS
+def meanbin(y,x,bins):
+    bins.sort() # sorting is needed for binning
+    meanvalues = []
+    for iii in range(len(bins)):
+        ## GET VALUES OVER RELEVANT INTERVALS
+        if iii == 0:
+          ind = x<bins[iii+1] ; ys = y[ind]
+        elif iii == len(bins)-1:
+          ind = x>=bins[iii-1] ; ys = y[ind]
+        else:
+          ind = x>=bins[iii-1] ; interm = x[ind] ; intermf = y[ind]
+          ind = interm<bins[iii+1] ; xs = interm[ind] ; ys = intermf[ind]
+        ## COMPUTE MEAN and TREAT NAN CASE
+        meanvalues.append(mean(ys))
+    ## RETURN A NUMPY ARRAY
+    meanvalues = np.array(meanvalues)
+    return meanvalues
 
 ################

trunk/UTIL/PYTHON/planetoplot_v2/ppplot.py

@@ -180 +180 @@
 # useful functions #
 ####################
+
+# continuity with matplotlib
+def close():
+    mpl.close()
+
+# a function for predefined figure sizes
+def figuref(x=16,y=9):
+    mpl.figure(figsize=(x,y))
 
 # a function to change color lines according to a color map (idea by A. Pottier)
@@ -649 +657 @@
                 objC = mpl.contour(x, y, what_I_contour, \
                             zelevelsc, colors = ccol, linewidths = cline)
-                #mpl.clabel(objC, inline=1, fontsize=10)
+                #mpl.clabel(objC, inline=1, fontsize="small",\
+                #             inline_spacing=1,fmt="%.0f")
             mpl.contourf(x, y, \
                          self.field, \
                          zelevels, cmap=palette)
+            #mpl.pcolor(x,y,\
+            #             self.field, \
+            #             cmap=palette)
             # make log axes and/or invert ordinate
             ax = mpl.gca()