source/wham/src-M/openunits.F

   1       subroutine openunits
   2 #ifdef WIN
   3       use dfport
   4 #endif
   5       implicit real*8 (a-h,o-z)
   6       include 'DIMENSIONS'
   7       include 'DIMENSIONS.ZSCOPT'
   8 #ifdef MPI
   9       include 'mpif.h'
  10       include 'COMMON.MPI'
  11       integer MyRank
  12       character*3 liczba
  13 #endif
  14       include 'COMMON.IOUNITS'
  15       integer lenpre,lenpot,ilen
  16       external ilen
  17
  18 #ifdef MPI
  19       MyRank=Me
  20 #endif
  21       call mygetenv('PREFIX',prefix)
  22       call mygetenv('SCRATCHDIR',scratchdir)
  23       call mygetenv('POT',pot)
  24       lenpre=ilen(prefix)
  25       lenpot=ilen(pot)
  26       call mygetenv('POT',pot)
  27       entname=prefix(:lenpre)//'_'//pot(:lenpot)//'.entr'
  28 C Get the names and open the input files
  29       open (1,file=prefix(:ilen(prefix))//'.inp',status='old')
  30 C Get parameter filenames and open the parameter files.
  31       call mygetenv('BONDPAR',bondname)
  32       open (ibond,file=bondname,status='old')
  33       call mygetenv('THETPAR',thetname)
  34       open (ithep,file=thetname,status='old')
  35       call mygetenv('ROTPAR',rotname)
  36       open (irotam,file=rotname,status='old')
  37       call mygetenv('TORPAR',torname)
  38       open (itorp,file=torname,status='old')
  39       call mygetenv('TORDPAR',tordname)
  40       open (itordp,file=tordname,status='old')
  41       call mygetenv('FOURIER',fouriername)
  42       open (ifourier,file=fouriername,status='old')
  43       call mygetenv('SCCORPAR',sccorname)
  44       open (isccor,file=sccorname,status='old')
  45       call mygetenv('ELEPAR',elename)
  46       open (ielep,file=elename,status='old')
  47       call mygetenv('SIDEPAR',sidename)
  48       open (isidep,file=sidename,status='old')
  49       call mygetenv('SIDEP',sidepname)
  50       open (isidep1,file=sidepname,status="old")
  51 #ifndef OLDSCP
  52 C
  53 C 8/9/01 In the newest version SCp interaction constants are read from a file
  54 C Use -DOLDSCP to use hard-coded constants instead.
  55 C
  56       call mygetenv('SCPPAR',scpname)
  57       open (iscpp,file=scpname,status='old')
  58 #endif
  59 #ifdef MPL
  60       if (MyID.eq.BossID) then
  61       MyRank = MyID/fgProcs
  62 #endif
  63 #ifdef MPI
  64       print *,'OpenUnits: processor',MyRank
  65       call numstr(MyRank,liczba)
  66       outname=prefix(:lenpre)//'.out_'//pot(:lenpot)//liczba
  67 #else
  68       outname=prefix(:lenpre)//'.out_'//pot(:lenpot)
  69 #endif
  70       open(iout,file=outname,status='unknown')
  71       write (iout,'(80(1h-))')
  72       write (iout,'(30x,a)') "FILE ASSIGNMENT"
  73       write (iout,'(80(1h-))')
  74       write (iout,*) "Input file                      : ",
  75      &  prefix(:ilen(prefix))//'.inp'
  76       write (iout,*) "Output file                     : ",
  77      &  outname(:ilen(outname))
  78       write (iout,*)
  79       write (iout,*) "Sidechain potential file        : ",
  80      &  sidename(:ilen(sidename))
  81 #ifndef OLDSCP
  82       write (iout,*) "SCp potential file              : ",
  83      &  scpname(:ilen(scpname))
  84 #endif
  85       write (iout,*) "Electrostatic potential file    : ",
  86      &  elename(:ilen(elename))
  87       write (iout,*) "Cumulant coefficient file       : ",
  88      &  fouriername(:ilen(fouriername))
  89       write (iout,*) "Torsional parameter file        : ",
  90      &  torname(:ilen(torname))
  91       write (iout,*) "Double torsional parameter file : ",
  92      &  tordname(:ilen(tordname))
  93       write (iout,*) "Backbone-rotamer parameter file : ",
  94      &  sccorname(:ilen(sccorname))
  95       write (iout,*) "Bond & inertia constant file    : ",
  96      &  bondname(:ilen(bondname))
  97       write (iout,*) "Bending parameter file          : ",
  98      &  thetname(:ilen(thetname))
  99       write (iout,*) "Rotamer parameter file          : ",
 100      &  rotname(:ilen(rotname))
 101       write (iout,'(80(1h-))')
 102       write (iout,*)
 103       return
 104       end
 105