Changeset 38d77eb for minim.f


Ignore:
Timestamp:
11/19/09 11:29:48 (14 years ago)
Author:
baerbaer <baerbaer@…>
Branches:
master
Children:
7137e5d
Parents:
32289cd
Message:

Redirected standard out to logString.

SMMP produced a lot of log messages. This became an issue when run in massively
parallel environments. I replaced all writes to standard out to a write to logString.
The next step is to pass this string to a function that writes the messages to a log
file according to their importance and the chosen log level.

git-svn-id: svn+ssh://svn.berlios.de/svnroot/repos/smmp/trunk@34 26dc1dd8-5c4e-0410-9ffe-d298b4865968

File:
1 edited

Legend:

Unmodified
Added
Removed
  • minim.f

    r32289cd r38d77eb  
    5858
    5959      if (msv.gt.msvmx) then
    60         write (*,'(a,i5)') ' minreg> parameter MSV > ',msvmx
     60        write (logString, '(a,i5)') ' minreg> parameter MSV > ',msvmx
    6161        stop
    6262      endif
     
    6666      call gradient()
    6767
    68       write(*,'(/,a,/)') ' Energy BEFORE minimization:'
     68      write (logString, '(/,a,/)') ' Energy BEFORE minimization:'
    6969
    7070      if (ireg.eq.0) then
    7171
    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,
    7474     &    '   Coulomb: ',eyel,' Lennard-Jones: ',eyvw,' HB: ',eyhb,
    7575     &    '   Variables: ',eyvr,'  Solvatation: ',eysl
     
    7878
    7979
    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,
    8282     &    '   Coulomb: ',eyel,' Lennard-Jones: ',eyvw,' HB: ',eyhb,
    8383     &    '   Variables: ',eyvr,'  Solvatation: ',eysl,
     
    148148!          gdrg2 = max(acc,gdrg2)
    149149!          wtrg = wtrg * sqrt(gdey2/gdrg2)
    150 !          write(*,*)  ' -->    Wt_energy = ',wtey,'  Wt_regul. = ',wtrg
    151 !          write(*,*)  '  '
     150!          write (logString, *)  ' -->    Wt_energy = ',wtey,'  Wt_regul. = ',wtrg
     151!          write (logString, *)  '  '
    152152!        endif
    153153
     
    184184
    185185      if (nop.lt.mxop) then
    186         write (*,'(a)') ' ---- CONVERGENCE ----'
     186        write (logString, '(a)') ' ---- CONVERGENCE ----'
    187187      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 ',
    193193     & '__________________________________________________'
    194194
     
    197197      if (ireg.eq.0) then
    198198
    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: ',eysl
     199        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
    202202
    203203      else
    204204
    205         write (*,'(a,e12.5,/,3(a,e11.4),/,3(a,e11.4))') ' Total: ',
    206      &      wtey*eysm + wtrg*eyrg,
     205        write (logString, '(a,e12.5,/,3(a,e11.4),/,3(a,e11.4))')
     206     &  ' Total: ', wtey*eysm + wtrg*eyrg,
    207207     &  '   Coulomb: ',eyel,' Lennard-Jones: ',eyvw,' HB: ',eyhb,
    208208     &  '   Variables: ',eyvr,'  Solvatation: ',eysl,
     
    211211      endif
    212212
    213       write (*,'(/,a,/)') ' Variables _________________'
     213      write (logString, '(/,a,/)') ' Variables _________________'
    214214
    215215      nv = 0
     
    221221          if (abs(vr).gt.pi) vr=vr-sign(pi2,vr)
    222222
    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,')'
    225225
    226226          vlvr(i) = vr
     
    231231      if (ireg.ne.0) then
    232232
    233         write (*,'(/,a,/)') ' Global Variables ___________'
     233        write (logString,'(/,a,/)') ' Global Variables ___________'
    234234
    235235        do i=1,ntlml
    236           write(*,*) ' Molecule #',i,' old          new'
     236          write (logString, *) ' Molecule #',i,' old          new'
    237237          do j=1,3
    238             write(*,*) gbpro(j,i),' ',gbpr(j,i)
     238            write (logString, *) gbpro(j,i),' ',gbpr(j,i)
    239239          enddo
    240240          do j=4,6
    241             write(*,*) gbpro(j,i)*crd,' ',gbpr(j,i)*crd
     241            write (logString, *) gbpro(j,i)*crd,' ',gbpr(j,i)*crd
    242242          enddo
    243243        enddo
     
    245245      endif
    246246
    247       write (*,'(/,2a)') ' Gradient ',
     247      write (logString,'(/,2a)') ' Gradient ',
    248248     & '______________________________________________________________'
    249249
    250       write (*,'(8(1x,f8.3))') (gdvr(i),i=1,nv)
     250      write (logString,'(8(1x,f8.3))') (gdvr(i),i=1,nv)
    251251
    252252      if (ireg.ne.0) then
    253253
    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)
    256257
    257258      endif
     
    336337        esm=eysm
    337338
    338         write (*,'(a,i5,a,2(e13.6,a))') ' Step ',nop,': energy ',esm
    339      &                                 ,'  (',gdsmey,' )'
     339        write (logString, '(a,i5,a,2(e13.6,a))') ' Step ',nop,
     340     &                    ': energy ', esm, '  (',gdsmey,' )'
    340341
    341342      else
     
    355356        esm=wtey*eysm+wtrg*eyrg
    356357
    357         write (*,'(a,i5,a,3(e13.6,a))') ' Step ',nop,': energy ',esm
    358      &                                 ,'  (',gdsmey,',',gdsmrg,' )'
     358        write (logString, '(a,i5,a,3(e13.6,a))') ' Step ',nop,
     359     &             ': energy ',esm,'  (',gdsmey,',',gdsmrg,' )'
    359360
    360361      endif
Note: See TracChangeset for help on using the changeset viewer.