source: trunk/UTIL/PYTHON/powerlaw/ddstat.py @ 1142

#! /usr/bin/env python
import ppplot
import numpy as np

##########################################################
namefile = "/home/aymeric/Big_Data/LES_dd/psfc_f18.ncm1_2.txt"
namefile = "/home/aymeric/Big_Data/LES_dd/psfc_f18.ncm2_2.txt"
##########################################################
namefile = "/home/aymeric/Big_Data/LES_dd/psfc.LMD_LES_MARS.160564.ncm1_2.txt"
#namefile = "/home/aymeric/Big_Data/LES_dd/sav/psfc.LMD_LES_MARS.160564.ncm2_2.txt"
##########################################################
#namefile = "/home/aymeric/Big_Data/LES_dd/psfc_oldinsight100m.ncm1_2.txt"
##########################################################


 

#case = "188324p"
#case = "191798"
#case = "160564p"
#case = "156487"
#case = "2007p"
#case = "13526p"
#case = "172097"


namefile = "/planeto/aslmd/LESdata/"+case+".ncm1_2.txt"





# load data
data = np.loadtxt(namefile,delimiter=";")
t = data[:,0] ; n = data[:,1] ; s = data[:,2] ; d = data[:,3]

# remove size and drop point when no vortex detected
d[np.where(n==0)] = np.nan
s[np.where(n==0)] = np.nan

## PLOTS

number = ppplot.plot1d()
number.f = n
number.x = t
number.linestyle = ''
number.marker = '.'
number.color = 'b'
number.xlabel = "Local time (hour)"
number.ylabel = "Detected vortices"
number.makeshow()

drop = ppplot.plot1d()
drop.f = d
drop.x = t
drop.linestyle = ''
drop.marker = '.'
drop.color = 'r'
drop.fmt = "%.1f"
drop.xlabel = "Local time (hour)"
drop.ylabel = "Maximum drop of detected vortices (Pa)"
drop.makeshow()

size = ppplot.plot1d()
size.f = s
size.x = t
size.linestyle = ''
size.marker = '.'
size.color = 'g'
size.xlabel = "Local time (hour)"
size.ylabel = "Maximum size of detected vortices (m)"
size.makeshow()