Changeset 38d77eb for mulcan_par_mod.f90
- Timestamp:
- 11/19/09 11:29:48 (14 years ago)
- Branches:
- master
- Children:
- 7137e5d
- Parents:
- 32289cd
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
mulcan_par_mod.f90
r32289cd r38d77eb 103 103 104 104 eol = energy() 105 write ( *,'(a,e12.5,/)') 'Energy of start configuration: ',eol106 write (*,*)105 write (logString, '(a,e12.5,/)') 'Energy of start configuration: ',eol 106 write (logString, *) 107 107 108 108 call outpdb(1, 'start.pdb') … … 158 158 ! Final output 159 159 acz = acz/dble(nsw*nvr) 160 write (*,*) 'last energy',eol161 write (*,*) 'acceptance rate:',acz162 163 write (*,*) 'Histogram:'160 write (logString, *) 'last energy',eol 161 write (logString, *) 'acceptance rate:',acz 162 163 write (logString, *) 'Histogram:' 164 164 do i=kmin,kmax 165 165 if(xhist(i).gt.0.0d0) then 166 write (*,*) i,xhist(i)166 write (logString, *) i,xhist(i) 167 167 end if 168 168 end do … … 266 266 end do 267 267 end do 268 write (*,*) 'Number of contacts in reference conformation:',nci268 write (logString, *) 'Number of contacts in reference conformation:',nci 269 269 270 270 ! Read in fields with multicanonical parameter … … 282 282 vlvr(iv) = x 283 283 end do 284 write (*,*) 'Last iteration, energy:',nswm,eol_old284 write (logString, *) 'Last iteration, energy:',nswm,eol_old 285 285 else 286 286 ! _________________________________ random start … … 294 294 295 295 eol = energy() 296 write ( *,'(e12.5,/)') eol296 write (logString, '(e12.5,/)') eol 297 297 call contacts(nhy,nhx,dham) 298 write (*,*) 'Number of contacts in start configuration:',nhy299 write (*,*) 'Number of native contacts in start configuration:', &298 write (logString, *) 'Number of contacts in start configuration:',nhy 299 write (logString, *) 'Number of native contacts in start configuration:', & 300 300 & nhx 301 301 do i=1,nresi 302 write (*,'(62I1)') (ijcont(i,j), j=1,nresi)303 end do 304 write (*,*)302 write (logString, '(62I1)') (ijcont(i,j), j=1,nresi) 303 end do 304 write (logString, *) 305 305 306 306 … … 355 355 356 356 acz = acz/dble(nsw*nvr) 357 write (*,*) 'last energy',eol358 write (*,*) 'acceptance rate:',acz357 write (logString, *) 'last energy',eol 358 write (logString, *) 'acceptance rate:',acz 359 359 360 360 ! WRITE DOWN (UN-REWEIGHTED) HISTOGRAM OF MULTICANONICAL SIMULATION 361 361 do i=kmin,kmax 362 362 if(xhist(i).gt.0.0d0) then 363 write (*,*) i,xhist(i)363 write (logString, *) i,xhist(i) 364 364 end if 365 365 end do
Note:
See TracChangeset
for help on using the changeset viewer.