1. 
   
   ok

2. boucle sur chaque sous ensemble : select='gas', disk=('bulge','gas','disk')
   
   ok

3. 
   
   a continuer

gextractall.py gplotall.py gplotvelocities.py gplotstability.py gplotsdensity.py --> a faire

./gextractall.py -t gadget --cmd='''nb.histocenter()''' --nmin=2 --nR=32 --Rmax=50 --eps=0.25 snap/Albanese0000.dat ./gextractall.py -t gadget --cmd='''nb.histocenter()''' --nmin=2 --nR=32 --Rmax=50 --eps=0.25 snap/Albanese0000.dat --components="('gas',)" --y Sden --mode=Sden ./gextractall.py -t gadget --cmd='''nb.histocenter()''' --nmin=2 --nR=32 --Rmax=50 --eps=0.25 snap/Albanese0000.dat --components="('gas',)" --y Q --mode=stability

./gplotall.py gas.h5 halo.h5 --x R --y vc --mode=sum2 ./gplotall.py gas.h5 halo.h5 --x R --y Sden --mode=sum --log y

./gextractall.py -t gadget --eps 0.1 --Rmax=20 --nR=64 ~/pgm/python/condinit/miyamoto.dat ./gplotvelocities.py gas.h5

./gplotstability.py gas.h5

---> legends ---> labels ---> choix des couleurs

options pour les figures...

pt.figure(figsize=(8*2,6*2)) taille des points... palette...

noplot

• calculer :
  
  vc, Phi0
  
  
  
  Sden,Sdend
  
  vm,sr,sz,st
  
  Q,X,A
  
  -> affichage -> sortie en txt ou h5

• vitesse de rotation par composant
  
  composants = ("gas",'disk','halo',...)