parameter (zero = 0.d0, half = 0.5d0, one = 1.d0, two = 2.d0
     &   , four = 4.d0, five = 5.d0, hc = 1.9876d-16
     &   , cs = 2.9979d10, rkb = 1.38d-16, amu = 1.66d-24
     &   , hck = 1.439565d0, rda = 1.745329d-2
     &   , pi = 3.141592653589793238462643d0
     &   , sqpi = 1.772453850905516027298167d0, nd6=3
     &   , nd1 = 101, nd2 = 3000, nd3=200, nd4=15000, nd5=100)