|
Last change
on this file since 1242 was
1070,
checked in by aslmd, 11 years ago
|
|
PYTHON. updated analysis tools for dust devils.
|
-
Property svn:executable set to
*
|
|
File size:
1.1 KB
|
| Line | |
|---|
| 1 | #! /usr/bin/env python |
|---|
| 2 | import matplotlib.pyplot as mpl |
|---|
| 3 | import ppplot |
|---|
| 4 | import numpy as np |
|---|
| 5 | |
|---|
| 6 | ########################################################## |
|---|
| 7 | namefile = "/home/aymeric/Big_Data/LES_dd/psfc_f18.ncm1_1.txt" |
|---|
| 8 | namefile = "/home/aymeric/Big_Data/LES_dd/psfc_f18.ncm2_1.txt" |
|---|
| 9 | #namefile = "/home/aymeric/Big_Data/LES_dd/psfc.LMD_LES_MARS.160564.ncm1_1.txt" |
|---|
| 10 | #namefile = "/home/aymeric/Big_Data/LES_dd/psfc.LMD_LES_MARS.160564.ncm2_1.txt" |
|---|
| 11 | #namefile = "/home/aymeric/Big_Data/LES_dd/press_ustm_exomars.nc1.txt" |
|---|
| 12 | #namefile = "/home/aymeric/Big_Data/LES_dd/psfc_oldinsight100m.ncm1_1.txt" |
|---|
| 13 | ########################################################## |
|---|
| 14 | |
|---|
| 15 | #case = "188324p" |
|---|
| 16 | #case = "191798" |
|---|
| 17 | #case = "160564p" |
|---|
| 18 | #case = "156487" |
|---|
| 19 | #case = "2007p" |
|---|
| 20 | #case = "13526p" |
|---|
| 21 | #case = "172097" |
|---|
| 22 | |
|---|
| 23 | |
|---|
| 24 | namefile = "/planeto/aslmd/LESdata/"+case+".ncm1_1.txt" |
|---|
| 25 | |
|---|
| 26 | # load data |
|---|
| 27 | data = np.loadtxt(namefile,delimiter=";") |
|---|
| 28 | t = data[:,0] ; s = data[:,1] ; d = data[:,2] |
|---|
| 29 | i = data[:,3] ; j = data[:,4] # int? |
|---|
| 30 | |
|---|
| 31 | drop = ppplot.plot1d() ; drop.f = d ; drop.x = s |
|---|
| 32 | drop.linestyle = '' ; drop.marker = '.' ; drop.color = 'r' ; drop.fmt = "%.1f" |
|---|
| 33 | drop.xlabel = "Vortex size (m)" ; drop.ylabel = "Pressure drop (Pa)" |
|---|
| 34 | drop.makeshow() |
|---|
| 35 | |
|---|
Note: See
TracBrowser
for help on using the repository browser.
