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

Пример на основе https://arxiv.org/pdf/1611.07672.pdf галактики ESO534-G001 и 2dFRS S833Z022


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>