write_tab.c File Reference

Write tabulated 1D or 2D ascii data files. More...

#include "pluto.h"

Include dependency graph for write_tab.c:

Go to the source code of this file.

Functions
void	WriteTabArray (Output *output, char *filename, Grid *grid)

Detailed Description

Write tabulated 1D or 2D ascii data files.

WriteTabArray() provides ascii output for 2-D and 1-D data in tabulated multi-column format:

    .     .         .         .             .
    .     .         .         .             .
    .     .         .         .             .
  x1(i) x2(j)   var_1(i,j)  var_2(i,j)   var_3(i,j) ...
    .     .         .         .             .
    .     .         .         .             .
    .     .         .         .             .

Blocks with different x2 coordinates are separated by a blank row (in more than 1D). One file (with all variables) per output is written to disk. This format does not work in parallel mode and can be used for simple serial runs.

Author

A. Mignone (migno.nosp@m.ne@p.nosp@m.h.uni.nosp@m.to.i.nosp@m.t)

Date

June 27, 2014

Definition in file write_tab.c.

Function Documentation

void WriteTabArray	(	Output *	output,
		char *	filename,
		Grid *	grid
	)

Write tabulated array.

Parameters

[in]	output	a pointer to the output structure corresponding to the tab format
[in]	filename	the output file name
[in]	grid	pointer to an array of Grid structures

Definition at line 31 of file write_tab.c.

{
  int  nv, i, j, k;
  FILE *fout;

  #ifdef PARALLEL
   print1 ("! WriteTabArray: tab output not supported in parallel\n");
   return;
  #endif

  #if DIMENSIONS == 3
   print1 ("! WriteTabArray: tab output not supported in 3D\n");
   return;
  #endif

/* ------------------------------------------
      dump arrays in ascii format to disk 
   ------------------------------------------ */

  fout = fopen (filename, "w");
  k = 0;
  IDOM_LOOP (i){
    JDOM_LOOP(j){
      fprintf (fout, "%f %f ", grid[IDIR].x[i], grid[JDIR].x[j]);
      for (nv = 0; nv < output->nvar; nv++) {
        if (output->dump_var[nv]) 
          fprintf (fout, "%12.6e ", output->V[nv][k][j][i]);
      }
      fprintf (fout, "\n");   /* newline */
    }
    #if DIMENSIONS > 1
     fprintf (fout, "\n");   /* skip one more empty line in 2D */
    #endif
  }
  fclose (fout);
}

Here is the call graph for this function:

Here is the caller graph for this function: