fnames = [os.path.join(tdir, 'x%02d.jpg' % i) for i in range(n_img) ] f = m.figure(1, size = (1000, 500), fgcolor = (0, 0, 0), bgcolor = (1., 1., 1.)) s = m.surf(n_e_arr[1], n_e_arr[2], n_mu_q_arr[0, :, :]) ms = s.mlab_source m.axes(s, color = (.7, .7, .7), extent = (-1, 1, 0, 1, 0, 1), ranges = (-0.21, 0.21, 0.1, 20, 0, max_mu_q), xlabel = 'x[mm]', ylabel = 'Lr[mm]', zlabel = 'f[N]',) m.view(-60.0, 70.0, focalpoint = [0., 0.45, 0.45]) m.savefig(fnames[0]) for i, fname in enumerate(fnames[1:]): ms.scalars = n_mu_q_arr[i, :, :] m.savefig(fname) images = string.join(fnames, ' ') destination = os.path.join('fig', 'fiber_cb_8p_anim.gif') import platform if platform.system() == 'Linux': os.system('convert ' + images + ' ' + destination) else: raise NotImplementedError, 'film production available only on linux' m.show()
m.axes( s, color=(.7, .7, .7), extent=(-1, 1, 0, 1, 0, 1), ranges=(-0.21, 0.21, 0.1, 20, 0, max_mu_q), xlabel='x', ylabel='Lr', zlabel='Force', ) m.view(-60.0, 70.0, focalpoint=[0., 0.45, 0.45]) # Store the information view = m.view() roll = m.roll() print(('view', view)) print(('roll', roll)) print((n_mu_q_arr.shape[2])) ms = s.mlab_source for i in range(1, n_mu_q_arr.shape[0]): ms.scalars = n_mu_q_arr[i, :, :] fname = 'x%02d.jpg' % i print(('saving', fname)) m.savefig(fname) sleep(0.1) # m.surf( n_e_arr[0], n_e_arr[1], n_mu_q_arr + n_std_q_arr ) # m.surf( n_e_arr[0], n_e_arr[1], n_mu_q_arr - n_std_q_arr ) m.show()
s = m.surf(n_e_arr[1], n_e_arr[2], n_mu_q_arr[0, :, :]) ms = s.mlab_source m.axes( s, color=(.7, .7, .7), extent=(-1, 1, 0, 1, 0, 1), ranges=(-0.21, 0.21, 0.1, 20, 0, max_mu_q), xlabel='x[mm]', ylabel='Lr[mm]', zlabel='f[N]', ) m.view(-60.0, 70.0, focalpoint=[0., 0.45, 0.45]) m.savefig(fnames[0]) for i, fname in enumerate(fnames[1:]): ms.scalars = n_mu_q_arr[i, :, :] m.savefig(fname) images = string.join(fnames, ' ') destination = os.path.join('fig', 'fiber_cb_8p_anim.gif') import platform if platform.system() == 'Linux': os.system('convert ' + images + ' ' + destination) else: raise NotImplementedError, 'film production available only on linux' m.show()
f = m.figure(1, size = (1000, 500), fgcolor = (0, 0, 0), bgcolor = (1., 1., 1.)) s = m.surf(n_e_arr[1], n_e_arr[2], n_mu_q_arr[0, :, :]) m.axes(s, color = (.7, .7, .7), extent = (-1, 1, 0, 1, 0, 1), ranges = (-0.21, 0.21, 0.1, 20, 0, max_mu_q), xlabel = 'x', ylabel = 'Lr', zlabel = 'Force',) m.view(-60.0, 70.0, focalpoint = [0., 0.45, 0.45]) # Store the information view = m.view() roll = m.roll() print 'view', view print 'roll', roll print n_mu_q_arr.shape[2] ms = s.mlab_source for i in range(1, n_mu_q_arr.shape[0]): ms.scalars = n_mu_q_arr[i, :, :] fname = 'x%02d.jpg' % i print 'saving', fname m.savefig(fname) sleep(0.1) # m.surf( n_e_arr[0], n_e_arr[1], n_mu_q_arr + n_std_q_arr ) # m.surf( n_e_arr[0], n_e_arr[1], n_mu_q_arr - n_std_q_arr ) m.show()