# File download/unzip written 2012 by Lenna X. Peterson (arklenna@gmail.com)
# Dictionary extraction written 2011 by Hongbo Zhu
#
# This code is part of the Biopython distribution and governed by its
# license.  Please see the LICENSE file that should have been included
# as part of this package.

"""Download PDB Chemical Component Dictionary and generate dict.

Download and parse PDB Chemical Component Dictionary,
then write out dict for to_one_letter_code.
"""

from __future__ import print_function

import gzip
import inspect
import os
import warnings

from Bio._py3k import urlopen as _urlopen

url = "ftp://ftp.wwpdb.org/pub/pdb/data/monomers/components.cif.gz"

# extract name of gzip file
gzname = os.path.basename(url)
# extract name of cif file (split by sep, remove last, rejoin)
cifname = os.extsep.join(gzname.split(os.extsep)[:-1])

url_handle = urlopen(url)

with open(gzname, 'wb') as gzh:
    print("Downloading file... (approx. 29 MB)")
    while True:
        data = url_handle.read(1024)
        if len(data) == 0:
            break
        gzh.write(data)

# size as of 13 April 2012
if os.path.getsize(gzname) < 29944258:
    warnings.warn("ERROR: Downloaded file is too small",
                  RuntimeWarning)

fh = gzip.open(gzname, 'rb')

# write extracted file to disk (not necessary)
#with open(cifname, 'wb') as cifh:
    #print("Extracting file...")
    #cifh.write(fh.read())

# The following code written by Hongbo Zhu
# generate three_to_one_dict
# two records in PDB Chemical Component Dictionary are parsed to
# generate the dictionary:
# _chem_comp.one_letter_code
# _chem_comp.three_letter_code

three_to_one_buf = []      # all three-letter codes
three_to_one_buf_noq = []  # only those with non-'?' one-letter codes

current_line = 'to_one_letter_code = {'
current_line_noq = 'to_one_letter_code = {'

found_one = False    # found one-letter code
found_three = False  # found three-letter code

counter = 0
counter_noq = 0

line = fh.readline()

while line:
    if line.startswith('_chem_comp.one_letter_code'):
        one = line.strip().split()[-1]
        found_one = True
    if line.startswith('_chem_comp.three_letter_code'):
        three = '%-3s' % (line.strip().split()[-1],)  # make it three-letter
        found_three = True

    if found_one and found_three:
        if counter % 5 == 0:
            three_to_one_buf.append('%s\n' % (current_line,))
            current_line = '    '

        current_line = '%s\'%s\':\'%s\',' % (current_line, three, one)
        counter += 1

        if one != '?':
            if counter_noq % 5 == 0:
                three_to_one_buf_noq.append('%s\n' % (current_line_noq,))
                current_line_noq = '    '

            current_line_noq = '%s\'%s\':\'%s\',' % (current_line_noq, three, one)
            counter_noq += 1

        found_one = False
        found_three = False

    line = fh.readline()

if len(current_line) < 5:
    three_to_one_buf[-1] = three_to_one_buf[:-1]  # remove the last comma
    three_to_one_buf.append('}')
else:
    three_to_one_buf.append('%s }' % (current_line[:-1]))

if len(current_line_noq) < 5:
    three_to_one_buf_noq[-1] = three_to_one_buf_noq[:-1]
    three_to_one_buf_noq.append('}')
else:
    three_to_one_buf_noq.append('%s }' % (current_line_noq[:-1]))

# Find path of current script
_scriptPath = os.path.abspath(os.path.split(inspect.getfile(inspect.currentframe()))[0])
# Path to SCOP module
_rafPath = os.path.normpath(os.path.join(_scriptPath, "..", "..", "Bio", "SCOP"))
_threeAllPath = os.path.join(_rafPath, 'three_to_one_all.py')
_threePath = os.path.join(_rafPath, 'three_to_one_dict.py')

#with open(_threeAllPath, 'w') as fh:
    #fh.writelines(three_to_one_buf)
with open(_threePath, 'w') as fh:
    fh.writelines(three_to_one_buf_noq)