#include <stdio.h>
#include <stddef.h>
#include <string.h>
#include <stdarg.h>
#include <stdlib.h>
#include "genlib.h"
#include "simpio.h"
#include "geostat.h"
#include "toolsIO.h"

void inputfiledata(string inputfile,long *seed,struct grid_mod *grid,string filename[7],struct vario_mod *variogram,struct statistic_mod *stat,struct pressure_mod *pression)
{
  FILE *fp;
  int i,j;

  fp=fopen(inputfile,"r");
  if(fp== NULL)
    {
      printf("Erreur d'ouverture du fichier\n");
      exit(0);
    }
  *seed=atoi(ReadLine(fp));
  (*grid).NX = atoi(ReadLine(fp));
  (*grid).NY = atoi(ReadLine(fp));
  (*grid).NZ = atoi(ReadLine(fp));
  (*grid).DX = atof(ReadLine(fp));
  (*grid).DY = atof(ReadLine(fp));
  (*grid).DZ = atof(ReadLine(fp));
  filename[0] = ReadLine(fp);
  (*variogram).Nvario = atoi(ReadLine(fp));

  (*variogram).vario = (int *) malloc((*variogram).Nvario * sizeof(int));
  (*variogram).alpha = (double *) malloc((*variogram).Nvario * sizeof(double));      
  (*variogram).ap = (double *) malloc(9*(*variogram).Nvario * sizeof(double));
  (*variogram).scf = (double *) malloc(3*(*variogram).Nvario * sizeof(double));
  (*variogram).var = (double *) malloc((*variogram).Nvario * sizeof(double));
  for (i= 0; i < (*variogram).Nvario; i++) {
    (*variogram).var[i] = atof(ReadLine(fp));
    (*variogram).vario[i] = atoi(ReadLine(fp));
    (*variogram).alpha[i] = atof(ReadLine(fp));
    for (j = 0; j < 3; j++) 
      (*variogram).scf[i*3+j]= atof(ReadLine(fp));
    for (j = 0; j < 6; j++) 
      (*variogram).ap[i*9+j] = atof(ReadLine(fp));
  }
  /*statistical data*/
  (*stat).nblock_mean = 1;
  (*stat).nblock_var = 1;
  (*stat).mean = (double *)malloc((*stat).nblock_mean * sizeof(double));
  if ((*stat).mean == NULL)
    Error("No memory available");
  (*stat).variance = (double *)malloc((*stat).nblock_var * sizeof(double));
  if ((*stat).variance == NULL)
    Error("No memory available");

  (*stat).mean[0] = atof(ReadLine(fp));
  (*stat).variance[0] = atof(ReadLine(fp));
  (*stat).type = atoi(ReadLine(fp));

  filename[1] = ReadLine(fp);
  (*pression).x = atof(ReadLine(fp));
  (*pression).y = atof(ReadLine(fp));
  (*pression).z = atof(ReadLine(fp));
  filename[2] = ReadLine(fp);
  filename[3] = ReadLine(fp);
  filename[4] = ReadLine(fp);
  filename[5] = ReadLine(fp);
  filename[6] = ReadLine(fp);
  fclose(fp);
  return;
}