% math symbols

\newcommand{\Fr}{F}                       %{F_{\scriptscriptstyle r}}
\newcommand{\Fc}{F_{\scriptscriptstyle c}}
\newcommand{\Ft}{\Psi}                    %{F_{\scriptscriptstyle t}}

\newcommand{\Er}{E_{\scriptscriptstyle r}}
\newcommand{\er}{E}                       %{e_{\scriptscriptstyle r}}
\newcommand{\eg}{e}                       %{e_{\scriptscriptstyle g}}
\newcommand{\et}{\varepsilon}             %{e_{\scriptscriptstyle t}}
\newcommand{\ei}{e_{\scriptscriptstyle ion}}

\newcommand{\Ne}{N_{\scriptscriptstyle e}}
\newcommand{\Pe}{P_{\scriptscriptstyle e}}
\newcommand{\Pg}{p}                       %{P_{\scriptscriptstyle g}}
\newcommand{\pr}{P}                       %{P_{\scriptscriptstyle r}}
\newcommand{\Pv}{P_{\scriptscriptstyle v}}
\newcommand{\Pt}{\pi}                     %{P_{\scriptscriptstyle t}}
\newcommand{\PQ}{P_{\scriptscriptstyle Q}}

\newcommand{\St}{\Sigma}                  %{S_{\scriptscriptstyle t}}
\newcommand{\Kt}{K}                       %{K_{\scriptscriptstyle t}}
\newcommand{\Mt}{M_{\scriptscriptstyle t}}
\newcommand{\Dt}{{D_t}}
\newcommand{\nt}{\nu_{\scriptscriptstyle t}}
\newcommand{\Ct}{C_{\scriptscriptstyle t}}
\newcommand{\Cr}{C_{\scriptscriptstyle r}}
\newcommand{\Teff}{T_{\scriptscriptstyle eff}}
\newcommand{\fG}{f_{\scriptscriptstyle G}}
\newcommand{\fr}{\zeta}                   %{f_r}}

\newcommand{\Me}{M_{\scriptscriptstyle ext}}
\newcommand{\sr}{\sigma_{\scriptscriptstyle \rho}}
\newcommand{\se}{\sigma_{\scriptscriptstyle e}}
\newcommand{\tr}{\tau_{\scriptscriptstyle r}}
\newcommand{\ellQ}{\ell_{\scriptscriptstyle Q}}
\newcommand{\nuQ}{\nu_{\scriptscriptstyle Q}}
\newcommand{\mq}{\mu_{\scriptscriptstyle Q}}
\newcommand{\CQ}{C_{\scriptscriptstyle Q}}
\newcommand{\lQ}{l_{\scriptscriptstyle Q}}
\newcommand{\cs}{c_{\scriptscriptstyle s}}
\newcommand{\cv}{c_{\scriptscriptstyle V}}
\newcommand{\cp}{c_{\scriptscriptstyle P}}
\newcommand{\ar}{a_{\scriptscriptstyle r}}
\newcommand{\ur}{u_{\scriptscriptstyle rel}}
\newcommand{\fe}{f_{\scriptscriptstyle E}}
\newcommand{\gE}{g_{\scriptscriptstyle E}}
\newcommand{\kb}{k_{\scriptscriptstyle B}}
\newcommand{\kP}{\kappa_{\scriptscriptstyle P}}
\newcommand{\kE}{\kappa_{\scriptscriptstyle E}}
\newcommand{\kF}{\chi_{\scriptscriptstyle F}}
\newcommand{\kR}{\chi_{\scriptscriptstyle Ross}}
\newcommand{\en}{\varepsilon_{\scriptscriptstyle N}}

\newcommand{\om}{\varpi}

\newcommand{\p}{\partial}
\newcommand{\overp}{\over\partial}

\newcommand{\eavh}{\hat{\mbox{\itbf e}}_{\alpha}}
\newcommand{\evh}{\hat{\mbox{\itbf e}}}

\newcommand{\eavhs}{\hat{\mbox{\itbfs e}}_{\alpha}}
\newcommand{\evhs}{\hat{\mbox{\itbfs e}}}

\newcommand{\faxt}{f_{\alpha}(\mbox{\itbf x}, \, t)}
\newcommand{\faxetl}{f_{\alpha}(\mbox{\itbf x} + 
            \hat{\mbox{\itbf e}}_{\alpha}, \, t + 1)}

\newcommand{\frgt}{f(\mbox{\itbf r},\, \Gamma,\, t)}

\newcommand{\Hkv}{H(\mbox{\itbf k})}

\newcommand{\kv}{\mbox{\itbf k}}
\newcommand{\kvh}{\hat{\mbox{\itbf k}}}
\newcommand{\kvs}{\mbox{\itbfs k}}
\newcommand{\kvhs}{\hat{\mbox{\itbfs k}}}

\newcommand{\Naxt}{N_{\alpha}(\mbox{\itbf x}, \, t) }
\newcommand{\Naxetl}{N_{\alpha}(\mbox{\itbf x} + 
            \hat{\mbox{\itbf e}}_{\alpha}, \, t + 1) }

\newcommand{\phiaxt}{\phi_{\alpha}(\mbox{\itbf x}, \, t)}
\newcommand{\phibxt}{\phi_{\beta}(\mbox{\itbf x}, \, t)}
\newcommand{\phiaxetl}{\phi_{\alpha}(\mbox{\itbf x} + 
            \hat{\mbox{\itbf e}}_{\alpha}, \, t + 1)}
\newcommand{\phiakt}{\phi_{\alpha}(\mbox{\itbf k}, \, t)}
\newcommand{\phibkt}{\phi_{\beta}(\mbox{\itbf k}, \, t)}
\newcommand{\phiaktl}{\phi_{\alpha}(\mbox{\itbf k}, \, t + 1)}
\newcommand{\phikt}{| \phi (\mbox{\itbf k}, \, t) \rangle}
\newcommand{\phiktl}{| \phi(\mbox{\itbf k}, \, t + 1) \rangle}
\newcommand{\phitkt}{\phi_{t}(\mbox{\itbf k}, \, t)}
\newcommand{\phitktl}{\phi_{t}(\mbox{\itbf k}, \, t + 1)}
\newcommand{\phipmkt}{\phi_{\pm}(\mbox{\itbf k}, \, t)}
\newcommand{\phipmktl}{\phi_{\pm}(\mbox{\itbf k}, \, t + 1)}

\newcommand{\pv}{\mbox{\itbf p}}
\newcommand{\pvp}{\mbox{\itbf p}^{\prime}}

\newcommand{\rhort}{\rho (\mbox{\itbf r},\, t)}

\newcommand{\rv}{\mbox{\itbf r}}

\newcommand{\sv}{\mbox{\itbf s}}
\newcommand{\svp}{\mbox{\itbf s}'}

\newcommand{\td}{\tilde{d}}

\newcommand{\temprt}{T(\mbox{\itbf r},\, t)}

\newcommand{\vvec}{\mbox{\itbf v}}
\newcommand{\vvecrt}{\mbox{\itbf v}(\mbox{\itbf r},\, t)}

\newcommand{\xiv}{\mbox{\boldmath $\xi$}}

\newcommand{\xv}{\mbox{\itbf x}}
\newcommand{\xvh}{\hat{\mbox{\itbf x}}}
\newcommand{\yvh}{\hat{\mbox{\itbf y}}}

% making plots for collision tables.

\newcommand{\stagmesh}{     
\bpt(00,80)
 \thinlines
 \put(060,15){\line( 0,2){60}}
 \put(105,41){X}
 \put(150,15){\line( 0,2){45}}
 \put(195,41){X}
 \put(240,15){\line( 0,2){60}}
 \put(295,41){X}
 
 \put(054,08){k-1}
 \put(101,30){k-1}
 \put(144,08){ k }
 \put(139,65){r,m,u}
 \put(192,30){ k }
 \put(192,55){e,P}
 \put(234,08){k+1}
 \put(291,30){k+1}
 \put(351,41){k}

 \thicklines
 \put(45,45){\vector( 1,0){300}}
 \put(45,45){\line  (-1,0){030}}
\ept}

\newcommand{\sodinit}{     
\bpt(00,135)
 \thinlines
 \put(030,15){\line( 1,0){300}}
 
 \put(090,090){P=1.0}
 \put(090,080){d=1.0}
 \put(270,050){P=0.1}
 \put(270,040){d=0.125}

 \put(024,03){0.0}
 \put(176,03){0.5}
 \put(324,03){1.0}

 \thicklines
 \put(030,015){\line( 0,2){100}}
 \put(180,025){\line( 0,2){090}}
 \put(330,015){\line( 0,2){010}}

 \put(030,115){\line( 1,0){150}}
 \put(180,025){\line( 1,0){150}}
\ept}

\newcommand{\woodinit}{     
\bpt(00,145)
 \thinlines
 \put(030,15){\line( 1,0){300}}
 
 \put(025,120){P=1000}
 \put(175,030){P=0}
 \put(175,020){d=1}
 \put(298,030){P=100}

 \put(024,03){0.0}
 \put(056,03){0.1}
 \put(296,03){0.9}
 \put(324,03){1.0}

 \thicklines
 \put(030,015){\line( 0,2){100}}
 \put(030,115){\line( 1,0){030}}
 \put(060,015){\line( 0,2){100}}

 \put(060,015){\line( 1,0){240}}

 \put(300,015){\line( 0,2){010}}
 \put(300,025){\line( 1,0){030}}
 \put(330,015){\line( 0,2){010}}
\ept}

\newcommand{\flowchart}{
\bpt(00,20)
 \thinlines
 \put(135,-010){\line(1, 1){010}}
 \put(135,-010){\line(1,-1){010}}
 \put(145,-020){\line(1, 0){071}}
 \put(145, 000){\line(1, 0){071}}
 \put(215,-020){\line(1, 1){010}}
 \put(215,-000){\line(1,-1){010}}
 \put(143,-014){initialize model}

 \put(045,-040){\line(1, 0){071}}
 \put(045,-060){\line(1, 0){071}}
 \put(045,-060){\line(0, 1){020}}
 \put(116,-060){\line(0, 1){020}}
 \put(053,-054){analytically}

 \put(240,-040){\line(1, 0){071}}
 \put(240,-060){\line(1, 0){071}}
 \put(240,-060){\line(0, 1){020}}
 \put(311,-060){\line(0, 1){020}}
 \put(248,-054){numerically}

 \put(045,-080){\line(1, 0){071}}
 \put(045,-120){\line(1, 0){071}}
 \put(045,-120){\line(0, 1){040}}
 \put(116,-120){\line(0, 1){040}}
 \put(047,-094){apply formula} 
 \put(062,-104){relax on}
 \put(060,-114){fixed grid}

 \put(240,-080){\line(1, 0){071}}
 \put(240,-120){\line(1, 0){071}}
 \put(240,-120){\line(0, 1){040}}
 \put(311,-120){\line(0, 1){040}}
 \put(244,-094){read from file}
 \put(257,-104){relax on}
 \put(255,-114){fixed grid}

 \put(135,-150){\line(1, 1){010}}
 \put(135,-150){\line(1,-1){010}}
 \put(145,-160){\line(1, 0){071}}
 \put(145,-140){\line(1, 0){071}}
 \put(215,-160){\line(1, 1){010}}
 \put(215,-140){\line(1,-1){010}}
 \put(145,-154){discontinuous?}

 \put(285,-154){yes}

 \put(240,-180){\line(1, 0){071}}
 \put(240,-210){\line(1, 0){071}}
 \put(240,-210){\line(0, 1){030}}
 \put(311,-210){\line(0, 1){030}}
 \put(255,-194){relax on}
 \put(245,-204){adaptive grid}

 \put(135,-260){\line(1, 0){091}}
 \put(135,-230){\line(1, 0){091}}
 \put(135,-230){\line(0,-1){030}}
 \put(226,-230){\line(0,-1){030}}
 \put(150,-244){model startup}
 \put(160,-254){complete}

 \put(135,-300){\line(1, 0){091}}
 \put(135,-280){\line(1, 0){091}}
 \put(135,-280){\line(0,-1){020}}
 \put(226,-280){\line(0,-1){020}}
 \put(150,-294){evolve model}

 \thicklines

 \put(080,-010){\line(1, 0){055}}
 \put(080,-010){\line(0,-1){030}}

 \put(225,-010){\line(1, 0){055}}
 \put(280,-010){\line(0,-1){030}}

 \put(080,-060){\line(0,-1){020}}

 \put(280,-060){\line(0,-1){020}}

 \put(080,-120){\line(0,-1){010}}
 \put(280,-120){\line(0,-1){010}}
 \put(080,-130){\line(1, 0){200}}
 \put(180,-130){\line(0,-1){010}}

 \put(180,-160){\line(0,-1){070}}

 \put(225,-150){\line(1, 0){055}}
 \put(280,-150){\line(0,-1){030}}

 \put(280,-210){\line(0,-1){010}}
 \put(180,-220){\line(1, 0){100}}

 \put(180,-260){\line(0,-1){020}}

\ept}

.