Open goncalojustino opened 4 years ago
Terei então de adicionar ai o seu resid Entao talvez seria mais facil escrever apenas basic = universe.select_atoms("not resname SOL or not name C or not name H") e mais a baixo, na seleção fazer como foi feito no trabalho anterior e por algo do estilo : acid[i].atoms.FE3.position ?
Ou entao algo como: resname FE3 and name FE que depois dava, nos outros casos: resname [V3/VO/VO2/VO4] and name [V3/V4/V]
https://github.com/goncalojustino/md_code/blob/a345179fdf8d88e8f52278482469d9368ebff550/distancia_fe3.py#L15
a ideia é calcular distâncias a UM atomo de metal especifico, neste caso vai ter que ser resname FE3 e atom name adequado.