Esempio n. 1
0
chunks = int(sys.argv[1])
chunkno = int(sys.argv[2])

suitelist = planetengine.utilities.suite_list(
    {
        'f': [x / 10. for x in range(1, 11)],
        'Ra': [10.**(x / 2.) for x in range(6, 15)],
        'eta0': [10.**(x / 2.) for x in range(11)],
    },
    shuffle=True,
    chunks=chunks)

localJobs = suitelist[chunkno]

planetengine.log("Starting chunk no# " + str(chunkno), 'logs')

for index, job in enumerate(localJobs):

    model = planetengine.frame.make_frame(
        modelscripts.arrhenius.build(**job, res=16),
        {'temperatureField': planetengine.initials.sinusoidal.IC(freq=1.)},
        outputPath)

    conditions = {
        'stopCondition':
        lambda: model.modeltime > 0.3,
        'collectConditions':
        lambda: model.step % 10 == 0,
        'checkpointCondition':
        lambda: any([
Esempio n. 2
0
chunks = int(sys.argv[1])
chunkno = int(sys.argv[2])
iterno = int(sys.argv[3])

suitelist = planetengine.utilities.suite_list(
    {
        'f': [0.2, 0.4, 0.6, 0.8, 1.],
        'eta': [1., 10., 100., 1000., 10000.],
        'Ra': [1e3, 1e4, 1e5, 1e6, 1e7],
    },
    shuffle=True,
    chunks=chunks)

localJob = suitelist[chunkno][iterno]

planetengine.log(
    "Starting chunk no# " + str(chunkno) + ", iter no# " + str(iterno), 'logs')

for index, job in enumerate(localJobs):

    model = planetengine.frame.make_frame(
        modelscripts.isovisc.build(**job, res=64),
        {'temperatureField': planetengine.initials.sinusoidal.IC(freq=1.)},
        outputPath)

    conditions = {
        'stopCondition':
        lambda: model.modeltime() > 0.3,
        'checkpointCondition':
        lambda: any([
            model.status == 'pre-traverse',
            model.step() % 1000 == 0,