User: knarf Date: 2011/10/17 12:21 AM

Added: /examples/collapse/ rho_max_conv.pdf, rho_max_conv.py

Modified: /examples/collapse/ make_plots.sh, radii.pdf, radii.py, rho_max.pdf

Log: one more matplotlib plot

Directory Changes:

Directory: /svn:executable/ ===========================

+ * Directory: /svn:mime-type/ ==========================

+ application/pdf

File Changes:

Directory: /examples/collapse/ ==============================

File [modified]: make_plots.sh Delta lines: +3 -3 =================================================================== --- examples/collapse/make_plots.sh 2011-10-17 04:58:40 UTC (rev 149) +++ examples/collapse/make_plots.sh 2011-10-17 05:21:57 UTC (rev 150) @@ -1,5 +1,5 @@ #!/bin/bash +#gnuplot *.plt

-# use find_edge.awk to find the outer edge of the star if needed - -gnuplot *.plt +echo ./rho_max.py ./radii.py ./rho_max_conv.py \ +| xargs -n 1 -P 4 bash -c

File [modified]: radii.pdf Delta lines: +0 -0 =================================================================== (Binary files differ)

File [modified]: radii.py Delta lines: +1 -2 =================================================================== --- examples/collapse/radii.py 2011-10-17 04:58:40 UTC (rev 149) +++ examples/collapse/radii.py 2011-10-17 05:21:57 UTC (rev 150) @@ -36,8 +36,7 @@ ax2.xaxis.set_major_locator(mticker.MaxNLocator(7)) ax2.xaxis.set_minor_locator(mticker.MaxNLocator(14)) ax.set_ylabel(r'$R/R_G$') -ax.yaxis.set_major_locator(mticker.MaxNLocator(5)) -ax.yaxis.set_minor_locator(mticker.MaxNLocator(10)) +ax.yaxis.set_major_locator(mticker.MaxNLocator(8)) ax.yaxis.grid(False) set_tick_sizes(ax, 8, 4)

File [modified]: rho_max.pdf Delta lines: +0 -0 =================================================================== (Binary files differ)

File [added]: rho_max_conv.pdf Delta lines: +0 -0 =================================================================== (Binary files differ)

Property changes on: examples/collapse/rho_max_conv.pdf ___________________________________________________________________

File [added]: rho_max_conv.py Delta lines: +46 -0 =================================================================== --- examples/collapse/rho_max_conv.py (rev 0) +++ examples/collapse/rho_max_conv.py 2011-10-17 05:21:57 UTC (rev 150) @@ -0,0 +1,46 @@ +#!/usr/bin/python + +# Plot evolution of convergence factor of hamiltonian constraint violation at x=0 +from plot_defaults import * +import math + +# stuff +xlim = (0,1400) +ylim = (1,5) + +# load data +Fx, Fy3, Fy5, Fy6 = np.loadtxt("rho_maximum_convergence.asc", comments="#", usecols=(1,2,4,5), unpack=True) + +# Plot basics +fig = plt.figure(figsize=(6, 3)) +fig.subplots_adjust(top=0.85, bottom=0.16, left=0.11,right=0.97) +ax = fig.add_subplot(1,1,1) + +for i in range(len(Fy3)): + Fy3[i] = math.log(max(1.e-10, (Fy3[i]-Fy5[i])/(Fy5[i]-Fy6[i])),2) +# plot +ax.plot(Fx, Fy3, linestyle='-' , color='black') + +# plot properties +#ax.set_xlim(xlim) +ax.set_ylim(ylim) + +#plt.title('PSD', fontsize=fontsize) + +ax.set_xlabel(r't [M]') +ax.xaxis.set_major_locator(mticker.MaxNLocator(7)) +ax.xaxis.set_minor_locator(mticker.MaxNLocator(14)) +ax.xaxis.grid(False) +ax2 = ax.twiny() +ax2.set_xlabel(r't [ms]') +ax2.set_xlim((ax.get_xlim()[0]/M_to_ms, ax.get_xlim()[1]/M_to_ms)) +ax2.xaxis.set_major_locator(mticker.MaxNLocator(7)) +ax2.xaxis.set_minor_locator(mticker.MaxNLocator(14)) +ax.set_ylabel(r'$C(\varrho_C)$') +ax.yaxis.set_major_locator(mticker.MaxNLocator(8)) +ax.yaxis.grid(False) +set_tick_sizes(ax, 8, 4) + +#plt.show() +plt.savefig('rho_max_conv.pdf') +

Property changes on: examples/collapse/rho_max_conv.py ___________________________________________________________________