if (bRead)
snew(pull->grp,pull->ngrp+1);
for(g=0; g<pull->ngrp+1; g++)
- do_pullgrp(&pull->grp[g],bRead,file_version);
+ do_pullgrp(fio,&pull->grp[g],bRead,file_version);
}
- static void do_rotgrp(t_rotgrp *rotg,bool bRead, int file_version)
++
++static void do_rotgrp(t_fileio *fio, t_rotgrp *rotg,bool bRead, int file_version)
+{
+ bool bDum=TRUE;
+ int i;
+
- do_int(rotg->eType);
- do_int(rotg->bMassW);
- do_int(rotg->nat);
++ gmx_fio_do_int(fio,rotg->eType);
++ gmx_fio_do_int(fio,rotg->bMassW);
++ gmx_fio_do_int(fio,rotg->nat);
+ if (bRead)
+ snew(rotg->ind,rotg->nat);
- ndo_int(rotg->ind,rotg->nat,bDum);
++ gmx_fio_ndo_int(fio,rotg->ind,rotg->nat);
+ if (bRead)
+ snew(rotg->x_ref,rotg->nat);
- ndo_rvec(rotg->x_ref,rotg->nat);
- do_rvec(rotg->vec);
- do_rvec(rotg->pivot);
- do_real(rotg->rate);
- do_real(rotg->k);
- do_real(rotg->slab_dist);
- do_real(rotg->min_gaussian);
- do_real(rotg->eps);
- do_int(rotg->eFittype);
++ gmx_fio_ndo_rvec(fio,rotg->x_ref,rotg->nat);
++ gmx_fio_do_rvec(fio,rotg->vec);
++ gmx_fio_do_rvec(fio,rotg->pivot);
++ gmx_fio_do_real(fio,rotg->rate);
++ gmx_fio_do_real(fio,rotg->k);
++ gmx_fio_do_real(fio,rotg->slab_dist);
++ gmx_fio_do_real(fio,rotg->min_gaussian);
++ gmx_fio_do_real(fio,rotg->eps);
++ gmx_fio_do_int(fio,rotg->eFittype);
+}
+
- static void do_rot(t_rot *rot,bool bRead, int file_version)
++static void do_rot(t_fileio *fio, t_rot *rot,bool bRead, int file_version)
+{
+ int g;
+
- do_int(rot->ngrp);
- do_int(rot->nstrout);
- do_int(rot->nsttout);
++ gmx_fio_do_int(fio,rot->ngrp);
++ gmx_fio_do_int(fio,rot->nstrout);
++ gmx_fio_do_int(fio,rot->nsttout);
+ if (bRead)
+ snew(rot->grp,rot->ngrp);
+ for(g=0; g<rot->ngrp; g++)
- do_rotgrp(&rot->grp[g],bRead,file_version);
++ do_rotgrp(fio, &rot->grp[g],bRead,file_version);
+}
+
- static void do_inputrec(t_inputrec *ir,bool bRead, int file_version,
- real *fudgeQQ)
++
+ static void do_inputrec(t_fileio *fio, t_inputrec *ir,bool bRead,
+ int file_version, real *fudgeQQ)
{
int i,j,k,*tmp,idum=0;
bool bDum=TRUE;
ir->ePull = epullNO;
}
- do_int(ir->bRot);
+ /* Enforced rotation */
+ if (file_version >= 71) {
- do_rot(ir->rot,bRead,file_version);
++ gmx_fio_do_int(fio,ir->bRot);
+ if (ir->bRot == TRUE) {
+ if (bRead)
+ snew(ir->rot,1);
++ do_rot(fio, ir->rot,bRead,file_version);
+ }
+ } else {
+ ir->bRot = FALSE;
+ }
+
/* grpopts stuff */
- do_int(ir->opts.ngtc);
+ gmx_fio_do_int(fio,ir->opts.ngtc);
if (file_version >= 69) {
- do_int(ir->opts.nhchainlength);
+ gmx_fio_do_int(fio,ir->opts.nhchainlength);
} else {
ir->opts.nhchainlength = 1;
}
{ efXVG, "-runav", "runaver", ffOPTWR },
{ efXVG, "-px", "pullx", ffOPTWR },
{ efXVG, "-pf", "pullf", ffOPTWR },
+ { efXVG, "-ro", "rotation", ffOPTWR },
+ { efLOG, "-ra", "rotangles",ffOPTWR },
+ { efLOG, "-rs", "rotslabs", ffOPTWR },
+ { efLOG, "-rt", "rottorque",ffOPTWR },
{ efMTX, "-mtx", "nm", ffOPTWR },
{ efNDX, "-dn", "dipole", ffOPTWR }
- #endif
};
#define NFILE asize(fnm)