A new procedure to design numerical schemes for coupled partial differential equations is proposed. The resulting schemes have discrete counterparts of conservative or dissipative quantity in original system. They also enjoy another welcome feature that they are constructed on staggered time meshes, by which each variables can be computed alternately with less computational costs than usual schemes. The procedure is demonstrated in the case of the coupled KdV equations.