module defs
 integer, parameter :: i4b = selected_int_kind(9)
 integer, parameter :: i2b = selected_int_kind(4)
 integer, parameter :: ip = kind(1)
 integer, parameter :: sp = kind(1.0)
 integer, parameter :: dp = kind(1.0d0)
 integer(i4b), dimension(8) :: scalefactors                        &
 = (/ 1, 10, 100, 0, 1, 1, 1, 1 /) 
 real(dp), parameter :: ncab = 1.45_dp                             & 
 , cacpb = 1.52_dp, cpnb = 1.33_dp, cpob = 1.23_dp, nhb = 1.0_dp    &
 , cacbb = 1.54_dp, cahb = 1.07_dp                                  &
! the hydrogen bond is something like 5.4 - 1.23 -1.0 = 3.17 A; b = 3.808
 , ancacp = 1.22173047639603070383_dp & ! 180 - 110 = 70 degrees
 , an_ca_cp = 1.91986217719376253460_dp & ! 110 degrees 
 , acacpn = 1.12399203828434824753_dp & ! 180 - 115.6 = 64.4 degrees
 , aca_cp_n = 2.01760061530544499090_dp & ! 115.6 degrees
 , acpnca = 1.01403629540870548418_dp & ! 180 - 121.9 = 58.1 degrees
 , acp_n_ca = 2.12755635818108775425_dp & ! 121.9 degrees
 , acacpo = 1.02799892942466012080_dp & ! 180 - 121.1 = 58.9 degrees
 , aca_cp_o = 2.11359372416513311763_dp & ! 121.1 degrees
 , acpnh = 1.05592419745656939403_dp &  ! 180 - 119.5 = 60.5 degrees
 , acp_n_h = 2.08566845613322384440_dp & ! 119.5 degrees
 , pi = 3.14159265358979323844_dp                                   &
 , twopi = 6.28318530717958647688_dp, eps = 0.01_dp                 &
 , twopiover3 = 2.09439510239319549229_dp                           &
 , sin37 = 0.59972214027798415249_dp &  ! sin(36.85)
 , cos37 = 0.80020831941463461536_dp &  ! cos(36.85)
 , cos106 = -0.28066670892078779654_dp  ! cos(106.3)
!  real(dp), dimension(3,3), parameter :: Id = (/ 1.0_dp, 0.0_dp, 0.0_dp, &
!     0.0_dp, 1.0_dp, 0.0_dp, 0.0_dp, 0.0_dp, 1.0_dp /) 
!  Why doesn't this work?  One must use reshape!
    ! -75 < phi=-62.5 < -50  GLY Zone Parameters
 real(dp), parameter :: phigly1b = -75*pi/180, phigly1 = -62.5*pi/180, &
  phigly1t = -50*pi/180
    ! -60 < psi=-40 < -20
 real(dp), parameter :: psigly1b = -60*pi/180, psigly1 = -40*pi/180, &
  psigly1t = -20*pi/180

    !  60 < phi=80 < 100     
 real(dp), parameter :: phigly2b = 60*pi/180, phigly2 = 80*pi/180, &
  phigly2t = 100*pi/180
    ! -30 < psi=7.5 < 45
 real(dp), parameter :: psigly2b = -30*pi/180, psigly2 = 7.5*pi/180, &
  psigly2t = 45*pi/180

   !   60 < phi=120 < 180    
 real(dp), parameter :: phigly3b = 60*pi/180, phigly3 = 120*pi/180, &
  phigly3t=pi
   !  120 < psi=150 < 180  
 real(dp), parameter :: psigly3b = 120*pi/180, psigly3 = 150*pi/180, &
  psigly3t = pi

   ! -180 < phi=-120 < -60   
 real(dp), parameter :: phigly4b = -pi, phigly4 = -120*pi/180, &
  phigly4t = -60*pi/180
   ! -180 < psi=-150 < -120  
 real(dp), parameter :: psigly4b = -pi, psigly4 = -150*pi/180, &
  psigly4t = -120*pi/180

   ! -80 < phi=-60 < -50    PRO Zone 1  
 real(dp), parameter :: phipro1b = -65*pi/180, phipro1 = -58*pi/180, &
  phipro1t = -50*pi/180
   ! -50 < psi=-30 < -10   
 real(dp), parameter :: psipro1b = -50*pi/180, psipro1 = -30*pi/180, &
  psipro1t = -10*pi/180
   ! -80 < phi=-60 < -50    PRO Zone 2
 real(dp), parameter :: phipro2b = -65*pi/180, phipro2 = -58*pi/180, &
  phipro2t = -50*pi/180 
   ! 125 < psi=148 < 171  
 real(dp), parameter :: psipro2b = 125*pi/180, psipro2 = 148*pi/180, &
  psipro2t = 171*pi/180

   !  -90 < phi=-65 < -40    ALA Zone 1 
 real(dp), parameter :: phiala1b = -94*pi/180, phiala1 = -65*pi/180, &
  phiala1t = -33*pi/180
   !  -60 < psi=-32 < -4     
 real(dp), parameter :: psiala1b = -81*pi/180, psiala1 = -32*pi/180, &
  psiala1t = -4*pi/180 

   ! -170 < phi=-110 < -50    ALA Zone 2 
 real(dp), parameter :: phiala2b = -170*pi/180, phiala2 = -110*pi/180, &
  phiala2t = -50*pi/180
   !  115 < psi=145 < 175     
 real(dp), parameter :: psiala2b = 115*pi/180, psiala2 = 145*pi/180, &
  psiala2t = 175*pi/180 

! external ranlux
 end module defs