Пример работы с двухжидкостной неустойчивостью для карт

In [1]:

    

from IPython.display import HTML
from IPython.display import Image
import os
from PIL import Image as ImagePIL

%pylab
%matplotlib inline
%run ../../utils/load_notebook.py


    

    




Using matplotlib backend: Qt5Agg
Populating the interactive namespace from numpy and matplotlib

In [2]:

    

from photometry import *


    

    




importing Jupyter notebook from photometry.ipynb
Using matplotlib backend: Qt5Agg
Populating the interactive namespace from numpy and matplotlib

In [3]:

    

from instabilities import *


    

    




importing Jupyter notebook from instabilities.ipynb
Using matplotlib backend: Qt5Agg
Populating the interactive namespace from numpy and matplotlib

In [4]:

    

from utils import *


    

    




importing Jupyter notebook from utils.ipynb

In [5]:

    

name = 'ESO534-G001'
gtype = 'SBm' #LEDA
incl = 74.0  #from arxiv
distance = 20.0 #Mpc, arxiv
scale = 0.097 #kpc/arcsec according to Noordermeer distance

data_path = '../../data/demo_data'
sin_i, cos_i = np.sin(incl*np.pi/180.), np.cos(incl*np.pi/180.)
print sin_i, cos_i


    

    




0.961261695938 0.275637355817

In [6]:

    

%%javascript 
$.getScript('https://kmahelona.github.io/ipython_notebook_goodies/ipython_notebook_toc.js')


    

In [7]:

    




    

In [8]:

    

basewidth, baseheight = 814, 814


    

In [ ]:

    




    

In [9]:

    

Image('full_HI_Sigma.png', width=700)


    

    
Out[9]:

In [10]:

    

Sigma_raw = ImagePIL.open('HI_Sigma.png')


    

In [11]:

    

plt.figure(figsize=[8, 8])
plt.imshow(Sigma_raw)


    

    
Out[11]:





<matplotlib.image.AxesImage at 0xcd6bef0>