#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 "pressure.h"
#include "toolsIO.h"

void  pgeneration2(int n[3],struct grid_mod grid,struct statistic_mod stat,struct vario_mod variogram,string filename[7],struct pressure_mod pression,struct realization_mod *Y,struct realization_mod *P,struct realization_mod *VX,struct realization_mod *VY,struct realization_mod *VZ, int solver, int format_file)
   {
     int i,ntot;
     struct realization_mod GP;

     //FFTPressure(n,grid,Y,stat,pression,P,VX,VY,VZ,solver);

     /*save the delta-pressure realization*/
     switch (format_file)
       {
       case 0:
	 writefile(filename[2],P);
	 break;
       case 1:
	 writefile_bin(filename[2],P);
	 break;
       }
     
     //total_pressure(grid,pression,P);
     total_velocity(grid,stat.mean[0],1,pression,VX);
     total_velocity(grid,stat.mean[0],2,pression,VY);
     total_velocity(grid,stat.mean[0],3,pression,VZ);

     /*save the total pressure realization*/
     switch (format_file)
       {
       case 0:
	 writefile(filename[3],P);
	 break;
       case 1:
	 writefile_bin(filename[3],P);
	 break;
       }

     /*save the x-velocity realization*/
     switch (format_file)
       {
       case 0:
	 writefile(filename[4],VX);
	 break;
       case 1:
	 writefile_bin(filename[4],VX);
	 break;
       }

     /*save the y-velocity realization*/
     switch (format_file)
       {
       case 0:
	 writefile(filename[5],VY);
	 break;
       case 1:
	 writefile_bin(filename[5],VY);
	 break;
       }
     /*save the z-velocity realization*/
     switch (format_file)
       {
       case 0:
	 writefile(filename[6],VZ);
	 break;
       case 1:
	 writefile_bin(filename[6],VZ);
	 break;
       }

     ntot=grid.NX*grid.NY*grid.NZ;
     GP.vector = (double *) malloc(ntot * sizeof(double));
     testmemory(GP.vector);
     GP.n= ntot-1;
     for (i = 0; i < ntot-1; i++) 
       {
	 GP.vector[i]=(*P).vector[i]-(*P).vector[i-1];
       }

     /*save the pressure gradient realization*/
      switch (format_file)
       {
       case 0:
	 writefile(filename[7],&GP);
	 break;
       case 1:
	 writefile_bin(filename[7],&GP);
	 break;
       }

      return;
   }