call xdrffloat(ixdrf, potE, iret)
call xdrffloat(ixdrf, uconst, iret)
call xdrffloat(ixdrf, t_bath, iret)
+ print *,"potE",potE," uconst",uconst," t_bath",t_bath
#ifdef NEWUNRES
call xdrffloat(ixdrf, uconst_back, iret)
#endif
+ print *,"uconst_back",uconst_back
call xdrfint(ixdrf, nss, iret)
do j=1,nss
call xdrfint(ixdrf, ihpb(j), iret)
prec=10000.0
isize=0
+ print *," call xdrf3coord"
call xdrf3dfcoord(ixdrf, coord, isize, prec, iret)