Needs["hirota`"]

n=3;

B[f_, g_]:=Dxt[3,1,x,t,f,g]-Dx[2,x,f,g]-Dxt[1,1,x,t,f,g];

name="Shallow_Water";

myfile="h_shwat.o";

OpenLog[myfile];

(* Do not run this case with the symbolic test for 4 solitons set to true *)
(* unless you have a lot of patience or a large computer *)

hirota[B,name,n,2,2,True,False];

CloseLog[];