Changeset 1053 for trunk/UTIL/PYTHON

Timestamp:

Sep 26, 2013, 1:29:22 PM (11 years ago)

Author:

aslmd

Message:

UTIL PYTHON. tools to study dust devils in LES.

Location:

trunk/UTIL/PYTHON/powerlaw

Files:

• ddhist.py (added)
• ddstat.py (added)
• dropsize.py (added)
• finddd.py (added)
• performdd.py (added)
• previous_tests (added)
• previous_tests/find_devils.py (moved) (moved from trunk/UTIL/PYTHON/powerlaw/find_devils.py) (6 diffs)
• previous_tests/plfit.py (moved) (moved from trunk/UTIL/PYTHON/powerlaw/plfit.py)
• previous_tests/plot_granulo.py (moved) (moved from trunk/UTIL/PYTHON/powerlaw/plot_granulo.py)
• previous_tests/plplot.py (moved) (moved from trunk/UTIL/PYTHON/powerlaw/plplot.py)
• previous_tests/plpva.py (moved) (moved from trunk/UTIL/PYTHON/powerlaw/plpva.py)
• previous_tests/plvar.py (moved) (moved from trunk/UTIL/PYTHON/powerlaw/plvar.py)
• previous_tests/randht.py (moved) (moved from trunk/UTIL/PYTHON/powerlaw/randht.py)
• showdd.py (added)

trunk/UTIL/PYTHON/powerlaw/previous_tests/find_devils.py

@@ -22 +22 @@
     import numpy as np
     from scipy.ndimage.measurements import minimum_position
+    from scipy import ndimage
     from netCDF4 import Dataset
     import matplotlib.pyplot as plt
@@ -29 +30 @@
     nc = Dataset(filename) 
     psfc = nc.variables["PSFC"]
+    print "yeah"
     
     ### LOOP on TIME
@@ -39 +41 @@
     stride = 20
     #stride = 50
-    #stride = 100
+    stride = 100
     start = 0
     start = stride
@@ -66 +68 @@
         depression = np.append(depression,np.ravel(psfc2d[where])-ave)
    
+        ## lab is 0 or 1
         lab = np.zeros(np.array(psfc2d).shape) ## points to be treated by the minimum_position routine
         lab[where] = 1.  ## do not treat points close to 'mean' (background) pressure
-
-        draw = False
-        #draw = True
-        if draw:
-        ##################################################################################
-            vmin = -0.3
-            vmax =  0.0
-            ndiv = 3
-            palette = plt.get_cmap(name="YlGnBu") 
-            what_I_plot = psfc2d-ave
-            zevmin, zevmax = myp.calculate_bounds(what_I_plot,vmin=vmin,vmax=vmax)
-            what_I_plot = myp.bounds(what_I_plot,zevmin,zevmax)
-            zelevels = np.linspace(zevmin,zevmax)
-            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(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() 
-        ##################################################################################
   
         xx = []
@@ -151 +127 @@
 #save = False
 pression = False
-#pression = True
+pression = True
 
 filename = "/home/aymeric/Big_Data/psfc_f18.nc"
 
 if save:
+    ### getsize
     allsizesx, allsizesy, depression = getsize(filename)
     ### sauvegarde texte pour inspection
@@ -162 +139 @@
     mym.writeascii(depression,'alldepression.txt')
     ### sauvegarde binaire pour utilisation python
-    myfile = open('allsizex.bin', 'wb')
-    pickle.dump(allsizesx, myfile)
-    myfile.close()
-    myfile = open('allsizey.bin', 'wb')
-    pickle.dump(allsizesy, myfile)
-    myfile.close()
-    myfile = open('alldepression.bin', 'wb')
-    pickle.dump(depression, myfile)
-    myfile.close()
+    myfile = open('allsizex.bin', 'wb') ; pickle.dump(allsizesx, myfile) ; myfile.close()
+    myfile = open('allsizey.bin', 'wb') ; pickle.dump(allsizesy, myfile) ; myfile.close()
+    myfile = open('alldepression.bin', 'wb') ; pickle.dump(depression, myfile) ; myfile.close()
 
 ### load files