In [1]:

    

import pynoddy.experiment
reload(pynoddy.experiment)


    

    
Out[1]:





<module 'pynoddy.experiment' from '/Users/flow/git/pynoddy/pynoddy/experiment/__init__.pyc'>

In [2]:

    

import os


    

In [3]:

    

from pynoddy.experiment import MonteCarlo


    

In [6]:

    

#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)

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)


    

    




Constructing out_0001... 






    




---------------------------------------------------------------------------
KeyError                                  Traceback (most recent call last)
<ipython-input-6-18c2c5242f83> in <module>()
     23 output_name = "mc_out"
     24 n = 10
---> 25 mc.generate_model_instances(output_name,n) # ,threads=4)
     26 

/Users/flow/git/pynoddy/pynoddy/experiment/MonteCarlo.py in generate_model_instances(self, path, count, **kwds)
    251 
    252                 #do random perturbation
--> 253                 self.random_perturbation(verbose=vb)
    254 
    255                 #save history

/Users/flow/git/pynoddy/pynoddy/experiment/__init__.pyc in random_perturbation(self, **kwds)
    185             model_type = kwds.get('model_type', 'base')
    186             if model_type == 'base':
--> 187                 ori_val = self.base_events[param['event']].properties[param['parameter']]
    188             else:
    189                 ori_val = self.events[param['event']].properties[param['parameter']]

KeyError: 'parameter'

In [31]:

    




    

In [ ]: