/*---------------------------------------------------------------------------*/
/*gccPSHD-gnuplot.c*/
/*---------------------------------------------------------------------------*/
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#define Nmax 100

void main_part(char fnameR[],char fnameW[]);
void Trim(char *s);
/*---------------------------------------------------------------------------*/

int main(int argc,char *argv[])
{
    char fnameR[50]; /*入力csvファイル名*/
    char fnameW[50]; /*出力txtファイル名*/

    strcpy(fnameR,argv[1]);
    strcpy(fnameW,argv[2]);

    main_part(fnameR,fnameW);

    return 0;
}
/*---------------------------------------------------------------------------*/

void main_part(char fnameR[],char fnameW[])
{
    FILE *fin,*fout;
    char dat[256];
    char str[100],str0[50],str1[50],str2[50],str3[50],str4[50],str5[50];

    char fnameF[50]; /*出力epsファイル名*/
    char sxjiku[100],syjiku[100];
    double xmin,xmax,dx,ymin,ymax,dy,xx,yy,dlx,dly;
    int nd,i,k;

    char ps_name[Nmax][50],ps_type[Nmax][50],st_class[Nmax][50],tbrl[Nmax][50];
    double co_YEAR[Nmax],ps_HD[Nmax];
    int mark[Nmax];

    /***********************************************/
    /*データ入力*/
    /***********************************************/
    fin=fopen(fnameR,"r");
    fgets(dat,sizeof dat,fin);
        strcpy(fnameF,strtok(dat,","));
    fgets(dat,sizeof dat,fin);
        strcpy(sxjiku,strtok(dat,","));
        strcpy(syjiku,strtok(NULL,","));
    fgets(dat,sizeof dat,fin);
        xmin=atof(strtok(dat,","));
        xmax=atof(strtok(NULL,","));
        dx  =atof(strtok(NULL,","));
    fgets(dat,sizeof dat,fin);
        ymin=atof(strtok(dat,","));
        ymax=atof(strtok(NULL,","));
        dy  =atof(strtok(NULL,","));
    fgets(dat,sizeof dat,fin);
        nd=atoi(strtok(dat,","));
    fgets(dat,sizeof dat,fin);
    k = 0;
    for (i=0;i<=nd-1;i++){
        fgets(dat,sizeof dat,fin);
            strcpy(str0,strtok(dat,",")) ;Trim(str0);
            strcpy(str1,strtok(NULL,","));Trim(str1);
            strcpy(str2,strtok(NULL,","));Trim(str2);
            strcpy(str3,strtok(NULL,","));Trim(str3);
            strcpy(str4,strtok(NULL,","));Trim(str4);
            strcpy(str5,strtok(NULL,","));Trim(str5);
        if(strstr(str5,"0")==NULL){
            strcpy(ps_name[k],str0);
            strcpy(ps_type[k],str1);
            co_YEAR[k]=atof(str2);
            ps_HD[k] =atof(str3);
            strcpy(st_class[k],str4);
            strcpy(tbrl[k],str5);
            if(0<strstr(st_class[k],"400MPa") && 0<strstr(ps_type[k],"C"))mark[k]=0;
            if(0<strstr(st_class[k],"490MPa") && 0<strstr(ps_type[k],"C"))mark[k]=1;
            if(0<strstr(st_class[k],"490MPa") && 0<strstr(ps_type[k],"P"))mark[k]=2;
            if(0<strstr(st_class[k],"570MPa") && 0<strstr(ps_type[k],"C"))mark[k]=3;
            if(0<strstr(st_class[k],"570MPa") && 0<strstr(ps_type[k],"P"))mark[k]=4;
            if(0<strstr(st_class[k],"690MPa") && 0<strstr(ps_type[k],"C"))mark[k]=5;
            if(0<strstr(st_class[k],"690MPa") && 0<strstr(ps_type[k],"P"))mark[k]=6;
            if(0<strstr(st_class[k],"780MPa") && 0<strstr(ps_type[k],"P"))mark[k]=7;
            if(0<strstr(st_class[k],"950MPa") && 0<strstr(ps_type[k],"P"))mark[k]=8;
            k=k+1;
        }
    }
    fclose(fin);
    nd = k;

    /***********************************************/
    /*gnuplot制御用データ書き込み*/
    /***********************************************/
    fout=fopen(fnameW,"w");
    fprintf(fout,"set terminal postscript eps font \"Arial Narrow\" 14\n");
    fprintf(fout,"set output \"%s\"\n",fnameF);
    fprintf(fout,"set key left top\n");
    fprintf(fout,"set key width -10\n");
    fprintf(fout,"set key box\n");
    fprintf(fout,"set key title \"Tensile strength and Type of P/S\"\n");
    fprintf(fout,"set key Left top\n");
    fprintf(fout,"set key spacing 1.2\n");
    fprintf(fout,"set xlabel \"%s\"\n",sxjiku);
    fprintf(fout,"set ylabel \"%s\"\n",syjiku);
    fprintf(fout,"set xrange [%.0f:%.0f]\n",xmin,xmax);
    fprintf(fout,"set xtics %.0f\n",dx);
    fprintf(fout,"set format x \"%s\"\n","%.0f");
    fprintf(fout,"set yrange [%.0f:%.0f]\n",ymin,ymax);
    fprintf(fout,"set ytics %.0f\n",dy);
    fprintf(fout,"set format y \"%s\"\n","%.0f");
    dlx=0.0;dly=0.0;
    strcpy(str,"");
    for(i=0;i<=nd-1;i++){
        if(0<strstr(tbrl[i],"t")){dlx= 0.0;dly= 150.0;strcpy(str,"center");}
        if(0<strstr(tbrl[i],"b")){dlx= 0.0;dly=-120.0;strcpy(str,"center");}
        if(0<strstr(tbrl[i],"r")){dlx=-0.7;dly=   0.0;strcpy(str,"right");}
        if(0<strstr(tbrl[i],"l")){dlx= 0.7;dly=   0.0;strcpy(str,"left");}
        xx=co_YEAR[i]+dlx;
        yy=ps_HD[i]+dly;
        fprintf(fout,"set label \"%s\" at %.0f,%.0f %s font \"Arial Narrow,12\"\n",ps_name[i],xx,yy,str);
    }
    /*軸作成*/
    fprintf(fout,"set grid xtics ytics mxtics mytics\n");
    fprintf(fout,"plot \\\n");
    fprintf(fout,"\"-\" with points pointtype 12 pointsize 1.0 title \" 400MPa: Conventional  \",\\\n");
    fprintf(fout,"\"-\" with points pointtype  6 pointsize 1.0 title \" 490MPa: Conventional  \",\\\n");
    fprintf(fout,"\"-\" with points pointtype  7 pointsize 1.0 title \" 490MPa: Pumped-storage\",\\\n");
    fprintf(fout,"\"-\" with points pointtype  4 pointsize 1.0 title \" 570MPa: Conventional  \",\\\n");
    fprintf(fout,"\"-\" with points pointtype  5 pointsize 1.0 title \" 570MPa: Pumped-storage\",\\\n");
    fprintf(fout,"\"-\" with points pointtype  8 pointsize 1.0 title \" 690MPa: Conventional  \",\\\n");
    fprintf(fout,"\"-\" with points pointtype  9 pointsize 1.0 title \" 690MPa: Pumped-storage\",\\\n");
    fprintf(fout,"\"-\" with points pointtype 21 pointsize 1.3 title \" 780MPa: Pumped-storage\",\\\n");
    fprintf(fout,"\"-\" with points pointtype 37 pointsize 1.3 title \" 950MPa: Pumped-storage\n");
    for(k=0;k<=8;k++){
        for(i=0;i<=nd-1;i++){
            if(mark[i]==k){
                fprintf(fout,"%.0f %.0f  #%s\n",co_YEAR[i],ps_HD[i],ps_name[i]);
            }
        }
        fprintf(fout,"e\n");
    }
    fclose(fout);
}
/*---------------------------------------------------------------------------*/
void Trim(char *s)
{
    int i;

    /*末尾から空白でない位置を捜し終端文字を付加*/
    /*strlenは終端文字は数えない*/
    /*配列s[]の添字の範囲は0からstrlen(s)-1*/
    for(i=strlen(s)-1;i>=0;i--){
        if(s[i]!=' ')break;
    }
    s[i+1]='\0';

    /*先頭から空白でない位置を捜しその位置以降をstrcpy*/
    for(i=0;i<=strlen(s);i++){
        if(s[i]!=' ')break;
    }
    strcpy(s,&s[i]);
}
/*---------------------------------------------------------------------------*/