wham and cluster_wham Adam's new constr_dist multichain
[unres.git] / source / unres / src_MD-M / restbin2asc.F
1       implicit real*8 (a-h,o-z)
2       include 'DIMENSIONS'
3 c      include 'mpif.h'
4       include 'COMMON.MD'
5       include 'COMMON.IOUNITS'
6       include 'COMMON.REMD'
7       include 'COMMON.SETUP'
8       include 'COMMON.CHAIN'
9       include 'COMMON.SBRIDGE'
10       include 'COMMON.INTERACT'
11       integer*2 i_index
12      &            (maxprocs/4,maxprocs/20,maxprocs/200,maxprocs/200)
13       character*128 restin,restout
14       integer idir,n
15       irest2=33
16       ixdrf=34
17       n=iargc()
18       print *,"n",n
19       if (n.lt.5) then
20        print'(2a)',"Usage: restbin 0/1 NTRAJ NRES INRESTFILE",
21      &   " OUTRESTFILE [REMD_M] [USAMPL]"
22        print'(a)',"0 - conversion from bin to ascii / from ascii to bin"
23        print'(a)',"NTRAJ: number of trajectories"
24        print'(a)',"NRES: number of residues"
25        print'(a)',"INRESTFILE: input UNRES/MREMD restart file"
26        print'(a)',"OUTRESTFILE: output UNRES/MREMD restart file"
27        print'(a)',"USAMPL [T/F] flag for umbrella-sampling"
28        stop
29       endif  
30       call getarg(1,liczba)
31       read (liczba,*) idir
32       if (idir. ne. 0 .and. idir.ne.1) stop "idir must be 0 or 1"
33       call getarg(2,liczba)
34       read (liczba,*) nodes
35       call getarg(3,liczba)
36       read (liczba,*) nres
37       call getarg (4,restin)
38       call getarg (5,restout)
39       print *,"idir",idir," nodes",nodes," nres",nres
40       print *,"restin ",restin
41       print *,"restout ",restout
42       if (n.eq.6) then
43         call getarg(6,liczba)
44         read(liczba,*) remd_m(1)
45         do i=2,nodes
46           remd_m(i)=remd_m(1)
47         enddo
48       else
49         do i=2,nodes
50           remd_m(i)=1
51         enddo
52       endif
53       if (n.eq.7) then
54         call getarg(7,liczba)
55         read(liczba,'(L1)') usampl
56       endif
57       print *,"usampl",usampl
58       print *,"remd_m",(remd_m(i),i=1,nodes)
59       if (idir.eq.0) then
60         mremd_rst_name=restin
61         call read1restart(i_index)
62         mremd_rst_name=restout
63         call write1rst_asc(i_index)
64       else
65         mremd_rst_name=restin
66         call read1restart_asc(i_index)
67         mremd_rst_name=restout
68         call write1rst(i_index)
69       endif
70       end
71 c------------------------------------------------------------------------------
72       subroutine read1restart(i_index)
73       implicit real*8 (a-h,o-z)
74       include 'DIMENSIONS'
75 c      include 'mpif.h'
76       include 'COMMON.MD'
77       include 'COMMON.IOUNITS'
78       include 'COMMON.REMD'
79       include 'COMMON.SETUP'
80       include 'COMMON.CHAIN'
81       include 'COMMON.SBRIDGE'
82       include 'COMMON.INTERACT'
83       real d_restart1(3,2*maxres*maxprocs),r_d(3,2*maxres),
84      &     d_restart2(3,2*maxres*maxprocs),t5_restart1(5)
85       integer*2 i_index
86      &            (maxprocs/4,maxprocs/20,maxprocs/200,maxprocs/200)
87       common /przechowalnia/ d_restart1,d_restart2
88       write (*,*) "Processor",me," called read1restart"
89
90 #ifdef AIX
91       call xdrfopen_(ixdrf,mremd_rst_name, "r", iret)
92       print *,"restart file ",mremd_rst_name
93       print *,"iret",iret
94
95       do i=0,nodes-1
96         call xdrfint_(ixdrf, i2rep(i), iret)
97       enddo
98       print *,"iret",iret
99       do i=1,remd_m(1)
100         call xdrfint_(ixdrf, ifirst(i), iret)
101       enddo
102       do il=1,nodes
103         call xdrfint_(ixdrf, nupa(0,il), iret)
104         print *,"nupa(0,",il,")=",nupa(0,il)
105         do i=1,nupa(0,il)
106           call xdrfint_(ixdrf, nupa(i,il), iret)
107         enddo
108         print *,"nupa",(nupa(i,il),i=1,nupa(0,il))
109
110         call xdrfint_(ixdrf, ndowna(0,il), iret)
111         print *,"ndowna(0",il,")=",ndowna(0,il)
112         do i=1,ndowna(0,il)
113           call xdrfint_(ixdrf, ndowna(i,il), iret)
114         enddo
115         print *,"ndowna",(ndowna(i,il),i=1,ndowna(0,il))
116       enddo
117       do il=1,nodes
118         do j=1,4
119           call xdrffloat_(ixdrf, t_restart1(j,il), iret)
120         enddo
121         print *,"il",il,"t_restart",(t_restart1(j,il),j=1,4)
122       enddo
123 #else
124       call xdrfopen(ixdrf,mremd_rst_name, "r", iret)
125
126       do i=0,nodes-1
127         call xdrfint(ixdrf, i2rep(i), iret)
128       enddo
129       do i=1,remd_m(1)
130         call xdrfint(ixdrf, ifirst(i), iret)
131       enddo
132       do il=1,nodes
133         call xdrfint(ixdrf, nupa(0,il), iret)
134         do i=1,nupa(0,il)
135           call xdrfint(ixdrf, nupa(i,il), iret)
136         enddo
137
138         call xdrfint(ixdrf, ndowna(0,il), iret)
139         do i=1,ndowna(0,il)
140           call xdrfint(ixdrf, ndowna(i,il), iret)
141         enddo
142       enddo
143       do il=1,nodes
144         do j=1,4
145           call xdrffloat(ixdrf, t_restart1(j,il), iret)
146         enddo
147       enddo
148 #endif
149       totT=t5_restart1(1)              
150       EK=t5_restart1(2)
151       potE=t5_restart1(3)
152       t_bath=t5_restart1(4)
153
154       do il=0,nodes-1
155         do i=1,2*nres
156           do j=1,3
157 #ifdef AIX
158             call xdrffloat_(ixdrf, d_restart1(j,i+2*nres*il), iret)
159 #else
160             call xdrffloat(ixdrf, d_restart1(j,i+2*nres*il), iret)
161 #endif
162           enddo
163         enddo
164       enddo
165
166       do il=0,nodes-1
167         do i=1,2*nres
168           do j=1,3
169 #ifdef AIX
170             call xdrffloat_(ixdrf, d_restart2(j,i+2*nres*il), iret)
171 #else
172             call xdrffloat(ixdrf, d_restart2(j,i+2*nres*il), iret)
173 #endif
174           enddo
175         enddo
176       enddo
177
178       if (usampl) then
179 #ifdef AIX
180         call xdrfint_(ixdrf, nset, iret)
181         do i=1,nset
182           call xdrfint_(ixdrf,mset(i), iret)
183         enddo
184         do i=0,nodes-1
185           call xdrfint_(ixdrf,i2set(i), iret)
186         enddo
187         do il=1,nset
188           do il1=1,mset(il)
189             do i=1,nrep
190               do j=1,remd_m(i)
191                 call xdrfint_(ixdrf,itmp, iret)
192                 i_index(i,j,il,il1)=itmp
193               enddo
194             enddo
195           enddo
196         enddo
197 #else
198         call xdrfint(ixdrf, nset, iret)
199         do i=1,nset
200           call xdrfint(ixdrf,mset(i), iret)
201         enddo
202         do i=0,nodes-1
203           call xdrfint(ixdrf,i2set(i), iret)
204         enddo
205         do il=1,nset
206           do il1=1,mset(il)
207             do i=1,nrep
208               do j=1,remd_m(i)
209                 call xdrfint(ixdrf,itmp, iret)
210                 i_index(i,j,il,il1)=itmp
211               enddo
212             enddo
213           enddo
214         enddo
215 #endif
216 c              call mpi_scatter(i2set,1,mpi_integer,
217 c     &           iset,1,mpi_integer,king,
218 c     &           CG_COMM,ierr) 
219 c
220       endif
221 #ifdef AIX
222       call xdrfclose_(ixdrf, iret)
223 #else
224       call xdrfclose(ixdrf, iret)
225 #endif
226       return
227       end
228 c------------------------------------------------------------------------------------
229       subroutine write1rst(i_index)
230       implicit real*8 (a-h,o-z)
231       include 'DIMENSIONS'
232 c      include 'mpif.h'
233       include 'COMMON.MD'
234       include 'COMMON.IOUNITS'
235       include 'COMMON.REMD'
236       include 'COMMON.SETUP'
237       include 'COMMON.CHAIN'
238       include 'COMMON.SBRIDGE'
239       include 'COMMON.INTERACT'
240                
241       real d_restart1(3,2*maxres*maxprocs),r_d(3,2*maxres),
242      &     d_restart2(3,2*maxres*maxprocs)
243       real t5_restart1(5)
244       integer iret,itmp
245       integer*2 i_index
246      &            (maxprocs/4,maxprocs/20,maxprocs/200,maxprocs/200)
247        common /przechowalnia/ d_restart1,d_restart2
248
249 #ifdef AIX
250          call xdrfopen_(ixdrf,mremd_rst_name, "w", iret)
251          do i=0,nodes-1
252           call xdrfint_(ixdrf, i2rep(i), iret)
253          enddo
254          do i=1,remd_m(1)
255           call xdrfint_(ixdrf, ifirst(i), iret)
256          enddo
257          do il=1,nodes
258               do i=0,nupa(0,il)
259                call xdrfint_(ixdrf, nupa(i,il), iret)
260               enddo
261
262               do i=0,ndowna(0,il)
263                call xdrfint_(ixdrf, ndowna(i,il), iret)
264               enddo
265          enddo
266
267          do il=1,nodes
268            do j=1,4
269             call xdrffloat_(ixdrf, t_restart1(j,il), iret)
270            enddo
271          enddo
272
273          do il=0,nodes-1
274            do i=1,2*nres
275             do j=1,3
276              call xdrffloat_(ixdrf, d_restart1(j,i+2*nres*il), iret)
277             enddo
278            enddo
279          enddo
280          do il=0,nodes-1
281            do i=1,2*nres
282             do j=1,3
283              call xdrffloat_(ixdrf, d_restart2(j,i+2*nres*il), iret)
284             enddo
285            enddo
286          enddo
287
288          if(usampl) then
289            call xdrfint_(ixdrf, nset, iret)
290            do i=1,nset
291              call xdrfint_(ixdrf,mset(i), iret)
292            enddo
293            do i=0,nodes-1
294              call xdrfint_(ixdrf,i2set(i), iret)
295            enddo
296            do il=1,nset
297              do il1=1,mset(il)
298                do i=1,nrep
299                  do j=1,remd_m(i)
300                    itmp=i_index(i,j,il,il1)
301                    call xdrfint_(ixdrf,itmp, iret)
302                  enddo
303                enddo
304              enddo
305            enddo
306            
307          endif
308          call xdrfclose_(ixdrf, iret)
309 #else
310          call xdrfopen(ixdrf,mremd_rst_name, "w", iret)
311          do i=0,nodes-1
312           call xdrfint(ixdrf, i2rep(i), iret)
313          enddo
314          do i=1,remd_m(1)
315           call xdrfint(ixdrf, ifirst(i), iret)
316          enddo
317          do il=1,nodes
318               do i=0,nupa(0,il)
319                call xdrfint(ixdrf, nupa(i,il), iret)
320               enddo
321
322               do i=0,ndowna(0,il)
323                call xdrfint(ixdrf, ndowna(i,il), iret)
324               enddo
325          enddo
326
327          do il=1,nodes
328            do j=1,4
329             call xdrffloat(ixdrf, t_restart1(j,il), iret)
330            enddo
331          enddo
332
333          do il=0,nodes-1
334            do i=1,2*nres
335             do j=1,3
336              call xdrffloat(ixdrf, d_restart1(j,i+2*nres*il), iret)
337             enddo
338            enddo
339          enddo
340          do il=0,nodes-1
341            do i=1,2*nres
342             do j=1,3
343              call xdrffloat(ixdrf, d_restart2(j,i+2*nres*il), iret)
344             enddo
345            enddo
346          enddo
347
348
349              if(usampl) then
350               call xdrfint(ixdrf, nset, iret)
351               do i=1,nset
352                 call xdrfint(ixdrf,mset(i), iret)
353               enddo
354               do i=0,nodes-1
355                 call xdrfint(ixdrf,i2set(i), iret)
356               enddo
357               do il=1,nset
358                do il1=1,mset(il)
359                 do i=1,nrep
360                  do j=1,remd_m(i)
361                    itmp=i_index(i,j,il,il1)
362                    call xdrfint(ixdrf,itmp, iret)
363                  enddo
364                 enddo
365                enddo
366               enddo
367            
368              endif
369          call xdrfclose(ixdrf, iret)
370 #endif
371       return
372       end
373 c----------------------------------------------------------------------------
374       subroutine read1restart_asc(i_index)
375       implicit real*8 (a-h,o-z)
376       include 'DIMENSIONS'
377       include 'COMMON.MD'
378       include 'COMMON.IOUNITS'
379       include 'COMMON.REMD'
380       include 'COMMON.SETUP'
381       include 'COMMON.CHAIN'
382       include 'COMMON.SBRIDGE'
383       include 'COMMON.INTERACT'
384                
385       real d_restart1(3,2*maxres*maxprocs),r_d(3,2*maxres),
386      &     d_restart2(3,2*maxres*maxprocs)
387       real t5_restart1(5)
388       integer iret,itmp
389       integer*2 i_index
390      &            (maxprocs/4,maxprocs/20,maxprocs/200,maxprocs/200)
391       common /przechowalnia/ d_restart1,d_restart2
392
393       open(irest2,file=mremd_rst_name,status="old")
394       read (irest2,*) (i2rep(i),i=0,nodes-1)
395       read (irest2,*) (ifirst(i),i=1,remd_m(1))
396       do il=1,nodes
397         read (irest2,*) nupa(0,il),(nupa(i,il),i=1,nupa(0,il))
398         read (irest2,*) ndowna(0,il),(ndowna(i,il),i=1,ndowna(0,il))
399       enddo
400       do il=1,nodes
401         read (irest2,*) (t_restart1(j,il),j=1,4)
402       enddo
403       do il=0,nodes-1
404         do i=1,2*nres
405           read (irest2,*) (d_restart1(j,i+2*nres*il),j=1,3)
406         enddo
407       enddo
408       do il=0,nodes-1
409         do i=1,2*nres
410           read (irest2,*) (d_restart2(j,i+2*nres*il),j=1,3)
411         enddo
412       enddo
413
414       if (usampl) then
415         read (irest2,*) nset
416         read (irest2,*) (mset(i),i=1,nset)
417         read (irest2,*) (i2set(i),i=1,nodes-1)
418         do il=1,nset
419           do il1=1,mset(il)
420             do i=1,nrep
421               read (irest2,*) (i_index(i,j,il,il1),j=1,remd_m(i))
422             enddo
423           enddo
424         enddo
425       endif
426       return
427       end
428 c---------------------------------------------------------------------------
429       subroutine write1rst_asc(i_index)
430       implicit real*8 (a-h,o-z)
431       include 'DIMENSIONS'
432       include 'COMMON.MD'
433       include 'COMMON.IOUNITS'
434       include 'COMMON.REMD'
435       include 'COMMON.SETUP'
436       include 'COMMON.CHAIN'
437       include 'COMMON.SBRIDGE'
438       include 'COMMON.INTERACT'
439                
440       real d_restart1(3,2*maxres*maxprocs),r_d(3,2*maxres),
441      &     d_restart2(3,2*maxres*maxprocs)
442       real t5_restart1(5)
443       integer iret,itmp
444       integer*2 i_index
445      &            (maxprocs/4,maxprocs/20,maxprocs/200,maxprocs/200)
446       common /przechowalnia/ d_restart1,d_restart2
447
448       open(irest2,file=mremd_rst_name,status="unknown")
449       write (irest2,*) (i2rep(i),i=0,nodes-1)
450       write (irest2,*) (ifirst(i),i=1,remd_m(1))
451       do il=1,nodes
452         write (irest2,*) nupa(0,il),(nupa(i,il),i=1,nupa(0,il))
453         write (irest2,*) ndowna(0,il),(ndowna(i,il),i=1,ndowna(0,il))
454       enddo
455       do il=1,nodes
456         write (irest2,*) (t_restart1(j,il),j=1,4)
457       enddo
458       do il=0,nodes-1
459         do i=1,2*nres
460           write (irest2,*) (d_restart1(j,i+2*nres*il),j=1,3)
461         enddo
462       enddo
463       do il=0,nodes-1
464         do i=1,2*nres
465           write (irest2,*) (d_restart2(j,i+2*nres*il),j=1,3)
466         enddo
467       enddo
468
469       if (usampl) then
470         write (irest2,*) nset
471         write (irest2,*) (mset(i),i=1,nset)
472         write (irest2,*) (i2set(i),i=1,nodes-1)
473         do il=1,nset
474           do il1=1,mset(il)
475             do i=1,nrep
476               write (irest2,*) (i_index(i,j,il,il1),j=1,remd_m(i))
477             enddo
478           enddo
479         enddo
480       endif
481       return
482       end