def Deal_pdb(fileA,fileB,index,parameterA,replaceParaA,line_paraA,parameterB,replaceParaB,line_paraB):
from_file=open(fileA,"r")
to_file=open(fileB,"a")
lines=from_file.readlines()
from_file.close()
for line in lines:
line_data=line.split()
print len(lines)
if line_data[index]==parameterA and int(line_data[1])>line_paraA:
line_data[index]=replaceParaA
line=line.replace(" %s "%parameterA," %s "%replaceParaA)
to_file.write(line)
elif line_data[index]==parameterB and int(line_data[1])>line_paraB:
line_data[index]=replaceParaB
line=line.replace(" %s "%parameterB," %s "%replaceParaB)
to_file.write(line)
else:
to_file.write(line)
to_file.close()
if __name__=="__main__":
Deal_pdb("3BEQ-tet.pdb","1.pdb",4,"A","C",2983,"B","D",9856)