2 # -*- coding: iso-8859-1 -*-
4 import matplotlib
as mpl
5 mpl
.rc('text', usetex
= True)
8 import matplotlib
.cm
as cm
9 import matplotlib
.mlab
as mlab
10 import matplotlib
.pyplot
as plt
11 import matplotlib
.patches
as patch
12 import matplotlib
.colors
as col
14 FOUT
= "../latex/images/periodadding.png"
16 FIN
.append("../rk_out0.dat")
17 #FIN.append("../rk_out1.dat")
18 #FIN.append("../rk_out175.dat")
19 FIN
.append("../rk_out2.dat")
20 #FIN.append("../rk_out225.dat")
21 #FIN.append("../rk_out23.dat")
22 FIN
.append("../rk_out25.dat")
23 #FIN.append("../rk_out26.dat")
24 #FIN.append("../rk_out27.dat")
25 #FIN.append("../rk_out35.dat")
26 #FIN.append("../rk_out4.dat")
27 FIN
.append("../rk_out5.dat")
28 FIN
.append("../rk_out6.dat")
35 res
= (1024/80,768/80)
37 #######################################################################
39 plots
= range(len(FIN
))
42 nPlots
= nCols
*len(plots
)
49 data
.append(mlab
.csv2rec(FIN
[i
], delimiter
='\t', comments
='#'))
60 phi
.append(data
[i
].phi
)
62 tstep
= t
[0][1] - t
[0][0]
67 fromt_eff
= fromt
/tstep
81 trunc_t
.append(tmpt
[fromt
*tau
/tstep
:tot
*tau
/tstep
])
82 trunc_x
.append(tmpx
[fromt
*tau
/tstep
:tot
*tau
/tstep
])
83 trunc_y
.append(tmpy
[fromt
*tau
/tstep
:tot
*tau
/tstep
])
84 trunc_phi
.append(tmpphi
[fromt
*tau
/tstep
:tot
*tau
/tstep
])
86 fig
= plt
.figure(figsize
=res
)
88 ax
= [ [] for DUMMYVAR
in range(nPlots
) ]
91 for i
in range(nPlots
):
92 ax
[i
] = fig
.add_subplot(nRows
,nCols
,i
+1)
94 tticks
= np
.arange(0*tau
,100*tau
,2*tau
)
98 for i
in range(len(tticks
)):
99 tticks_eff
.append(int(tticks
[i
]))
100 xticks
= np
.arange(-100*L
,100*L
,2*L
)
101 yticks
= np
.arange(0,10*L
,L
)
102 for i
in range(len(xticks
)):
103 xticks_eff
.append(int(xticks
[i
]))
104 for i
in range(len(yticks
)):
105 yticks_eff
.append(int(yticks
[i
]))
106 phiticks
= np
.arange(0,360+1,360/4)
121 ax
[j
].plot(t
[i
],x
[i
], color
='black')
122 ax
[j
+1].plot(x
[i
],y
[i
], color
='black')
123 ax
[j
+2].plot(t
[i
],phi
[i
], color
='black')
130 ax
[j
].set_xticks(tticks_eff
)
132 ax
[j
].set_yticks(yticks_eff
)
134 ax
[j
].set_yticks(xticks_eff
)
135 ax
[j
].set_ylabel(r
'$x$')
137 ax
[j
].set_xlabel(r
'$t$')
138 ax
[j
].set_xlim(fromt
*tau
-1,tot
*tau
+1)
139 ax
[j
].set_ylim(min(trunc_x
[l
])-1,max(trunc_x
[l
])+1)
147 ax
[j
].set_xticks(yticks_eff
)
149 ax
[j
].set_xticks(xticks_eff
)
150 ax
[j
].set_yticks(yticks_eff
)
152 ax
[j
].set_xlabel(r
'$x$')
153 ax
[j
].set_ylabel(r
'$y$')
154 ax
[j
].set_xlim(min(trunc_x
[l
])-1,max(trunc_x
[l
])+1)
155 ax
[j
].set_ylim(min(trunc_y
[l
])-1,max(trunc_y
[l
])+1)
162 ax
[j
].set_xticks(tticks_eff
)
163 ax
[j
].set_yticks(phiticks
)
164 ax
[j
].set_xlim(min(trunc_t
[l
])-1,max(trunc_t
[l
])+1)
165 ax
[j
].set_ylim(135,315)
167 ax
[j
].set_xlabel(r
'$t$')
168 ax
[j
].set_ylabel(r
'$\varphi$')
172 for i
in range(nPlots
/nCols
):
173 fig
.text(0.01,1-0.18-i
*0.17,eta_label
[i
], fontsize
=fs_label
)