660d4c8afba2fa28fe4826f6382b9108ec86a021
[unres.git] / source / unres / src_MD-M / xdrf2pdb / xdrf2x.F
1       real coord(3,1000)
2       real prec,time,potE,uconst,t_bath,qfrag(100)
3       integer ihpb(100),jhpb(100),iset
4       character*80 arg
5       
6       print *,"start"
7       call getarg(1,arg)
8
9       print *,arg
10 #ifdef BLUEGENE
11       call xdrfopen_(ixdrf,arg, "r", iret)
12        
13       do while(.true.) 
14        call xdrffloat_(ixdrf, time, iret)
15        if(iret.eq.0) goto 10
16        call xdrffloat_(ixdrf, potE, iret)
17        call xdrffloat_(ixdrf, uconst, iret)
18        call xdrffloat_(ixdrf, t_bath, iret)
19        call xdrfint_(ixdrf, nss, iret) 
20        do j=1,nss
21         call xdrfint_(ixdrf, ihpb(j), iret)
22         call xdrfint_(ixdrf, jhpb(j), iret)
23        enddo
24        call xdrfint_(ixdrf, nfrag, iret)
25        call xdrfint_(ixdrf, iset, iret)
26        do i=1,nfrag
27         call xdrffloat_(ixdrf, qfrag(i), iret)
28        enddo
29        prec=10000.0
30
31        isize=0
32        call xdrf3dfcoord_(ixdrf, coord, isize, prec, iret)
33 #else
34       call xdrfopen(ixdrf,arg, "r", iret)
35        
36       do while(.true.) 
37        call xdrffloat(ixdrf, time, iret)
38        if(iret.eq.0) goto 10
39        call xdrffloat(ixdrf, potE, iret)
40        call xdrffloat(ixdrf, uconst, iret)
41        call xdrffloat(ixdrf, t_bath, iret)
42        call xdrfint(ixdrf, nss, iret) 
43        do j=1,nss
44         call xdrfint(ixdrf, ihpb(j), iret)
45         call xdrfint(ixdrf, jhpb(j), iret)
46        enddo
47        call xdrfint(ixdrf, nfrag, iret)
48        call xdrfint(ixdrf, iset, iret)
49        do i=1,nfrag
50         call xdrffloat(ixdrf, qfrag(i), iret)
51        enddo
52        prec=10000.0
53
54        isize=0
55        call xdrf3dfcoord(ixdrf, coord, isize, prec, iret)
56 #endif
57
58        write (*,'(e15.8,2e15.5,f12.5,$)') time,potE,uconst,t_bath
59        write (*,'(i4,$)') nss,(ihpb(j),jhpb(j),j=1,nss)
60        write (*,'(i4)') iset
61        write (*,'(i4,20f7.4)') nfrag,(qfrag(i),i=1,nfrag)
62        write (*,'(8f10.5)') ((coord(k,j),k=1,3),j=1,isize)
63       enddo
64   10  continue
65       end