def METHOD_NAME(nx, ny, name="Unknown"):
    data = np.zeros([nx, ny])
    print("Reading 2D variable " + name)
    for i in np.arange(nx):
        data[i, :] = read_array(ny, name + "[" + str(i) + "]")
    return data