Hello all, I'm trying to get surface patches of PDB structures in VTK format. So I want to define a patch as the set of atoms within a given radius of a specific residue, build the surface corresponding to that patch and export the result in VTK format. Additionally, I want to do all this with a python script. Here the code I'm trying but doesn't work : *import pychimerapychimera.patch_environ()pychimera.enable_chimera()import osimport chimerafrom chimera import runCommand as rc* *list_resid = [31, 64]* *radius = 3list_pdbs = ['ensayo.pdb']for resid in list_resid: for line in list_pdbs: list_atoms = [] rc('open '+line) rc('sel :{} za<{}'.format(resid, radius)) for a in chimera.selection.currentAtoms(): list_atoms.append('serialNumber=' + str(a.serialNumber)) sel_atoms = ' or '.join(list_atoms) rc('surfcat patch @*/{}'.format(sel_atoms)) rc('surface patch') rc("export format VTK {}_resid_{}_radius_{}.vtk".format(line[:-4], resid, radius)) rc("close all")* Thanks in advance for your help. -- Diego A. Amaya Ramírez -- *Aviso legal:* El contenido de este mensaje y los archivos adjuntos son confidenciales y de uso exclusivo de la Universidad Nacional de Colombia. Se encuentran dirigidos sólo para el uso del destinatario al cual van enviados. La reproducción, lectura y/o copia se encuentran prohibidas a cualquier persona diferente a este y puede ser ilegal. Si usted lo ha recibido por error, infórmenos y elimínelo de su correo. Los Datos Personales serán tratados conforme a la Ley 1581 de 2012 y a nuestra Política de Datos Personales que podrá consultar en la página web www.unal.edu.co <http://www.unal.edu.co/>.* *Las opiniones, informaciones, conclusiones y cualquier otro tipo de dato contenido en este correo electrónico, no relacionados con la actividad de la Universidad Nacional de Colombia, se entenderá como personales y de ninguna manera son avaladas por la Universidad.