From: Adam Kazimierz Sieradzan <adasko@sun1.chem.univ.gda.pl>
Date: Thu, 6 Dec 2012 09:31:09 +0000 (-0500)
Subject: zmiany w xdrfie
X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=commitdiff_plain;h=e04fbf112bb74a38b66d5c0db6cf0fe216b5885e;p=unres.git

zmiany w xdrfie
---

diff --git a/source/xdrfpdb/src-M/COMMON.NAMES b/source/xdrfpdb/src-M/COMMON.NAMES
index e6f926b..1c85c4f 100644
--- a/source/xdrfpdb/src-M/COMMON.NAMES
+++ b/source/xdrfpdb/src-M/COMMON.NAMES
@@ -1,6 +1,6 @@
       character*3 restyp
       character*1 onelet
-      common /names/ restyp(ntyp+1),onelet(ntyp+1)
+      common /names/ restyp(-ntyp1:ntyp1),onelet(-ntyp1:ntyp1)
       character*10 ename,wname
       integer nprint_ene,print_order
       common /namterm/ ename(n_ene),wname(n_ene),nprint_ene,
diff --git a/source/xdrfpdb/src-M/DIMENSIONS b/source/xdrfpdb/src-M/DIMENSIONS
index 32a8ae5..43e77a2 100644
--- a/source/xdrfpdb/src-M/DIMENSIONS
+++ b/source/xdrfpdb/src-M/DIMENSIONS
@@ -38,7 +38,7 @@ C Max. number of contacts per residue
       parameter (maxconts=maxres)
 C Number of AA types (at present only natural AA's will be handled
       integer ntyp,ntyp1
-      parameter (ntyp=20,ntyp1=ntyp+1)
+      parameter (ntyp=24,ntyp1=ntyp+1)
 C Max. number of types of dihedral angles & multiplicity of torsional barriers
 C and the number of terms in double torsionals
       integer maxtor,maxterm,maxlor,maxtermd_1,maxtermd_2
diff --git a/source/xdrfpdb/src-M/nazwy.f b/source/xdrfpdb/src-M/nazwy.f
index 9f61cfb..d495120 100644
--- a/source/xdrfpdb/src-M/nazwy.f
+++ b/source/xdrfpdb/src-M/nazwy.f
@@ -4,10 +4,17 @@
       include 'COMMON.NAMES'
       include 'COMMON.FFIELD'
       data restyp /
+     &'DD','DAU','DAI','DDB','DSM','DPR','DLY','DAR','DHI','DAS','DGL',
+     & 'DSG','DGN','DSN','DTH',
+     &'DYY','DAL','DTY','DTR','DVA','DLE','DIL','DPN','MED','DCY','ZER',
      &'CYS','MET','PHE','ILE','LEU','VAL','TRP','TYR','ALA','GLY','THR',
-     &'SER','GLN','ASN','GLU','ASP','HIS','ARG','LYS','PRO','D'/
+     &'SER','GLN','ASN','GLU','ASP','HIS','ARG','LYS','PRO','SME','DBZ',
+     &'AIB','ABU','D'/
       data onelet /
+     &'z','z','z','z','z','p','k','r','h','d','e','n','q','s','t','g',
+     &'a','y','w','v','l','i','f','m','c','x',
      &'C','M','F','I','L','V','W','Y','A','G','T',
-     &'S','Q','N','E','D','H','R','K','P','X'/
+     &'S','Q','N','E','D','H','R','K','P','z','z','z','z','X'/
+      data potname /'LJ','LJK','BP','GB','GBV'/
       end 
 
diff --git a/source/xdrfpdb/src-M/rescode.f b/source/xdrfpdb/src-M/rescode.f
index 1fff9ec..efa6c5c 100644
--- a/source/xdrfpdb/src-M/rescode.f
+++ b/source/xdrfpdb/src-M/rescode.f
@@ -7,7 +7,7 @@
 
       if (ione.eq.0) then
 
-      do i=1,ntyp1
+      do i=-ntyp1,ntyp1
         if (ucase(nam).eq.restyp(i)) then
           rescode=i
           return
@@ -15,7 +15,7 @@
       enddo
       else
 
-      do i=1,ntyp1
+      do i=-ntyp1,ntyp1
         if (nam(1:1).eq.onelet(i)) then
           rescode=i
           return