Skip to content
Snippets Groups Projects
Commit a6326fac authored by Mohcine Chraibi's avatar Mohcine Chraibi
Browse files

Upload New File

parent 83437183
Branches
No related tags found
No related merge requests found
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Tue Apr 9 12:29:26 2019
@author: administer
"""
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Mon Apr 8 20:36:43 2019
@author: administer
"""
import subprocess
import os
import glob
import shutil
import numpy as np
import sys
#path to jpscore
executable = "/Users/administer/Workspace/jpscore/bin/jpscore"
#width here represents only a list ,its not a directory here
Model = glob.glob("m*")
for model in Model:
os.chdir(model)
Inflowrates=glob.glob('inflowrate=*')
for inflowrate in Inflowrates:
os.chdir(inflowrate)
subprocess.call(["python", "makeini.py", "-f", "master_ini.xml"])
os.chdir('..')
os.chdir('..')
print("inifiles are finished")
print('from here begins the simulation ')
Model = glob.glob("m*")
for model in Model:
os.chdir(model)
Inflowrates=glob.glob('inflowrate=*')
for inflowrate in Inflowrates:
os.chdir(inflowrate)
inifiles = glob.glob("inifiles/*.xml")
for inifile in inifiles:
subprocess.call([executable, "--inifile=%s"%inifile])
for filename in os.listdir(os.path.abspath("inifiles")):
if filename.endswith('.txt'):
shutil.move(os.path.abspath("inifiles")+"/"+filename,os.path.abspath("clogginglog")+'/')
os.path.abspath("clogginglog")
os.chdir("clogginglog")
print(os.getcwd())
os.rename(filename,inflowrate.split('=')[1]+filename)
os.chdir("..")
os.chdir("..")
os.chdir('..')
os.chdir('..')
import matplotlib.pyplot as plt
import os
import glob
import numpy as np
j=0
Model = glob.glob("m*")
BIG_meanNmax=[[],[],[],[],[]]
BIG_meanTimespan=[[],[],[],[],[]]
BIG_C=[[],[],[],[],[]]
W=['inflowrate=2','inflowrate=4','inflowrate=8','inflowrate=10','inflowrate=16']
Nmax_w1=[]
Nmax_w2=[]
Nmax_w3=[]
Nmax_w4=[]
Nmax_w5=[]
Timespan_w1=[]
Timespan_w2=[]
Timespan_w3=[]
Timespan_w4=[]
Timespan_w5=[]
for model in Model:
MeanNmax=[0,0,0,0,0]
meanTimespan=[0,0,0,0,0]
os.chdir(model)
Inflowrates=glob.glob('=*')
for inflowrate in Inflowrates:
os.chdir(inflowrate+'/clogginglog/')
Logs=glob.glob('*.txt')
for log in Logs:
if log.split('Clo')[0] == '2':
if np.loadtxt(log).shape ==(3,):
Nmax_w1.append(np.loadtxt(log)[-1])
MeanNmax[0]=sum(Nmax_w1)/len(Nmax_w1)
Timespan_w1.append(np.loadtxt(log)[1])
meanTimespan[0]=sum(Timespan_w1)/len(Timespan_w1)
else :
Nmax_w1.append(np.loadtxt(log)[:,2][-1])
MeanNmax[0]=sum(Nmax_w1)/len(Nmax_w1)
Timespan_w1.append(np.loadtxt(log)[-1][1]-np.loadtxt(log)[0][1])
meanTimespan[0]=sum(Timespan_w1)/len(Timespan_w1)
if log.split('Clo')[0] == '4':
if np.loadtxt(log).shape ==(3,):
Nmax_w2.append(np.loadtxt(log)[-1])
MeanNmax[1]=sum(Nmax_w2)/len(Nmax_w2)
Timespan_w2.append(np.loadtxt(log)[1])
meanTimespan[1]=sum(Timespan_w2)/len(Timespan_w2)
else :
Nmax_w2.append(np.loadtxt(log)[:,2][-1])
MeanNmax[1]=sum(Nmax_w2)/len(Nmax_w2)
Timespan_w2.append(np.loadtxt(log)[-1][1]-np.loadtxt(log)[0][1])
meanTimespan[1]=sum(Timespan_w2)/len(Timespan_w2)
if log.split('Clo')[0] == '8':
if np.loadtxt(log).shape ==(3,):
Nmax_w3.append(np.loadtxt(log)[-1])
MeanNmax[2]=sum(Nmax_w3)/len(Nmax_w3)
Timespan_w3.append(np.loadtxt(log)[1])
meanTimespan[2]=sum(Timespan_w3)/len(Timespan_w3)
else :
Nmax_w3.append(np.loadtxt(log)[:,2][-1])
MeanNmax[2]=sum(Nmax_w3)/len(Nmax_w3)
Timespan_w3.append(np.loadtxt(log)[-1][1]-np.loadtxt(log)[0][1])
meanTimespan[2]=sum(Timespan_w3)/len(Timespan_w3)
if log.split('Clo')[0] == '10':
if np.loadtxt(log).shape ==(3,):
Nmax_w4.append(np.loadtxt(log)[-1])
MeanNmax[3]=sum(Nmax_w4)/len(Nmax_w4)
Timespan_w4.append(np.loadtxt(log)[1])
meanTimespan[3]=sum(Timespan_w4)/len(Timespan_w4)
else :
Nmax_w4.append(np.loadtxt(log)[:,2][-1])
MeanNmax[3]=sum(Nmax_w4)/len(Nmax_w4)
Timespan_w4.append(np.loadtxt(log)[-1][1]-np.loadtxt(log)[0][1])
meanTimespan[3]=sum(Timespan_w4)/len(Timespan_w4)
if log.split('Clo')[0] == '16':
if np.loadtxt(log).shape ==(3,):
Nmax_w5.append(np.loadtxt(log)[-1])
MeanNmax[4]=sum(Nmax_w5)/len(Nmax_w5)
Timespan_w5.append(np.loadtxt(log)[1])
meanTimespan[4]=sum(Timespan_w5)/len(Timespan_w5)
else:
Nmax_w5.append(np.loadtxt(log)[:,2][-1])
MeanNmax[4]=sum(Nmax_w5)/len(Nmax_w5)
Timespan_w5.append(np.loadtxt(log)[-1][1]-np.loadtxt(log)[0][1])
meanTimespan[4]=sum(Timespan_w5)/len(Timespan_w5)
os.chdir('../..')
print(MeanNmax)
print(meanTimespan)
BIG_meanNmax[j]= MeanNmax
BIG_meanTimespan[j]=meanTimespan
C=[0,0,0,0,0]
for i in range(len(MeanNmax)):
if meanTimespan[i]==0:
pass
else:
C[i]=MeanNmax[i]/meanTimespan[i]
print(C)
BIG_C[j]=C
j=j+1
plt.figure()
plt.plot(W,MeanNmax)
plt.savefig('inflowrate-n')
plt.figure()
plt.plot(W,C,'-or')
plt.savefig('inflowrate-nt')
print('now leave'+model)
os.chdir('..')
std1=np.std(BIG_meanNmax)
plt.figure()
plt.grid()
plt.title('standard deviation is %s'%std1)
plt.plot(W,BIG_meanNmax[0],label='M1')
plt.plot(W,BIG_meanNmax[1],label='M0')
plt.plot(W,BIG_meanNmax[2],label='M2')
plt.plot(W,BIG_meanNmax[3],label='M3')
plt.plot(W,BIG_meanNmax[4],label='M4')
plt.xlabel('inflowrate in (m)')
plt.ylabel('mean clogging times')
plt.legend()
plt.savefig('figure/multiM-inflowrate-n',dpi=300)
std2=np.std(BIG_C)
plt.figure()
plt.grid()
plt.title('standard deviation is %s'%std2)
plt.plot(W,BIG_C[0],label='M1')
plt.plot(W,BIG_C[1],label='M0')
plt.plot(W,BIG_C[2],label='M2')
plt.plot(W,BIG_C[3],label='M3')
plt.plot(W,BIG_C[4],label='M4')
plt.xlabel('inflowrate in (m)')
plt.ylabel('mean clogging times per seconds')
plt.legend()
plt.savefig('figure/multiM-inflowrate-nt',dpi=300)
plt.xlabel('inflowrate in (m)')
plt.ylabel('mean clogging times per seconds')
plt.legend()
plt.savefig('figure/multiM-inflowrate-nt',dpi=300)
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment