tla.xmds

<?xml version="1.0"?>
<!--Two Level Atom Example simulation to illustrate a
cross propagating field-->

<simulation>

  <prop_dim>z</prop_dim>

  <globals>
    <![CDATA[
    const double g  = 1;
    const double t0 = 1;
    ]]>
  </globals>

  <field>
    <dimensions>   t    </dimensions>
    <lattice>     100   </lattice>
    <domains>  (-10,15) </domains>
    <samples>1 0</samples>

    <vector>
      <name>main</name>
      <type>double</type>
      <components>E</components>
      <![CDATA[
      E = 2/t0/cosh(t/t0);
      ]]>
    </vector>

    <vector>
      <name>cross</name>
      <type>double</type>
      <components>P N</components>
      <![CDATA[
      P =  0;
      N = -1;
      ]]>
    </vector>
  </field>

  <sequence>
    <integrate>
      <algorithm>RK4EX</algorithm>
      <interval>4</interval>
      <lattice>50</lattice>
      <samples>50 50</samples>
      <vectors>main cross</vectors>
      <![CDATA[
      dE_dz =  g*P;
      ]]>
      <cross_propagation>
	<vectors>cross</vectors>
	<prop_dim>t</prop_dim>
	<![CDATA[
	dP_dt =  E*N;
	dN_dt = -E*P;
	]]>
      </cross_propagation>
    </integrate>
  </sequence>

  <output>

    <group>
      <sampling>
	<lattice>50</lattice>
	<moments>pow_dens</moments>
	<![CDATA[
	pow_dens = E*E;
	]]>
      </sampling>
    </group>

    <group>
      <sampling>
	<vectors>main cross</vectors>
	<lattice>50</lattice>
	<moments>P_out N_out</moments>
	<![CDATA[
	P_out = P;
	N_out = N;
	]]>
      </sampling>
    </group>
  </output>
</simulation>