- Timestamp:
- 11/19/09 11:29:48 (14 years ago)
- Branches:
- master
- Children:
- 7137e5d
- Parents:
- 32289cd
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
minim.f
r32289cd r38d77eb 58 58 59 59 if (msv.gt.msvmx) then 60 write ( *,'(a,i5)') ' minreg> parameter MSV > ',msvmx60 write (logString, '(a,i5)') ' minreg> parameter MSV > ',msvmx 61 61 stop 62 62 endif … … 66 66 call gradient() 67 67 68 write (*,'(/,a,/)') ' Energy BEFORE minimization:'68 write (logString, '(/,a,/)') ' Energy BEFORE minimization:' 69 69 70 70 if (ireg.eq.0) then 71 71 72 write ( *,'(a,e12.5,/,3(a,e11.4),/,2(a,e11.4),/)') ' Total: ',73 & eysm,72 write (logString, '(a,e12.5,/,3(a,e11.4),/,2(a,e11.4),/)') 73 & ' Total: ', eysm, 74 74 & ' Coulomb: ',eyel,' Lennard-Jones: ',eyvw,' HB: ',eyhb, 75 75 & ' Variables: ',eyvr,' Solvatation: ',eysl … … 78 78 79 79 80 write ( *,'(a,e12.5,/,3(a,e11.4),/,3(a,e11.4),/)') ' Total: ',81 & wtey*eysm + wtrg*eyrg,80 write (logString, '(a,e12.5,/,3(a,e11.4),/,3(a,e11.4),/)') 81 & ' Total: ', wtey*eysm + wtrg*eyrg, 82 82 & ' Coulomb: ',eyel,' Lennard-Jones: ',eyvw,' HB: ',eyhb, 83 83 & ' Variables: ',eyvr,' Solvatation: ',eysl, … … 148 148 ! gdrg2 = max(acc,gdrg2) 149 149 ! wtrg = wtrg * sqrt(gdey2/gdrg2) 150 ! write (*,*) ' --> Wt_energy = ',wtey,' Wt_regul. = ',wtrg151 ! write (*,*) ' '150 ! write (logString, *) ' --> Wt_energy = ',wtey,' Wt_regul. = ',wtrg 151 ! write (logString, *) ' ' 152 152 ! endif 153 153 … … 184 184 185 185 if (nop.lt.mxop) then 186 write ( *,'(a)') ' ---- CONVERGENCE ----'186 write (logString, '(a)') ' ---- CONVERGENCE ----' 187 187 else 188 write ( *,'(a)') '---- STEP LIMIT ----'189 endif 190 191 192 write ( *,'(/,2a,/)') ' Final energies ',188 write (logString, '(a)') '---- STEP LIMIT ----' 189 endif 190 191 192 write (logString, '(/,2a,/)') ' Final energies ', 193 193 & '__________________________________________________' 194 194 … … 197 197 if (ireg.eq.0) then 198 198 199 write ( *,'(a,e12.5,/,3(a,e11.4),/,2(a,e11.4))') ' Total: ',eysm,200 & ' Coulomb: ',eyel,' Lennard-Jones: ',eyvw,' HB: ',eyhb,201 & ' Variables: ',eyvr,' Solvatation: ',eysl199 write (logString, '(a,e12.5,/,3(a,e11.4),/,2(a,e11.4))') 200 & ' Total: ',eysm,' Coulomb: ',eyel,' Lennard-Jones: ',eyvw, 201 & ' HB: ',eyhb, ' Variables: ',eyvr,' Solvatation: ',eysl 202 202 203 203 else 204 204 205 write ( *,'(a,e12.5,/,3(a,e11.4),/,3(a,e11.4))') ' Total: ',206 & 205 write (logString, '(a,e12.5,/,3(a,e11.4),/,3(a,e11.4))') 206 & ' Total: ', wtey*eysm + wtrg*eyrg, 207 207 & ' Coulomb: ',eyel,' Lennard-Jones: ',eyvw,' HB: ',eyhb, 208 208 & ' Variables: ',eyvr,' Solvatation: ',eysl, … … 211 211 endif 212 212 213 write ( *,'(/,a,/)') ' Variables _________________'213 write (logString, '(/,a,/)') ' Variables _________________' 214 214 215 215 nv = 0 … … 221 221 if (abs(vr).gt.pi) vr=vr-sign(pi2,vr) 222 222 223 write ( *,'(1x,a,1x,i4,f8.1,a,f5.1,a)') nmvr(i),nursvr(i),224 & vr*crd,' (',abs(difang(vr,vlvro(i)))*crd,')'223 write (logString, '(1x,a,1x,i4,f8.1,a,f5.1,a)') nmvr(i), 224 & nursvr(i), vr*crd,' (',abs(difang(vr,vlvro(i)))*crd,')' 225 225 226 226 vlvr(i) = vr … … 231 231 if (ireg.ne.0) then 232 232 233 write ( *,'(/,a,/)') ' Global Variables ___________'233 write (logString,'(/,a,/)') ' Global Variables ___________' 234 234 235 235 do i=1,ntlml 236 write (*,*) ' Molecule #',i,' old new'236 write (logString, *) ' Molecule #',i,' old new' 237 237 do j=1,3 238 write (*,*) gbpro(j,i),' ',gbpr(j,i)238 write (logString, *) gbpro(j,i),' ',gbpr(j,i) 239 239 enddo 240 240 do j=4,6 241 write (*,*) gbpro(j,i)*crd,' ',gbpr(j,i)*crd241 write (logString, *) gbpro(j,i)*crd,' ',gbpr(j,i)*crd 242 242 enddo 243 243 enddo … … 245 245 endif 246 246 247 write ( *,'(/,2a)') ' Gradient ',247 write (logString,'(/,2a)') ' Gradient ', 248 248 & '______________________________________________________________' 249 249 250 write ( *,'(8(1x,f8.3))') (gdvr(i),i=1,nv)250 write (logString,'(8(1x,f8.3))') (gdvr(i),i=1,nv) 251 251 252 252 if (ireg.ne.0) then 253 253 254 write (*,*) ' -------------- global variables ------------' 255 write (*,'(6(1x,f8.3))') (gdvr(i+nv),i=1,ngbvr) 254 write (logString, *) 255 & ' -------------- global variables ------------' 256 write (logString, '(6(1x,f8.3))') (gdvr(i+nv),i=1,ngbvr) 256 257 257 258 endif … … 336 337 esm=eysm 337 338 338 write ( *,'(a,i5,a,2(e13.6,a))') ' Step ',nop,': energy ',esm339 & ,' (',gdsmey,' )'339 write (logString, '(a,i5,a,2(e13.6,a))') ' Step ',nop, 340 & ': energy ', esm, ' (',gdsmey,' )' 340 341 341 342 else … … 355 356 esm=wtey*eysm+wtrg*eyrg 356 357 357 write ( *,'(a,i5,a,3(e13.6,a))') ' Step ',nop,': energy ',esm358 & 358 write (logString, '(a,i5,a,3(e13.6,a))') ' Step ',nop, 359 & ': energy ',esm,' (',gdsmey,',',gdsmrg,' )' 359 360 360 361 endif
Note:
See TracChangeset
for help on using the changeset viewer.