-
Notifications
You must be signed in to change notification settings - Fork 0
/
compute_rmsd_2.py
104 lines (81 loc) · 2.85 KB
/
compute_rmsd_2.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
#!/usr/bin/python
import sys
import numpy as np
## THIS .py scans 14 .rmsd.dat (each replica) and returns a file with the Temperature(T) LowestValue(LV) LowerQuartile(Q1) MiddleQuartile(Q2) HigherQuartile(Q3) e HighestValue(HV)
## usage: python compute_AMM.py output input1 input2 input3 ... inputn
#in_file = open(sys.argv[1], "rw") ........
#IMPORTANTE, O C-MD eh o 4 na numeracao das pastas, ele continua o dados4.dat, e isso eh solucionado na hora de criar o grafico yerrorbars para os graficos com quartis, eu modifiquei antes, ou seja, o dados7 virou o 8, arrumando a ordem e assim acabando com a necessidade de tratar isso na hora de criar o grafico de candlebars
text = ''
x = 2
while (x < 16):
in_file = sys.argv[x]
method = in_file.split('-')
if method[0] == "rmsd" :
with open(in_file, 'r') as f:
next(f)
#print "pulei"
#print method[0]
col=[float(row.split(" ")[-1]) for row in f]
else:
with open(in_file, 'r') as f:
#print "nao pulei, gdt" NAO ENTENDI PQ O RMSD TA ENTRANDO AQUI!
#print method[0]
next(f)
col=[float(row.split(" ")[-1]) for row in f]
words = in_file.split('.')
#aqui eu tenho o col contendo todos os rmsds daquela temperatura
#print col
#print "teste"
LV = min(col) # LV WRONG?
HV = max(col) #HV WRONG?
#print words
if words[1] == "001" or words[1] == "269" :
saida = "001"
if words[1] == "002" or words[1] == "285" :
saida = "003"
if words[1] == "003" or words[1] == "301" :
saida = "005"
if words[1] == "004" or words[1] == "318" :
saida = "007"
if words[1] == "005" or words[1] == "336" :
saida = "009"
if words[1] == "006" or words[1] == "355" :
saida = "011"
if words[1] == "007" or words[1] == "374" :
saida = "013"
if words[1] == "008" or words[1] == "395" :
saida = "015"
if words[1] == "009" or words[1] == "416" :
saida = "017"
if words[1] == "010" or words[1] == "438" :
saida = "019"
if words[1] == "011" or words[1] == "461" :
saida = "021"
if words[1] == "012" or words[1] == "485" :
saida = "023"
if words[1] == "013" or words[1] == "511" :
saida = "025"
if words[1] == "014" or words[1] == "537" :
saida = "027"
line = saida + ', {}, {},'.format(LV,HV) + '\n' #ultima virgula necessaria pq uso o paste dps em mais de um arquivo
#print line
text += line
x = x + 1
#print text
out = open(sys.argv[1], 'w')
out.write(text)
out.close()
# "269.50" 1 ,
# "285.22" 3 ,
# "301.66" 5 ,
# "318.76" 7 ,
# "336.61" 9 ,
# "355.26" 11 ,
# "374.74" 13 ,
# "395.07" 15 ,
# "416.32" 17 ,
# "438.50" 19 ,
# "461.67" 21 ,
# "485.87" 23 ,
# "511.14" 25 ,
# "537.54" 27 )