#include "plotting.h"
#include "analysis.h"


 // *****************************************************************
 // *                                                               *
 // *                           MAIN PROGRAM                        *
 // *                                                               *
 // *****************************************************************

const int A1g   =   1;
const int A1gp  =   2;
const int A1gpp =   3;
const int A1u   =   4;
const int A2g   =   5;
const int A2u   =   6;
const int B1g   =   7;
const int B1gp  =   8;
const int B1u   =   9;
const int B2g   =  10;
const int B2gp  =  11;
const int B2u   =  12;
const int Eg    =  13;
const int Egp   =  14;
const int Egpp  =  15;
const int Eu    =  16;
const int Eup   =  17;
const int Eupp  =  18;
const int Euppp =  19;


int main()
{
int run_choice,string_level,R_minor_tick_num=0,Rinv_minor_tick_num=0,
    E_minor_tick_num=0,gap_minor_tick_num=0;
double Rmin=0,Rmax=0,Rinvmin=0,Rinvmax=0,Emin=0,Emax=0,gapmax=0,gapmin=0,
       Zxi=0,xi=0,Zxi_err=0,string_dev_max=1.0,string_dev_min=-1.0,
       R_major_tick_spacing=0,Rinv_major_tick_spacing=0,E_major_tick_spacing=0,
       gap_major_tick_spacing=0;
charString beta,plotdir,datafile,gapdatafile,outfile,
           rundir,htmlfile,footerfile;
int max_line_length=8;

Nchannels=19;
channel_match.reset(Nchannels);
channel_label.reset(Nchannels);

channel_match(A1g  )="A1g";     channel_label(A1g  )="A\\s1g";
channel_match(A1gp )="A1g'";    channel_label(A1gp )="A\\m{1}\\s1g\\N\\M{1}\\x\\c\\\"\\C";
channel_match(A1gpp)="A1g''";   channel_label(A1gpp)="A\\m{1}\\s1g\\N\\M{1}\\x\\c\\\"\\\"\\C";
channel_match(A1u  )="A1u";     channel_label(A1u  )="A\\s1u";
channel_match(A2g  )="A2g";     channel_label(A2g  )="A\\s2g";
channel_match(A2u  )="A2u";     channel_label(A2u  )="A\\s2u";
channel_match(B1g  )="B1g";     channel_label(B1g  )="B\\s1g";
channel_match(B1gp )="B1g'";    channel_label(B1gp )="B\\m{1}\\s1g\\N\\M{1}\\x\\c\\\"\\C";
channel_match(B1u  )="B1u";     channel_label(B1u  )="B\\s1u";
channel_match(B2g  )="B2g";     channel_label(B2g  )="B\\s2g";
channel_match(B2gp )="B2g'";    channel_label(B2gp )="B\\m{1}\\s2g\\N\\M{1}\\x\\c\\\"\\C";
channel_match(B2u  )="B2u";     channel_label(B2u  )="B\\s2u";
channel_match(Eg   )="Eg";      channel_label(Eg   )="E\\sg";
channel_match(Egp  )="Eg'";     channel_label(Egp  )="E\\m{1}\\sg\\N\\M{1}\\x\\c\\\"\\C";
channel_match(Egpp )="Eg''";    channel_label(Egpp )="E\\m{1}\\sg\\N\\M{1}\\x\\c\\\"\\\"\\C";
channel_match(Eu   )="Eu";      channel_label(Eu   )="E\\su";
channel_match(Eup  )="Eu'";     channel_label(Eup  )="E\\m{1}\\su\\N\\M{1}\\x\\c\\\"\\C";
channel_match(Eupp )="Eu''";    channel_label(Eupp )="E\\m{1}\\su\\N\\M{1}\\x\\c\\\"\\\"\\C";
channel_match(Euppp)="Eu'''";   channel_label(Euppp)="E\\m{1}\\su\\N\\M{1}\\x\\c\\\"\\\"\\\"\\C";

subtitle="SU(3) in 4 dimensions";
plotdir="plots";


run_choice=0;
while ((run_choice<1)||(run_choice>3)){
   cout << "  Which run?"<<endl;
   cout << "   [1] Run A:  beta = 3.0, xi=3 "<<endl;
   cout << "   [2] Run B:  beta = 2.5, xi=5 "<<endl;
   cout << "   [3] Run C:  beta = 2.2, xi=5 "<<endl;
   read_int("  Run? ",run_choice,max_line_length);
   }

 // ****************************************************************

if (run_choice==1){

   rundir="/home1/STATIC_POTENTIAL_RESULTS/SU3_4D/B30_AR3/";
   datafile=rundir+"all_results.dat";
   gapdatafile=rundir+"all_gap_results.dat";
   footerfile="../B30_AR3.html";
   xi=3.0;
   Zxi=0.943;
   Zxi_err=0.010;
   beta="\\xb\\f{normal}=3.0;  \\xx\\f{}=3.0; ";
   beta=beta+"(18\\S2\\N\\x\\c4\\C24)\\x\\c4\\C54\\N\\f{} Lattice";
   Rmin=0.0;
   Rmax=10.0;
   R_major_tick_spacing=2.0;     R_minor_tick_num=3;  
   Rinvmin=0.0;
   Rinvmax=0.4;
   Rinv_major_tick_spacing=0.1;     Rinv_minor_tick_num=4;
   Emin=0.0;
   Emax=1.4;
   E_major_tick_spacing=0.5;     E_minor_tick_num=4;
   gapmin=0.0;
   gapmax=1.2;
   gap_major_tick_spacing=0.2;     gap_minor_tick_num=4;
   string_dev_min=-1.0;
   string_dev_max=0.3;
   
   }

 // ****************************************************************
 // *                                                              *
 // *                             RUN B                            *
 // *                                                              *
 // ****************************************************************

else if (run_choice==2){

   rundir="/home1/STATIC_POTENTIAL_RESULTS/SU3_4D/B25_AR5/";
   datafile=rundir+"all_results.dat";
   gapdatafile=rundir+"all_gap_results.dat";
   footerfile="../B25_AR5.html";
   xi=5.0;
   Zxi=0.964;
   Zxi_err=0.010;
   beta="\\xb\\f{normal}=2.5;  \\xx\\f{}=5.0; ";
   beta=beta+"(16\\S2\\N\\x\\c4\\C20)\\x\\c4\\C80\\N\\f{} Lattice";
   Rmin=0.0;
   Rmax=11.0;
   R_major_tick_spacing=4.0;     R_minor_tick_num=3;  
   Rinvmin=0.0;
   Rinvmax=0.22;
   Rinv_major_tick_spacing=0.1;     Rinv_minor_tick_num=4;
   Emin=0.1;
   Emax=1.0;
   E_major_tick_spacing=0.5;     E_minor_tick_num=4;
   gapmin=0.0;
   gapmax=0.6;
   gap_major_tick_spacing=0.1;     gap_minor_tick_num=4;
   string_dev_min=-1.0;
   string_dev_max=0.8;

   }
   
 // ****************************************************************
 // *                                                              *
 // *                             RUN C                            *
 // *                                                              *
 // ****************************************************************

else if (run_choice==3){

   rundir="/home1/STATIC_POTENTIAL_RESULTS/SU3_4D/B22_AR5/";
   datafile=rundir+"all_results.dat";
   gapdatafile=rundir+"all_gap_results.dat";
   footerfile="../B22_AR5.html";
   xi=5.0;
   Zxi=0.9541;
   Zxi_err=0.0077;
   beta="\\xb\\f{normal}=2.2;  \\xx\\f{}=5.0; ";
   beta=beta+"12\\S3\\N\\x\\c4\\C48\\N\\f{} Lattice";
   Rmin=0.0;
   Rmax=6.0;
   R_major_tick_spacing=2.0;     R_minor_tick_num=3;  
   Rinvmin=0.0;
   Rinvmax=0.4;
   Rinv_major_tick_spacing=0.1;     Rinv_minor_tick_num=4;
   Emin=0.0;
   Emax=1.2;
   E_major_tick_spacing=0.5;     E_minor_tick_num=4;
   gapmin=0.0;
   gapmax=0.8;
   gap_major_tick_spacing=0.1;     gap_minor_tick_num=4;
   string_dev_min=-1.0;
   string_dev_max=0.5;

   }
   
 
 // ***********************************************************


     // preliminaries over....now do the plots!!

label_strings.assign(1,beta,0.065,0.915);
in.open(convert_to_charstar(datafile));
if (!in){
    cerr << " could not open data file";
    exit(1);}

   // do the R plots of the levels first
   
Rplot_type=1;
title="Gauge field energy with static quark-antiquark pair";
ylabel="a\\st\\NE(R)";
xlabel="R/a\\ss";
xmin=Rmin;  xmax=Rmax;  ymin=Emin;  ymax=Emax;
x_major_tick_spacing=R_major_tick_spacing;
x_minor_tick_num=R_minor_tick_num;
y_major_tick_spacing=E_major_tick_spacing;
y_minor_tick_num=E_minor_tick_num;
html_header(rundir+plotdir+"/Rplots.html",convert_to_charstar(footerfile));

outfile=rundir+plotdir+"/A1g_plot";   // A1g_plot  (A1g, A1g', A1g'')
pump_out_plot(outfile,A1g,0.20,-1,   A1gp,0.03,-1,    A1gpp,0.1,1,  0);
outfile=rundir+plotdir+"/A1u_plot";   // A1u_plot  (A1u)
pump_out_plot(outfile,A1u,0.3,-1);
outfile=rundir+plotdir+"/A2g_plot";   // A2g_plot  (A2g)
pump_out_plot(outfile,A2g,0.3,-1);
outfile=rundir+plotdir+"/A2u_plot";   // A2u_plot  (A2u)
pump_out_plot(outfile,A2u,0.3,-1);
outfile=rundir+plotdir+"/B1g_plot";   // B1g_plot  (B1g, B1g')
pump_out_plot(outfile,B1g,0.1,-1,   B1gp,0.1,1,   0);
outfile=rundir+plotdir+"/B1u_plot";   // B1u_plot  (B1u)
pump_out_plot(outfile,B1u,0.3,-1);
outfile=rundir+plotdir+"/B2g_plot";   // B2g_plot  (B2g, B2g')
pump_out_plot(outfile,B2g,0.1,-1,   B2gp,0.1,1,   0);
outfile=rundir+plotdir+"/B2u_plot";   // B2u_plot  (B2u)
pump_out_plot(outfile,B2u,0.3,-1);
outfile=rundir+plotdir+"/Eg_plot1";   // Eg_plot  (Eg, Eg')
pump_out_plot(outfile,Eg,0.1,-1,   Egp,0.1,1,   0);
outfile=rundir+plotdir+"/Eg_plot2";   // Eg_plot  (Eg, Eg'')
pump_out_plot(outfile,Eg,0.1,-1,   Egpp,0.1,1,   1);
outfile=rundir+plotdir+"/Eu_plot1";   // Eu_plot  (Eu, Eu')
pump_out_plot(outfile,Eu,0.1,-1,   Eup,0.1,1,   0);
outfile=rundir+plotdir+"/Eu_plot2";   // Eu_plot  (Eu, Eu'')
pump_out_plot(outfile,Eu,0.1,-1,   Eupp,0.1,1,   1);
outfile=rundir+plotdir+"/Eu_plot3";   // Eu_plot  (Eu, Eu''')
pump_out_plot(outfile,Eu,0.1,-1,   Euppp,0.1,1,   1);

in.close();
html_footer(convert_to_charstar(footerfile));


      // do the R plots of the gaps above the ground state next

in.open(convert_to_charstar(gapdatafile));
if (!in){
    cerr << " could not open data file";
    exit(1);}
Rplot_type=1;
title="Gauge field energy gap with static quark-antiquark pair";
ylabel="a\\st\\N(E(R)-E\\sA\\s1g\\N(R))";
xlabel="R/a\\ss";
xmin=Rmin;  xmax=Rmax;  ymin=gapmin;  ymax=gapmax;
x_major_tick_spacing=R_major_tick_spacing;
x_minor_tick_num=R_minor_tick_num;
y_major_tick_spacing=gap_major_tick_spacing;
y_minor_tick_num=gap_minor_tick_num;
setup_goldstone(Rplot_type,xmin+0.55*(xmax-xmin),xmax,xi,Zxi);
html_header(rundir+plotdir+"/gap_Rplots.html",convert_to_charstar(footerfile));
           
outfile=rundir+plotdir+"/A1g_gap_plot";   // A1g_gap_plot  (A1g', A1g'')
pump_out_plot(outfile,A1gp,0.1,-1,   A1gpp,0.1,1,   0);
outfile=rundir+plotdir+"/A1u_gap_plot";   // A1u_gap_plot  (A1u)
pump_out_plot(outfile,A1u,0.1,-1);
outfile=rundir+plotdir+"/A2g_gap_plot";   // A2g_gap_plot  (A2g)
pump_out_plot(outfile,A2g,0.1,-1);
outfile=rundir+plotdir+"/A2u_gap_plot";   // A2u_gap_plot  (A2u)
pump_out_plot(outfile,A2u,0.1,-1);
outfile=rundir+plotdir+"/B1g_gap_plot";   // B1g_gap_plot  (B1g, B1g')
pump_out_plot(outfile,B1g,0.1,-1,   B1gp,0.1,1,   0);
outfile=rundir+plotdir+"/B1u_gap_plot";   // B1u_gap_plot  (B1u)
pump_out_plot(outfile,B1u,0.1,-1);
outfile=rundir+plotdir+"/B2g_gap_plot";   // B2g_gap_plot  (B2g, B2g')
pump_out_plot(outfile,B2g,0.1,-1,   B2gp,0.1,1,   0);
outfile=rundir+plotdir+"/B2u_gap_plot";   // B2u_gap_plot  (B2u)
pump_out_plot(outfile,B2u,0.1,-1);
outfile=rundir+plotdir+"/Bg_gap_plot";    // Bg_gap_plot  (B1g, B1g', B2g, B2g')
pump_out_plot(outfile,B1g,0.1,-1,   B1gp,0.3,-1,   B2g,0.15,1,   B2gp,0.35,1,   0);
outfile=rundir+plotdir+"/Bu_gap_plot";    // Bu_gap_plot  (B1u, B2u)
pump_out_plot(outfile,B1u,0.2,-1,   B2u,0.2,1,   0);
outfile=rundir+plotdir+"/Eg_gap_plot1";   // Eg_gap_plot  (Eg, Eg')
pump_out_plot(outfile,Eg,0.1,-1,   Egp,0.1,1,   0);
outfile=rundir+plotdir+"/Eg_gap_plot2";   // Eg_gap_plot  (Eg, Eg'')
pump_out_plot(outfile,Eg,0.1,-1,   Egpp,0.1,1,   1);
outfile=rundir+plotdir+"/Eg_gap_plot";    // Eg_gap_plot  (Eg, Eg', Eg'')
pump_out_plot(outfile,Eg,0.1,-1,   Egp,0.01,1,   Egpp,0.1,1,   1);
outfile=rundir+plotdir+"/Eu_gap_plot1";   // Eu_gap_plot  (Eu, Eu')
pump_out_plot(outfile,Eu,0.1,-1,   Eup,0.1,1,   0);
outfile=rundir+plotdir+"/Eu_gap_plot2";   // Eu_gap_plot  (Eu, Eu'')
pump_out_plot(outfile,Eu,0.1,-1,   Eupp,0.1,1,   1);
outfile=rundir+plotdir+"/Eu_gap_plot3";   // Eu_gap_plot  (Eu, Eu''')
pump_out_plot(outfile,Eu,0.1,-1,   Euppp,0.1,1,   1);
outfile=rundir+plotdir+"/Eu_gap_plot";    // Eu_gap_plot  (Eu, Eu', Eu'', Eu''')
pump_out_plot(outfile,Eu,0.1,-1,   Eup,0.2,-1,   Eupp,0.1,1,   Euppp,0.8,1,   1);

html_footer(convert_to_charstar(footerfile));


   // do the 1/R plots of the gaps above the ground state next
 
Rplot_type=0;
title="Gauge field energy gap with static quark-antiquark pair";
ylabel="a\\st\\N(E(R)-E\\sA\\s1g\\N(R))";
xlabel="a\\ss\\N/R";
xmin=Rinvmin;  xmax=Rinvmax;  ymin=gapmin;  ymax=gapmax;
x_major_tick_spacing=Rinv_major_tick_spacing;
x_minor_tick_num=Rinv_minor_tick_num;
y_major_tick_spacing=gap_major_tick_spacing;
y_minor_tick_num=gap_minor_tick_num;
setup_goldstone(Rplot_type,xmin,0.6*(xmax-xmin),xi,Zxi);
html_header(rundir+plotdir+"/gap_Rinv_plots.html",convert_to_charstar(footerfile));
           
outfile=rundir+plotdir+"/A1g_gap_Rinv_plot";   // A1g_gap_plot  (A1g', A1g'')
pump_out_plot(outfile,A1gp,0.7,-1,   A1gpp,0.7,-1,   0);
outfile=rundir+plotdir+"/A1u_gap_Rinv_plot";   // A1u_gap_plot  (A1u)
pump_out_plot(outfile,A1u,0.7,-1);
outfile=rundir+plotdir+"/A2g_gap_Rinv_plot";   // A2g_gap_plot  (A2g)
pump_out_plot(outfile,A2g,0.7,-1);
outfile=rundir+plotdir+"/A2u_gap_Rinv_plot";   // A2u_gap_plot  (A2u)
pump_out_plot(outfile,A2u,0.7,-1);
outfile=rundir+plotdir+"/B1g_gap_Rinv_plot";   // B1g_gap_plot  (B1g, B1g')
pump_out_plot(outfile,B1g,0.7,-1,   B1gp,0.7,-1,   0);
outfile=rundir+plotdir+"/B1u_gap_Rinv_plot";   // B1u_gap_plot  (B1u)
pump_out_plot(outfile,B1u,0.7,-1);
outfile=rundir+plotdir+"/B2g_gap_Rinv_plot";   // B2g_gap_plot  (B2g, B2g')
pump_out_plot(outfile,B2g,0.7,-1,   B2gp,0.7,-1,   0);
outfile=rundir+plotdir+"/B2u_gap_Rinv_plot";   // B2u_gap_plot  (B2u)
pump_out_plot(outfile,B2u,0.7,-1);
outfile=rundir+plotdir+"/Bg_gap_Rinv_plot";    // Bg_gap_plot  (B1g, B1g', B2g, B2g')
pump_out_plot(outfile,B1g,0.7,-1,   B1gp,0.7,-1,   B2g,0.7,1,   B2gp,0.7,1,   0);
outfile=rundir+plotdir+"/Bu_gap_Rinv_plot";    // Bu_gap_plot  (B1u, B2u)
pump_out_plot(outfile,B1u,0.7,-1,   B2u,0.7,1,  0);
outfile=rundir+plotdir+"/Eg_gap_Rinv_plot1";   // Eg_gap_plot  (Eg, Eg')
pump_out_plot(outfile,Eg,0.7,-1,   Egp,0.7,1,   0);
outfile=rundir+plotdir+"/Eg_gap_Rinv_plot2";   // Eg_gap_plot  (Eg, Eg'')
pump_out_plot(outfile,Eg,0.7,-1,   Egpp,0.7,-1,   1);
outfile=rundir+plotdir+"/Eg_gap_Rinv_plot";    // Eg_gap_plot  (Eg, Eg', Eg'')
pump_out_plot(outfile,Eg,0.7,-1,   Egp,0.7,-1,   Egpp,0.7,1,   1);
outfile=rundir+plotdir+"/Eu_gap_Rinv_plot1";   // Eu_gap_plot  (Eu, Eu')
pump_out_plot(outfile,Eu,0.7,-1,   Eup,0.7,-1,   0);
outfile=rundir+plotdir+"/Eu_gap_Rinv_plot2";   // Eu_gap_plot  (Eu, Eu'')
pump_out_plot(outfile,Eu,0.7,-1,   Eupp,0.7,-1,   1);
outfile=rundir+plotdir+"/Eu_gap_Rinv_plot3";   // Eu_gap_plot  (Eu, Eu''')
pump_out_plot(outfile,Eu,0.7,-1,   Euppp,0.7,-1,   1);
outfile=rundir+plotdir+"/Eu_gap_Rinv_plot";    // Eu_gap_plot  (Eu, Eu', Eu'', Eu''')
pump_out_plot(outfile,Eu,0.7,-1,   Eup,0.7,-1,   Eupp,-0.1,0,   Euppp,0.7,1,   1);

html_footer(convert_to_charstar(footerfile));


   // do the R plots of the gaps above the ground state compared to N*Pi/R

Rplot_type=1;
label_strings.reset(2);
label_strings.assign(1,beta,0.065,0.915);
title="Gauge field energy gap with static quark-antiquark pair";
ylabel="\\xD\\f{}E(R) / (N\\xp\\f{}/R) - 1";
xlabel="R/a\\ss";
xmin=Rmin;  xmax=Rmax;  ymin=string_dev_min;  ymax=string_dev_max;
x_major_tick_spacing=R_major_tick_spacing;
x_minor_tick_num=R_minor_tick_num;
y_major_tick_spacing=0.2;
y_minor_tick_num=3;
html_header(rundir+plotdir+"/string_plots.html",convert_to_charstar(footerfile));

string_level=1;
label_strings.assign(2,"N="+int_to_string(string_level),0.15,0.85);
outfile=rundir+plotdir+"/Eu_string_plot"; 
pump_out_string_plot(outfile,  Eu, 0.4, -1,   string_level,xi,Zxi,Zxi_err);
string_level=2;
label_strings.assign(2,"N="+int_to_string(string_level),0.15,0.85);
outfile=rundir+plotdir+"/A1gp_string_plot"; 
pump_out_string_plot(outfile,  A1gp, 0.4, -1,   string_level,xi,Zxi,Zxi_err);
outfile=rundir+plotdir+"/B1g_string_plot"; 
pump_out_string_plot(outfile,  B1g, 0.4, -1,   string_level,xi,Zxi,Zxi_err);
outfile=rundir+plotdir+"/B2g_string_plot"; 
pump_out_string_plot(outfile,  B2g, 0.4, -1,   string_level,xi,Zxi,Zxi_err);
outfile=rundir+plotdir+"/Eg_string_plot"; 
pump_out_string_plot(outfile,  Eg, 0.4, -1,   string_level,xi,Zxi,Zxi_err);
string_level=3;
label_strings.assign(2,"N="+int_to_string(string_level),0.15,0.85);
outfile=rundir+plotdir+"/A1u_string_plot"; 
pump_out_string_plot(outfile,  A1u, 0.4, -1,   string_level,xi,Zxi,Zxi_err);
outfile=rundir+plotdir+"/A2u_string_plot"; 
pump_out_string_plot(outfile,  A2u, 0.4, -1,   string_level,xi,Zxi,Zxi_err);
outfile=rundir+plotdir+"/B1u_string_plot"; 
pump_out_string_plot(outfile,  B1u, 0.4, -1,   string_level,xi,Zxi,Zxi_err);
outfile=rundir+plotdir+"/B2u_string_plot"; 
pump_out_string_plot(outfile,  B2u, 0.4, -1,   string_level,xi,Zxi,Zxi_err);
outfile=rundir+plotdir+"/Eup_string_plot"; 
pump_out_string_plot(outfile,  Eup, 0.4, -1,   string_level,xi,Zxi,Zxi_err);
outfile=rundir+plotdir+"/Eupp_string_plot"; 
pump_out_string_plot(outfile,  Eupp, 0.4, -1,   string_level,xi,Zxi,Zxi_err);
outfile=rundir+plotdir+"/Euppp_string_plot"; 
pump_out_string_plot(outfile,  Euppp, 0.4, -1,   string_level,xi,Zxi,Zxi_err);
string_level=4;
label_strings.assign(2,"N="+int_to_string(string_level),0.15,0.85);
outfile=rundir+plotdir+"/A1gpp_string_plot"; 
pump_out_string_plot(outfile,  A1gpp, 0.4, -1,   string_level,xi,Zxi,Zxi_err);
outfile=rundir+plotdir+"/A2g_string_plot"; 
pump_out_string_plot(outfile,  A2g, 0.4, -1,   string_level,xi,Zxi,Zxi_err);
outfile=rundir+plotdir+"/B1gp_string_plot"; 
pump_out_string_plot(outfile,  B1gp, 0.4, -1,   string_level,xi,Zxi,Zxi_err);
outfile=rundir+plotdir+"/B2gp_string_plot"; 
pump_out_string_plot(outfile,  B2gp, 0.4, -1,   string_level,xi,Zxi,Zxi_err);
outfile=rundir+plotdir+"/Egp_string_plot"; 
pump_out_string_plot(outfile,  Egp, 0.4, -1,   string_level,xi,Zxi,Zxi_err);
outfile=rundir+plotdir+"/Egpp_string_plot"; 
pump_out_string_plot(outfile,  Egpp, 0.4, -1,   string_level,xi,Zxi,Zxi_err);


html_footer(convert_to_charstar(footerfile));
in.close(); 


 return 0;
}