def METHOD_NAME(ds: xr.Dataset, nudging_variables: list):
    for var in nudging_variables:
        ds[f"tendency_of_{var}_due_to_scream_physics"] = (
            ds[f"physics_{var}_tend"] - ds[f"nudging_{var}_tend"]
        )
    return ds