include 'COMMON.IOUNITS'
include 'COMMON.CHAIN'
include 'COMMON.DFA'
+ include 'COMMON.FFIELD'
C NOTE THAT FILENAMES are FIXED, CURRENTLY!!
open(iodfa, file = 'dist_dfa.dat', status = 'old', err=33)
goto 34
33 write(iout,'(a)') 'Error opening dist_dfa.dat file'
+ print *,'Error opening dist_dfa.dat file'
stop
34 continue
write(iout,'(a)') 'dist_dfa.dat is opened!'
C PHI RESTRAINTS
open(iodfa, file='phi_dfa.dat',status='old',err=35)
goto 36
- 35 write(iout,'(a)') 'Error opening dist_dfa.dat file'
+ 35 write(iout,'(a)') 'Error opening phi_dfa.dat file'
+ print *, 'Error opening phi_dfa.dat file'
stop
36 continue
C THETA RESTRAINTS
open(iodfa, file='theta_dfa.dat',status='old',err=41)
goto 42
- 41 write(iout,'(a)') 'Error opening dist_dfa.dat file'
+ 41 write(iout,'(a)') 'Error opening theta_dfa.dat file'
+ print *,'Error opening theta_dfa.dat file'
stop
42 continue
write(iout,'(a)') 'theta_dfa.dat is opened!'
open(iodfa,file='nei_dfa.dat',status='old',err=37)
goto 38
37 write(iout,'(a)') 'Error opening nei_dfa.dat file'
+ print *,'Error opening nei_dfa.dat file'
stop
38 continue
write(iout,'(a)') 'nei_dfa.dat is opened!'
C END OF NEIGHBORING CA
C READ BETA RESTRAINT
+ if (wdfa_beta.eq.0.0) return
open(iodfa, file='beta_dfa.dat',status='old',err=39)
goto 40
39 write(iout,'(a)') 'Error opening beta_dfa.dat file'
+ print *,'Error opening beta_dfa.dat file'
stop
40 continue
write(iout,'(a)') 'beta_dfa.dat is opened!'