how combined contour plot and 3d plo

figure;
                             figure
restart;
K := -1;
                               -1
mu := 1;
                               1
a := 1;
                               1
h := 1;
                               1
delta := 1;
                               1
r := 1;
                               1
m := 0;
                               0
l := 1;
                               1

                               3
omega := -a*r^2-delta*l*r-l^2*mu+2*K*a+2*K*delta+2*K*mu;
                               -9
xi := x+y-(-2*a*r-delta*l-delta*r-2*l*mu)*t;
                         x + y + 6 t
eta := -l*y+omega*t-r*x;
                        -9 t- x - y

u[0] :=sqrt(2)*sqrt(h*(a+delta+mu))/(h*tanh(xi));
                          (1/2)  (1/2)  
                         2      3       
                       -----------------
                       tanh(x + y + 6 t)

                         -10
Q := exp(sqrt(-1)*eta)*u[0];
                                     (1/2)  (1/2)
              exp(I (-9 t - x - y)) 2      3     
              -----------------------------------
                      tanh(x + y + 6 t)         
y := 0;
                            0
plot3d(Re(Q), x = -10 .. 10, t = -10 .. 10);

plot3d(Im(Q), x = -10 .. 10, t = -10 .. 10);

plot3d(abs(Q), x = -10 .. 10, t = -10 .. 10);

t := 0;
                            0
plot(abs(Q), x = -10 .. 10);

set2;
                             set2
restart;
{m = m, omega = -a*r^2-delta*l*r-l^2*mu+2*K*a+2*K*delta+2*K*mu, alpha[0] = -2*m*(a+delta+mu)/(h*RootOf(_Z^2*h-2*a-2*delta-2*mu)), alpha[1] = RootOf(_Z^2*h-2*a-2*delta-2*mu), beta[1] = 0};
 /       
 |m = m, 
<        
 |       
 \       

              2                2                                  
  omega = -a r  - delta l r - l  mu + 2 K a + 2 K delta + 2 K mu, 

                        2 m (a + delta + mu)           
  alpha[0] = - --------------------------------------, 
                       /  2                         \  
               h RootOf\_Z  h - 2 a - 2 delta - 2 mu/  

                   /  2                         \             \ 
  alpha[1] = RootOf\_Z  h - 2 a - 2 delta - 2 mu/, beta[1] = 0| 
                                                           >
                                                              | 
                                                           / 

figure;
                            figure
restart;
K := -1;
                              -1
mu := 1;
                              1
a := 1;
                              1
h := 1;
                               1
delta := 1;
                               1
r := 1;
                               1
m := 0;
                               0
l := 1;
                             1

                               3
omega := -a*r^2-delta*l*r-l^2*mu+2*K*a+2*K*delta+2*K*mu;
                               -9
xi := x+y-(-2*a*r-delta*l-delta*r-2*l*mu)*t;
                          x + y + 6 t
eta := -l*y+omega*t-r*x;
                          -9 t - x - y

u[0] := -m*(a+delta+mu)*sqrt(2)/sqrt(h*(a+delta+mu))+sqrt(2)*sqrt(h*(a+delta+mu))*(m-tanh(xi))/h;
                  (1/2)  (1/2)                  
                -2     3      tanh(x + y + 6 t)

                            -10
Q := exp(sqrt(-1)*eta)*u[0];
                          (1/2)  (1/2)                
    -exp(I (-9 t - x - y)) 2      3      tanh(x + y + 6 t)
y := 0;
                            0
plot3d(Re(Q), x = -10 .. 10, t = -10 .. 10);

plot3d(Im(Q), x = -10 .. 10, t = -10 .. 10);

plot3d(abs(Q), x = -10 .. 10, t = -10 .. 10);

t := 0;
                               0
plot(abs(Q), x = -10 .. 10);