Dear all
Im trying to use the file viewstar.py at the postporcessing folder with one of my 2d models but all that I get if a black square. I can't find what is wrong. This is what im doing
rom mpl_toolkits.mplot3d import Axes3D
from matplotlib import cm
from matplotlib.colors import LogNorm
import matplotlib.pyplot as plt
import matplotlib.colors as pltc
import numpy as np
import string
from ester import *
a=star2d('MI_FILE2d')
nth=a.nth+2 # we take care that the plotted fields as read by "a" have
# nth+2 points in latitude
u = np.linspace(0, 2 * np.pi, 100) # generate angle phi grid
v = a.th[-1,:] # angle theta
w=v+np.pi/2
v=np.concatenate([w[0:nth-1],v]) # angle theta pi (South pole) to 0
r_sud=a.r[-1,:]
r_nord=a.r[-1,nth:0:-1]
r=np.concatenate([r_nord,r_sud])/a.Re
teff_rsh=np.reshape(a.Teff,nth)
teff_ext=np.concatenate([teff_rsh[nth:0:-1],teff_rsh])
x = np.outer(np.cos(u), rnp.sin(v))
y = np.outer(np.sin(u), rnp.sin(v))
z = np.outer(np.ones(np.size(u)), r*np.cos(v))
Teff = np.outer(np.ones(np.size(u)), teff_ext)
N=1-(Teff/Teff.max())**4 # so that white of "Blues" shows the max of
Teff
m=cm.ScalarMappable(cmap=cm.Blues, norm = pltc.Normalize(vmin=0,vmax=1.,clip=False))
m.set_array(N)
i=112.
print ('inclination angle i = ',i)
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.plot_surface(x, y, z, rstride=1, cstride=1, facecolors=m.to_rgba(N),linewidth=0, antialiased=True, shade=True, alpha=0.9) # alpha = opacity of the surface
ax.set_xlim(-1,1)
ax.set_ylim(-1,1)
ax.set_zlim(-1,1)
ax._axis3don=False
ax.patch.set_facecolor('black')
ax.auto_scale_xyz([-1,1],[-1,1],[-1,1]) # true axis ratio
incl=90-i
ax.view_init(elev=incl,azim=60) # view angle, here azim is arbitrary
# the star is axisymmetric!
plt.axis('image')
srtnum='{0:03}'.format(i)
plt.show()
Any ideas???
THANKS!!!!
Dear all
Im trying to use the file viewstar.py at the postporcessing folder with one of my 2d models but all that I get if a black square. I can't find what is wrong. This is what im doing
rom mpl_toolkits.mplot3d import Axes3D
from matplotlib import cm
from matplotlib.colors import LogNorm
import matplotlib.pyplot as plt
import matplotlib.colors as pltc
import numpy as np
import string
from ester import *
a=star2d('MI_FILE2d')
nth=a.nth+2 # we take care that the plotted fields as read by "a" have
# nth+2 points in latitude
u = np.linspace(0, 2 * np.pi, 100) # generate angle phi grid
v = a.th[-1,:] # angle theta
w=v+np.pi/2
v=np.concatenate([w[0:nth-1],v]) # angle theta pi (South pole) to 0
r_sud=a.r[-1,:]
r_nord=a.r[-1,nth:0:-1]
r=np.concatenate([r_nord,r_sud])/a.Re
teff_rsh=np.reshape(a.Teff,nth)
teff_ext=np.concatenate([teff_rsh[nth:0:-1],teff_rsh])
x = np.outer(np.cos(u), rnp.sin(v))
y = np.outer(np.sin(u), rnp.sin(v))
z = np.outer(np.ones(np.size(u)), r*np.cos(v))
Teff = np.outer(np.ones(np.size(u)), teff_ext)
N=1-(Teff/Teff.max())**4 # so that white of "Blues" shows the max of
Teff
m=cm.ScalarMappable(cmap=cm.Blues, norm = pltc.Normalize(vmin=0,vmax=1.,clip=False))
m.set_array(N)
i=112.
print ('inclination angle i = ',i)
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.plot_surface(x, y, z, rstride=1, cstride=1, facecolors=m.to_rgba(N),linewidth=0, antialiased=True, shade=True, alpha=0.9) # alpha = opacity of the surface
ax.set_xlim(-1,1)
ax.set_ylim(-1,1)
ax.set_zlim(-1,1)
ax._axis3don=False
ax.patch.set_facecolor('black')
ax.auto_scale_xyz([-1,1],[-1,1],[-1,1]) # true axis ratio
incl=90-i
ax.view_init(elev=incl,azim=60) # view angle, here azim is arbitrary
# the star is axisymmetric!
plt.axis('image')
srtnum='{0:03}'.format(i)
plt.show()
Any ideas???
THANKS!!!!