def METHOD_NAME(sele,preserve=0):
    global state
    global model
    state = State()
    model = cmd.get_model(sele)
    sm = 0
    for a in model.atom:
        a.resi = str(a.resi_number)
        sm = sm + a.partial_charge
    print(" lig: net charge on ligand  is %8.4f\n" % sm)