<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def vecs_to_datmesh(x, y):
 """ Converts input arguments x and y to a 2d meshgrid, suitable for calling Means, Covariances and Realizations. """

    x, y = meshgrid(x, y)
    out = zeros(x.shape + (2,), dtype=float)
    out[:, :, 0] = x
    out[:, :, 1] = y
    return out