added count.seqs command and made some modifcations to the uchime code to allow it...

[mothur.git] / parsimonycommand.cpp

diff --git a/parsimonycommand.cpp b/parsimonycommand.cpp
index be64ec273ec8c62fcba9a7cd206955878f7e21ed..403b2b86d9addf1f9291e7f87fac63ef1755c002 100644 (file)
--- a/parsimonycommand.cpp
+++ b/parsimonycommand.cpp
@@ -74,6 +74,7 @@ ParsimonyCommand::ParsimonyCommand(string option)  {
                        
                //allow user to run help
                if(option == "help") { help(); abort = true; calledHelp = true; }
+               else if(option == "citation") { citation(); abort = true; calledHelp = true;}
                
                else {
                        vector<string> myArray = setParameters();
@@ -125,6 +126,10 @@ ParsimonyCommand::ParsimonyCommand(string option)  {
                        }
                        
                        m->runParse = true;
+                       m->Groups.clear();
+                       m->namesOfGroups.clear();
+                       m->Treenames.clear();
+                       m->names.clear();
                        
                        outputDir = validParameter.validFile(parameters, "outputdir", false);           if (outputDir == "not found"){  outputDir = ""; }
                        
@@ -188,6 +193,9 @@ int ParsimonyCommand::execute() {
                //randomtree will tell us if user had their own treefile or if they just want the random distribution
                //user has entered their own tree
                if (randomtree == "") { 
+                       
+                       m->setTreeFile(treefile);
+                       
                        if (groupfile != "") {
                                //read in group map info.
                                tmap = new TreeMap(groupfile);