integral += x1;
}
integral *= dt;
- ffclose(fp);
+ gmx_ffclose(fp);
printf("%s lifetime = %.2f ps\n", bContact ? "Contact" : "HB", integral);
printf("Note that the lifetime obtained in this manner is close to useless\n");
printf("Use the -ac option instead and check the Forward lifetime\n");
{
return;
}
- fp = ffopen("debug-ac.xvg", "w");
+ fp = gmx_ffopen("debug-ac.xvg", "w");
for (j = 0; (j < nframes); j++)
{
fprintf(fp, "%10.3f", hb->time[j]);
}
fprintf(fp, "\n");
}
- ffclose(fp);
+ gmx_ffclose(fp);
}
static real calc_dg(real tau, real temp)
fprintf(fp, "%10g %10g %10g %10g %10g\n",
hb->time[j]-hb->time[0], ct[j], cct[j], ght[j], kt[j]);
}
- ffclose(fp);
+ gmx_ffclose(fp);
analyse_corr(nn, hb->time, ct, ght, kt, NULL, NULL, NULL,
fit_start, temp, smooth_tail_start, oenv);
fp = opt2FILE("-hbn", nfile, fnm, "w");
if (opt2bSet("-g", nfile, fnm))
{
- fplog = ffopen(opt2fn("-g", nfile, fnm), "w");
+ fplog = gmx_ffopen(opt2fn("-g", nfile, fnm), "w");
fprintf(fplog, "# %10s %12s %12s\n", "Donor", "Hydrogen", "Acceptor");
}
else
}
}
}
- ffclose(fp);
+ gmx_ffclose(fp);
if (fplog)
{
- ffclose(fplog);
+ gmx_ffclose(fplog);
}
}
close_trj(status);
if (fpnhb)
{
- ffclose(fpnhb);
+ gmx_ffclose(fpnhb);
}
/* Compute maximum possible number of different hbonds */
aver_nhb += hb->nhb[i];
aver_dist += hb->ndist[i];
}
- ffclose(fp);
+ gmx_ffclose(fp);
aver_nhb /= nframes;
aver_dist /= nframes;
/* Print HB distance distribution */
{
fprintf(fp, "%10g %10g\n", (i+0.5)*rbin, rdist[i]/(rbin*(real)sum));
}
- ffclose(fp);
+ gmx_ffclose(fp);
}
/* Print HB angle distribution */
{
fprintf(fp, "%10g %10g\n", (i+0.5)*abin, adist[i]/(abin*(real)sum));
}
- ffclose(fp);
+ gmx_ffclose(fp);
}
/* Print HB in alpha-helix */
}
fprintf(fp, "\n");
}
- ffclose(fp);
+ gmx_ffclose(fp);
}
if (!bNN)
{
}
fp = opt2FILE("-hbm", NFILE, fnm, "w");
write_xpm_m(fp, mat);
- ffclose(fp);
+ gmx_ffclose(fp);
for (x = 0; x < mat.nx; x++)
{
sfree(mat.matrix[x]);
}
fprintf(fp, "\n");
}
- ffclose(fp);
+ gmx_ffclose(fp);
}
return 0;