if(iret.eq.0) exit
kk = kk + 1
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
if(iret.eq.0) exit
kk = kk + 1
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