do while (iret.gt.0)
#if (defined(AIX) && !defined(JUBL))
+#ifdef DEBUG
+ write (iout,*) "ii",ii," itraj",itraj," it",it
+#endif
call xdrffloat_(ixdrf, rtime, iret)
-c print *,"rtime",rtime," iret",iret
call xdrffloat_(ixdrf, rpotE, iret)
-c write (iout,*) "rpotE",rpotE," iret",iret
+#ifdef DEBUG
+ write (iout,*) "rtime",rtime," rpotE",rpotE," iret",iret
+#endif
call flush(iout)
call xdrffloat_(ixdrf, ruconst, iret)
call xdrffloat_(ixdrf, rt_bath, iret)
call xdrfint_(ixdrf, nss, iret)
+#ifdef DEBUG
+ write (iout,*) "ruconst",ruconst," rt_bath",rt_bath," nss",nss
+#endif
do j=1,nss
call xdrfint_(ixdrf, ihpb(j), iret)
call xdrfint_(ixdrf, jhpb(j), iret)
call xdrffloat_(ixdrf, rprop(i), iret)
enddo
#else
+#ifdef DEBUG
+ write (iout,*) "ii",ii," itraj",itraj," it",it
+#endif
call xdrffloat(ixdrf, rtime, iret)
call xdrffloat(ixdrf, rpotE, iret)
-c write (iout,*) "rpotE",rpotE," iret",iret
+#ifdef DEBUG
+ write (iout,*) "rtime",rtime," rpotE",rpotE," iret",iret
+#endif
call flush(iout)
call xdrffloat(ixdrf, ruconst, iret)
call xdrffloat(ixdrf, rt_bath, iret)
call xdrfint(ixdrf, nss, iret)
+#ifdef DEBUG
+ write (iout,*) "ruconst",ruconst," rt_bath",rt_bath," nss",nss
+#endif
do j=1,nss
call xdrfint(ixdrf, ihpb(j), iret)
call xdrfint(ixdrf, jhpb(j), iret)
#endif
if (iret.eq.0) exit
itraj=mod(it,totraj(iR,iparm))
-#ifdef DEBUG
- write (iout,*) "ii",ii," itraj",itraj," it",it
-#endif
if (iset.eq.0) iset = 1
call flush(iout)
it=it+1