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       call mygetenv('LIPTRANPAR',liptranname)
  52       open (iliptranpar,file=liptranname,status='old',action='read')
  53 #ifndef OLDSCP
  54 C
  55 C 8/9/01 In the newest version SCp interaction constants are read from a file
  56 C Use -DOLDSCP to use hard-coded constants instead.
  57 C
  58       call mygetenv('SCPPAR',scpname)
  59       open (iscpp,file=scpname,status='old')
  60 #endif
  61 #ifdef MPL
  62       if (MyID.eq.BossID) then
  63       MyRank = MyID/fgProcs
  64 #endif
  65 #ifdef MPI
  66       print *,'OpenUnits: processor',MyRank
  67       call numstr(MyRank,liczba)
  68       outname=prefix(:lenpre)//'.out_'//pot(:lenpot)//liczba
  69 #else
  70       outname=prefix(:lenpre)//'.out_'//pot(:lenpot)
  71 #endif
  72       open(iout,file=outname,status='unknown')
  73       write (iout,'(80(1h-))')
  74       write (iout,'(30x,a)') "FILE ASSIGNMENT"
  75       write (iout,'(80(1h-))')
  76       write (iout,*) "Input file                      : ",
  77      &  prefix(:ilen(prefix))//'.inp'
  78       write (iout,*) "Output file                     : ",
  79      &  outname(:ilen(outname))
  80       write (iout,*)
  81       write (iout,*) "Sidechain potential file        : ",
  82      &  sidename(:ilen(sidename))
  83 #ifndef OLDSCP
  84       write (iout,*) "SCp potential file              : ",
  85      &  scpname(:ilen(scpname))
  86 #endif
  87       write (iout,*) "Electrostatic potential file    : ",
  88      &  elename(:ilen(elename))
  89       write (iout,*) "Cumulant coefficient file       : ",
  90      &  fouriername(:ilen(fouriername))
  91       write (iout,*) "Torsional parameter file        : ",
  92      &  torname(:ilen(torname))
  93       write (iout,*) "Double torsional parameter file : ",
  94      &  tordname(:ilen(tordname))
  95       write (iout,*) "Backbone-rotamer parameter file : ",
  96      &  sccorname(:ilen(sccorname))
  97       write (iout,*) "Bond & inertia constant file    : ",
  98      &  bondname(:ilen(bondname))
  99       write (iout,*) "Bending parameter file          : ",
 100      &  thetname(:ilen(thetname))
 101       write (iout,*) "Rotamer parameter file          : ",
 102      &  rotname(:ilen(rotname))
 103       write (iout,'(80(1h-))')
 104       write (iout,*)
 105       return
 106       end
 107