Handling/DataAccess/SWtools/wget/10104/shot/analysis/Basics/0612AdvancedAnalysis.ON/main.py














<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN"
  "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">

<html xmlns="http://www.w3.org/1999/xhtml">
  <head>
    <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />

    <title>Golem No: 10104 - the only fully opensource tokamak </title>
    <link rel="shortcut icon" href="/favicon.ico" type="image/x-icon" />
    <link rel="stylesheet" href="/_static/default.css" type="text/css" />
    <link rel="stylesheet" href="/_static/pygments.css" type="text/css" />
    <link rel="stylesheet" href="/_static/lightbox.css" type="text/css" media="screen" />

    <script type="text/javascript">
      var DOCUMENTATION_OPTIONS = {
	URL_ROOT:    '',
	VERSION:     '0.2',
	COLLAPSE_INDEX: false,
	FILE_SUFFIX: '.html',
	HAS_SOURCE:  true
      };
    </script>
    <script type="text/javascript" src="/_static/jquery.js"></script>
    <script type="text/javascript" src="/_static/underscore.js"></script>
    <script type="text/javascript" src="/_static/doctools.js"></script>
    <script type="text/javascript" src="/_static/sidebar.js"></script>

    <script src="/_static/lightbox_plus.js"></script>

    <script type="text/javascript">
	function getURL(val){
	location = '/shots/' + val + '/' +  "analysis/Basics/0612AdvancedAnalysis.ON/";
	return false;
	}
    </script>

  </head>
  <body>
    <div class="related">
      <h3>Navigation</h3>

	<ul>
	<li class="right" ><a href= /shots/0/analysis/Basics/0612AdvancedAnalysis.ON/ title="current shot" >current</a> </li>
	<li class="right" ><a href= /shots/10105/analysis/Basics/0612AdvancedAnalysis.ON/ title="next shot" >next</a> |</li>
	<li class="right" ><a href= /shots/10103/analysis/Basics/0612AdvancedAnalysis.ON/ title="previous shot" >previous</a> | </li>
	</ul>
	<ul>
	<li><a href="/" >  GOLEM  </a> &raquo;
	<a href= /shots/10104 >Shot  &#35;10104</a>  &raquo;   </li>
	    <li><a href=/shots/10104/analysis/Basics/0612AdvancedAnalysis.ON/ >  main  </a>	    &raquo;
	</ul>
    </div>

    <div class="document">
      <div class="documentwrapper">
	<div class="bodywrapper">
	  <div class="body">
	    <div class="section" >

		<h1>Source code :: main</a></h1>


		<pre><font face="Lucida,Courier New"><a class="viewcode-back" href="javascript:history.go(-1)">[Return]</a><br/><a class="viewcode-back" href="main.py">[Download]</a><font color="#008000">#!/usr/bin/python2</font>
<font color="#008000"># -*- coding: utf-8 -*-</font>


<font color="#004080">""" CREATED: 7/2012
    AUTHOR: MICHAL ODSTRČIL
"""</font>


<font color="#C00000">import</font> <font color="#000000">matplotlib</font>
<font color="#000000">matplotlib</font><font color="#0000C0">.</font><font color="#000000">rcParams</font><font color="#0000C0">[</font><font color="#004080">'backend'</font><font color="#0000C0">]</font> <font color="#0000C0">=</font> <font color="#004080">'Agg'</font>
<font color="#000000">matplotlib</font><font color="#0000C0">.</font><font color="#000000">rc</font><font color="#0000C0">(</font><font color="#004080">'font'</font><font color="#0000C0">,</font>  <font color="#000000">size</font><font color="#0000C0">=</font><font color="#004080">'10'</font><font color="#0000C0">)</font>
<font color="#000000">matplotlib</font><font color="#0000C0">.</font><font color="#000000">rc</font><font color="#0000C0">(</font><font color="#004080">'text'</font><font color="#0000C0">,</font> <font color="#000000">usetex</font><font color="#0000C0">=</font><font color="#000000">True</font><font color="#0000C0">)</font>  <font color="#008000"># FIXME !! nicer but slower !!!</font>

<font color="#C00000">import</font> <font color="#000000">pygolem_lite</font>

<font color="#C00000">from</font> <font color="#000000">pygolem_lite</font> <font color="#C00000">import</font> <font color="#000000">Shot</font>


<font color="#C00000">from</font> <font color="#000000">numpy</font> <font color="#C00000">import</font> <font color="#0000C0">*</font>
<font color="#C00000">from</font> <font color="#000000">pygolem_lite</font><font color="#0000C0">.</font><font color="#000000">config</font> <font color="#C00000">import</font> <font color="#0000C0">*</font>
<font color="#C00000">from</font> <font color="#000000">pygolem_lite</font><font color="#0000C0">.</font><font color="#000000">modules</font> <font color="#C00000">import</font> <font color="#0000C0">*</font>
<font color="#C00000">from</font> <font color="#000000">matplotlib</font><font color="#0000C0">.</font><font color="#000000">pyplot</font> <font color="#C00000">import</font> <font color="#0000C0">*</font>
<font color="#C00000">import</font> <font color="#000000">time</font>
<font color="#C00000">from</font> <font color="#000000">shutil</font> <font color="#C00000">import</font> <font color="#000000">copy</font><font color="#0000C0">,</font> <font color="#000000">move</font>
<font color="#C00000">import</font> <font color="#000000">os</font>
<font color="#C00000">from</font> <font color="#000000">scipy</font><font color="#0000C0">.</font><font color="#000000">interpolate</font> <font color="#C00000">import</font> <font color="#000000">interp1d</font>


<font color="#C00000">def</font> <font color="#000000">analysis</font><font color="#0000C0">(</font><font color="#0000C0">)</font><font color="#0000C0">:</font>
  <font color="#C00000">try</font><font color="#0000C0">:</font>
    <font color="#000000">Data</font> <font color="#0000C0">=</font> <font color="#000000">Shot</font><font color="#0000C0">(</font><font color="#0000C0">)</font>


    <font color="#000000">I</font> <font color="#0000C0">=</font> <font color="#000000">Data</font><font color="#0000C0">[</font><font color="#004080">'plasma_current_mean'</font><font color="#0000C0">]</font> <font color="#0000C0">*</font> <font color="#0080C0">1e-6</font>
    <font color="#000000">Bt</font> <font color="#0000C0">=</font> <font color="#000000">Data</font><font color="#0000C0">[</font><font color="#004080">'toroidal_field_mean'</font><font color="#0000C0">]</font>
    <font color="#000000">P</font> <font color="#0000C0">=</font> <font color="#000000">Data</font><font color="#0000C0">[</font><font color="#004080">'input_power_mean'</font><font color="#0000C0">]</font><font color="#0000C0">*</font><font color="#0080C0">1e-3</font>
    <font color="#000000">n</font> <font color="#0000C0">=</font> <font color="#000000">Data</font><font color="#0000C0">[</font><font color="#004080">'electron_density_mean'</font><font color="#0000C0">]</font> <font color="#0000C0">*</font> <font color="#0080C0">1e-19</font>
    <font color="#000000">M</font>  <font color="#0000C0">=</font> <font color="#0080C0">1</font> <font color="#008000"># hydrogen !!</font>
    <font color="#000000">R</font> <font color="#0000C0">=</font> <font color="#000000">MajorRadius</font>
    <font color="#000000">eps</font> <font color="#0000C0">=</font> <font color="#000000">Aspect</font>
    <font color="#000000">kappa</font> <font color="#0000C0">=</font> <font color="#0080C0">1</font>
    <font color="#000000">const</font> <font color="#0000C0">=</font> <font color="#0080C0">5.62e-2</font>
    <font color="#C00000">print</font> <font color="#004080">"I"</font><font color="#0000C0">,</font><font color="#000000">I</font>
    <font color="#C00000">print</font> <font color="#004080">"B"</font><font color="#0000C0">,</font> <font color="#000000">Bt</font>
    <font color="#C00000">print</font> <font color="#004080">"P"</font><font color="#0000C0">,</font> <font color="#000000">P</font>
    <font color="#C00000">print</font> <font color="#004080">"n"</font><font color="#0000C0">,</font> <font color="#000000">n</font>
    <font color="#C00000">print</font> <font color="#004080">"R"</font><font color="#0000C0">,</font> <font color="#000000">R</font>
    <font color="#C00000">print</font> <font color="#004080">"Aspect"</font><font color="#0000C0">,</font> <font color="#000000">Aspect</font>

    <font color="#000000">t_98</font> <font color="#0000C0">=</font> <font color="#000000">const</font><font color="#0000C0">*</font><font color="#000000">I</font><font color="#0000C0">**</font><font color="#0080C0">0.93</font><font color="#0000C0">*</font><font color="#000000">Bt</font><font color="#0000C0">**</font><font color="#0080C0">0.15</font><font color="#0000C0">*</font><font color="#000000">P</font><font color="#0000C0">**</font><font color="#0000C0">-</font><font color="#0080C0">0.69</font><font color="#0000C0">*</font><font color="#000000">n</font><font color="#0000C0">**</font><font color="#0080C0">0.41</font><font color="#0000C0">*</font><font color="#000000">M</font><font color="#0000C0">**</font><font color="#0080C0">0.18</font><font color="#0000C0">*</font><font color="#000000">R</font><font color="#0000C0">**</font><font color="#0080C0">1.97</font><font color="#0000C0">*</font><font color="#000000">eps</font><font color="#0000C0">**</font><font color="#0080C0">0.58</font><font color="#0000C0">*</font><font color="#000000">kappa</font><font color="#0000C0">**</font><font color="#0080C0">0.78</font>

    <font color="#C00000">print</font> <font color="#004080">"confinent time = "</font><font color="#0000C0">,</font> <font color="#000000">t_98</font> <font color="#0000C0">*</font> <font color="#0080C0">1e6</font>

    <font color="#000000">saveconst</font><font color="#0000C0">(</font><font color="#004080">'t_98'</font><font color="#0000C0">,</font> <font color="#000000">t_98</font><font color="#0000C0">)</font>
  <font color="#C00000">except</font> <font color="#000000">Exception</font><font color="#0000C0">,</font> <font color="#000000">e</font><font color="#0000C0">:</font>
    <font color="#C00000">print</font> <font color="#004080">" confinent time failed "</font><font color="#0000C0">,</font> <font color="#000000">e</font><font color="#0000C0">.</font><font color="#000000">__str__</font><font color="#0000C0">(</font><font color="#0000C0">)</font>


<font color="#C00000">def</font> <font color="#000000">plot_data</font><font color="#0000C0">(</font><font color="#000000">file_type</font><font color="#0000C0">)</font><font color="#0000C0">:</font>


    <font color="#000000">data</font> <font color="#0000C0">=</font> <font color="#000000">Shot</font><font color="#0000C0">(</font><font color="#0000C0">)</font>
    <font color="#000000">plasma</font> <font color="#0000C0">=</font> <font color="#000000">data</font><font color="#0000C0">[</font><font color="#004080">'plasma'</font><font color="#0000C0">]</font>
    <font color="#000000">start</font> <font color="#0000C0">=</font> <font color="#000000">data</font><font color="#0000C0">[</font><font color="#004080">'plasma_start'</font><font color="#0000C0">]</font><font color="#0000C0">*</font><font color="#0080C0">1e3</font>
    <font color="#000000">end</font> <font color="#0000C0">=</font> <font color="#000000">data</font><font color="#0000C0">[</font><font color="#004080">'plasma_end'</font><font color="#0000C0">]</font><font color="#0000C0">*</font><font color="#0080C0">1e3</font>

    <font color="#000000">nGW_limit</font> <font color="#0000C0">=</font> <font color="#0080C0">3</font>  <font color="#008000"># n/nGW &lt; 3 </font>
    <font color="#000000">calb_photo</font> <font color="#0000C0">=</font> <font color="#0080C0">10</font> <font color="#008000"># guess of absolute calbration of photodiode !!</font>

    <font color="#C00000">if</font> <font color="#C00000">not</font> <font color="#000000">plasma</font><font color="#0000C0">:</font>
        <font color="#000000">t_cd</font><font color="#0000C0">,</font> <font color="#000000">I_cd</font> <font color="#0000C0">=</font> <font color="#000000">Shot</font><font color="#0000C0">(</font><font color="#0000C0">)</font><font color="#0000C0">[</font><font color="#004080">'current_cd_coils'</font><font color="#0000C0">]</font>   <font color="#008000"># current drive coils current</font>
        <font color="#008000">#plot(t_cd, I_cd)</font>
        <font color="#008000">#savefig('I_cd.png')</font>
        <font color="#008000">#clf()</font>
        <font color="#000000">_</font><font color="#0000C0">,</font> <font color="#000000">I_cd</font><font color="#0000C0">,</font> <font color="#000000">_</font> <font color="#0000C0">=</font> <font color="#000000">DiffFilter</font><font color="#0000C0">(</font><font color="#000000">I_cd</font><font color="#0000C0">,</font> <font color="#000000">mean</font><font color="#0000C0">(</font><font color="#000000">diff</font><font color="#0000C0">(</font><font color="#000000">t_cd</font><font color="#0000C0">)</font><font color="#0000C0">)</font><font color="#0000C0">,</font> <font color="#0080C0">1000</font><font color="#0000C0">,</font> <font color="#0080C0">1e2</font><font color="#0000C0">)</font>  <font color="#008000"># smooth the data </font>
        <font color="#000000">t_PHI</font><font color="#0000C0">,</font> <font color="#000000">PHI</font> <font color="#0000C0">=</font> <font color="#000000">Shot</font><font color="#0000C0">(</font><font color="#0000C0">)</font><font color="#0000C0">[</font><font color="#004080">'total_magnetic_flux'</font><font color="#0000C0">]</font>


        <font color="#000000">I_cd</font> <font color="#0000C0">=</font> <font color="#000000">interp</font><font color="#0000C0">(</font><font color="#000000">t_PHI</font><font color="#0000C0">,</font> <font color="#000000">t_cd</font><font color="#0000C0">,</font> <font color="#000000">squeeze</font><font color="#0000C0">(</font><font color="#000000">I_cd</font><font color="#0000C0">)</font><font color="#0000C0">,</font> <font color="#000000">left</font><font color="#0000C0">=</font><font color="#0080C0">0</font><font color="#0000C0">,</font> <font color="#000000">right</font><font color="#0000C0">=</font><font color="#0080C0">0</font><font color="#0000C0">)</font>

        <font color="#000000">data</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0000C0">[</font><font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#004080">'total_magnetic_flux'</font><font color="#0000C0">,</font> <font color="#004080">'Total magnetic flux  $\Phi$'</font><font color="#0000C0">,</font> <font color="#004080">'$\Psi$ [Vs]'</font><font color="#0000C0">,</font>  <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font> <font color="#0080C0">40</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">)</font><font color="#0000C0">,</font>
                <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#000000">array</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font> <font color="#0080C0">40</font><font color="#0000C0">]</font><font color="#0000C0">)</font><font color="#0000C0">*</font><font color="#0080C0">1e-3</font><font color="#0000C0">,</font> <font color="#000000">ones</font><font color="#0000C0">(</font><font color="#0080C0">2</font><font color="#0000C0">)</font><font color="#0000C0">*</font><font color="#000000">MaxTransformatorSaturation</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#004080">'Max saturation'</font><font color="#0000C0">,</font> <font color="#004080">'$\Psi$ [Vs]'</font> <font color="#0000C0">,</font> <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#0080C0">0.18</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font> <font color="#0080C0">40</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">)</font><font color="#0000C0">]</font><font color="#0000C0">,</font>
                <font color="#000000">get_data</font><font color="#0000C0">(</font> <font color="#0000C0">[</font><font color="#000000">t_PHI</font><font color="#0000C0">,</font> <font color="#000000">I_cd</font><font color="#0000C0">]</font> <font color="#0000C0">,</font> <font color="#004080">'Current drive I'</font><font color="#0000C0">,</font> <font color="#004080">'I [A]'</font><font color="#0000C0">,</font> <font color="#000000">xlabel</font> <font color="#0000C0">=</font> <font color="#004080">'Time [ms]'</font><font color="#0000C0">,</font>  <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font> <font color="#0080C0">40</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">tvec_rescale</font><font color="#0000C0">=</font><font color="#0080C0">1e3</font><font color="#0000C0">)</font><font color="#0000C0">,</font>
                <font color="#0000C0">]</font>


        <font color="#000000">paralel_multiplot</font><font color="#0000C0">(</font><font color="#000000">data</font><font color="#0000C0">,</font> <font color="#004080">'Total magnetic flux'</font><font color="#0000C0">,</font> <font color="#004080">'magnetic_flux'</font><font color="#0000C0">,</font> <font color="#0000C0">(</font><font color="#0080C0">9</font><font color="#0000C0">,</font><font color="#0080C0">6</font><font color="#0000C0">)</font><font color="#0000C0">,</font> <font color="#0080C0">100</font><font color="#0000C0">,</font> <font color="#004080">'vertical'</font><font color="#0000C0">,</font> <font color="#000000">file_type</font><font color="#0000C0">)</font>

        <font color="#000000">plot</font><font color="#0000C0">(</font> <font color="#000000">I_cd</font><font color="#0000C0">,</font><font color="#000000">PHI</font><font color="#0000C0">)</font>
        <font color="#008000">#plot( PHI)</font>

        <font color="#000000">savefig</font><font color="#0000C0">(</font><font color="#004080">'I.png'</font><font color="#0000C0">)</font>
        <font color="#000000">clf</font><font color="#0000C0">(</font><font color="#0000C0">)</font>
        <font color="#000000">data</font> <font color="#0000C0">=</font> <font color="#000000">get_data</font><font color="#0000C0">(</font> <font color="#0000C0">[</font><font color="#000000">I_cd</font><font color="#0000C0">,</font> <font color="#000000">PHI</font> <font color="#0000C0">]</font> <font color="#0000C0">,</font> <font color="#004080">'Hysteresis curve'</font><font color="#0000C0">,</font> <font color="#004080">'$\Phi$ [Vs]'</font><font color="#0000C0">,</font> <font color="#000000">xlabel</font> <font color="#0000C0">=</font> <font color="#004080">'I [A]'</font><font color="#0000C0">,</font><font color="#000000">xlim</font><font color="#0000C0">=</font><font color="#0000C0">[</font><font color="#000000">amin</font><font color="#0000C0">(</font><font color="#000000">I_cd</font><font color="#0000C0">)</font><font color="#0000C0">,</font> <font color="#000000">amax</font><font color="#0000C0">(</font><font color="#000000">I_cd</font><font color="#0000C0">)</font><font color="#0000C0">]</font><font color="#0000C0">,</font>  <font color="#000000">plot_limits</font> <font color="#0000C0">=</font> <font color="#000000">False</font><font color="#0000C0">,</font> <font color="#000000">tvec_rescale</font><font color="#0000C0">=</font><font color="#0080C0">1</font><font color="#0000C0">)</font>
        <font color="#000000">paralel_multiplot</font><font color="#0000C0">(</font><font color="#000000">data</font><font color="#0000C0">,</font> <font color="#004080">'Hysteresis curve'</font><font color="#0000C0">,</font> <font color="#004080">'hysteresis'</font><font color="#0000C0">,</font> <font color="#0000C0">(</font><font color="#0080C0">9</font><font color="#0000C0">,</font><font color="#0080C0">3</font><font color="#0000C0">)</font><font color="#0000C0">,</font> <font color="#0080C0">100</font><font color="#0000C0">,</font> <font color="#004080">'vertical'</font><font color="#0000C0">,</font> <font color="#000000">file_type</font><font color="#0000C0">)</font>


        <font color="#000000">paralel_multiplot</font><font color="#0000C0">(</font><font color="#000000">data</font><font color="#0000C0">,</font> <font color="#004080">""</font> <font color="#0000C0">,</font> <font color="#004080">'icon'</font><font color="#0000C0">,</font> <font color="#0000C0">(</font><font color="#0080C0">4</font><font color="#0000C0">,</font><font color="#0080C0">3</font><font color="#0000C0">)</font><font color="#0000C0">,</font> <font color="#0080C0">40</font><font color="#0000C0">)</font>


        <font color="#C00000">return</font>



    <font color="#000000">t0</font> <font color="#0000C0">=</font> <font color="#000000">time</font><font color="#0000C0">.</font><font color="#000000">time</font><font color="#0000C0">(</font><font color="#0000C0">)</font>

    <font color="#000000">tvec</font><font color="#0000C0">,</font> <font color="#000000">Q95</font> <font color="#0000C0">=</font> <font color="#000000">Shot</font><font color="#0000C0">(</font><font color="#0000C0">)</font><font color="#0000C0">[</font><font color="#004080">'Q95'</font><font color="#0000C0">]</font>

    <font color="#C00000">try</font><font color="#0000C0">:</font>
        <font color="#000000">ne_corr</font> <font color="#0000C0">=</font> <font color="#000000">r0</font><font color="#0000C0">/</font><font color="#000000">a</font><font color="#0000C0">*</font><font color="#0000C0">(</font><font color="#000000">pos_part</font><font color="#0000C0">(</font><font color="#0080C0">1</font><font color="#0000C0">-</font><font color="#000000">R</font><font color="#0000C0">/</font><font color="#000000">r0</font><font color="#0000C0">)</font><font color="#0000C0">)</font><font color="#0000C0">**</font><font color="#0080C0">1.5</font>
    <font color="#C00000">except</font><font color="#0000C0">:</font>
        <font color="#000000">tvec_n_e</font> <font color="#0000C0">=</font> <font color="#000000">tvec</font>


    <font color="#C00000">try</font><font color="#0000C0">:</font>
        <font color="#000000">tvec_r</font><font color="#0000C0">,</font> <font color="#000000">plasma_radius</font> <font color="#0000C0">=</font>  <font color="#000000">Shot</font><font color="#0000C0">(</font><font color="#0000C0">)</font><font color="#0000C0">[</font><font color="#004080">'plasma_radius'</font><font color="#0000C0">]</font>
        <font color="#000000">plasma_radius_f</font> <font color="#0000C0">=</font> <font color="#000000">interp1d</font><font color="#0000C0">(</font><font color="#000000">tvec_r</font><font color="#0000C0">,</font> <font color="#000000">plasma_radius</font><font color="#0000C0">,</font> <font color="#000000">fill_value</font> <font color="#0000C0">=</font> <font color="#0080C0">0</font><font color="#0000C0">,</font> <font color="#000000">bounds_error</font>  <font color="#0000C0">=</font>  <font color="#000000">False</font><font color="#0000C0">)</font>

        <font color="#000000">corr</font> <font color="#0000C0">=</font> <font color="#0000C0">(</font><font color="#000000">plasma_radius_f</font><font color="#0000C0">(</font><font color="#000000">tvec</font><font color="#0000C0">)</font> <font color="#0000C0">/</font> <font color="#000000">MeanPlasmaRadius</font> <font color="#0000C0">)</font> <font color="#0000C0">**</font> <font color="#0080C0">2</font>  <font color="#008000"># corrected for plasma readius</font>
        <font color="#000000">corr</font><font color="#0000C0">[</font> <font color="#000000">corr</font> <font color="#0000C0">&lt;</font> <font color="#0080C0">0.3</font> <font color="#0000C0">]</font> <font color="#0000C0">=</font> <font color="#000000">nan</font>

        <font color="#C00000">if</font> <font color="#000000">all</font><font color="#0000C0">(</font><font color="#000000">isnan</font><font color="#0000C0">(</font><font color="#000000">corr</font><font color="#0000C0">)</font><font color="#0000C0">)</font><font color="#0000C0">:</font>
          <font color="#C00000">del</font> <font color="#000000">corr</font>


    <font color="#C00000">except</font><font color="#0000C0">:</font>
        <font color="#C00000">print</font> <font color="#004080">"correction on plasma position failed"</font>

    <font color="#000000">tvec_gw</font><font color="#0000C0">,</font> <font color="#000000">nGW</font> <font color="#0000C0">=</font>  <font color="#000000">Shot</font><font color="#0000C0">(</font><font color="#0000C0">)</font><font color="#0000C0">[</font><font color="#004080">'greenwald_density'</font><font color="#0000C0">]</font>
    <font color="#000000">nGW_med</font> <font color="#0000C0">=</font> <font color="#000000">max</font><font color="#0000C0">(</font><font color="#000000">nGW</font><font color="#0000C0">[</font><font color="#0000C0">(</font><font color="#000000">tvec_gw</font> <font color="#0000C0">&gt;</font> <font color="#000000">start</font><font color="#0000C0">*</font><font color="#0080C0">1e-3</font><font color="#0000C0">)</font> <font color="#0000C0">&amp;</font> <font color="#0000C0">(</font><font color="#000000">tvec_gw</font> <font color="#0000C0">&lt;</font> <font color="#000000">end</font><font color="#0000C0">*</font><font color="#0080C0">1e-3</font><font color="#0000C0">)</font><font color="#0000C0">]</font><font color="#0000C0">)</font>


    <font color="#C00000">try</font><font color="#0000C0">:</font>
        <font color="#000000">tvec_n_e</font><font color="#0000C0">,</font> <font color="#000000">n_e</font> <font color="#0000C0">=</font>  <font color="#000000">Shot</font><font color="#0000C0">(</font><font color="#0000C0">)</font><font color="#0000C0">[</font><font color="#004080">'electron_density'</font><font color="#0000C0">]</font>
    <font color="#C00000">except</font><font color="#0000C0">:</font>
        <font color="#C00000">print</font> <font color="#004080">"electron_density failed"</font>

    <font color="#C00000">try</font><font color="#0000C0">:</font>
        <font color="#000000">tvec_r</font><font color="#0000C0">,</font> <font color="#000000">R_position</font> <font color="#0000C0">=</font>  <font color="#000000">Shot</font><font color="#0000C0">(</font><font color="#0000C0">)</font><font color="#0000C0">[</font><font color="#004080">'plasma_position_r'</font><font color="#0000C0">]</font>
        <font color="#000000">plasma_R_pos_f</font> <font color="#0000C0">=</font> <font color="#000000">interp1d</font><font color="#0000C0">(</font><font color="#000000">tvec_r</font><font color="#0000C0">,</font> <font color="#000000">R_position</font><font color="#0000C0">,</font> <font color="#000000">fill_value</font> <font color="#0000C0">=</font> <font color="#0080C0">0</font><font color="#0000C0">,</font> <font color="#000000">bounds_error</font>  <font color="#0000C0">=</font>  <font color="#000000">False</font><font color="#0000C0">)</font>


        <font color="#000000">a</font>  <font color="#0000C0">=</font> <font color="#0080C0">0.1</font> <font color="#008000">#[m]</font>
        <font color="#000000">r0</font> <font color="#0000C0">=</font> <font color="#000000">plasma_radius_f</font><font color="#0000C0">(</font><font color="#000000">tvec_n_e</font><font color="#0000C0">)</font>
        <font color="#000000">R</font> <font color="#0000C0">=</font> <font color="#000000">plasma_R_pos_f</font><font color="#0000C0">(</font><font color="#000000">tvec_n_e</font><font color="#0000C0">)</font><font color="#0000C0">-</font><font color="#0080C0">0.4</font>
        <font color="#000000">pos_part</font> <font color="#0000C0">=</font> <font color="#C00000">lambda</font> <font color="#000000">x</font><font color="#0000C0">:</font> <font color="#0000C0">(</font><font color="#000000">x</font><font color="#0000C0">+</font><font color="#000000">abs</font><font color="#0000C0">(</font><font color="#000000">x</font><font color="#0000C0">)</font><font color="#0000C0">)</font><font color="#0000C0">/</font><font color="#0080C0">2</font>

        <font color="#000000">ne_corr_nu1</font>   <font color="#0000C0">=</font> <font color="#000000">r0</font><font color="#0000C0">/</font><font color="#000000">a</font><font color="#0000C0">*</font><font color="#0000C0">(</font><font color="#000000">pos_part</font><font color="#0000C0">(</font><font color="#0080C0">1</font><font color="#0000C0">-</font><font color="#000000">R</font><font color="#0000C0">/</font><font color="#000000">r0</font><font color="#0000C0">)</font><font color="#0000C0">)</font><font color="#0000C0">**</font><font color="#0080C0">1.5</font>  <font color="#008000">#peaking factor 1; profile = (1-(r/a)^2)</font>
        <font color="#000000">ne_corr_nu0_5</font> <font color="#0000C0">=</font> <font color="#000000">r0</font><font color="#0000C0">/</font><font color="#000000">a</font><font color="#0000C0">*</font><font color="#0000C0">(</font><font color="#000000">pos_part</font><font color="#0000C0">(</font><font color="#0080C0">1</font><font color="#0000C0">-</font><font color="#000000">R</font><font color="#0000C0">/</font><font color="#000000">r0</font><font color="#0000C0">)</font><font color="#0000C0">)</font>  <font color="#008000">#peaking factor 0.5; profile = sqrt(1-(r/a)^2)</font>

        <font color="#000000">ne_corr_nu1</font><font color="#0000C0">[</font><font color="#000000">ne_corr_nu1</font><font color="#0000C0">&lt;</font><font color="#0080C0">0.2</font><font color="#0000C0">]</font> <font color="#0000C0">=</font> <font color="#000000">nan</font>
        <font color="#000000">ne_corr_nu0_5</font><font color="#0000C0">[</font><font color="#000000">ne_corr_nu0_5</font><font color="#0000C0">&lt;</font><font color="#0080C0">0.2</font><font color="#0000C0">]</font> <font color="#0000C0">=</font> <font color="#000000">nan</font>

        <font color="#C00000">if</font> <font color="#000000">all</font><font color="#0000C0">(</font><font color="#000000">isnan</font><font color="#0000C0">(</font><font color="#000000">ne_corr_nu0_5</font><font color="#0000C0">)</font><font color="#0000C0">)</font><font color="#0000C0">:</font>
          <font color="#C00000">del</font> <font color="#000000">ne_corr_nu0_5</font>
          <font color="#C00000">del</font> <font color="#000000">ne_corr_nu1</font>

    <font color="#C00000">except</font><font color="#0000C0">:</font>
      <font color="#C00000">print</font> <font color="#004080">"electron_density correction failed"</font>

    <font color="#000000">Q95</font><font color="#0000C0">[</font><font color="#000000">Q95</font> <font color="#0000C0">&gt;</font>  <font color="#0080C0">2</font><font color="#0000C0">*</font> <font color="#000000">median</font><font color="#0000C0">(</font><font color="#000000">Q95</font><font color="#0000C0">[</font><font color="#0000C0">(</font><font color="#000000">tvec</font> <font color="#0000C0">&gt;</font> <font color="#000000">start</font><font color="#0000C0">*</font><font color="#0080C0">1e-3</font><font color="#0000C0">)</font> <font color="#0000C0">&amp;</font> <font color="#0000C0">(</font><font color="#000000">tvec</font> <font color="#0000C0">&lt;</font> <font color="#000000">end</font><font color="#0000C0">*</font><font color="#0080C0">1e-3</font><font color="#0000C0">)</font><font color="#0000C0">]</font><font color="#0000C0">)</font><font color="#0000C0">]</font> <font color="#0000C0">=</font> <font color="#000000">nan</font>

    <font color="#C00000">try</font><font color="#0000C0">:</font>
      <font color="#000000">temp_spec</font> <font color="#0000C0">=</font>  <font color="#000000">Shot</font><font color="#0000C0">(</font><font color="#0000C0">)</font><font color="#0000C0">[</font><font color="#004080">'spectrometr\temperature'</font><font color="#0000C0">]</font>
      <font color="#000000">power_spec</font> <font color="#0000C0">=</font>  <font color="#000000">Shot</font><font color="#0000C0">(</font><font color="#0000C0">)</font><font color="#0000C0">[</font><font color="#004080">'spectrometr\total_radiated_power'</font><font color="#0000C0">]</font>
    <font color="#C00000">except</font><font color="#0000C0">:</font>
      <font color="#C00000">print</font> <font color="#004080">"missing data spectrometer"</font>


    <font color="#000000">data_1</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font>
        <font color="#0000C0">[</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#000000">tvec</font><font color="#0000C0">,</font><font color="#000000">Q95</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#004080">'Safety factor edge (Q$_{edge}$)'</font><font color="#0000C0">,</font> <font color="#004080">'Q [-]'</font> <font color="#0000C0">,</font> <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#000000">None</font><font color="#0000C0">]</font> <font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">)</font><font color="#0000C0">,</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#000000">tvec</font><font color="#0000C0">,</font><font color="#000000">Q95</font><font color="#0000C0">/</font><font color="#0080C0">3</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#004080">'(Q$_{center}$) (peaking factor = 2)'</font><font color="#0000C0">,</font> <font color="#004080">'Q [-]'</font> <font color="#0000C0">,</font> <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#000000">None</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">)</font><font color="#0000C0">,</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#000000">tvec</font><font color="#0000C0">,</font><font color="#000000">corr</font><font color="#0000C0">*</font><font color="#000000">Q95</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#004080">'Q edge with radius'</font><font color="#0000C0">,</font> <font color="#004080">'Q [-]'</font> <font color="#0000C0">,</font> <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#000000">None</font><font color="#0000C0">]</font> <font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">,</font> <font color="#000000">line_format</font><font color="#0000C0">=</font><font color="#004080">'b--'</font><font color="#0000C0">)</font> <font color="#C00000">if</font> <font color="#004080">'corr'</font> <font color="#C00000">in</font> <font color="#000000">locals</font><font color="#0000C0">(</font><font color="#0000C0">)</font>   <font color="#C00000">else</font> <font color="#0000C0">[</font><font color="#0000C0">]</font><font color="#0000C0">,</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#000000">tvec</font><font color="#0000C0">,</font><font color="#000000">corr</font><font color="#0000C0">*</font><font color="#000000">Q95</font><font color="#0000C0">/</font><font color="#0080C0">3</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#004080">'Q center with radius'</font><font color="#0000C0">,</font> <font color="#004080">'Q [-]'</font> <font color="#0000C0">,</font> <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#000000">None</font><font color="#0000C0">]</font> <font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">,</font> <font color="#000000">line_format</font><font color="#0000C0">=</font><font color="#004080">'r--'</font><font color="#0000C0">)</font> <font color="#C00000">if</font> <font color="#004080">'corr'</font> <font color="#C00000">in</font> <font color="#000000">locals</font><font color="#0000C0">(</font><font color="#0000C0">)</font>   <font color="#C00000">else</font> <font color="#0000C0">[</font><font color="#0000C0">]</font>
        <font color="#0000C0">]</font><font color="#0000C0">,</font>
        <font color="#0000C0">[</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#000000">tvec</font><font color="#0000C0">,</font><font color="#0080C0">1</font><font color="#0000C0">/</font><font color="#000000">Q95</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#004080">'Inverse Safety factor $\iota_{edge}$'</font><font color="#0000C0">,</font> <font color="#004080">'$\iota$ [-]'</font> <font color="#0000C0">,</font> <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#000000">None</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">)</font><font color="#0000C0">,</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#000000">tvec</font><font color="#0000C0">,</font><font color="#0080C0">3</font><font color="#0000C0">/</font><font color="#000000">Q95</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#004080">' $\iota_{center}$ (peaking factor = 2)'</font><font color="#0000C0">,</font> <font color="#004080">'$\iota$ [-]'</font> <font color="#0000C0">,</font> <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#000000">None</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">)</font><font color="#0000C0">,</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#000000">tvec</font><font color="#0000C0">,</font><font color="#0080C0">1</font><font color="#0000C0">/</font><font color="#0000C0">(</font><font color="#000000">Q95</font><font color="#0000C0">*</font><font color="#000000">corr</font><font color="#0000C0">)</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#004080">' $\iota_{edge}$ with radius'</font><font color="#0000C0">,</font> <font color="#004080">'$\iota$ [-]'</font> <font color="#0000C0">,</font> <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#000000">None</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">line_format</font><font color="#0000C0">=</font><font color="#004080">'b--'</font><font color="#0000C0">,</font>  <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">)</font> <font color="#C00000">if</font> <font color="#004080">'corr'</font> <font color="#C00000">in</font> <font color="#000000">locals</font><font color="#0000C0">(</font><font color="#0000C0">)</font>   <font color="#C00000">else</font> <font color="#0000C0">[</font><font color="#0000C0">]</font><font color="#0000C0">,</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#000000">tvec</font><font color="#0000C0">,</font><font color="#0080C0">3</font><font color="#0000C0">/</font><font color="#0000C0">(</font><font color="#000000">Q95</font><font color="#0000C0">*</font><font color="#000000">corr</font><font color="#0000C0">)</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#004080">' $\iota_{center}$ with radius'</font><font color="#0000C0">,</font> <font color="#004080">'$\iota$ [-]'</font> <font color="#0000C0">,</font> <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#000000">None</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">line_format</font><font color="#0000C0">=</font><font color="#004080">'r--'</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">)</font> <font color="#C00000">if</font> <font color="#004080">'corr'</font> <font color="#C00000">in</font> <font color="#000000">locals</font><font color="#0000C0">(</font><font color="#0000C0">)</font>   <font color="#C00000">else</font> <font color="#0000C0">[</font><font color="#0000C0">]</font><font color="#0000C0">,</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#000000">array</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font><font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">)</font><font color="#0000C0">*</font><font color="#0080C0">1e-3</font><font color="#0000C0">,</font> <font color="#000000">array</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#0080C0">0.5</font><font color="#0000C0">,</font> <font color="#0080C0">0.5</font><font color="#0000C0">]</font><font color="#0000C0">)</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#004080">'Q=2'</font><font color="#0000C0">,</font> <font color="#004080">'$\iota$ [-]'</font> <font color="#0000C0">,</font> <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#0080C0">0.55</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">,</font> <font color="#000000">line_format</font> <font color="#0000C0">=</font> <font color="#004080">"k:"</font><font color="#0000C0">)</font>
        <font color="#0000C0">]</font>
        <font color="#0000C0">]</font>

    <font color="#000000">multiplot</font><font color="#0000C0">(</font><font color="#000000">data_1</font><font color="#0000C0">,</font> <font color="#004080">'Safety factor'</font><font color="#0000C0">,</font> <font color="#004080">'plot_1'</font><font color="#0000C0">,</font> <font color="#0000C0">(</font><font color="#0080C0">9</font><font color="#0000C0">,</font><font color="#0080C0">5</font><font color="#0000C0">)</font><font color="#0000C0">,</font> <font color="#0080C0">100</font><font color="#0000C0">,</font> <font color="#004080">'vertical'</font><font color="#0000C0">,</font> <font color="#000000">file_type</font><font color="#0000C0">)</font>

    <font color="#000000">data_2</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0000C0">[</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#004080">'electron_temperature_median'</font><font color="#0000C0">,</font> <font color="#004080">'Electron temperature'</font><font color="#0000C0">,</font> <font color="#004080">'T [eV]'</font> <font color="#0000C0">,</font> <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#000000">None</font><font color="#0000C0">]</font> <font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">)</font><font color="#0000C0">,</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#004080">'electron_temperature_median'</font><font color="#0000C0">,</font> <font color="#004080">'$T_e$ corrected by plasma radius'</font><font color="#0000C0">,</font> <font color="#004080">'T [eV]'</font> <font color="#0000C0">,</font> <font color="#000000">data_rescale</font> <font color="#0000C0">=</font> <font color="#0080C0">1</font><font color="#0000C0">/</font><font color="#000000">corr</font><font color="#0000C0">,</font>  <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#000000">None</font><font color="#0000C0">]</font> <font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">line_format</font><font color="#0000C0">=</font><font color="#004080">'r--'</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">)</font> <font color="#C00000">if</font> <font color="#004080">'corr'</font> <font color="#C00000">in</font> <font color="#000000">locals</font><font color="#0000C0">(</font><font color="#0000C0">)</font>   <font color="#C00000">else</font> <font color="#0000C0">[</font><font color="#0000C0">]</font><font color="#0000C0">,</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#000000">temp_spec</font><font color="#0000C0">,</font> <font color="#004080">'$T_e$ spectrometer'</font><font color="#0000C0">,</font> <font color="#004080">'T [eV]'</font> <font color="#0000C0">,</font>  <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#000000">None</font><font color="#0000C0">]</font> <font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">line_format</font><font color="#0000C0">=</font><font color="#004080">'r--'</font><font color="#0000C0">)</font> <font color="#C00000">if</font> <font color="#004080">'temp_spec'</font> <font color="#C00000">in</font> <font color="#000000">locals</font><font color="#0000C0">(</font><font color="#0000C0">)</font>   <font color="#C00000">else</font> <font color="#0000C0">[</font><font color="#0000C0">]</font><font color="#0000C0">,</font>
        <font color="#008000">#get_data('electron_temperature', 'Electron temperature', 'T [eV]' , ylim = [0,None], xlim = [start, end])</font>
        <font color="#0000C0">]</font><font color="#0000C0">]</font>
    <font color="#C00000">try</font><font color="#0000C0">:</font>
      <font color="#000000">multiplot</font><font color="#0000C0">(</font><font color="#000000">data_2</font><font color="#0000C0">,</font> <font color="#004080">'Electron temperature'</font><font color="#0000C0">,</font> <font color="#004080">'plot_2'</font><font color="#0000C0">,</font> <font color="#0000C0">(</font><font color="#0080C0">9</font><font color="#0000C0">,</font><font color="#0080C0">4</font><font color="#0000C0">)</font><font color="#0000C0">,</font> <font color="#0080C0">100</font><font color="#0000C0">,</font> <font color="#004080">'vertical'</font><font color="#0000C0">,</font> <font color="#000000">file_type</font><font color="#0000C0">)</font>
    <font color="#C00000">except</font> <font color="#000000">Exception</font><font color="#0000C0">,</font> <font color="#000000">e</font><font color="#0000C0">:</font>
      <font color="#C00000">print</font> <font color="#004080">'Electron temperature failed'</font><font color="#0000C0">,</font> <font color="#000000">e</font><font color="#0000C0">.</font><font color="#000000">__str__</font><font color="#0000C0">(</font><font color="#0000C0">)</font>

    <font color="#000000">data_3</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0000C0">[</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#004080">'input_power_total'</font><font color="#0000C0">,</font> <font color="#004080">'Ohmic heating - total'</font><font color="#0000C0">,</font> <font color="#004080">'P [kW]'</font> <font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">data_rescale</font><font color="#0000C0">=</font><font color="#0080C0">1e-3</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">)</font><font color="#0000C0">,</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#004080">'input_power_magnetic'</font><font color="#0000C0">,</font> <font color="#004080">'Ohmic heating - magnetic'</font><font color="#0000C0">,</font> <font color="#004080">'P [kW]'</font> <font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">data_rescale</font><font color="#0000C0">=</font><font color="#0000C0">-</font><font color="#0080C0">1e-3</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">)</font><font color="#0000C0">,</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#004080">'input_power_plasma'</font><font color="#0000C0">,</font> <font color="#004080">'Ohmic heating - plasma'</font><font color="#0000C0">,</font> <font color="#004080">'P [kW]'</font> <font color="#0000C0">,</font>  <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">data_rescale</font><font color="#0000C0">=</font><font color="#0000C0">-</font><font color="#0080C0">1e-3</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">)</font><font color="#0000C0">,</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#004080">'input_power_chamber'</font><font color="#0000C0">,</font> <font color="#004080">'Ohmic heating - chamber'</font><font color="#0000C0">,</font> <font color="#004080">'P [kW]'</font> <font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">data_rescale</font><font color="#0000C0">=</font><font color="#0000C0">-</font><font color="#0080C0">1e-3</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">)</font><font color="#0000C0">,</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#004080">'photodiode'</font><font color="#0000C0">,</font> <font color="#004080">'Output power (diod)'</font><font color="#0000C0">,</font> <font color="#004080">'P [kW]'</font> <font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">data_rescale</font><font color="#0000C0">=</font> <font color="#000000">calb_photo</font> <font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">,</font> <font color="#000000">line_format</font> <font color="#0000C0">=</font> <font color="#004080">"k--"</font><font color="#0000C0">)</font><font color="#0000C0">,</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#000000">power_spec</font><font color="#0000C0">,</font> <font color="#004080">'Output power (spectrometer)'</font><font color="#0000C0">,</font> <font color="#004080">'P [kW]'</font> <font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">data_rescale</font><font color="#0000C0">=</font><font color="#0080C0">1e-3</font> <font color="#0000C0">,</font>  <font color="#000000">line_format</font> <font color="#0000C0">=</font> <font color="#004080">"r--"</font><font color="#0000C0">)</font> <font color="#C00000">if</font> <font color="#004080">'power_spec'</font> <font color="#C00000">in</font> <font color="#000000">locals</font><font color="#0000C0">(</font><font color="#0000C0">)</font> <font color="#C00000">else</font> <font color="#0000C0">[</font><font color="#0000C0">]</font><font color="#0000C0">,</font>
        <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#000000">array</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font><font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">)</font><font color="#0000C0">*</font><font color="#0080C0">1e-3</font><font color="#0000C0">,</font> <font color="#000000">array</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font> <font color="#0080C0">0</font><font color="#0000C0">]</font><font color="#0000C0">)</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#004080">''</font><font color="#0000C0">,</font> <font color="#004080">'P [kW]'</font> <font color="#0000C0">,</font>  <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">,</font> <font color="#000000">line_format</font> <font color="#0000C0">=</font> <font color="#004080">"k:"</font><font color="#0000C0">)</font>
        <font color="#0000C0">]</font><font color="#0000C0">]</font>
    <font color="#C00000">try</font><font color="#0000C0">:</font>
      <font color="#000000">multiplot</font><font color="#0000C0">(</font><font color="#000000">data_3</font><font color="#0000C0">,</font> <font color="#004080">'Input / Output power'</font><font color="#0000C0">,</font> <font color="#004080">'plot_3'</font><font color="#0000C0">,</font> <font color="#0000C0">(</font><font color="#0080C0">9</font><font color="#0000C0">,</font><font color="#0080C0">4</font><font color="#0000C0">)</font><font color="#0000C0">,</font> <font color="#0080C0">100</font><font color="#0000C0">,</font> <font color="#004080">'vertical'</font><font color="#0000C0">,</font> <font color="#000000">file_type</font><font color="#0000C0">)</font>
    <font color="#C00000">except</font> <font color="#000000">Exception</font><font color="#0000C0">,</font> <font color="#000000">e</font><font color="#0000C0">:</font>
      <font color="#C00000">print</font> <font color="#004080">'Input / Output power'</font><font color="#0000C0">,</font> <font color="#000000">e</font><font color="#0000C0">.</font><font color="#000000">__str__</font><font color="#0000C0">(</font><font color="#0000C0">)</font>

    <font color="#C00000">if</font> <font color="#004080">'n_e'</font> <font color="#C00000">in</font> <font color="#000000">locals</font><font color="#0000C0">(</font><font color="#0000C0">)</font><font color="#0000C0">:</font>
      <font color="#000000">data_4</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0000C0">[</font>
          <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#004080">'electron_density'</font><font color="#0000C0">,</font> <font color="#004080">'Electron density ($n_e$)'</font><font color="#0000C0">,</font> <font color="#004080">'n$_{e}$ [10$^{19}\cdot$m$^{-3}$]'</font> <font color="#0000C0">,</font> <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#000000">None</font><font color="#0000C0">]</font> <font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">data_rescale</font><font color="#0000C0">=</font> <font color="#0080C0">1e-19</font><font color="#0000C0">)</font> <font color="#C00000">if</font> <font color="#004080">'n_e'</font> <font color="#C00000">in</font> <font color="#000000">locals</font><font color="#0000C0">(</font><font color="#0000C0">)</font>   <font color="#C00000">else</font> <font color="#0000C0">[</font><font color="#0000C0">]</font><font color="#0000C0">,</font>
          <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#004080">'electron_density'</font><font color="#0000C0">,</font> <font color="#004080">'$n_e$ - position corrected, $\\nu = 1$'</font><font color="#0000C0">,</font> <font color="#004080">'n$_{e}$ [10$^{19}\cdot$m$^{-3}$]'</font> <font color="#0000C0">,</font> <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#000000">None</font><font color="#0000C0">]</font> <font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">data_rescale</font><font color="#0000C0">=</font> <font color="#0080C0">1e-19</font><font color="#0000C0">/</font><font color="#000000">ne_corr_nu1</font><font color="#0000C0">,</font> <font color="#000000">line_format</font><font color="#0000C0">=</font><font color="#004080">'k--'</font><font color="#0000C0">)</font> <font color="#C00000">if</font> <font color="#0000C0">(</font><font color="#004080">'ne_corr_nu1'</font> <font color="#C00000">in</font> <font color="#000000">locals</font><font color="#0000C0">(</font><font color="#0000C0">)</font><font color="#0000C0">)</font> <font color="#C00000">and</font> <font color="#0000C0">(</font><font color="#004080">'n_e'</font> <font color="#C00000">in</font> <font color="#000000">locals</font><font color="#0000C0">(</font><font color="#0000C0">)</font> <font color="#0000C0">)</font>  <font color="#C00000">else</font> <font color="#0000C0">[</font><font color="#0000C0">]</font><font color="#0000C0">,</font>
          <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#004080">'electron_density'</font><font color="#0000C0">,</font> <font color="#004080">'$n_e$ - position corrected, $\\nu = 0.5$'</font><font color="#0000C0">,</font> <font color="#004080">'n$_{e}$ [10$^{19}\cdot$m$^{-3}$]'</font> <font color="#0000C0">,</font> <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#000000">None</font><font color="#0000C0">]</font> <font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">data_rescale</font><font color="#0000C0">=</font> <font color="#0080C0">1e-19</font><font color="#0000C0">/</font><font color="#000000">ne_corr_nu0_5</font><font color="#0000C0">,</font> <font color="#000000">line_format</font><font color="#0000C0">=</font><font color="#004080">'k:'</font><font color="#0000C0">)</font> <font color="#C00000">if</font> <font color="#0000C0">(</font><font color="#004080">'ne_corr_nu0_5'</font> <font color="#C00000">in</font> <font color="#000000">locals</font><font color="#0000C0">(</font><font color="#0000C0">)</font><font color="#0000C0">)</font> <font color="#C00000">and</font>  <font color="#0000C0">(</font><font color="#004080">'n_e'</font> <font color="#C00000">in</font> <font color="#000000">locals</font><font color="#0000C0">(</font><font color="#0000C0">)</font><font color="#0000C0">)</font> <font color="#C00000">else</font> <font color="#0000C0">[</font><font color="#0000C0">]</font><font color="#0000C0">,</font>
          <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#004080">'greenwald_density'</font><font color="#0000C0">,</font> <font color="#004080">'Greenwald density ($n_{GW}$)'</font><font color="#0000C0">,</font> <font color="#004080">'n$_{GW}$ [10$^{19}\cdot$m$^{-3}$]'</font> <font color="#0000C0">,</font> <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#000000">None</font><font color="#0000C0">]</font> <font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">data_rescale</font><font color="#0000C0">=</font> <font color="#0080C0">1e-19</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">)</font><font color="#0000C0">,</font>
          <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#004080">'greenwald_density'</font><font color="#0000C0">,</font> <font color="#004080">'Greenwald density with radius ($n_{GW}$)'</font><font color="#0000C0">,</font> <font color="#004080">'n$_{GW}$ [10$^{19}\cdot$m$^{-3}$]'</font> <font color="#0000C0">,</font> <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#000000">None</font><font color="#0000C0">]</font> <font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#000000">data_rescale</font><font color="#0000C0">=</font><font color="#0080C0">1e-19</font><font color="#0000C0">/</font><font color="#000000">corr</font><font color="#0000C0">,</font> <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">,</font> <font color="#000000">line_format</font><font color="#0000C0">=</font><font color="#004080">'r--'</font><font color="#0000C0">)</font> <font color="#C00000">if</font> <font color="#004080">'corr'</font> <font color="#C00000">in</font> <font color="#000000">locals</font><font color="#0000C0">(</font><font color="#0000C0">)</font>   <font color="#C00000">else</font> <font color="#0000C0">[</font><font color="#0000C0">]</font><font color="#0000C0">,</font>
          <font color="#000000">get_data</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#000000">array</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font><font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">)</font><font color="#0000C0">*</font><font color="#0080C0">1e-3</font><font color="#0000C0">,</font> <font color="#000000">array</font><font color="#0000C0">(</font><font color="#0000C0">[</font><font color="#000000">nGW_med</font><font color="#0000C0">,</font> <font color="#000000">nGW_med</font><font color="#0000C0">]</font><font color="#0000C0">)</font><font color="#0000C0">*</font><font color="#000000">nGW_limit</font><font color="#0000C0">]</font><font color="#0000C0">,</font> <font color="#004080">'~GW limit'</font><font color="#0000C0">,</font> <font color="#004080">'n$_{e}$ [10$^{19}\cdot$m$^{-3}$]'</font> <font color="#0000C0">,</font> <font color="#000000">ylim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#0080C0">0</font><font color="#0000C0">,</font><font color="#000000">None</font><font color="#0000C0">]</font> <font color="#0000C0">,</font> <font color="#000000">xlim</font> <font color="#0000C0">=</font> <font color="#0000C0">[</font><font color="#000000">start</font><font color="#0000C0">,</font> <font color="#000000">end</font><font color="#0000C0">]</font><font color="#0000C0">,</font><font color="#000000">data_rescale</font><font color="#0000C0">=</font> <font color="#0080C0">1e-19</font> <font color="#0000C0">,</font>  <font color="#000000">reduction</font> <font color="#0000C0">=</font> <font color="#000000">True</font><font color="#0000C0">,</font> <font color="#000000">line_format</font> <font color="#0000C0">=</font> <font color="#004080">"k:"</font><font color="#0000C0">)</font> <font color="#C00000">if</font> <font color="#004080">'n_e'</font> <font color="#C00000">in</font> <font color="#000000">locals</font><font color="#0000C0">(</font><font color="#0000C0">)</font> <font color="#C00000">else</font> <font color="#0000C0">[</font><font color="#0000C0">]</font><font color="#0000C0">,</font>
          <font color="#0000C0">]</font><font color="#0000C0">]</font>


    <font color="#C00000">try</font><font color="#0000C0">:</font>
      <font color="#000000">multiplot</font><font color="#0000C0">(</font><font color="#000000">data_4</font><font color="#0000C0">,</font> <font color="#004080">'Density'</font><font color="#0000C0">,</font> <font color="#004080">'plot_4'</font><font color="#0000C0">,</font> <font color="#0000C0">(</font><font color="#0080C0">9</font><font color="#0000C0">,</font><font color="#0080C0">4</font><font color="#0000C0">)</font><font color="#0000C0">,</font> <font color="#0080C0">100</font><font color="#0000C0">,</font> <font color="#004080">'vertical'</font><font color="#0000C0">,</font> <font color="#000000">file_type</font><font color="#0000C0">)</font>
    <font color="#C00000">except</font> <font color="#000000">Exception</font><font color="#0000C0">,</font> <font color="#000000">e</font><font color="#0000C0">:</font>
      <font color="#C00000">print</font> <font color="#004080">'Density'</font><font color="#0000C0">,</font> <font color="#000000">e</font><font color="#0000C0">.</font><font color="#000000">__str__</font><font color="#0000C0">(</font><font color="#0000C0">)</font>


    <font color="#008000">#os.system('convert -resize 150x120\! total_plot.png icon.png')</font>

    <font color="#C00000">print</font> <font color="#004080">"===== plotting done ===== "</font>
    <font color="#000000">data</font> <font color="#0000C0">=</font> <font color="#000000">data_1</font> <font color="#0000C0">+</font> <font color="#000000">data_2</font> <font color="#0000C0">+</font><font color="#000000">data_3</font> <font color="#0000C0">+</font> <font color="#0000C0">(</font><font color="#000000">data_4</font> <font color="#C00000">if</font> <font color="#004080">'n_e'</font> <font color="#C00000">in</font> <font color="#000000">locals</font><font color="#0000C0">(</font><font color="#0000C0">)</font> <font color="#C00000">else</font> <font color="#0000C0">[</font><font color="#0000C0">]</font><font color="#0000C0">)</font>

    <font color="#000000">paralel_multiplot</font><font color="#0000C0">(</font><font color="#000000">data</font><font color="#0000C0">,</font> <font color="#004080">""</font> <font color="#0000C0">,</font> <font color="#004080">'icon'</font><font color="#0000C0">,</font> <font color="#0000C0">(</font><font color="#0080C0">4</font><font color="#0000C0">,</font><font color="#0080C0">3</font><font color="#0000C0">)</font><font color="#0000C0">,</font> <font color="#0080C0">40</font><font color="#0000C0">)</font>


<font color="#C00000">def</font> <font color="#000000">main</font><font color="#0000C0">(</font><font color="#0000C0">)</font><font color="#0000C0">:</font>
    <font color="#C00000">if</font> <font color="#000000">sys</font><font color="#0000C0">.</font><font color="#000000">argv</font><font color="#0000C0">[</font><font color="#0080C0">1</font><font color="#0000C0">]</font> <font color="#0000C0">==</font>  <font color="#004080">"analysis"</font><font color="#0000C0">:</font>
        <font color="#000000">analysis</font><font color="#0000C0">(</font><font color="#0000C0">)</font>

    <font color="#C00000">if</font> <font color="#000000">sys</font><font color="#0000C0">.</font><font color="#000000">argv</font><font color="#0000C0">[</font><font color="#0080C0">1</font><font color="#0000C0">]</font> <font color="#0000C0">==</font>  <font color="#004080">"plots"</font><font color="#0000C0">:</font>
        <font color="#000000">plot_data</font><font color="#0000C0">(</font><font color="#004080">'png'</font><font color="#0000C0">)</font>
        <font color="#000000">saveconst</font><font color="#0000C0">(</font><font color="#004080">'status'</font><font color="#0000C0">,</font> <font color="#0080C0">0</font><font color="#0000C0">)</font>

    <font color="#C00000">if</font> <font color="#000000">sys</font><font color="#0000C0">.</font><font color="#000000">argv</font><font color="#0000C0">[</font><font color="#0080C0">1</font><font color="#0000C0">]</font> <font color="#0000C0">==</font>  <font color="#004080">"postanalysis"</font><font color="#0000C0">:</font>
        <font color="#000000">plot_data</font><font color="#0000C0">(</font><font color="#004080">'png'</font><font color="#0000C0">)</font>
        <font color="#000000">saveconst</font><font color="#0000C0">(</font><font color="#004080">'status'</font><font color="#0000C0">,</font> <font color="#0080C0">0</font><font color="#0000C0">)</font>


<font color="#C00000">if</font> <font color="#000000">__name__</font> <font color="#0000C0">==</font> <font color="#004080">"__main__"</font><font color="#0000C0">:</font>
    <font color="#000000">main</font><font color="#0000C0">(</font><font color="#0000C0">)</font><font color="#000000"></font><a class="viewcode-back" href="javascript:history.go(-1)">[Return]</a></font></pre>
	    </div>
	  </div>
	</div>
      </div>

      <div class="sphinxsidebar">
	<div class="sphinxsidebarwrapper">


	    <a href= /shots/10104 ><img src="/_static/logos/golem.svg" alt="" width=80 /></a>



		<h4>Navigation</h4>
		<ul class="this-page-menu">
		    
    

     <li>
    
    <a  class="reference internal" href="javascript:history.go(-1)"><b>Return</b></a>
    </li>

		</ul>

	    <div id="searchbox" style="display: none">
		<h3>Go to shot</h3>
		    <form  style='display:inline;' method="post" action="" onsubmit="return getURL(this.shot.value)">
		    <input size="5" type="text" name="shot" value="10104"  />
		    <input  size="5" type="submit" name="Go" value="Go" />
		    </form>
	    </div>
	    <script type="text/javascript">$('#searchbox').show(0);</script>
	</div>
	</div>
      <div class="clearer"></div>
    </div>

    
    <div class="related">
	<span class="base_template"><a  href="/shots/10104/base.php.template" title="main template">[Main template]</a> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Disk usage: 492KB </span>

	
    <h3>Navigation</h3>
    <ul>
	<li class="right" ><a href= /shots/0/analysis/Basics/0612AdvancedAnalysis.ON/ title="current shot" >current</a> </li>
	<li class="right" ><a href= /shots/10105/analysis/Basics/0612AdvancedAnalysis.ON/ title="next shot" >next</a> |</li>
	<li class="right" ><a href= /shots/10103/analysis/Basics/0612AdvancedAnalysis.ON/ title="previous shot" >previous</a> | </li>
    </ul>
    </div>

    <div class="footer">
	&copy; Copyright 2012, <a href="mailto:svoboda@fjfi.cvut.cz">GOLEM Team</a>
    </div>

  </body>
</html>