/*-----------------------------*/
/* プログラム：gccGRPARA-png.c */
/*-----------------------------*/
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#define MM 10
/*gnuplot用構造体宣言*/
/*scom,sxjiku,syjiku :タイトル(凡例),x軸ラベル,ｙ軸ラベル*/
/*slogx,xmin,xmax,dx,sxfmt :x軸種,x軸最小値,x軸最大値,x軸増分,x軸書式*/
/*slogy,ymin,ymax,dy,syfmt :y軸種,y軸最小値,y軸最大値,y軸増分,y軸書式*/
struct GPdata{
    char scom[100];char sxjiku[100];char syjiku[100];
    char slogx[3];double xmin;double xmax;double dx;char sxfmt[20];
    char slogy[3];double ymin;double ymax;double dy;char syfmt[20];
};
/*関数プロトタイプ*/
void main_part(char fnameR[50],char fnameW[50]);
void GPLOT(char fnameW[50],int mcol,struct GPdata gp,char legend[MM][20]);
/*---------------------------------------------------------------------------*/

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

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

    main_part(fnameR,fnameW);

    printf("リターンキーを押してください");
    getchar();
    return 0;
}
/*---------------------------------------------------------------------------*/

void main_part(char fnameR[50],char fnameW[50])
{
    char dat[256]="",sv[100]="";
    int i,j,ndata,mcol;
    double datax,datay[MM];
    struct GPdata gp={
        "","","",
        "",0.0,0.0,0.0,"",
        "",0.0,0.0,0.0,""
    };
    char legend[MM][20]={"","","","","","","","","",""};/*凡例*/
    FILE *fin,*fout,*ftemp;

    /*データ入力*/
    fin=fopen(fnameR,"r");
    fgets(dat,sizeof dat,fin);
        strcpy(sv,strtok(dat,","));
        strncat(gp.scom,sv,strlen(sv)-1);
    fgets(dat,sizeof dat,fin);
        strcpy(gp.sxjiku,strtok(dat,","));
        strcpy(sv,strtok(NULL,","));
        strncpy(gp.syjiku,sv,strlen(sv)-1);
    fgets(dat,sizeof dat,fin);
        strcpy(gp.slogx,strtok(dat,","));
        gp.xmin=atof(strtok(NULL,","));
        gp.xmax=atof(strtok(NULL,","));
        if(strncmp(gp.slogx,"N",1)==0){
            gp.dx=atof(strtok(NULL,","));
            strcpy(sv,strtok(NULL,","));
            strncpy(gp.sxfmt,sv,strlen(sv)-1);
        }
    fgets(dat,sizeof dat,fin);
        strcpy(gp.slogy,strtok(dat,","));
        gp.ymin=atof(strtok(NULL,","));
        gp.ymax=atof(strtok(NULL,","));
        if(strncmp(gp.slogy,"N",1)==0){
            gp.dy=atof(strtok(NULL,","));
            strcpy(sv,strtok(NULL,","));
            strncpy(gp.syfmt,sv,strlen(sv)-1);
        }
    fgets(dat,sizeof dat,fin);
        ndata=atoi(strtok(dat,","))-1;
        mcol=atoi(strtok(NULL,","));
    fgets(dat,sizeof dat,fin);
        strtok(dat,",");
        for(j=1;j<=mcol;j++){
            strcpy(sv,strtok(NULL,","));
            if(j==mcol){
                strncpy(legend[j-1],sv,strlen(sv)-1);
            }else{
                strcpy(legend[j-1],sv);
            }
        }
    /*数値データの読み込みと一時ファイル書き出し*/
    ftemp=fopen("temp.prn","w");
    for(i=0;i<=ndata;i++){
        fgets(dat,sizeof dat,fin);
            datax=atof(strtok(dat,","));
            if(strncmp(gp.slogx,"L",1)==0){
                if(datax<gp.xmin)datax=gp.xmin;
            }
            for(j=1;j<=mcol;j++){
                datay[j-1]=atof(strtok(NULL,","));
                if(strncmp(gp.slogy,"L",1)==0){
                    if(datay[j-1]<gp.ymin)datay[j-1]=gp.ymin;
                }
            }
        sprintf(sv,"%g",datax);strcpy(dat,sv);
        for(j=1;j<=mcol;j++){
            sprintf(sv,"%g",datay[j-1]);strcat(dat," ");strcat(dat,sv);
        }
        strcat(dat,"\n");
        fputs(dat,ftemp);
    }
    fclose(ftemp);
    fclose(fin);


    GPLOT(fnameW,mcol,gp,legend);

}
/*---------------------------------------------------------------------------*/

void GPLOT(char fnameW[50],int mcol,struct GPdata gp,char legend[MM][20])
{
    FILE *fp;
    int j;

    /****************************/
    /*gnuplot制御用ファイル作成*/
    /****************************/
    fp=fopen("gpl_temp.txt","w");
        /*出力ファイル定義*/
        fprintf(fp,"set terminal png\n");
        fprintf(fp,"set output \"%s\"\n",fnameW);
        /*グラフタイトル定義*/
        fprintf(fp,"set title \"%s\"\n",gp.scom);
        /*軸ラベル定義*/
        fprintf(fp,"set xlabel \"%s\"\n",gp.sxjiku);
        fprintf(fp,"set ylabel \"%s\"\n",gp.syjiku);
        /*ｘ軸定義*/
        if(strncmp(gp.slogx,"L",1)==0){
            fprintf(fp,"set logscale x\n");
            fprintf(fp,"set mxtics 10\n");
        }
        if(strncmp(gp.slogx,"N",1)==0){
            fprintf(fp,"set xtics %f\n",gp.dx);
            fprintf(fp,"set format x %s\n",gp.sxfmt);
        }
        /*ｙ軸定義*/
        if(strncmp(gp.slogy,"L",1)==0){
            fprintf(fp,"set logscale y\n");
            fprintf(fp,"set mytics 10\n");
        }
        if(strncmp(gp.slogy,"N",1)==0){
            fprintf(fp,"set ytics %f\n",gp.dy);
            fprintf(fp,"set format y %s\n",gp.syfmt);
        }
        /*軸範囲定義*/
        fprintf(fp,"set xrange [%f:%f]\n",gp.xmin,gp.xmax);
        fprintf(fp,"set yrange [%f:%f]\n",gp.ymin,gp.ymax);
        /*軸グリッド描画設定*/
        fprintf(fp,"set grid xtics ytics mxtics mytics\n");
        /*プロット実行（プロットは各種定義が終わって最後に実行）*/
        fprintf(fp,"plot \\\n");
        for(j=1;j<=mcol-1;j++){
            fprintf(fp,"\"temp.prn\" using 1:%d with lines linewidth 2 title \"%s\", \\\n",j+1,legend[j-1]);
        }
        fprintf(fp,"\"temp.prn\" using 1:%d with lines linewidth 2 title \"%s\"\n",mcol+1,legend[mcol-1]);
    fclose(fp);

    /****************************/
    /*gnuplot制御*/
    /****************************/
    fp=popen("gnuplot.exe","w");
        fprintf(fp,"load \"gpl_temp.txt\"\n");
        fflush(fp);
    pclose(fp);
}
/*---------------------------------------------------------------------------*/