numerical solution for coupled pde system with periodic boundary conditions

I would like to solve:

s_t+[S(s,q,s_x,q_x,h_xxx)]_x=0

q_t+[B(s,q,s_x,q_x,h_xxx)]_x=G[s,q]

S and B include functions of s and q and their derivatives...