Question: how to solve coupled pde to using pdsolve get plot and table values

Variation of (a) Skin friction ∂W/∂Z​, (b) Heat Transfer ∂θ/∂Z​, and (c) Mass Transfer ∂ϕ/∂Z​ for γ=10.0, Pr=7.0, ε=1.0, Nt=0.4, and Nb=0.2.

how to get this values by solving the PDE by using the pdsolve method 

 # Equation 37 (continuity): Simplified

    pde1 := diff(U(X,Z), X) + diff(W(X,Z), Z) = 0;

    # Equation 38 (momentum): Simplified

    pde2 := U(X,Z)*diff(U(X,Z), X) + W(X,Z)*diff(U(X,Z), Z) = 

           diff(U(X,Z), Z, Z) - theta(X,Z) - phi(X,Z);

    # Equation 39 (energy): Simplified

    pde3 := U(X,Z)*diff(theta(X,Z), X) + W(X,Z)*diff(theta(X,Z), Z) = 

           (1/Pr)*diff(theta(X,Z), Z, Z) + 

           Nb*diff(phi(X,Z), Z)*diff(theta(X,Z), Z) + 

           Nt*diff(theta(X,Z), Z)^2;

    # Equation 40 (concentration): Simplified

    pde4 := U(X,Z)*diff(phi(X,Z), X) + W(X,Z)*diff(phi(X,Z), Z) = 

           (1/Sc_val)*diff(phi(X,Z), Z, Z) + 

           (Nt/Nb)*diff(theta(X,Z), Z, Z);

    # Boundary conditions

    bcs := {

        # At Z = 0

        U(X,0) = 0,

        W(X,0) = 0,

        theta(X,0) = 1,

        phi(X,0) = 1,

        # At Z = Z_max

        U(X,Z_max) = 0,

        theta(X,Z_max) = 0,

        phi(X,Z_max) = 0,

    };

Variation of (a)W(b) θ and (c) ϕ for different value of γ when ε = 10.0, Nb = 0.4, ε = 10.0, Sc = 0.5 and Pr = 7.0.
consider X as 0.1