#!/usr/bin/env python


import glob
import os

cMs = [0.0, 0.0001, 0.0002, 0.0004, 0.0006, 0.0008, 0.001, 0.0012, 0.0015, 0.002, 0.0025, 0.003, 0.004, 0.005, 0.01, 0.02, 0.03]

for cM in cMs:

    print ('#### cM%6.4f ####' % (cM,))
 
    cg = glob.glob('cM%6.4f_cg/*pdb' % (cM,))
    aa = glob.glob('cM%6.4f_aa/*pdb' % (cM,))

    aabase = []
    for f in aa:
        aabase.append(os.path.basename(f))
    
    for f in cg:
        base = os.path.basename(f)
        if base in aabase:
            pass
        else:
            print base
