 \begin{pspicture}(0,0)(12,8)
  % nodes
  \rput(4.5,6.5){\rnode{init}{\psframebox{\parbox{7cm}{
   \begin{itemize}
    \item Initial configuration: $31\times31\times31$ unit cells Si
    \item Periodic boundary conditions
    \item $T=450\, ^{\circ}C$
    \item Equilibration of $E_{kin}$ and $E_{pot}$ for $600\, fs$
   \end{itemize}
  }}}}
  \rput(4.5,4.5){\rnode{tc1}{\psframebox[fillstyle=solid,fillcolor=red]{
                             $T=450\pm 1\, ^{\circ}C$}}}
  \rput(7,3.5){\rnode{insert}{\psframebox[fillstyle=solid,fillcolor=red]{
                              \parbox{3cm}{
                              Insertion of 10 atoms\\
                              at random positions}}}}
  \rput(2,3.5){\rnode{adj1}{\psframebox[fillstyle=solid,fillcolor=red]{
                            \parbox{3.5cm}{
                            Adjusting temperature\\
                            for another $100\, fs$}}}}
  \rput(7,2.5){\rnode{nc}{\psframebox[fillstyle=solid,fillcolor=red]{
                          $N_{atoms}=6000$}}}
  \rput(4.5,2){\rnode{tc2}{\psframebox[fillstyle=solid,fillcolor=cyan]{
                           $T=T_{set}$}}}
  \rput(7,1){\rnode{td}{\psframebox[fillstyle=solid,fillcolor=cyan]{
                        $T_{set}:=T_{set}-1\, ^{\circ}C$}}}
  \rput(2,1){\rnode{adj2}{\psframebox[fillstyle=solid,fillcolor=cyan]{
                          \parbox{3.5cm}{
                          Adjusting temperature\\
                          for another $50\, fs$}}}}
  \rput(7,0){\rnode{tc3}{\psframebox[fillstyle=solid,fillcolor=cyan]{
                         $T_{set}=20\, ^{\circ}C$}}}
  \rput(10,0){\rnode{end}{\psframebox{End}}}
  % help nodes
  \rput(7,4.5){\pnode{tc1-h}}
  \rput(2,4.5){\pnode{tc1-hh}}
  \rput(4.5,2.5){\pnode{nc-h}}
  \rput(9,2.5){\pnode{nc-hh}}
  \rput(9,2){\pnode{tc2-h}}
  \rput(2,2){\pnode{tc2-hh}}
  \rput(4.5,0){\pnode{tc3-h}}
  % direct lines
  \ncline[]{->}{init}{tc1}
  \ncline[]{->}{adj1}{tc1}
  \ncline[]{->}{insert}{nc}
  \ncline[]{->}{adj2}{tc2}
  \ncline[]{->}{tc2}{td}
  \lput*{0}{yes}
  \ncline[]{->}{td}{tc3}
  \ncline[]{->}{tc3}{end}
  \lput*{0}{yes}
  % lines using help nodes
  \ncline[]{tc1}{tc1-h}
  \lput*{0}{yes}
  \ncline[]{->}{tc1-h}{insert}
  \ncline[]{tc1}{tc1-hh}
  \lput*{0}{no}
  \ncline[]{->}{tc1-hh}{adj1}
  \ncline[]{nc}{nc-h}
  \lput*{0}{no}
  \ncline[]{->}{nc-h}{tc1}
  \ncline[]{nc}{nc-hh}
  \ncline[]{-}{nc-hh}{tc2-h}
  \ncline[]{->}{tc2-h}{tc2}
  \lput*{0}{yes, {\footnotesize $T_{set}:=450\, ^{\circ}C$}}
  \ncline[]{tc2}{tc2-hh}
  \lput*{0}{no}
  \ncline[]{->}{tc2-hh}{adj2}
  \ncline[]{tc3}{tc3-h}
  \lput*{0}{no}
  \ncline[]{->}{tc3-h}{tc2}
  % insertion volumes
  \psframe[fillstyle=solid,fillcolor=white](9.5,1.3)(13.5,5.3)
  \psframe[fillstyle=solid,fillcolor=lightgray](10,1.8)(13,4.8)
  \psframe[fillstyle=solid,fillcolor=gray](10.5,2.3)(12.5,4.3)
  \rput(9.75,3){\pnode{ins1}}
  \rput(10.25,3.3){\pnode{ins2}}
  \rput(10.75,3.6){\pnode{ins3}}
  \ncline[]{-}{insert}{ins1}
  \ncline[]{-}{insert}{ins2}
  \ncline[]{-}{insert}{ins3}
  \psframe[fillstyle=solid,fillcolor=white](9.5,7.6)(13.5,8.1)
  \psframe[fillstyle=solid,fillcolor=lightgray](9.5,6.8)(13.5,7.3)
  \psframe[fillstyle=solid,fillcolor=gray](9.5,6)(13.5,6.5)
  \rput(11.5,7.85){{\tiny Simulation volume:
                          $31\times31\times31\, a^3_{Si}$}}
  \rput(11.5,7.05){{\tiny Volume of minimal SiC precipitation}}
  \rput(11.5,6.25){{\tiny Volume of necessary amount of Si}}
 \end{pspicture}
