In [7]:
import pynoddy.experiment
reload(pynoddy.experiment)
Out[7]:
In [8]:
import os
In [9]:
from pynoddy.experiment import MonteCarlo
In [13]:
#setup
pynoddy.ensure_discrete_volumes = True
###################################################
#MONTE CARLO PERTURBATION OF HIS FILE EXAMPLE
###################################################
#setup working directory
# os.chdir(r'C:\Users\Sam\OneDrive\Documents\Masters\Models\Primitive\monte carlo test')
os.chdir("/Users/flow/git/pynoddy/sandbox")
# his_file = "foldUC.his"
his_file = "simple_two_faults_no_gps.his"
# params_file = "foldUC_params.csv"
params_file = "params.csv"
#create new MonteCarlo experiment
mc = MonteCarlo.MonteCarlo(his_file,params_file)
param_stats = [{'event' : 3,
'parameter': 'Slip',
'stdev': 1000.0,
'type': 'normal'}]
mc.set_parameter_statistics(param_stats)
mc.freeze()
#generate 100 random perturbations using 4 separate threads (in TOPOLOGY mode)
output_name = "mc_out"
n = 10
mc.generate_model_instances(output_name,n) # ,threads=4)
In [31]:
In [ ]: