Ignore:
Timestamp:
Mar 4, 2014, 12:00:10 PM (11 years ago)
Author:
aslmd
Message:

PYTHON. updated analysis tools for dust devil analysis.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/UTIL/PYTHON/powerlaw/ddstat.py

    r1070 r1197  
    33import numpy as np
    44
    5 ##########################################################
    6 namefile = "/home/aymeric/Big_Data/LES_dd/psfc_f18.ncm1_2.txt"
    7 namefile = "/home/aymeric/Big_Data/LES_dd/psfc_f18.ncm2_2.txt"
    8 ##########################################################
    9 namefile = "/home/aymeric/Big_Data/LES_dd/psfc.LMD_LES_MARS.160564.ncm1_2.txt"
    10 #namefile = "/home/aymeric/Big_Data/LES_dd/sav/psfc.LMD_LES_MARS.160564.ncm2_2.txt"
    11 ##########################################################
    12 #namefile = "/home/aymeric/Big_Data/LES_dd/psfc_oldinsight100m.ncm1_2.txt"
    13 ##########################################################
     5### STATDD
     6def statdd(namefile):
    147
    15 
    16  
    17 
    18 #case = "188324p"
    19 #case = "191798"
    20 #case = "160564p"
    21 #case = "156487"
    22 #case = "2007p"
    23 #case = "13526p"
    24 #case = "172097"
    25 
    26 
    27 namefile = "/planeto/aslmd/LESdata/"+case+".ncm1_2.txt"
    28 
    29 
    30 
    31 
    32 
    33 # load data
    34 data = np.loadtxt(namefile,delimiter=";")
    35 t = data[:,0] ; n = data[:,1] ; s = data[:,2] ; d = data[:,3]
    36 
    37 # remove size and drop point when no vortex detected
    38 d[np.where(n==0)] = np.nan
    39 s[np.where(n==0)] = np.nan
    40 
    41 ## PLOTS
    42 
    43 number = ppplot.plot1d()
    44 number.f = n
    45 number.x = t
    46 number.linestyle = ''
    47 number.marker = '.'
    48 number.color = 'b'
    49 number.xlabel = "Local time (hour)"
    50 number.ylabel = "Detected vortices"
    51 number.makeshow()
    52 
    53 drop = ppplot.plot1d()
    54 drop.f = d
    55 drop.x = t
    56 drop.linestyle = ''
    57 drop.marker = '.'
    58 drop.color = 'r'
    59 drop.fmt = "%.1f"
    60 drop.xlabel = "Local time (hour)"
    61 drop.ylabel = "Maximum drop of detected vortices (Pa)"
    62 drop.makeshow()
    63 
    64 size = ppplot.plot1d()
    65 size.f = s
    66 size.x = t
    67 size.linestyle = ''
    68 size.marker = '.'
    69 size.color = 'g'
    70 size.xlabel = "Local time (hour)"
    71 size.ylabel = "Maximum size of detected vortices (m)"
    72 size.makeshow()
     8    # load data
     9    data = np.loadtxt(namefile,delimiter=";")
     10    t = data[:,0] ; n = data[:,1] ; s = data[:,2] ; d = data[:,3]
     11   
     12    # remove size and drop point when no vortex detected
     13    d[np.where(n==0)] = np.nan ; s[np.where(n==0)] = np.nan
     14   
     15    ## PLOTS
     16    number = ppplot.plot1d()
     17    number.f = n
     18    number.x = t
     19    number.linestyle = ''
     20    number.marker = '.'
     21    number.color = 'b'
     22    number.xlabel = "Local time (hour)"
     23    number.ylabel = "Detected vortices"
     24    number.makeshow()
     25   
     26    drop = ppplot.plot1d()
     27    drop.f = d
     28    drop.x = t
     29    drop.linestyle = ''
     30    drop.marker = '.'
     31    drop.color = 'r'
     32    drop.fmt = "%.1f"
     33    drop.xlabel = "Local time (hour)"
     34    drop.ylabel = "Maximum drop of detected vortices (Pa)"
     35    drop.makeshow()
     36   
     37    size = ppplot.plot1d()
     38    size.f = s
     39    size.x = t
     40    size.linestyle = ''
     41    size.marker = '.'
     42    size.color = 'g'
     43    size.xlabel = "Local time (hour)"
     44    size.ylabel = "Maximum size of detected vortices (m)"
     45    size.makeshow()
Note: See TracChangeset for help on using the changeset viewer.