Changeset 943 for trunk/UTIL/PYTHON/powerlaw

Timestamp:

Apr 27, 2013, 2:49:02 PM (12 years ago)

Author:

aslmd

Message:

clean and organized UTIL/PYTHON folder

Location:

trunk/UTIL/PYTHON/powerlaw

Files:

• . (added)
• find_devils.py (copied) (copied from trunk/UTIL/PYTHON/examples/scripts/find_devils.py) (9 diffs)
• plfit.py (added)
• plplot.py (added)
• plpva.py (added)
• plvar.py (added)
• randht.py (added)

trunk/UTIL/PYTHON/powerlaw/find_devils.py

@@ -27 +27 @@
    
     ### LOAD NETCDF DATA
-    #filename = "/home/aymeric/Big_Data/psfc.nc"
     nc = Dataset(filename) 
     psfc = nc.variables["PSFC"]
@@ -38 +37 @@
     depression = []
     stride = 1 #5
-    for i in range(0,shape[0],stride):
+    stride = 20
+    #stride = 50
+    #stride = 100
+    start = 0
+    start = stride
+    for i in range(start,shape[0],stride):
 
         psfc2d = np.array ( psfc [ i, : , : ] )
@@ -52 +56 @@
                    ## 2.5 clearly too low, 3.5 not too bad, 4 probably good
         fac = 3.5
+        fac = 3.2
+        ##fac = 2.5
+        #fac = 3. ## final choice
         #fac = 2.5
-        fac = 3.
         lim = ave - fac*std
         where = np.where(psfc2d < lim)
@@ -64 +70 @@
 
         draw = False
+        #draw = True
         if draw:
         ##################################################################################
@@ -74 +81 @@
             what_I_plot = myp.bounds(what_I_plot,zevmin,zevmax)
             zelevels = np.linspace(zevmin,zevmax)
-            fig = plt.figure() 
-            sub = myp.definesubplot(2,fig)
-            plt.subplot(sub) 
+            fig = plt.figure(figsize=(16,8)) 
+            subv,subh = myp.definesubplot(2,fig)
+            plt.subplot(subv,subh,1) 
             plt.contourf(what_I_plot,zelevels,cmap=palette) 
             plt.colorbar(fraction=0.05,pad=0.03,format="%.1f",\
                            ticks=np.linspace(zevmin,zevmax,ndiv+1),\
                            extend='both',spacing='proportional')
-            plt.subplot(sub+1)
-            palette = plt.get_cmap(name="binary")
+            plt.subplot(subv,subh,2)
+            palette = plt.get_cmap(name="hot") #"binary")
+            #lab[np.where(lab > 0.99)] = np.NaN
+            #plt.pcolor(lab,cmap=palette)
             plt.contourf(lab,2,cmap=palette)
             plt.show() 
@@ -104 +113 @@
         #sizex = detsize( xx, res = 10, loga=False, thres=3 )
         #sizey = detsize( yy, res = 10, loga=False, thres=3 )
+        sizex = detsize( xx, res = 15, loga=False, thres=2 )
+        sizey = detsize( yy, res = 15, loga=False, thres=2 )
         ###
+        print sizex, sizey
         #if ( mym.max(sizex) > mym.max(sizey) ): sizey = sizex  ### un peu limite dans certains cas
         if (len(sizex) > len(sizey))     : sizey = sizex        ### plus fidele mais petit souci lorsque PBC
@@ -137 +149 @@
 
 save = True
-save = False
+#save = False
 pression = False
-pression = True
+#pression = True
+
+filename = "/home/aymeric/Big_Data/psfc_f18.nc"
 
 if save:
-    allsizesx, allsizesy, depression = getsize("/home/aymeric/Big_Data/psfc.nc")
+    allsizesx, allsizesy, depression = getsize(filename)
     ### sauvegarde texte pour inspection
     mym.writeascii(allsizesx,'allsizex.txt')
@@ -198 +212 @@
 if not pression:
     zebins = zebins [ zebins > 15. ] 
+    #zebins = zebins [ zebins > 20. ]
+    zebins = zebins [ zebins > 25. ]
     zebins = zebins [ zebins < 1000. ]
 else:
@@ -288 +304 @@
 #plothist = np.append(plothist,400.*np.ones(2)  )
 
+exit()
+
 plt.figure(3)
 [alpha, xmin, L] = plfit.plfit(plothist,'xmin',30)#50.)