now running simulations that will not finish

def get_parametercollection(seed,crack_length, offset):

    par = mn.ParameterCollection()

    nx = 8 # needs to == el**h

    nx = 16 # needs to == el**h

    ny = nx

    nz = 5 + offset # need to be >=  h +2 + o_set (+2 bc bound cond)    

    nz = 6 + offset # need to be >=  h +2 + o_set (+2 bc bound cond)    

    n = nx*ny*nz

    par.set_N(n)

    par.set_Nx(nx)