
c This function returns the uppercase atomic symbol of the atomic number N.
C It returns '  ' if the atomic number is not recognized.

      character*2 function AtomSymb(n)
      implicit none

      integer n

      character*2 SymbTable(111)
      data SymbTable /
     & 'H ','HE',
     & 'LI','BE',
     & 'B ','C ','N ','O ','F ','NE',
     & 'NA','MG',
     & 'AL','SI','P ','S ','CL','AR',
     & 'K ','CA',
     & 'SC','TI','V ','CR','MN','FE','CO','NI','CU','ZN',
     & 'GA','GE','AS','SE','BR','KR',
     & 'RB','SR',
     & 'Y ','ZR','NB','MO','TC','RU','RH','PD','AG','CD',
     & 'IN','SN','SB','TE','I ','XE',
     & 'CS','BA',
     & 'LA','CE','PR','ND','PM','SM','EU',
     &                               'GD','TB','DY','HO','ER','TM','YB',
     & 'LU','HF','TA','W ','RE','OS','IR','PT','AU','HG',
     & 'TL','PB','BI','PO','AT','RN',
     & 'FR','RA',
     & 'AC','TH','PA','U ','NP','PU','AM',
     &                               'CM','BK','CF','ES','FM','MD','NO',
     & 'LR','RF','DB','SG','BH','HS','MT','DS','RG'
     &               /

      AtomSymb = SymbTable(n)
      if (n.lt.1.or.111.lt.n) AtomSymb = '  '

      return
c     end function AtomSymb
      end

