<tool id="gmx_setup" name="GROMACS initial setup" version="@VERSION@">
    <description>of topology and GRO structure file</description>
    <macros>
        <import>macros.xml</import>
    </macros>

    <expand macro="requirements" />

    <command detect_errors="exit_code"><![CDATA[

        ln -s '$pdb_input' ./pdb_input.pdb &&
        gmx pdb2gmx -f ./pdb_input.pdb -o processed.gro -p topol.top -i posres.itp -water $water -ff $ff -${ignore_h}ignh &>> verbose.txt &&
        gmx editconf -f processed.gro -o newbox.gro -c -d $box_d -bt $box_type &>> verbose.txt

    ]]></command>

    <inputs>
        <!-- PDB2GMX VARS -->
        <param argument="pdb_input" type="data" format='pdb' label="PDB input file" help="Input structure"/>
        <param argument="water" type="select" label="Water model" help="Model for simulating water molecules in the solvent">
            <option value="tip3p" selected="true">TIP3P</option>
            <option value="tip4p">TIP4P</option>
            <option value="tips3p">TIPS3P</option>
            <option value="tip5p">TIP5P</option>
            <option value="spc">SPC</option>
            <option value="spce">SPC/E</option>
            <option value="none">None</option>
        </param>

        <param argument="ff" type="select" label="Force field" help="Specify a molecular force field">
            <option value="oplsaa" selected="true">OPLS/AA</option>
            <option value="gromos43a1">GROMOS43A1</option>
            <option value="amber96">AMBER96</option>
            <option value="gromos53a6">GROMOS53A6</option>
            <option value="amber99sb-ildn">AMBER99SB-ILDN</option>
            <option value="amber99sb">AMBER99SB</option>
            <option value="gromos53a5">GROMOS53A5</option>
            <option value="gromos43a2">GROMOS43A2</option>
            <option value="amberGS">AMBERGS</option>
            <option value="charmm27">CHARMM27</option>
            <option value="amber03">AMBER03</option>
            <option value="gromos54a7">GROMOS54A7</option>
            <option value="gromos45a3">GROMOS45A3</option>
            <option value="amber99">AMBER99</option>
            <option value="amber94">AMBER94</option>
        </param>

        <param argument="ignore_h" type="boolean" label="Ignore hydrogens" truevalue="" falsevalue="no" help="Ignore hydrogen atoms that are in the coordinate file." />
        
        <!-- EDITCONF VARS -->
        <param argument="box_d" type="float" label="Box dimensions in nanometers" value="0.0" min="0.0" max="10.0" help="Dimensions of the box in which simulation will take place" />
        <param argument="box_type" type="select" label="Box type"
            help="triclinic is a triclinic box, cubic is a rectangular box with all sides equal dodecahedron represents a rhombic dodecahedron and octahedron is a truncated octahedron. The last two are special cases of a triclinic box. The length of the three box vectors of the truncated octahedron is the shortest distance between two opposite hexagons. Relative to a cubic box with some periodic image distance, the volume of a dodecahedron with this same periodic distance is 0.71 times that of the cube, and that of a truncated octahedron is 0.77 times.">
            <option value="cubic">Rectangular box with all sides equal</option>
            <option value="triclinic">Triclinic</option>
            <option value="dodecahedron">represents a rhombic dodecahedron</option>
            <option value="octahedron">truncated octahedron</option>
        </param>

        <param name="capture_log" type="boolean" value="false" label="Generate detailed log" help="Generate detailed log information that can be summarized with ParseLog."/>

    </inputs>
    <outputs>
        <data name="output1" format="top" from_work_dir="topol.top"/>
        <data name="output2" format="gro" from_work_dir="newbox.gro"/>
        <data name="output3" format="itp" from_work_dir="posres.itp"/>
        <data name="report" format="txt" from_work_dir="verbose.txt">
            <filter>capture_log</filter>
        </data>
    </outputs>
    <tests>
        <test>
            <param name="pdb_input" value="1AKI.pdb" />
            <param name="water" value="spce" />
            <param name="ff" value="oplsaa" />
            <param name="ignore_h" value="no" />
            <param name="box_d" value="1.0" />
            <param name="box_type" value="cubic" />
            <output name="output1" file="topol.top" ftype="top" compare="diff" lines_diff="20"/>
            <output name="output2" file="newbox.gro" ftype="gro"/>
            <output name="output3" file="posres.itp" ftype="itp"/>
        </test>
    </tests>
    <help><![CDATA[

.. class:: infomark

**What it does**

This tool performs the initial setup prior to a GROMACS simulation. This entails producing a topology from an input structure using the pdb2gmx command and also defining a simulation box with editconf.

Please note that the tool will only successfully generate a topology for residues it recognizes (i.e. standard amino acids). If the structure contains other components such as ligands, these should be separately parameterized with the AmberTools or ACPYPE tools. In this case, make sure to also use an AMBER forcefield for this tool to ensure compatability with the ligand topology. 

_____

.. class:: infomark

**Input**

       - PDB file.
       - Water model and forcefield must be specified.
       - Parameters for the simulation box (dimensions and shape).
     
_____

        
.. class:: infomark

**Output**

       - GROMACS topology (TOP) file.
       - Position restraint (itp) file, which may be useful for system equilibration.
       - GRO structure file.

    ]]></help>

    <expand macro="citations" />
</tool>