- Timestamp:
- 11/19/09 11:29:48 (14 years ago)
- Branches:
- master
- Children:
- 7137e5d
- Parents:
- 32289cd
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
regul.f
r32289cd r38d77eb 44 44 wtey = 0.d0 45 45 46 write (*,'(/,a,2(a,f4.2),/)')46 write (logString, '(/,a,2(a,f4.2),/)') 47 47 & ' ====================== Regularization only', 48 48 & ' Wt(energy) = ',wtey,' Wt(regul.) = ',wtrg … … 50 50 call minim(2, nsteps, acc) 51 51 52 write (*,*) ' '53 write (*,*) ' ---------- contacts after 1st regularization'54 write (*,*) ' '52 write (logString, *) ' ' 53 write (logString, *) ' -------- contacts after 1st regularization' 54 write (logString, *) ' ' 55 55 call cnteny(nml) 56 write (*,*) ' '56 write (logString, *) ' ' 57 57 58 58 nrs = irsml2(nml)-irsml1(nml)+1 59 59 call rmsdopt(nml,1,nrs,ixatp,xatp,yatp,zatp,0,rm,av1,av2,rmsd) 60 60 61 write (*,*) ' RMSD = ',rmsd61 write (logString, *) ' RMSD = ',rmsd 62 62 63 63 ! --------------------------------------- fix vars. defined in PDB … … 70 70 ireg = 0 71 71 72 write (*,'(/,a,2(a,f4.2),/)')72 write (logString, '(/,a,2(a,f4.2),/)') 73 73 & ' ====================== Internal Energy for Hydrogens only', 74 74 & ' Wt(energy) = ',wtey,' Wt(regul.) = ',wtrg … … 76 76 call minim(1, nsteps, acc) 77 77 78 write (*,*) ' '79 write (*,*) ' ---------- contacs after Emin. for hydrogens'80 write (*,*) ' '78 write (logString, *) ' ' 79 write (logString, *) ' -------- contacs after Emin. for hydrogens' 80 write (logString, *) ' ' 81 81 call cnteny(nml) 82 82 … … 98 98 wtey = 1.d0 - wtrg 99 99 100 write (*,'(/,a,i2,2(a,e11.3),/)')100 write (logString, '(/,a,i2,2(a,e11.3),/)') 101 101 & ' ================ Minimization #',it, 102 102 & ' Wt(energy) = ',wtey,' Wt(regul.) = ',wtrg … … 107 107 call rmsdopt(nml,1,nrs,ixatp,xatp,yatp,zatp,0,rm,av1,av2,rmsd) 108 108 109 write (*,*) ' '110 write (*,*) ' RMSD = ',rmsd109 write (logString, *) ' ' 110 write (logString, *) ' RMSD = ',rmsd 111 111 112 112 enddo 113 113 114 write (*,*) ' '115 write (*,*) ' ---------- contacts after full regularization'116 write (*,*) ' '114 write (logString, *) ' ' 115 write (logString, *) ' ------- contacts after full regularization' 116 write (logString, *) ' ' 117 117 call cnteny(nml) 118 118 … … 120 120 121 121 ! Output of dihedral angles of the regularized structure 122 write(*,*) 'Dihedral angles of the regularized structure;' 122 write (logString, *) 123 & 'Dihedral angles of the regularized structure;' 123 124 call outvar(nml, 'regd.var') 124 125
Note:
See TracChangeset
for help on using the changeset viewer.