/*===========================*/
/* gccONDO4-splot.c            */
/*===========================*/
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <time.h>
#define NDmax 2000
#define NEmax 2000
#define MAmax 10
#define NTOUTmax 10

int main_part(char fnameR[50],char fnameF[50]);
void GPL_mesh(char fnameF[50],int NODT,int NELT,int kakom[NEmax+1][5],double x[NDmax+1],double y[NDmax+1]);
void GPL_splot(char fnameF[50],int NODT,int NELT,int kakom[NEmax+1][5],double x[NDmax+1],double y[NDmax+1],double tempand[NDmax+1][NTOUTmax],int ntout);
/*------------------------------------------------------------------------------*/

int main(int argc, char* argv[])
{
    char fnameR[50];
    char fnameF[50];

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

    main_part(fnameR,fnameF);

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

int main_part(char fnameR[],char fnameF[])
{
    /*２次元四角形要素熱伝導解析 */
    FILE *fin,*fout;
    char dat[1000]="";
    char sv[20]="";
    int i,j,k,l,kk,ne,n,ia,ja,mn;
    char strcom[256]="";  /* 書出用コメント */
    int NODT;             /*節点総数*/
    int NELT;             /*要素総数*/
    int MATEL;            /*材料種類数*/
    int KOT;              /*温度指定節点数*/
    int KOC;              /*熱伝達指定辺数*/
    /*--------------------------------------------------------------------------------------*/
    /*配列寸法宣言*/
    /*--------------------------------------------------------------------------------------*/
    int    kakom[NEmax+1][5];         /* 要素構成節点番号 */
    int    matno[NEmax+1];            /* 材料種別No */
    double x[NDmax+1];                /* 節点x座標 */
    double y[NDmax+1];                /* 節点y座標 */
    double wk;
    int itmax,iw,iw1,iw2,ntout;
    double tempand[NDmax+1][NTOUTmax];

    /******************************************************/
    /*データ入力*/
    /******************************************************/
    fin=fopen(fnameR,"r");
    fgets(dat,sizeof dat,fin);strncpy(strcom,dat,strlen(dat)-1);
    fgets(dat,sizeof dat,fin);
    fgets(dat,sizeof dat,fin);
    NODT  =atoi(strtok(dat,","));
    NELT  =atoi(strtok(NULL,","));
    MATEL =atoi(strtok(NULL,","));
    KOT   =atoi(strtok(NULL,","));
    KOC   =atoi(strtok(NULL,","));
    fgets(dat,sizeof dat,fin);
    fgets(dat,sizeof dat,fin);
    for(i=1;i<=NODT;i++){
        fgets(dat,sizeof dat,fin);
        iw  =atoi(strtok(dat,","));
        x[i]=atof(strtok(NULL,","));
        y[i]=atof(strtok(NULL,","));
    }
    fgets(dat,sizeof dat,fin);
    for(ne=1;ne<=NELT;ne++){
        fgets(dat,sizeof dat,fin);
        iw=atoi(strtok(dat,","));
        kakom[ne][1]=atoi(strtok(NULL,","));
        kakom[ne][2]=atoi(strtok(NULL,","));
        kakom[ne][3]=atoi(strtok(NULL,","));
        kakom[ne][4]=atoi(strtok(NULL,","));
        wk=atof(strtok(NULL,","));
        wk=atof(strtok(NULL,","));
        wk=atof(strtok(NULL,","));
        wk=atof(strtok(NULL,","));
        wk=atof(strtok(NULL,","));
        wk=atof(strtok(NULL,","));
        wk=atof(strtok(NULL,","));
        matno[ne]=atoi(strtok(NULL,","));
    }
    fgets(dat,sizeof dat,fin);
    fgets(dat,sizeof dat,fin);
        strtok(dat,",");
        itmax=atoi(strtok(NULL,","));
        printf("itmax=%d\n",itmax);
    fgets(dat,sizeof dat,fin);
    fgets(dat,sizeof dat,fin);
    for(i=0;i<=itmax;i++){
        fgets(dat,sizeof dat,fin);
    }
    fgets(dat,sizeof dat,fin);
    fgets(dat,sizeof dat,fin);
        strtok(dat,",");
        ntout=atoi(strtok(NULL,","));
        printf("ntout=%d\n",ntout);
    fgets(dat,sizeof dat,fin);
    for(i=1;i<=NODT;i++){
        fgets(dat,sizeof dat,fin);
            iw=atoi(strtok(dat,","));
            for(j=0;j<=ntout;j++){
                tempand[i][j]=atof(strtok(NULL,","));
            }
    }
    fclose(fin);
/*
    GPL_mesh(fnameF,NODT,NELT,kakom,x,y);
*/
    GPL_splot(fnameF,NODT,NELT,kakom,x,y,tempand,ntout);
}
/*------------------------------------------------------------------------------*/
void GPL_mesh(char fnameF[50],int NODT,int NELT,int kakom[NEmax+1][5],double x[NDmax+1],double y[NDmax+1])
{
    FILE *fp;
    char fnameW[50]="gpl_temp_ondo.txt";
    double x0,y0, x1,y1,x2,y2,x3,y3,ang;
    char str[50],sv[50],cn[10];
    int i,j,k,l,ne;
    double xmin,xmax,ymin,ymax;

    strcpy(cn,"#FFFFE0");
    xmin=x[1];
    xmax=x[1];
    ymin=y[1];
    ymax=y[1];
    for(i=1;i<=NODT;i++){
        if(x[i]<xmin)xmin=x[i];
        if(xmax<x[i])xmax=x[i];
        if(y[i]<ymin)ymin=y[i];
        if(ymax<y[i])ymax=y[i];
    }

    fp=fopen(fnameW,"w");
    fprintf(fp,"reset\n");
    fprintf(fp,"set terminal postscript eps enhanced color font \"Arial\" 12\n");
    fprintf(fp,"set output \"%s\"\n",fnameF);
    fprintf(fp,"#\n");
    fprintf(fp,"set size ratio -1\n");
    fprintf(fp,"#\n");
    fprintf(fp,"set xrange [%.3f:%.3f]\n",xmin,xmax);
    fprintf(fp,"set yrange [%.3f:%.3f]\n",ymin,ymax);
    fprintf(fp,"set border 0\n");
    fprintf(fp,"set xtics 0,0,0\n");
    fprintf(fp,"set ytics 0,0,0\n");
    fprintf(fp,"set xtics nomirror\n");
    fprintf(fp,"set ytics nomirror\n");
    fprintf(fp,"#\n");

    fprintf(fp,"#\n");

    for(ne=1;ne<=NELT;ne++){
        i=kakom[ne][1];
        j=kakom[ne][2];
        k=kakom[ne][3];
        l=kakom[ne][4];
        fprintf(fp,"set object %d polygon from %.3f,%.3f to %.3f,%.3f to %.3f,%.3f to %.3f,%.3f to %.3f,%.3f\n",
            ne,x[i],y[i],x[j],y[j],x[k],y[k],x[l],y[l],x[i],y[i]);
        fprintf(fp,"set object %d fillcolor rgb \"%s\" fillstyle solid 1.0 border lt -1\n",ne,cn);
    }
    fprintf(fp,"#\n");
    fprintf(fp,"plot \"-\" with point notitle\n");
    fprintf(fp,"%.3f %.3f\n",x[1],y[1]);
    fprintf(fp,"e\n");

    fprintf(fp,"#\n");
    fprintf(fp,"quit\n");
    fclose(fp);

    fp=popen("gnuplot.exe","w");
        fprintf(fp,"load \"%s\"\n",fnameW);
        fflush(fp);
    pclose(fp);

}
/*---------------------------------------------------------------------------*/
void GPL_splot(char fnameF[50],int NODT,int NELT,int kakom[NEmax+1][5],double x[NDmax+1],double y[NDmax+1],double tempand[NDmax+1][NTOUTmax],int ntout)
{
    FILE *fp;
    char fnameW[50]="gpl_temp_ondo.txt";
    char fnameT[50]="data_temp_ondo.prn";
    int i,j,k,l,ne,jj;
    double xmin,xmax,ymin,ymax;
    double tmax,tmin,tt;

    fp=fopen(fnameT,"w");
    jj=7;
    for(ne=1;ne<=NELT;ne++){
        i=kakom[ne][1];
        j=kakom[ne][2];
        k=kakom[ne][3];
        l=kakom[ne][4];
        fprintf(fp,"%.3f %.3f %.3f\n",x[i],y[i],tempand[i][jj]);
        fprintf(fp,"%.3f %.3f %.3f\n",x[j],y[j],tempand[j][jj]);
        fprintf(fp,"%.3f %.3f %.3f\n",x[k],y[k],tempand[k][jj]);
        fprintf(fp,"%.3f %.3f %.3f\n",x[l],y[l],tempand[l][jj]);
        fprintf(fp,"%.3f %.3f %.3f\n",x[i],y[i],tempand[i][jj]);
        fprintf(fp,"\n");
        fprintf(fp,"\n");
    }
    fclose(fp);

    tmax=tempand[1][0];
    tmin=tempand[1][0];
    for(i=1;i<=NODT;i++){
        for(j=0;j<=ntout;j++){
            if(tempand[i][j]<tmin)tmin=tempand[i][j];
            if(tmax<tempand[i][j])tmax=tempand[i][j];
    }}

    xmin=x[1];
    xmax=x[1];
    ymin=y[1];
    ymax=y[1];
    for(i=1;i<=NODT;i++){
        if(x[i]<xmin)xmin=x[i];
        if(xmax<x[i])xmax=x[i];
        if(y[i]<ymin)ymin=y[i];
        if(ymax<y[i])ymax=y[i];
    }

    printf("xmin=%.3f xmax=%.3f\n",xmin,xmax);
    printf("ymin=%.3f ymax=%.3f\n",ymin,ymax);
    printf("tmin=%.3f tmax=%.3f\n",tmin,tmax);

    fp=fopen(fnameW,"w");
    fprintf(fp,"reset\n");
    fprintf(fp,"set terminal postscript eps enhanced color font \"Arial\" 12\n");
    fprintf(fp,"set output \"%s\"\n",fnameF);
    fprintf(fp,"#\n");
    fprintf(fp,"set size ratio -1\n");
    fprintf(fp,"#\n");
    fprintf(fp,"set xrange [%.3f:%.3f]\n",xmin,xmax);
    fprintf(fp,"set yrange [%.3f:%.3f]\n",ymin,ymax);
    fprintf(fp,"set zrange [%.3f:%.3f]\n",tmin,tmax);
    fprintf(fp,"set ticslevel 0\n");
    fprintf(fp,"#\n");
    fprintf(fp,"splot \"%s\" with lines notitle\n",fnameT);
    fprintf(fp,"#\n");
    fprintf(fp,"quit\n");
    fclose(fp);

    fp=popen("gnuplot.exe","w");
        fprintf(fp,"load \"%s\"\n",fnameW);
        fflush(fp);
    pclose(fp);

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