#endif
}
-void list_tng_for_gmx_dump(const char *fn, gmx_bool bXVG)
+void list_tng_for_gmx_dump(const char *fn)
{
#ifdef GMX_USE_TNG
tng_trajectory_t tng;
gmx_int64_t nframe = 0;
gmx_int64_t i, *block_ids = NULL, step, ndatablocks;
- gmx_int64_t pos_block_id = TNG_TRAJ_POSITIONS;
gmx_bool bOK;
gmx_tng_open(fn, 'r', &tng);
gmx_print_tng_molecule_system(tng, stdout);
bOK = gmx_get_tng_data_block_types_of_next_frame(tng, -1,
- bXVG ? 1 : 0,
- bXVG ? &pos_block_id : NULL,
+ 0,
+ NULL,
&step, &ndatablocks,
&block_ids);
do
/* Can't write any output because we don't know what
arrays are valid. */
fprintf(stderr, "\nWARNING: Incomplete frame at time %g, will not write output\n", frame_time);
- list_tng_inner(fn, (0 == i), bXVG, values, step, frame_time,
+ list_tng_inner(fn, (0 == i), values, step, frame_time,
n_values_per_frame, n_atoms, prec, nframe, block_name);
}
}
nframe++;
}
while (gmx_get_tng_data_block_types_of_next_frame(tng, step,
- bXVG ? 1 : 0,
- bXVG ? &pos_block_id : NULL,
+ 0,
+ NULL,
&step,
&ndatablocks,
&block_ids));
gmx_tng_close(&tng);
#else
GMX_UNUSED_VALUE(fn);
- GMX_UNUSED_VALUE(bXVG);
#endif
}
close_trn(fpread);
}
-void list_xtc(const char *fn, gmx_bool bXVG)
+void list_xtc(const char *fn)
{
t_fileio *xd;
int indent;
nframe = 0;
do
{
- if (bXVG)
- {
- int i, d;
-
- fprintf(stdout, "%g", time);
- for (i = 0; i < natoms; i++)
- {
- for (d = 0; d < DIM; d++)
- {
- fprintf(stdout, " %g", x[i][d]);
- }
- }
- fprintf(stdout, "\n");
- }
- else
- {
- sprintf(buf, "%s frame %d", fn, nframe);
- indent = 0;
- indent = pr_title(stdout, indent, buf);
- pr_indent(stdout, indent);
- fprintf(stdout, "natoms=%10d step=%10d time=%12.7e prec=%10g\n",
- natoms, step, time, prec);
- pr_rvecs(stdout, indent, "box", box, DIM);
- pr_rvecs(stdout, indent, "x", x, natoms);
- }
+ sprintf(buf, "%s frame %d", fn, nframe);
+ indent = 0;
+ indent = pr_title(stdout, indent, buf);
+ pr_indent(stdout, indent);
+ fprintf(stdout, "natoms=%10d step=%10d time=%12.7e prec=%10g\n",
+ natoms, step, time, prec);
+ pr_rvecs(stdout, indent, "box", box, DIM);
+ pr_rvecs(stdout, indent, "x", x, natoms);
nframe++;
}
while (read_next_xtc(xd, natoms, &step, &time, box, x, &prec, &bOK));
/* Callback used by list_tng_for_gmx_dump. */
void list_tng_inner(const char *fn,
gmx_bool bFirstFrame,
- gmx_bool bXVG,
real *values,
gmx_int64_t step,
double frame_time,
char buf[256];
int indent = 0;
- if (bXVG)
+ if (bFirstFrame)
{
- gmx_int64_t j;
- int d;
-
- if (bFirstFrame)
+ sprintf(buf, "%s frame %" GMX_PRId64, fn, nframe);
+ indent = 0;
+ indent = pr_title(stdout, indent, buf);
+ pr_indent(stdout, indent);
+ fprintf(stdout, "natoms=%10" GMX_PRId64 " step=%10" GMX_PRId64 " time=%12.7e",
+ n_atoms, step, frame_time);
+ if (prec > 0)
{
- fprintf(stdout, "%g", (real)frame_time);
- }
- for (j = 0; j < n_atoms; j++)
- {
- for (d = 0; d < DIM; d++)
- {
- fprintf(stdout, " %g", values[j * DIM + d]);
- }
+ fprintf(stdout, " prec=%10g", prec);
}
fprintf(stdout, "\n");
}
- else
- {
- if (bFirstFrame)
- {
- sprintf(buf, "%s frame %" GMX_PRId64, fn, nframe);
- indent = 0;
- indent = pr_title(stdout, indent, buf);
- pr_indent(stdout, indent);
- fprintf(stdout, "natoms=%10" GMX_PRId64 " step=%10" GMX_PRId64 " time=%12.7e",
- n_atoms, step, frame_time);
- if (prec > 0)
- {
- fprintf(stdout, " prec=%10g", prec);
- }
- fprintf(stdout, "\n");
- }
- pr_reals_of_dim(stdout, indent, block_name, values, n_atoms, n_values_per_frame);
- }
+ pr_reals_of_dim(stdout, indent, block_name, values, n_atoms, n_values_per_frame);
}
-void list_trx(const char *fn, gmx_bool bXVG)
+void list_trx(const char *fn)
{
int ftp;
ftp = fn2ftp(fn);
if (ftp == efXTC)
{
- list_xtc(fn, bXVG);
+ list_xtc(fn);
}
else if ((ftp == efTRR) || (ftp == efTRJ))
{
}
else if (ftp == efTNG)
{
- list_tng_for_gmx_dump(fn, bXVG);
+ list_tng_for_gmx_dump(fn);
}
else
{
output_env_t oenv;
/* Command line options */
- static gmx_bool bXVG = FALSE;
static gmx_bool bShowNumbers = TRUE;
static gmx_bool bSysTop = FALSE;
t_pargs pa[] = {
- { "-xvg", FALSE, etBOOL, {&bXVG}, "HIDDENXVG layout for xtc" },
{ "-nr", FALSE, etBOOL, {&bShowNumbers}, "Show index numbers in output (leaving them out makes comparison easier, but creates a useless topology)" },
{ "-sys", FALSE, etBOOL, {&bSysTop}, "List the atoms and bonded interactions for the whole system instead of for each molecule type" }
};
}
else if (ftp2bSet(efTRX, NFILE, fnm))
{
- list_trx(ftp2fn(efTRX, NFILE, fnm), bXVG);
+ list_trx(ftp2fn(efTRX, NFILE, fnm));
}
else if (ftp2bSet(efEDR, NFILE, fnm))
{