Writing calculated script data back to the block model

Hi Ronald,

I'm not a python expert but this might help.

In order to add a field to a file I normally use the following code. I altered this code from the code produced by recording the File-> Create or Modify File process. You could copy this function and put it at the beginning or end of your script. I tend to keep my functions in a separate python file and import as I would do any module. By making it a function you can call it repeatedly.

def modify_add_column(file_name, field_name, mmpy, field_type= "s",  width = "", precision= "", start_val= ""):
    field_types = ['n', 'c', 'r', 'f', 'l', 's']# The function takes these field types.
    dict_field_types = {ftype: str(i) for i, ftype in enumerate(field_types)}
    Modparms_FormSet1= mmpy.FormSet("MODPARMS","15.0.5.697")
    Modparms_FormSet1.set_field("TRUNCATE","0")
    Modparms_FormSet1.set_field("ROUNDOFF","1")
    Modparms_FormSet1.set_field("NUMERIC","0")
    Modparms_FormSet1.set_field("JUSTIFY","0")
    Modparms_FormSet1.set_field("TYPE1",dict_field_types[field_type])#Changes field type to number
    #Modparms_FormSet1.set_field("INC1","")
    if start_val is not "":
        Modparms_FormSet1.set_field("START1",str(start_val))
    Modparms_FormSet1.set_field("FUNC1","0")
    if precision is not "":
        Modparms_FormSet1.set_field("PREC1",str(precision))
    if width is not "":
        Modparms_FormSet1.set_field("WIDTH1",str(width))
    Modparms_FormSet1.set_field("VAR1",field_name)
    Modparms_FormSet1.set_field("FILETYPE","0")
    Modparms_FormSet1.set_field("FILE",file_name)
    try:
        Modparms_FormSet1.run()
    except:
        Exception #Python won't fall over if the field already exists 

The file_name does not include ".DAT". so if you want to use it with blkmdlout you should write blkmdlout[:-4]

You could use the function by calling it using the following line to add a column of Short type:

modify_add_column(blkmdlout[:-4], "orsource ", MMpy, field_type= "s")

You can add this line after you have created the file i.e. after you close the first block model file and before you open the blkmdlout file in Python.

It looks like the rest of your code is fine if you have the indentations correct except that you may need to replace

orsource = 0 #with something like:

orsource = Myblkmdl_new.get_field_id("orsource")

Just to be sure, the indentations at the end should look like this:

orsource = Myblkmdl_new.get_field_id("orsource")

for record_id in range(1, NumRecords+1):
    Currentzone = Myblkmdl_new.get_num_block_value(BlockModelZoneField, record_id)     
    Currentoxstate = Myblkmdl_new.get_num_block_value(BlockmodeloxstateField, record_id)    
    CurrentRockType = Myblkmdl_new.get_num_block_value(BlockmodelRockTypeField, record_id)

    if Currentzone == 10 and Currentoxstate == 1:
        Myblkmdl_new.set_block_value(orsource,record_id,11)
    elif Currentzone == 10 and Currentoxstate == 2:
        Myblkmdl_new.set_block_value(orsource,record_id,12)

I hope that helps and gets things working. Apologies if I've totally missed the point!

Rupert

Hi Ronald,

For the sort of control that you want, you need to use MMpy.File rather than MMpy.BlockModel. The script snippet below adds 3 new fields to a COLLAR file, and then populates those fields.

I have 3 screenshots that show the original file, the file after the fields were added and then finally after the fields are populated.

# Screenshot 1

Collar_File_Name = "Collar"
Collar_Path_Name = MMpy.Project.path() + Collar_File_Name + '.DAT'

Collar_file = MMpy.File()if Collar_file.open(Collar_Path_Name) == False:
    sys.exit("Can't open Collarfile")

# Add Azimuth, Dip and Tenement Fields

struct = Collar_file.structure

struct.add_field("AZIMUTH", MMpy.FieldType.real,0,2)
struct.add_field("DIP", MMpy.FieldType.real,0,2)
struct.add_field("TENEMENT", MMpy.FieldType.character,20,0)

Collar_file.structure = struct

# Screenshot 2

field_names = ["AZIMUTH", "DIP", "TENEMENT"]
field_id = []

for i in range(0,len(field_names)):
    field_id.append(Collar_file.get_field_id(field_names[i]))

# Populate Azimuth, Dip and Tenement Fields

for i in range(1, Collar_file.records_count+1):
    Collar_file.set_field_value(field_id[0], i, 45.0)
    Collar_file.set_field_value(field_id[1], i, -60.0)
    Collar_file.set_field_value(field_id[2], i, "ABC")

Collar_file.close()

# Screenshot 3

Thanks guys,

Paul, thanks for your help Monday. I got it all sorted Monday night and yesterday and ran the script and worked as it should - will need to throw in a check for the fields before trying to add them though - only throws up an error and then still runs but it is annoying all the same.