runscript for running reference density 40, 50

import sys

import os

import numpy as np

import mechnet as mn

from multiprocessing import Pool

from lib.push_data import push_data_to, copy_runfile_to 

from lib.my_plot_edges import plot_edges_mechnet_network 

os.environ["OMP_NUM_THREADS"] = "1"

os.environ["OPENBLAS_NUM_THREADS"] = "1"

os.environ["MKL_NUM_THREADS"] = "1"

def get_parametercollection(seed,crack_length, offset):

    par = mn.ParameterCollection()

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

    simulation = mn.sim.Simulation(builder, solver, intrescal, outgen, applier, checker)

    dict_of_crs_edge_data = bounded_network.get_network_data() 

    simdata = simulation.start_simulation(dict_of_crs_edge_data, dict_of_simulation_parameters)

    simdata = simulation.start_simulation(dict_of_crs_edge_data, dict_of_simulation_parameters, maximum_runtime=2*8580)

    simdata_dict = simdata.get_simdata_dict()

    network_descriptor_dict = bounded_network.get_construction_data()

if __name__=="__main__":

    list_of_seeds =  []

    offset = 5

    crack_lengths = [20, 30, 40, 50]

    crack_lengths = [40, 50]

    for seed in range(1000, 1010):

            list_of_seeds.append(seed) 

    for crack_length in crack_lengths: