- character(len=1),dimension(10) :: chainid= (/'A','B','C','D','E','F','G','H','I','J'/)
+ character(len=1),dimension(58) :: chainid= (/'A','B','C','D','E','F','G','H','I','J',&
+ 'K','L','M','O','Q','P','R','S','T','U','W','X','Y','Z','a','b','c','d','e','f',&
+ 'g','h','i','j','k','l','m','n','o','p','r','s','t','u','w','x','y','z',&
+ '0','1','2','3','4','5','6','7','8','9'/)