use names
!-----------------------------------------------------------------------------
! Max. number of energy intervals
- integer,parameter :: max_ene=21 !10
+ integer,parameter :: max_ene=49 !10
!-----------------------------------------------------------------------------
! Maximum number of terms in SC bond-stretching potential
integer,parameter :: maxbondterm=3
integer,dimension(:),allocatable :: idih_constr,idih_nconstr,itheta_constr !(maxdih_constr)
integer :: idihconstr_start,idihconstr_end, &
ithetaconstr_start,ithetaconstr_end
- real(kind=8) :: ftors
+! real(kind=8) :: ftors
real(kind=8),dimension(:),allocatable :: drange,theta_constr0,theta_drange !(maxdih_constr)
- real(kind=8),dimension(:),allocatable :: phi0 !(maxdih_constr)
+ real(kind=8),dimension(:),allocatable :: phi0,ftors !(maxdih_constr)
real(kind=8),dimension(:),allocatable :: for_thet_constr !(maxdih_constr)
!-----------------------------------------------------------------------------
sigma_scpho, &
sigmap1_scpho,sigmap2_scpho,&
aa_scpho,bb_scpho,wqq_scpho,epsintab_scpho,alphapol_scpho,&
- rborn_scphoi,rborn_scphoj,dhead_scphoi
+ rborn_scphoi,rborn_scphoj,dhead_scphoi,alphi_scpho
real(kind=8),dimension(:,:), allocatable :: alphasur_scpho, &
chi_scpho,chipp_scpho,chis_scpho, &