using System;
using System.IO;

public class relaxxml
{
    //Extract model free parameters from NMR relax xml files
    //V.V. 03/01/2011. vvostri@gmail.com
    public static void Main()
    {
        string inputfile = "results";
        string outputfile = "variable initialized";
        StreamReader reader = new StreamReader(inputfile);
        string lineinfile;
        string[] splitline;

        lineinfile = reader.ReadLine();
        int counter = 0;

        while (lineinfile != null)
        {
            if (lineinfile.Contains("pipe desc"))
            {
                splitline = lineinfile.Split(new string[] {"\""}, StringSplitOptions.None);
                outputfile = splitline[3];
            }
            
            if (lineinfile.Contains("Residue"))
            {
                counter++;
            }
            lineinfile = reader.ReadLine();
        }
        reader.Close();

        //AA, Model, Tm, Tm_err, S2, S2_err, Te, Te_err, S2s, S2s_err, Ts, Ts_err, S2f, S2f_err, Tf, Tf_err, Rex, Rex_err, Chi2
        double[,] results = new double[counter, 19];

        reader = File.OpenText(inputfile);
        lineinfile = reader.ReadLine();
        counter = -1;
        
        while (lineinfile != null)
        {
            if (lineinfile.Contains("Residue"))
            {
                splitline = lineinfile.Split(new string[] {"\""}, StringSplitOptions.None);
                counter++;
                results[counter, 0] = Convert.ToDouble(splitline[5]);
            }

            if (lineinfile.Contains("<model desc=\"The model\""))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 1] = 0;
                }
                else
                {
                    splitline = lineinfile.Split(new string[] { "\'", "m"}, StringSplitOptions.None);
                    results[counter, 1] = Convert.ToDouble(splitline[2]);
                }
            }

            if (lineinfile.Contains("<local_tm desc"))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 2] = 0;
                }
                else
                {
                    results[counter, 2] = Math.Round(Convert.ToDouble(lineinfile) * 1e9, 3);
                }
            }

            if (lineinfile.Contains("<local_tm_err"))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 3] = 0;
                }
                else
                {
                    results[counter, 3] = Math.Round(Convert.ToDouble(lineinfile) * 1e9, 3);
                }
            }

            if (lineinfile.Contains("<s2 desc"))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 4] = 0;
                }
                else
                {
                    results[counter, 4] = Math.Round(Convert.ToDouble(lineinfile), 3);
                }
            }

            if (lineinfile.Contains("<s2_err"))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 5] = 0;
                }
                else
                {
                    results[counter, 5] = Math.Round(Convert.ToDouble(lineinfile), 3);
                }
            }

            if (lineinfile.Contains("te desc"))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 6] = 0;
                }
                else
                {
                    results[counter, 6] = Math.Round(Convert.ToDouble(lineinfile) * 1e12, 3);
                }
            }

            if (lineinfile.Contains("<te_err"))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 7] = 0;
                }
                else
                {
                    results[counter, 7] = Math.Round(Convert.ToDouble(lineinfile) * 1e12, 3);
                }
            }

            if (lineinfile.Contains("<s2s desc"))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 8] = 0;
                }
                else
                {
                    results[counter, 8] = Math.Round(Convert.ToDouble(lineinfile), 3);
                }
            }

            if (lineinfile.Contains("<s2s_err"))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 9] = 0;
                }
                else
                {
                    results[counter, 9] = Math.Round(Convert.ToDouble(lineinfile), 3);
                }
            }


            if (lineinfile.Contains("ts desc"))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 10] = 0;
                }
                else
                {
                    results[counter, 10] = Math.Round(Convert.ToDouble(lineinfile) * 1e9, 3);
                }
            }

            if (lineinfile.Contains("<ts_err"))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 11] = 0;
                }
                else
                {
                    results[counter, 11] = Math.Round(Convert.ToDouble(lineinfile) * 1e9, 3);
                }
            }

            if (lineinfile.Contains("<s2f desc"))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 12] = 0;
                }
                else
                {
                    results[counter, 12] = Math.Round(Convert.ToDouble(lineinfile), 3);
                }
            }

            if (lineinfile.Contains("<s2f_err"))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 13] = 0;
                }
                else
                {
                    results[counter, 13] = Math.Round(Convert.ToDouble(lineinfile), 3);
                }
            }

            if (lineinfile.Contains("tf desc"))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 14] = 0;
                }
                else
                {
                    results[counter, 14] = Math.Round(Convert.ToDouble(lineinfile) * 1e12, 3);
                }
            }

            if (lineinfile.Contains("<tf_err"))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 15] = 0;
                }
                else
                {
                    results[counter, 15] = Math.Round(Convert.ToDouble(lineinfile) * 1e12, 3);
                }
            }

            if (lineinfile.Contains("rex desc"))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 16] = 0;
                }
                else
                {
                    results[counter, 16] = Math.Round(Convert.ToDouble(lineinfile) * Math.Pow(2 * Math.PI * 500e6, 2),3);
                }
            }

            //Conversion to field-dependent value
            if (lineinfile.Contains("<rex_err"))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 17] = 0;
                }
                else
                {
                    results[counter, 17] = Math.Round(Convert.ToDouble(lineinfile) * Math.Pow(2 * Math.PI * 500e6, 2), 3);
                }
            }

            if (lineinfile.Contains("chi2 desc"))
            {
                lineinfile = reader.ReadLine();
                if (lineinfile.Contains("None"))
                {
                    results[counter, 18] = 0;
                }
                else
                {
                    results[counter, 18] = Convert.ToDouble(lineinfile);
                }
            }

            lineinfile = reader.ReadLine();
        }
        reader.Close();

        FileStream file = new FileStream(outputfile + ".csv", FileMode.Create, FileAccess.Write);
        StreamWriter sw = new StreamWriter(file);

        //AA, Model, Tm, Tm_err, S2, S2_err, Te, Te_err, S2s, S2s_err, Ts, Ts_err, S2f, S2f_err, Tf, Tf_err, Rex, Rex_err, Chi2
        sw.WriteLine("Rex values are for 500 MHz. Conversion to a different field is as below:");
        sw.WriteLine("Rex_1 / ((2 * Pi * Freq_1)^2) = Rex_2 / ((2 * Pi * Freq_2)^2)\n");
        sw.WriteLine("AA,Model,Tm,dTm,S2,dS2,Te,dTe,S2s,dS2s,Ts,dTs,S2f,dS2f,Tf,dTf,Rex,dRex,Chi2");
        sw.WriteLine("#,#,ns,ns,n/a,n/a,ps,ps,n/a,n/a,ns,ns,n/a,n/a,ps,ps,1/s,1/s,n/a");

        for (int index = 0; index != counter + 1; index++)
        {
            for (int index2 = 0; index2 != 19; index2++)
            {
                sw.Write(results[index, index2] + ",");                
            }
            sw.Write("\r\n");
        }
        sw.Close();
        file.Close();

    }
}