ó
óQc           @   s  d  Z  d d l Z d e j d <d d l Td d l j Z d d l m Z d d l	 m
 Z
 m Z m Z m Z m Z m Z d d l m Z m Z d d	 l m Z d d
 l m Z d d l m Z d d l m Z d d l m Z m Z d d l m Z d d d d d d d d d d d d g Z d d d d d d d g Z  e! e ƒ Z" e! e  ƒ Z# d  Z$ d! g d" d  d# g d$ g d% d& g d' d( d) g d* g d+ g g Z% d! d  d% d& g Z& d" d# d$ g Z' e( d" g ƒ Z) e( d" d d! g d! d" d g d" d! d g d" d! d! g g ƒ Z* e( d" g ƒ Z+ e( d" d g d d" g d" d! g g ƒ Z, e( d" d d! g d! d" d g d" d! d g d! d" d! g g ƒ Z- e( d" g ƒ Z. e( d" g ƒ Z/ e e) e* e+ e, e- e. e/ f ƒ Z0 e0 j1 ƒ  Z0 d" d  d# d% d& d' d( d) g Z2 d" d  d# d& d' d) d* d+ g Z3 e4 d, d- d. ƒZ5 e6 e5 e3 <e0 d d … e2 f e3 d d … f Z7 e8 j e0 ƒ Z9 d/ „  Z: d0 „  Z; d1 „  Z< d2 „  Z= d3 „  Z> d4 „  Z? d5 „  Z@ d S(6   sŽ  
 Estimate electron temperature and ion concentration based on their the spectral lines intensity. 
 
 Ions intensity is descrimed by a very simple model with minimal number of parameters:
 Expectations:
  1) that line intesity is proportional to impurity density and fnction of the temperature
  2) ratio of lines of the different ionization states of the impuritie as 
    function of the temperature is polynome of order k (k = 1 now)
  3) shape of f(T_e) of the first ionozation stage is gaussian with width a and center T0
    (for k = 1  => all other are also gaussians)
  4) this model can not exactly estimate the temperature, only something which is 
    monotonous function of the T_E => caibration from the spitzer resistivity was necassary. 
  5) it is expected that during the shot is impurity density as constant as possible. 
  
  Advantages:
  simple model => robust
  most of the spectra are explained very well even by this simple model 

  
  
  Disadvantages:
  
   OI and NI can be observed in low temperatures a  deviation from model, f(T) if this
  ions don't have a gaussian shape, it radiates strongyl alte low temperatures and the the radiation
  is alomost constant over large range of temperatures. 

 In helium shots is probably a different physic or a week he lines mess estimate of the NI density. 
 NI radiation is significantly underestimated. 
 
 Autor :Tomas Odstrcil
 Date 12.12.2012


iÿÿÿÿNt   Aggt   backend(   t   *(   t   norm(   t   qrt   invt   pinvt   eigt   choleskyt   solve_triangular(   t   pearsonrt   hmean(   t
   block_diag(   t   fmin_powell(   t   repmat(   t   time(   t   save_advt   load_adv(   t	   saveconstt   HIt   OIt   OIIt   OIIIt   HeIt   CIIt   CIIIt   NIt   NIIt   NIIIt   M1t   MIVt   Ht   Ot   Het   Ct   Ni   i    i   i   i   i   i   i   i   i	   i
   i   i   t   dtypet   boolc            sŽ   t  d ƒ } | d ‰  | d ‰ | d ‰ | d ‰ ‡  ‡ ‡ ‡ f d †  } | |  ƒ } | |  | ƒ | } | |  | ƒ | } | | | f S(   Ns)   model_constants/nonlin_tranform_coeff.txti    i   i   i   c            s¸   t  |  ƒ }  t |  ˆ d k ƒ } ˆ d ˆ d d t |  | ˆ d ƒ |  | <t |  d k ƒ } t |  ˆ ˆ ƒ ˆ ˆ  } t |  | d ƒ t d ˆ ˆ ƒ ˆ ˆ  | | <| S(   Ngffffffî?gš™™™™™©?i   g333333@g333333À(   t   copyt   wheret   expt   tan(   t   Tt   indt   T_(   t   at   ct   bt   d(    sC   ../includes/analysis/Radiation/1212Impurities_TO.ON/CalcDensTemp.pyt   T2T_eVd   s    03(   t   loadtxt(   R*   t   Terrt   nl_coeffR1   R,   t   Terr_ut   Terr_d(    (   R-   R.   R0   R/   sC   ../includes/analysis/Radiation/1212Impurities_TO.ON/CalcDensTemp.pyt   CalibrateTemperature\   s    



c         C   s  | } t  |  d ƒ } t  |  d ƒ } t | ƒ } t t t d  d  … f ƒ } t | t |  d ƒ j ƒ } t  | d ƒ } t t | d | |  d d j ƒ ƒ } | t	 | d  d  … d f | ƒ 8} | d  d  … d f }	 t
 |	 | d f ƒ j }	 t | |	 d d d |	 | d ƒ}
 t |	 | | d d d ƒt |	 | d d d ƒ}
 |	 |
 t d  d  … f 9}	 t | |	 | ƒ d t  | ƒ t  |
 ƒ } d G| GH|
 S(	   Ni    i   g-Cëâ6?i   gñhãˆµøä>t   axist   weightss
   chi2 temp:(   t   sizet   onest   arrayt   Mt   bound_inter_jt   dott   logR*   t   sqrtt   outert   tilet   averaget   sumt   newaxisR   (   t   pt   pEt   coefft   nt   mt   It   cutMt   Rt   Ret   BR*   t   chi2(    (    sC   ../includes/analysis/Radiation/1212Impurities_TO.ON/CalcDensTemp.pyt	   EstimateTy   s"    (#$6,	c
         C   s“  t  t d f ƒ }
 t d d d ƒ } t t  d ƒ | f ƒ j } t | |
 | |  | | | t |	 ƒ t } t | | | |  | | | t |	 ƒ } | d  d  … t |	 f | | d  d  … t |	 f } t j	 d d d ƒ } t
 t | d  d  … t |	 f ƒ t | ƒ ƒ \ } } t | ƒ d t | ƒ } t t t | ƒ ƒ ƒ } | | :} | d C} | t | t d | d | ƒ ƒ t | ƒ } d d d d	 g } x*t t t |	 ƒ ƒ D]} | } t | d  d  … t |	 | f t | |	 d  d  … f ƒ ƒ } t | d  d  … t |	 | f t | |	 d  d  … f ƒ ƒ } | | d k } t | | | | | | d
 | | d d d d d d d ƒt j | | d  d  … | f | | d t t |	 | ƒq¬Wt j d ƒ t j d ƒ t j d d d t ƒ } | j ƒ  j d ƒ t j  t! |	 ƒ | j" d d d | | | f d d d d ƒt j# d d ƒ t j$ d  d ƒ t j% d! t! |	 d" t& d# ƒ t' j( ƒ  d  S($   Nid   iýÿÿÿi   i   i   R/   t   rt   kt   yt   fmtt   .t   capsizei    t	   linewidthgš™™™™™¹?t
   markersizeg      à?t   labels   f(T) [-]s   intensity [a.u.]t   locs
   upper leftt   fancyboxgffffffæ?gffffffþ¿s7   $\chi^2$ = %2.1f 
$\chi^2$ robust = %2.1f 
$r$ =  %1.2ft   horizontalalignmentt   leftt   verticalalignmentt   bottomg      ø?iþÿÿÿs   grafy_advance/_t    s   .png()   R;   t	   n_elemetst   linspacet   vstackR*   t   CalcIntensityt
   ions_indext   et   pltt   subplotR
   t   ravelR   R:   t   doublet   amaxt   absRE   RA   t   ranget   lent   squeezeR(   t   errorbart   plott	   ion_namest   xlabelt   ylabelt   legendt   Truet	   get_framet	   set_alphat   titlet   element_namest   textt   ylimt   xlimt   savefigt   namet   figt   clf(   t   T0R-   t   densityR#   R*   RP   RI   RG   RH   t   it   density_modelt   T_modelt   N_modelt   c2RL   t   Residt   axt   r_tott   _RQ   t	   max_Residt   chi2_robustt   colourRT   t   xRU   t   y_errt   noisyt   leg(    (    sC   ../includes/analysis/Radiation/1212Impurities_TO.ON/CalcDensTemp.pyt   plotTprofiles•   s@    )%65

0>>?>!c         C   s  d | | d k <t  | d ƒ } t d | f d d ƒ}	 t | | j ƒ |	 t d  d  … f <t |  t d ƒ j }  |  | 8}  |  d C}  |  | 9}  |  j }  |  | | d  d  … t f 7}  |  |	 t d  d  … f <t t	 | d  d  … f |	 d ƒ j }
 t
 |
 d t ƒS(   Nid   i   i   R$   Rl   i   R&   (   R:   t   emptyR?   R*   t   bound_interR   Rc   RF   R(   t   iMR<   t   False(   R*   R…   R#   R„   R-   RP   RI   t   indexRJ   t   FRL   (    (    sC   ../includes/analysis/Radiation/1212Impurities_TO.ON/CalcDensTemp.pyRf   Â   s    "


	'c      
   C   sk  t  |  d ƒ } t g  t t ƒ D] } | | ^ q ƒ j }	 t t | f ƒ }
 t | |
 |	 | | | | t d ƒ ƒ } d | | d k  <t |  | d ƒ } | |  d } xb t	 t
 ƒ D]T \ } } t | d  d  … | f d d d | d  d  … | f ƒ|
 | d  d  … f <q³ Wt | |
 |	 | | | | t d ƒ ƒ } t |  | | ƒ d t  |  ƒ t  |
 ƒ } d G| GH|
 S(	   Ni    i   gü©ñÒMbP?R8   i   R9   i   s
   chi2 dens:(   R:   Re   Ro   RT   R*   t   zerosRc   Rf   R@   t	   enumerateRg   RD   R   (   RG   RH   R*   R„   R-   RP   RI   RJ   t   jR#   R…   RL   RS   t   reR†   R+   RQ   (    (    sC   ../includes/analysis/Radiation/1212Impurities_TO.ON/CalcDensTemp.pyt   EstimateDensityÖ   s    ,'L',	c    3         s
  t  d ƒ ‰ t  d ƒ ‰ t  d ƒ ‰	 t  d ƒ ‰ t d ƒ \ }  } | ‰  | j ‰ t ˆ  d  d  … t f d d ƒd k } t d	 d t ˆ  d	 ƒ d
 t ƒ} ˆ  | d  d  … f ‰  ˆ | d  d  … f ‰ |  | }  | | } t ˆ  d	 ƒ ‰ ˆ d	 k rù d  St ˆ d t	 f ƒ } t
 t ˆ  d d	 ƒƒ | d d … d  d  … f <t ˆ ƒ } xf t t	 ƒ D]X } t ˆ ƒ d } t ˆ d ƒ } t | | | d  d  … | f ƒ | d  d  … | f <qVWˆ d | 7‰ d ˆ ˆ d k <t ˆ t k ƒ ‰ ‡ ‡ ‡ ‡ ‡ ‡ ‡	 ‡ ‡  f	 d †  }	 t ˆ  ˆ ˆ ƒ }
 t ˆ  ˆ |
 ˆ	 ˆ ˆ ˆ ƒ } t |
 t | j ƒ f ƒ } d ‰ d GHt |	 | d d d d d d ƒ} d GH| ˆ  } t g  t t ƒ D] } | | ^ qŸƒ j ‰ t | ˆ t ˆ t f ƒ j } t | | ˆ ˆ	 ˆ ˆ ˆ t d ƒ ƒ } ‡ ‡ ‡ ‡	 ‡ ‡ f d †  } ‡  f d †  } | | | ƒ } t d d d d d d d d d d d d g ƒ } t ˆ ƒ } | t | ƒ 9} t | ƒ } | | :} t | | j ƒ } | j  d  d  t | d	 ƒ d … c d 7<t! | d t d t" ƒ} t# | j$ d	 ƒ } t% | | d d	 d t" d  t" ƒ} t& | d! | ƒ} | d" C} t t | d	 ƒ ƒ } t' d# ƒ } ˆ  | } t | d ƒ } | ˆ  } t | ˆ ˆ t f ƒ j } t( | | ƒ \ }  }! }" t) | ƒ d$ }# t* | d  d  … |# f d d d% | d  d  … |# f ƒ}$ t+ |$ ƒ }$ x= t, |$ t- ƒ D], \ }% }& d& |& |% f GHt. d' |& |% ƒ q Wt* |  |# d% d |" |# d" ƒ}' d( |' GHt. d) |' ƒ t t ƒ } x× t t ƒ D]É } t | | ˆ ˆ	 ˆ ˆ ˆ t/ | ƒ } ˆ  d  d  … t/ | f | ˆ d  d  … t/ | f }( t0 t1 t
 |( ƒ ƒ ƒ }) |( |) :}( |( d" C}( |) t |( t d |) d" |( ƒ ƒ t |( ƒ ˆ d | | <q¡Wt2 d* |  | d+ | d, d ƒt2 d- |  |  d+ |" |! g d, d ƒt2 d. t t ƒ | ƒ y6 t' d/ ƒ }* t3 t
 t4 |* ƒ t5 ƒ d k  ƒ d	 }+ Wn n Xt6 j7 d0 d  d1 dT d4 d5 d6 d7 d8 d9 ƒ }, |, j9 d: ƒ }- t | | ˆ ˆ	 ˆ ˆ ˆ t d ƒ ƒ } t6 j: t t ˆ  ƒ ƒ t ˆ  j ƒ d; t ˆ j ƒ ƒt6 j; t | j ƒ ƒ t< ˆ  | ˆ ƒ d" t ˆ  t | ƒ ƒ }. |- j= d< ƒ t ˆ  ˆ ˆ ƒ }/ t ˆ  ˆ |/ ˆ	 ˆ ˆ ˆ ƒ }0 t |/ |0 ˆ ˆ	 ˆ ˆ ˆ t d ƒ ƒ } t6 j; t | j ƒ d= ƒ |- j> d> ƒ |- j? d? ƒ |- j@ d	 t1 ˆ  ƒ d@ ƒ |- jA d" d" t1 ˆ  ƒ d dA |. dB dC dD dE ƒxK tB tC ƒ D]= \ } }& |- jA | ˆ d t1 ˆ  ƒ dF |& dB dC dD dE ƒqÇW|, j9 dG ƒ }- t6 j: |  dH | | dI d ƒt6 j; |  dH |/ ˆ  dJ ƒ y0 t' dK ƒ }1 t6 j; |* |+ t5 |1 |+ dL ƒ Wn n X|- j@ dM d ƒ |- j? dN ƒ |, j9 dO ƒ }- t6 jD |  dH t+ | j ƒ ƒ y^ t' dP ƒ }2 |2 ˆ d  d  … tE f 7}2 t6 jD |* |+ t5 t+ |2 d  d  … |+ f j ƒ dL ƒ Wn n X|- j@ d d ƒ |- j? dQ ƒ |- jF dR ƒ |, jG dS ƒ |, jH ƒ  d  S(U   Ns   model_constants/C.txts   model_constants/A.txts   model_constants/T0.txts   model_constants/B.txts   ./data/projectionR8   i   g{®Gáz„?i    t   endpointiÿÿÿÿg      à?gš™™™™™¹?gü©ñÒMbP?c      
      sì  t  |  ƒ }  |  ˆ  } t |  ˆ t ˆ t f ƒ j } t g  t t ƒ D] } | | ^ qG ƒ j } t | | | ˆ ˆ ˆ ˆ t d ƒ ƒ } ˆ | ˆ d } t  t t	 | ƒ ƒ d ƒ } | | :} | d C} | t
 | t d | d | ƒ ƒ t | ƒ t |  ƒ ˆ  } t | | d k ƒ }	 |	 t | | d k  ƒ 7}	 t | d d d ƒ}
 |
 d c d :<|
 d c d :<t |
 ƒ t |
 ƒ }
 t | d ƒ } d | | d k <t | ƒ } t t t | ƒ ƒ ƒ } | | d	 ˆ |
 d
 | d	 |	 S(   Ni   gñhãˆµøä>i   i   i   iýÿÿÿR8   i    i
   g333333ã?(   Rl   t   reshapeRc   R*   Re   Ro   RT   Rf   Rm   Rn   RE   RA   R:   R   t   diffRp   Rk   R(   (   R’   R*   R…   RŸ   R#   RL   R‹   R   RQ   R/   R0   t   Tdt   spars(	   t   n_inftyRH   t   At   lamR"   RP   R„   RJ   RG   (    sC   ../includes/analysis/Radiation/1212Impurities_TO.ON/CalcDensTemp.pyt   f_optim&  s*    
!,'

>i
   R   t   xtolgñhãˆµøä>t   ftolt   maxfung    €„.As   fmin_powell finishi   c      
      sW   |  ˆ  } t  |  ˆ ˆ t f ƒ j } t | | ˆ ˆ ˆ  ˆ ˆ t d ƒ ƒ } t | ƒ S(   Ni   (   R£   Rc   R*   Rf   Ro   Rk   (   R’   t   Tempt   DensRL   (   R¨   R"   RP   R„   RJ   R#   (    sC   ../includes/analysis/Radiation/1212Impurities_TO.ON/CalcDensTemp.pyt   RetrofitFunctionk  s    
'c            s·   t  ˆ  ƒ } d } t t |  ƒ | f d d ƒ} x€ t t |  ƒ ƒ D]l } t |  ƒ } | | c | 8<t |  ƒ } | | c | 7<| | ƒ | | ƒ d | | | d  d  … f <qC W| S(   Ngíµ ÷Æ°>R$   t   float32i   (   R:   R   Rp   Ro   R&   (   t   x0t   funt   n2t   dxt   JR†   t   xlt   xr(   RG   (    sC   ../includes/analysis/Radiation/1212Impurities_TO.ON/CalcDensTemp.pyt
   CalcJacobiq  s    2i   i   t   lowert   overwrite_at   transt   overwrite_bt   outi   s   ./data/energy_constants.npyiþÿÿÿR9   s   ion %2s concentration = %1.2fs   ./HistoricalAnalysis/s   flattop temperatute %1.1fs    ./HistoricalAnalysis/temperatures   ./data/densityt   data_errt   tvec_errs   ./data/temperatures   ./data/TempDensChi2s	   shots.npyt   numt   figsizei   i   t   dpiiP   t	   facecolort   wt	   edgecolorRT   i7  t   yerrs   Spectral data retrofits   r:t   tightt	   intensityg333333ó?s   $\chi^2$/doF = %.2fR^   R_   R`   Ra   gÍÌÌÌÌÌð?i8  g     @@RZ   t   :s   temperature.npys   --iüÿÿÿs   pseudo temperaturei9  s   density.npyt   Densitys   t [ms]s   ./graphs/proj_retrofit.png(   i   i   (I   R2   R   R¿   RE   t   bound_inter_iRd   R:   Rš   R   t   n_ionsRn   R¤   t
   zeros_likeRo   t   aranget   interpt   inftyRR   R¡   t   hstackRk   R*   R   Re   RT   R£   Rc   Rf   R<   R&   RA   R?   t   flatR   Rx   t   identityt   shapeR	   t   absolutet   loadR7   t   argsortRD   R(   t   zipR|   R   Rg   Rl   Rm   R   R'   t   floort   shot_numberRi   t   figuret   Nonet   add_subplotRr   Rs   R   t	   set_titleR8   t
   set_ylabelt   set_ylimR}   Rž   Rt   t   semilogyRF   t
   set_xlabelR€   Rƒ   (3   t   tvect	   proj_dictt   st
   shots_partt   dpt   dp_R†   R’   t   xpRª   t   temp0t   dens0R²   RU   R*   RŸ   R…   RL   R°   R¹   R¶   RQ   t   correctedPEt   errt   JJt   Lt   iLt   energy_constantst   p_calibt   P_totalt   temp_errt   dens_errt   T_calibt   T_calib_err_ut   T_calib_err_dR+   t   max_densR0   R   t   flt_tempR‹   R   t   shotst
   shot_indexR‚   RŒ   t   chi2doFt   T_estt   d_estt   temperature_fullt   density_full(    (
   RG   R"   R§   RH   RP   R#   RJ   R¨   R©   R„   sC   ../includes/analysis/Radiation/1212Impurities_TO.ON/CalcDensTemp.pyt   CalcDensTempñ   sü    	
(!

.6'%!
,!'0
)!


8"	%6

@"**'4,'15$8
c          C   sg  t  d ƒ \ }  } | j } t  d ƒ \ }  } | j } | d  d  … d f | d  d  … d f g } t  d ƒ \ } } d t j f d „  ƒ  Y} t j d d  d	 dA d d d d d d ƒ } t j d d d d ƒ xèt d t ƒ D]×}	 | j	 t d |	 d ƒ }
 |	 d k r|
 j
 d ƒ n  |
 j j t j ƒ  ƒ |
 j j | ƒ  ƒ t | |	 d  d  … f ƒ } t | |	 d  d  … f ƒ d | } t | |	 d  d  … f ƒ d | } t | d t | ƒ ƒ} t t | | ƒ d ƒ t t | ƒ ƒ d } t j |  d | d | | g d d d t |	 ƒ|
 j j t ƒ |
 j d d d | |	 | | f d d  d! d" d# |
 j ƒ|
 j d t d$ t | ƒ ƒ d% ƒ t j d& d' d( t ƒ } | j ƒ  j d ƒ qÝ W|
 j j t j ƒ  ƒ |
 j d) ƒ | j  d* d+ d, ƒ| j! ƒ  t j d d  d	 dB d d d d d d ƒ } | j	 d. ƒ }
 |
 j
 d/ ƒ t j |  d | d | d d ƒ|
 j" d0 ƒ |
 j d) ƒ |
 j d d ƒ |
 j j t ƒ | j  d1 d+ d, ƒt j# ƒ  yŸt j d d  d	 dC d d d d d d ƒ } | j	 d. ƒ }
 t$ d2 t% t& ƒ d3 ƒ } | d4 | d5 } t' | d6 | d7 t( | ƒ ƒ } y8 t) d2 t% t& ƒ d8 ƒ } t) d2 t% t& ƒ d9 ƒ } Wn n Xt& d: k r˜| d 9} n( | d :} | d 8} | d; :} | d; :} |
 j
 d< ƒ t j* | d | ƒ |
 j" d0 ƒ |
 j d) ƒ |
 j d d ƒ |
 j+ | d; | d; ƒ | j  d= t% t& ƒ d> d+ d, ƒt j# ƒ  d? GHWn d@ GHn Xd  S(D   Ns   ./data/densitys   ./data/temperaturei    i   s   ./data/TempDensChi2t   MyFormatterc           B   s   e  Z d d  „ Z RS(   c         S   s*   | d k r d St  j j |  | | ƒ Sd  S(   Ni    t    (   Ri   t   ScalarFormattert   __call__(   t   selfR’   t   pos(    (    sC   ../includes/analysis/Radiation/1212Impurities_TO.ON/CalcDensTemp.pyR    s    N(   t   __name__t
   __module__RÝ   R  (    (    (    sC   ../includes/analysis/Radiation/1212Impurities_TO.ON/CalcDensTemp.pyR    s   RÁ   RÂ   i
   RÃ   iP   RÄ   RÅ   RÆ   RT   t   hspacet   wspaces   Relative ions concentrationR9   gü©ñÒMbP?i   g     @@RÇ   RZ   g      à?R[   gš™™™™™©?gš™™™™™É?s*   $\chi^2/doF = %2.1f$ 
 $n = %1.2f\pm%1.1g$R^   R_   R`   Ra   t	   transformgš™™™™™¹?g      ø?R\   s
   upper leftR]   s   t [ms]s   ./graphs/density.pngt   bbox_inchesRÈ   i   io   s   Estimated temperatures
   $T_e$ [eV]s   ./graphs/temperature.pngs   ./electronTemperatures/s   /electron_temperature.npzt   datat   scalet   t_startt   t_ends   /PlasmaStarts
   /PlasmaEndiT$  iè  s   Spitzer temperatures   ./graphs/temperatures
   _spitz.pngs
   Te plotteds   no electron temperature(   i
   i
   (   i
   i   (   i
   i   (,   R   R¿   Ri   R  RÜ   RÝ   t   subplots_adjustRo   Rc   RÞ   Rß   t   xaxist   set_major_formattert   NullFormattert   yaxisR(   RD   Rn   R   Rr   R|   t   gridRx   R}   t	   transAxesRá   t   maxRm   Rw   Ry   Rz   Rã   R€   Rƒ   Rà   t   closeR×   t   strRÛ   Rd   Rp   R2   Rs   t   set_xlim(   Rä   R…   Rö   R*   Rõ   RŽ   Rþ   R  R‚   R†   RŒ   RU   t   y_err_upt
   y_err_downt   y_medt	   y_err_medR•   R  t   tempt
   tvec_spitzt   plasma_startt
   plasma_end(    (    sC   ../includes/analysis/Radiation/1212Impurities_TO.ON/CalcDensTemp.pyt   plotDensTemp  sŠ    		,	*$&.09#
* 
* 



!
	(A   t   __doc__t
   matplotlibt   rcParamst   numpyt   matplotlib.pyplott   pyplotRi   t   numpy.linalgR   t   scipy.linalgR   R   R   R   R   R	   t   scipy.stats.statsR
   R   t   scipy.sparseR   t   scipy.optimizeR   t   numpy.matlibR   R   t   pygolem_lite.modulesR   R   t   pygolem_liteR   Rt   R|   Rp   RÍ   Rc   RT   Rg   t
   single_iont	   multi_iont   matrixR   t   M2t   M3t   M4t   M5t   M6t   M7R=   t   todenseRÌ   R>   R   R˜   Rx   t	   M_limitedt   linalgR™   R7   RR   R–   Rf   R¡   R  R'  (    (    (    sC   ../includes/analysis/Radiation/1212Impurities_TO.ON/CalcDensTemp.pyt   <module>%   sX   
.*?<'<!
&			-			ÿ 