Changeset 38d77eb for mulcan_sim.f
- Timestamp:
- 11/19/09 11:29:48 (14 years ago)
- Branches:
- master
- Children:
- 7137e5d
- Parents:
- 32289cd
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
mulcan_sim.f
r32289cd r38d77eb 72 72 end do 73 73 end do 74 write (*,*) 'Number of contacts in reference conformation:',nci74 write (logString, *) 'Number of contacts in reference conformation:',nci 75 75 76 76 ! READ IN FIELDS WITH MULTICANONICAL PARAMETER … … 88 88 vlvr(iv) = x 89 89 end do 90 write (*,*) 'Last iteration, energy:',nswm,eol_old90 write (logString, *) 'Last iteration, energy:',nswm,eol_old 91 91 else 92 92 ! _________________________________ random start … … 100 100 ! 101 101 eol = energy() 102 write ( *,'(e12.5,/)') eol102 write (logString, '(e12.5,/)') eol 103 103 call contacts(nhy,nhx,dham) 104 write (*,*) 'Number of contacts in start configuration:',nhy105 write (*,*) 'Number of native contacts in start configuration:',104 write (logString, *) 'Number of contacts in start configuration:',nhy 105 write (logString, *) 'Number of native contacts in start configuration:', 106 106 & nhx 107 107 do i=1,nresi 108 write (*,'(62I1)') (ijcont(i,j), j=1,nresi)109 end do 110 write (*,*)108 write (logString, '(62I1)') (ijcont(i,j), j=1,nresi) 109 end do 110 write (logString, *) 111 111 ! 112 112 … … 165 165 ! FINAL OUTPUT: 166 166 acz = acz/dble(nsw*nvr) 167 write (*,*) 'last energy',eol168 write (*,*) 'aczeptance rate:',acz167 write (logString, *) 'last energy',eol 168 write (logString, *) 'aczeptance rate:',acz 169 169 170 170 ! WRITE DOWN (UN-REWEIGHTED) HISTOGRAM OF MULTICANONICAL SIMULATION 171 171 do i=kmin,kmax 172 172 if(xhist(i).gt.0.0d0) then 173 write (*,*) i,xhist(i)173 write (logString, *) i,xhist(i) 174 174 end if 175 175 end do
Note:
See TracChangeset
for help on using the changeset viewer.