##!/usr/bin/python2
## -*- coding: utf-8 -*-

<%page args="shot, page" />


<%!
    from pygolem_lite.modules import  cData
    from pygolem_lite  import Shot
    import random
    rand = random.randint(1,1e6)
    get = cData.get_nice
    
%>


<h3> <a href="../0411ShotHomepage.ONN/basic_diagn.py"> [Source - computations] </a> <a href="main.py"> [Source - plotting] </a> </h3>

%if cData['plasma']:

 <a name="safety_factor"><img src=plot_1.png?${rand}></a> 
 <a name="electron_temperature"><img src=plot_2.png?${rand}></a> 
 <a name="power_balance"><img src=plot_3.png?${rand}></a> 
 <a name="density"><img src=plot_4.png?${rand}></a> 

  
  <h4> Plasma position is not included in the safety factor !!! </h4>
  
  <h3>Other</h2>

First approximation of electron confinement time  t<sub>e</sub> = ${get('electron_confinement_time', "%.3f", 1e6) } [us] <br>
Electron confinement time predicted from t_98  t<sub>98</sub> = ${get('electron_confinement_t98', "%.3f", 1e6) } [us] <br>
(t<sub>IPB98</sub>(y,2) = 0.0562 I<sup>0.93</sup> B <sup>0.15</sup> n <sup>0.41</sup> P <sup>-0.69</sup> R <sup>1.39</sup> ka <sup>0.78</sup> a <sup>0.58</sup> M <sup>0.19</sup>) <br>
Approximation of electron density through equilibrium state  n<sub>e</sub> = ${get('electron_density_equilibrium', "%.2f" , 1e-19) } [10^19.m^-3] <br>
Mean electron density n<sub>e</sub> = ${get('electron_density_mean', "%.2f" , 1e-19) } [10^19.m^-3] <br>

% else:
    <img src=magnetic_flux.png?${rand}>
    <img src=hysteresis.png?${rand}>
% endif

${Shot(shot).get_pygolem_list(['total_magnetic_flux', 'safety_factor', 'electron_confinement_time', 'electron_density_equilibrium', 'greenwald_density', 'input_power_total','input_power_magnetic','input_power_plasma','input_power_chamber','electron_temperature', 'electron_temperature_max', 'temperature_mean', 'electron_density', 'electron_density_mean' , 'electron_confinement_t98', 'spectrometr:temperature', 'spectrometr:total_radiated_power' ])}