M500
M500
MineSight 500-Series Programs There are two primary MineSight programs used to calculate
drillhole composites:
The MineSight 500-series programs are used to create, update,
manipulate and display the drillhole composite File 8 (sorted M501V1 - composites within regular elevations
composites) and File 9 (unsorted composites).
M501SM - composites within variable elevations (seams)
A composite is the weighted average of a set of samples that fall for each drillhole
within a defined boundary. This boundary may be a bench or
a seam. The weighting factor is usually the sample length, but M501V1 uses the table of elevations stored in the PCF to
in some cases it may also include the sample specific gravity. determine which portions of each drillhole fall within each
vertical bench. The averages are computed using the actual
Composites are used in the interpolation of the deposit model length between the elevations of a bench as the usual weighting
instead of samples because they factor. This table of elevations is entered with and maintained
by M101V1. An option within M501V1 allows you to select a
Provide a mining basis for modeling separate item stored in File 11 (e.g., specific gravity) to be used
as an additional weighting factor. M501V1 requires that the
Reduce the amount of data used composite type BENCH is specified when initializing File 9.
Provide uniform support for geostatistics M501SM uses a seam code (integer values) to control the
compositing process rather than a table of elevations. All
The mining basis for an open pit operation is the bench height drillhole data with the same code is averaged into a single
which is defined by a set of elevations at the toe of each bench. composite using sample lengths as the weighting factor. An
The mining basis for an underground operation is the stoping optional geologic item may be used to select specified types
height or height to be drawn in relatively short term mining of material for the composite averages. M501SM requires that
operations. These heights are also defined by a set of elevations. the composite type SEAM is specified when initializing File 9.
The mining basis for a stratiform deposit is usually the height of
the seam which is defined by elevations at the top and bottom You also have the option to enter composites which have been
of the mining unit at each point in the deposit. calculated either manually or with another program. This data
is loaded into the composite Data Files with M500V1.
The general functions of the 500-series programs are:
When composites are calculated with M501V1 or M501SM,
Calculate composites from intervals an explicit linkage is created between File 12 and File 9. The
drillhole identifier and survey data are stored in File 12 but not
Load existing composites in File 9. Some 500-series programs refer to this data so you
should be careful that the relationship between Files 9 and 12 is
Edit and list composite values
not altered once M501V1 (or M501SM) has been run.
Interactive update of selected composite data
When composites are loaded with M500V1, a reference to File
Standard and user specified reports for file data 12 is not normally maintained. Therefore, you will not be able
to run some 500-series programs on the loaded File 9.
Preparation of composite maps in plan or section using
a plotter or a printer
M500V1 - Load Composites into MineSight Data Files This program accesses a user-coded subroutine for special
calculations and reports from the composite Data File. Calculated
This program is used to load composite or blasthole data into a values may be stored into the composite data file.
MineSight Data File. The composites may be computed outside
MineSight or created by program M507V1. M516V1 - Cross-Section of Composite Data (M122)
M501SM - Compute Seam Composites from Drillhole Data This program creates an M122 plot file for 2-D composite cross-
sections defined by two end points.
This program is used to compute weighted averages for grades
and quality parameters within predefined seams. M516V2 - Cross-Section of Composite Data (M122, without
dh traces)
M501V1 - Compute Bench Composites from Drillhole Data
This program creates an M122 plot file for a 2-D cross-section
This program is used to compute weighted average grades for displaying composite intervals without drillhole traces and
fixed vertical intervals (benches). labels.
M502V1 - Update Composites (Batch Run) M517V1 - Load VBM Codes to Assay or Composite Files
This program updates composites in a batch run by matching This program loads VBM feature codes to files 8, 9, or 11.
the elevation of the input with the composite file to find which
interval to update within a drillhole. M522V1 - Analyze Drillhole Spacing
M503V1 - List Composite Files This prgram accesses the composite data and computes average
drillhole spacing by bench based on a given search radius.
This program lists the composite data files. The data selection
and the order of data output from M503V1 may be specified. M523V1 - Decluster Composite Data
M504V1 - Bench Plan (M122) of Composite Data This program accesses the composite Data File and computes
the declustered composite values using the cell declustering
This program is used to prepare an M122 plot file from the technique.
composite data. Multiple items of data may be selected for
output and placed at any desired position around the location M524V1 - Point Validation for Interpolation Techniques
of each composite.
This program interpolates each composite from its neighbors
M505V1 - Addition of Codes to Composite Files by kriging and distance weighting, and then summarizes the
differences from the known values.
This program allows the addition of geologic or mining codes
to the composite Data File from ASCII input or an existing M525TS - Interactive Evaluation of Point Validation
assay file.
This program is used to evaluate point validation results from
M506V1 - Sort Composite Data File M524V1. The Program outputs a statistical analysis of the data
and histograms for specified items.
This program sorts the composite data from M501V1 and
M501SM by increasing northings (or eastings) and creates the M500-Technical Section
sorted composite data file, File 8.
Calculation of Composites
M507V1 - Extract Composite Data
Notes for M508V1 Driver Program
This program creates a sequential ASCII file from the composite
data file. Access to MSTorque Composite Set Data
Procedure Summary To Create File 9 interpolation as visually determined from the plan and
Composites cross-section maps?
1. List the PCF with M101V1 and check the toe elevations 10. Have the plan and cross-section maps from the M504
if calculating bench composites or check the seam names and M516 programs been compared with any manually
and codes if calculating seam composites. prepared maps?
2. Initialize File 9 using M102TS. Select the correct type of 11. Have engineers/geologists closely associated with the
compositeBENCH if M501V1 is to be used or SEAM project checked the composite data?
if M501SM is to be used.
12. Have engineers/geologists closely associated with the
3. Set up the appropriate run file for M501V1 or M501SM. project been consulted for their views and opinions
concerning:
4. Test the composite calculations by running the program
for just a few drillholes. any unreliable data or drillhole data not to be used?
5. Make a full production run with the program. any portions of the data that are different from the
majority?
6. Add additional data to the composites (i.e. geology
codes) with M505V1, M508V1 or M508RP. any differences in sample size, type or recovery?
7. List the composite file with M503V1. any special characteristics not indicated by grades?
5. Are there lengths within benches that are not equal to the
bench height other than at the top and bottom of the hole?
Suggestions For Bench Composite Items Other items of use to certain projects:
Each M501V1 drillhole composite to be stored in File 9 must SPGR - the specific gravity of rock
have the following items of information with the stated labels
after a regular M102TS run: REC% - the percent recovery of core
REF# - the drillhole reference number assigned by M201V1 TRUE - the true interval length corrected for dip
- Composite labels The composite data that is loaded should have the following
- Area selection variables present:
- Composites (optional)
DH-IDENT - drillhole identifier. Must be no more than 10
characters. System defined format is 10A1.
Input data fields that are blank will be read by the program as
0.0. Input data that is less than the minimum value for the item
will be stored by the program as missing values. Input data that
is greater than the maximum value for the item will be truncated
to the maximum value.
A single line of text that will be printed at the top of each page of printer output from the run.
NOTE: If IOP3=1, the specified boundaries can be in a separate file. Use one of the following command
lines:
CMD = INCL LOCPT FILE filename (to include points within the boundary) OR
CMD = EXCL LOCPT FILE filename (to exclude points within the boundary)
FMT1 = FORMAT FOR INPUT DATA (FOR FREE FORMAT, SET FMT1= (*))
NOTE: The first field for the input data must be DH-IDENT and must be specified as 10A1 in FMT1. If
IOP6 = 1, then the first field must be 10X.
NOTE: Optional items to override EAST and NORTH items for composites eastings and northings.
NOTE: Optional items to override default DHID1 and DHID2 labels. Used only if IOP13 = 1.
END
If the specified boundaries are to be added at the end of the run file rather than a separate file, then use the following format:
where 0 0 ends the input or input ends on end-of-file. Boundary points must close. The maximum number of boundary points
allowed is 2000.
corresponding to the format specified by the FMT1 line. Data with a completely blank DH-IDENT will be ignored. The M500V1
run ends on end-of-file.
USR = ABC
PUT09 = EAST NORTH ELEV. LNGTH T-CU MOLY EQCU / ITEMS TO STORE
END
will appear if you have not specified the standard coordinate items EAST, NORTH and ELEV. on the PUT09 or PUT08
line.
3. You should take care when loading File 9 and File 12 simultaneously. If you have already loaded File 12 (with
M201V1), you may either overwrite File 12 (if IOP2=1 and IOP4=0) or append to it (if IOP2=0 and IOP4=0).
4. If you are loading File 8, you should sort the input file on one of the coordinates (northing, easting or elevation) prior to
running M500V1. The file can be sorted with your computer systems sort routine.
PROGRAM SUMMARY The items to be composited are specified using ITMn commands
in the run file. The correspondence between the File 11 items
This program is used to compute weighted averages for grades (assays) and the File 9 items (composites) is also specified on
and quality parameters within predefined seams. the ITMn lines.
PROGRAM FLOW DIAGRAM Specific drillholes can be selected by using File 12 data unit
numbers. By default, all drillholes in File 12 will be composited.
RUN FILE: A maximum of 14,999 drillholes can be composited in a single
run. Multiple runs must be used when the drillhole database
- Quality & composite
labels exceeds this limit.
- Compositing parameters
- Seam table Any drillhole interval that has been deleted from File 11 (with
M202TS) will be ignored by M501SM.
- PCF - Composite quality M501SM establishes a linkage between File 9 and File 12. The
- Files 11 & 12 parameters over seams number of composite data units for a drillhole is stored with the
- Limit by litho code
collar data in File 12 (in items #CMPS and LOC09). There is
also a linkage between File 11 and File 12. The number of assay
data units for a drillhole are stored in File 12 (in items #INTS
and LOC11).
Note: The seam code on the waste interval must match that of
the coal. LITHO is a File 11 item.
Both Type 1 and Type 2 can exist in the same deposit. Use the
following commands:
IOP7 = 1
IOP8 = n
ITMn = WASTE PART PART
ITMn = LIMIT LITHO
A single line of text information that will be printed at the top of each page of printer output from the run.
NOTE: Use only with ITMn = LIMIT label-11 or ITMn = LITHO label-11.
The ITM commands below allow you to control the compositing calculations in several different ways. The items for SEAM,
LNGTH and ZTOP are required for all M501SM runs.
ITMn = SPLIT label-09 label-11 Optional split code; store split code in File 9
ITMn = SPLIT DUMMY label-11 Optional split code; do not store split code in File 9
ITMn = LNGTH label-09 label-11 Store dh length (LNGTH) or THICK vertical thickness
(THICK); required, label-11 is optional. If specified, it is
usually - AI -. If omitted, uses FROM and -TO- from File 11.
ITMn = ZTOP label-09 Store top elevation of seam; required.
Optional coordinates
ITMn = XBOT label-09 Store easting & northing of bottom of seam. Must specify
ITMn = YBOT label-09 both together.
ITMn = WASTE label-09 label-11 Store waste thickness; optional. Label-09 may be
or PART used alone.
ITMn = LIMIT label-11 Item from File 11 to limit intervals used to compute
or LITHO quality parameter averages; optional. Also use IOP7.
ITMn = READ SEAM FROM TO Read seam definitions from File 19, optional.
ITMn = #INTS label-09 Q.P. VALUE Store the number of Q.P./ore intervals within each seam;
optional.
END
Seam codes will be input at 10 per line in freefield format. Enter multiple lines, as necessary, to match the number of seams
specified by IOP6.
USR = ABC
COM
COM DH COMPOSITE FUNCTION
COM
COM* COMPS SEAM THICKNESS & QUALITY
I-O = 0
END
10 20 30
USR = ABC
I-O = 0
END
10 20 30
If your sample data resembles the example above, the following options are available to you.
HOLE# EAST NORTH ZTOP ZBOT BTU SEAM PART LNGTH FROM -TO-
77-012 17026.0 14999.0 925.6 905.5 12399 10. 0.00 20.1 134.0 154.1
77-012 17026.0 14999.0 881.8 875.0 10499 20. 0.00 6.8 177.8 184.6
77-012 17026.0 14999.0 856.6 851.2 9384 30. 0.00 5.4 203.0 208.4
To average the quality parameters for each seam and store the parting thickness
HOLE# EAST NORTH ZTOP ZBOT BTU SEAM PART LNGTH FROM -TO-
77-012 17026.0 14999.0 925.6 905.5 12399 10. 0.00 20.1 134.0 154.1
77-012 17026.0 14999.0 881.8 875.0 10499 20. 0.12 6.8 177.8 184.6
77-012 17026.0 14999.0 856.6 851.2 9384 30. 0.00 5.4 203.0 208.4
To average the quality parameters only for certain lithologies in each seamand treat the remainder of each seam as parting
HOLE# EAST NORTH ZTOP ZBOT BTU SEAM PART LNGTH FROM -TO-
77-012 17026.0 14999.0 925.6 905.5 12399 10. 0.00 20.1 134.0 154.1
77-012 17026.0 14999.0 881.8 875.0 10499 20. 0.00 6.8 177.8 184.6
77-012 17026.0 14999.0 856.6 851.2 9885 30. 0.40 5.4 203.0 208.4
To read the top and bottom of each seam from another file
COAL-501SM 19=CTSEAM.MOD
FILE CTSEAM.MOD
77-012 134.10 154.50 10
77-012 177.70 184.50 20
77-012 203.00 208.40 30
HOLE# EAST NORTH ZTOP ZBOT BTU SEAM PART LNGTH FROM -TO-
77-012 17026.0 14999.0 925.5 905.1 12399 10. 0.00 20.4 134.1 154.5
77-012 17026.0 14999.0 881.9 875.1 10499 20. 0.12 6.8 177.7 184.5
77-012 17026.0 14999.0 856.6 851.2 9384 30. 0.40 5.4 203.0 208.4
If your sample data resembles the example above, the following options are available to you.
To store partings inside each assay interval and store any portion of a seamnot assayed as parting
HOLE# EAST NORTH ZTOP ZBOT BTU SEAM PART LNGTH FROM -TO-
77-012 17026.0 14999.0 925.6 905.5 12399 10. 0.00 20.1 134.0 154.1
77-012 17026.0 14999.0 881.8 875.0 10499 20. 0.12 6.8 177.8 184.6
77-012 17026.0 14999.0 856.6 851.2 9885 30. 0.40 5.4 203.0 208.4
HOLE# EAST NORTH ZTOP ZBOT BTU SEAM PART LNGTH FROM -TO-
77-012 17026.0 14999.0 925.6 905.5 12399 10. 0.00 20.1 134.0 154.1
77-012 17026.0 14999.0 881.8 875.0 10499 20. 0.00 6.8 177.8 184.6
77-012 17026.0 14999.0 856.6 851.2 9885 30. 0.40 5.4 203.0 208.4
USR = ABC
COM
COM SEAM COMPOSITE FUNCTION
COM
I-O = 1
END
This error message will appear if there is some problem with File 12.
This error message will appear if there is some problem with File 11.
will appear if more than 1000 composites are assigned for one hole.
will appear if a composite control code was not specified (ITMn = SEAM label-09 label-11).
will appear if the line ITMn = ZTOP label-09 is missing from your run file.
9. If IOP7 and IOP8 are specified (lithology limit), then also use ITMn = LITHO label-11.
10. The MAXIMUM number of drillholes (14999) has been exceeded.All subsequent drillholes will be ignored.
Remaining drillholes may be composited in subsequent runs by setting IOP1 to 15000 (or 30000, 45000, etc.) and
setting IOP3=0.
FIXED LENGTH COMPOSITES NOTE: To add geologic codes to the intervals, include these two
lines in M505V1 run-file.
MineSight can calculate composites in several different ways.
You can composite to equal-height bench intervals, variable- PAR2 = ** / Downhole composite length
height bench intervals, between top and bottom of a coal seam, (must match PAR10 in M501V1)
between hangingwall and footwall of a vein, or by fixed length
from the collar of the drillhole. ITMn = DZBOT DEPTH
Program M501SM calculates seam composites. Vein composites / The keyword is DZBOT and the
are calculated usually by a special routine or by M501SM. label can be -to-, depth, from, etc.,
M501V1 calculates bench composites and fixed-length BUT IT MUST MATCH the label
composites. used in M501V1.
(for this example it is -TO-)
If you are calculating bench composites, horizontal drillholes
present a problem because they do not cross the bench toe or To store the coordinates of the composite endpoints, add the
crest. In this case, the drillhole is composited in fixed lengths following line to the run file:
from the collar. The length is set with PAR10.
CMD = COORD xtop ytop ztop xbot ybot zbot
To determine whether a drillhole is horizontal or vertical (for
compositing purposes), use the dip angle as the limiting angle Where:
(actually PAR4 is the sine of the dip angle). Any drillhole with a
xtop, etc are the appropriate File 9 item labels (in this specified
dip less than the limiting angle is treated as a horizontal drillhole.
order.) In this case the midpoint elevation is stored in the File 9
If you use fixed-length composites, make sure you use the proper item elev..
options in M505V1 and M516V1. For M501V1 and M620V1,
the programs will work without the ZMID options but may give
slightly different results. The options for M505V1 and M516V1 You can optionally merge small composite intervals to the
are needed for the programs to work. previous composite. The options include merging all small
composites to the previous composite, or merging only small
Please note that all the options for fixed length composites are
composites at the end of the drillhole to the previous composite.
in the menu system.
These options are specified using IOP14=1 or 2, together with a
NOTE: To compute just fixed length composites with M501V1. PAR12 threshold value.
Use these four lines in the run-file.
A single line of text that will be printed at the top of each page of printer output from the run.
NOTE: When rewriting File 9 (IOP3=1), recalculate all composites (IOP1; IOP2=0).
NOTE: If IOP4 = 1, the specified boundaries can be in a separate file. Use one of the following command
lines:
CMD = INCL LOCPT FILE filename (to include points within the boundary) -OR-
CMD = EXCL LOCPT FILE filename (to exclude points within the boundary)
NOTE: By default, DHID items are assumed to be: DHID1 and DHID2. To override the default
DHID labels, use ITMnn lines.
NOTE: PAR1 and PAR2 are used only if IOP5 is not equal to 0.
NOTE: If you chose a non-default limiting angle, then you must remember to include this value in your
M505V1 (PAR1) and M516V1 (PAR13) run files. See the Technical Section for more information.
NOTE: Use PAR5-8 only if IOP4 = 0. If PAR5-8 are not specified (i.e., PAR5 = PAR6 = 0.0 and/or PAR7
= PAR8 = 0.0), then values are retrieved from the PCF by default.
PAR11 = DEFAULT VALUE OF WEIGHTING FACTOR (used if weighting item is used, default=0).
NOTE: If PAR12 > 0, then the value of PAR12 will the the threshold used to determine whether to merge
a composite to a previous composite.
ITMn = THICK label-09 (of first grade) Store vertical composite length (or weight if WT501 command is used)
NOTE: For 3-D rotated models, vertical length is counted in the model Z direction.
ITMn = SLOPE label-09 (of first grade) Store slope composite length (or weight if WT501 command is used)
ITMn = SLOPE DUMMY (of first grade) Print but do not store composite length
NOTE: Only one (THICK label-09 , SLOPE label-09 OR SLOPE DUMMY) may be selected per run.
ITMn = DIST label-09 Store vertical (if THICK) or slope (if SLOPE) composite length of the
first grade item regardless of whether the WT501 command is specified.
NOTE: Normally the SLOPE or THICK command stores the length of the first grade item. However, a
WT501 command overrides this and instead stores the weight. The DIST command acts
independently of WT501 and always stores the length of the first grade item.
NOTE: Reference will be to the toe of the bench in which the drillhole terminates.
ITMn = GEOL label11 Geology code item to control compositing of fixed length
composites; a new composite is started each time the
geology item changes value.
Optional labels for DHID items (if different from DHID1 and DHID2):
END
If the specified boundaries are to be added at the end of the run file rather than a separate file, then use the following format:
where 0 0 ends the input or input ends on end-of-file. Boundary points must close. The maximum number of boundary points allowed
is 2000.
The number of values per line and their format are specified by the FMT1 line.
USR = ABC
IOP1 = 10 / M12 & N12 LIMITS
IOP3 = 1 / START LOCATION FOR COMPOSITES IN FILE 9
IOP4 = -1 / 0=COMPOSITE DH IN PCF LIMIT; -1=ALL DHS; 1=SPECIAL BOUND
IOP5 = 0 / 0=PCF TOES;N=READ TABLE;-1=CALC. A TABLE
IOP6 = 0 / 0=STD. OUTPUT; 1=PRINT DETAILS OF COMPOSITE CALCULATION
IOP8 = 1 / 1=DO NOT OUTPUT 0 LENGTH COMPS
IOP11 = 1 / 1=CALCULATE FIXED LENGTH
END
USR = ABC
USR = ABC
The PCF thinks your file 12 is empty. If this is not true, use M104TS to see how many data units are stored in file 12.
You may have to rebuild your file 12.
Theres some problem with your file 12 at survey 13. Use M202FS to view the offending drillhole.
4. ?M501V1: WARNING. DH SM-020 extends below PCF minimum. Ignoring rest of hole.
To composite entire hole, use IOP5 = -1 and PARs 1-3
Drillhole SM-202 will only be composited to the bottom of the PCF. If you want to composite the entire hole, you
should redefine the elevation limits of your program run by using PAR1 and PAR2.
Drillhole SM-202 will only be composited to the bottom of the PCF>. If you want to composite the entire hole, you
should redefine the elevatiopn limits of your program fun by using PAR1 and PAR2.
6. M501V1: WARNING. DH SM-020 extends outside the limits defined by PAR1 and PAR2.
Part of this hole will not be composited.
You have redefined the z limits of your program run using PAR1 and PAR2 but drillhole SM-020 still exceeds these
limits. The above message only occurs when I-O is greater than zero.
If IOP5 is greater than zero you need to specify a FORMAT for reading the array of toe elevations.
You have not specified a THICK or SLOPE length item in which to store the composite length.
10. You may check the internal calculations in M501V1 to verify results or to determine what assay intervals are being used
(IOP6=1). The following is a listing for vertical composites:
where
FROM & -TO- are the File 11 drillhole distances
ZTOP = elevation at FROM distance
ZBOT = elevation at -TO- distance
Whenever there are two lines printed for the same drillhole interval (i.e., 80.0 - 100.0), it means that this interval has been split
between two benches as shown under the bench # column.
The actual vertical interval which is used as the weighting factor is the difference ZTOP - ZBOT. Slope distance = TO - FROM.
where
FROM & -TO- are the File 11 drillhole distances
DIST1 = horizontal distance at FROM
DIST2 = horizontal distance at -TO-
INPUT: CALCS: For OCS projects, there probably is only one record per blasthole
so a very large tolerance can be used to allow the ELEV to get
- PCF - Select data and generate
- File 9 composite report
updated. For example, the bh data gets loaded without the
- File 12 coordinates as they are not yet available or the coordinates were
initially wrong. Later on, when the proper coordinates are
available the east, north, and elev can all be updated.
OUTPUT:
NOTE: File 29 will contain the names of any selected drillholes that are not found in File 9. The missing
drillhole names are always listed on the screen and in the report file.
A single line of text that will be printed at the top of each page of printer output from the run.
NOTE: Use IOP5=1 to update the value of the toe item. For example, in OCS projects, it is common to have
only one record per blasthole, and a very large tolerance can be used to allow the ELEV (toe item) to be
updated. If IOP5=1 is used, the toe item has to be first on the ITM1=VALUE label-09 list.
NOTE: All grade items must have real format. If FMT1 is not specified, free-field format is assumed.
END
4. INPUT DATA FILE (if IOP4=0)
A maximum of nine grade items is allowed. The grade items must be in real numbers.
USR = ABC
IOP1 = 100 / REF. #S FOR FILE 12 DHS TO BE LISTED
IOP4 = 0 / 0=DATA FOLLOWS END LINE; 1=DATA IN FILE 19
I-O = 2
END
OUTPUT:
A single line of text that will be printed at the top of each page of printer output from the run.
NOTE: If IOP3=1, the specified boundaries can be in a separate file. Use one of the following
command lines:
CMD = INCL LOCPT FILE filename (to include points within the boundary) or
CMD = EXCL LOCPT FILE filename (to exclude points within the boundary)
NOTE: Do not use IOP4 = 0 or 2 if composites were loaded with M500V1 and File 12 was not
created.
NOTE:
(1) By default, the first coordinate triplet is in items EAST, NORTH and ELEV. (not the dot), the second in XTOP,
YTOP, and ZTOP, and the third in XBOT, YBOT and ZBOT. To override the items use CMD statements
(below).
(2) For non-rotated models IOP9 is ignored. For horizontally rotated models only X and Y coordinates are used to
convert to model coordinates (the elevation in model and project coordinates is the same in this case).
GETn = LABELS OF ITEMS TO BE PRINTED (e.g., GET08 = EAST NORTH ELEV. T-CU MOLY)
NOTE: If you are using a FMT1 statement and IOP4=0 or 2, you must specify a field for DH-IDENT and
unit number. If IOP4=1, you must specify a field for REF# and unit number.
Optional labels for DHID items (if different from DHID1 and DHID2) (if IOP13 = 1):
END
If the specified boundaries are to be added at the end of the run file rather than a separate file, use the following format:
where 0 0 ends the input or input ends on end-of-file. Boundary points must close. The maximum number of boundary points
allowed is 2000.
DH-IDENT
USR = ABC
END
USR = ABC
END
USR = ABC
GET09 = EAST NORTH ELEV. LNGTH T-CU MOLY EQCU DHTOE ROCK
END
2. If you limit the composite data by ITMn = label-09 RANGE, OMIT or TABLE, the label-09 (or 08) must be specified
on the GET line.
OUTPUT:
A single line of text that will be printed at the top of each page of printer output from the run.
IOP1 = M08 FIRST SORTED COMPOSITE DATA UNIT TO ACCESS (if IOP6=1) OR
= M09 FIRST COMPOSITE DATA UNIT TO ACCESS (if IOP6=0) OR
= M12 FIRST SURVEY DATA UNIT TO ACCESS (if IOP6=1)
(default = 1)
IOP2 = N08 FIRST SORTED COMPOSITE DATA UNIT TO ACCESS (if IOP6=1) OR
= N09 LAST COMPOSITE DATA UNIT TO ACCESS (if IOP6=0) OR
= N12 LAST SURVEY DATA UNIT TO ACCESS (if IOP6=1)
(default = last data unit stored)
NOTE: If IOP3=1, the specified boundaries can be in a separate file. Use one of the following command
lines:
CMD = INCL LOCPT FILE filename (to include points within the boundary) -OR-
CMD = EXCL LOCPT FILE filename (to exclude points within the boundary)
NOTE: Do not use IOP6=1 if composites were stored with M500V1 and File 12 was not created.
NOTE: To use IOP7 = 1 for an inclined GSM project, you must specify a composite item for elevation
(e.g., ZTOP). Use CMD = ZCRD label-nn.
NOTE: If IOP6 = 1, drill-hole IDs will be used from file 12 and IOP8 will be ignored.
IOP16 = -1 SELECT PEN# BASED ON FIRST VALUE, BUT DO NOT PLOT FIRST VALUE
= 0 PLOT GRADES
= 1 DO NOT PLOT GRADES. USE GRADES TO SELECT PEN# FOR SYMBOL.
PAR5 = n PLOT ALL COMPOSITES FROM ELEVATION (PAR1) TO n UNITS ABOVE PAR1
Value Color
0.0 x vmax1 p1
vmax1 x vmax2 p2
vmax2 x vmax3 p3
vamx3 x vmax4 p4
x max4 p1 - Set vmax4 to 999 (or higher than any grade) to avoid reverting to pen 1.
NOTE: P1, p2, p3 and p4 are optional pen numbers. Default is 1, 2, 3, 4. Cutoffs are applied to ITMn item
and should be stated in the same precision as the data.
NOTE: The symbol numbers depend on your plotter, refer to the M122V1 documentation for symbol n
umbers and valid symbol sizes.
NOTE: Always list one more pen number, symbol number and symbol size than you have cutoffs. The
additional pen is for any material greater than the last cutoff.
Items to plot
NOTE: Use label-08 for the above ITM commands when IOP6 = 1.
To override default item labels for DHID items (if IOP8 = 1), use commands:
END
If the specified boundaries are to be added at the end of the run file rather than a separate file, then use the following format:
where 0 0 ends the input or input ends on end-of-file. Boundary points must close. The maximum number of boundary points
allowed is 2000.
filenametoe_elev.
where
filename= name of output M122 file
toe_elev.= elevation of composite to be selected or seam number
USR = ABC
IOP3 = 0 / USE COMPS WITHIN PCF BOUNDARY
IOP4 = 0 / OUTPUT M122 PLOT COMMANDS
END
USR = ABC
END
BORDR 1 10000 15000 10000 15000 0
BORDR 1 10015 14985 10000 14985 0
GRID 1 10000 15000 10000 15000 1000
TICKS 1 10000 10000 10000 15000 1000 .08
TICKS 1 10000 15000 15000 15000 1000 .08
COORD 1 9700 10025 9700 15025 1000 .05 0 0 10000 1000
COORD 1 9975 15025 14975 15025 1000 .05 90 0 10000 1000
BORDR 1 10000 15000 9500 10000 0BORDR 1 10015 14985 9515 9985 0
LINE 1 13500 9500 13500 10000 0
LINE 1 13500 9800 15000 9800 0
LINE 1 15000 9600 13500 9600 0
LABEL 1 13910 9880 0 .04 17/WIDGET MINING CO.
LABEL 1 13850 9680 0 .04 20/BENCH 22 T-CU VALUES
LABEL 1 13910 9540 0 .04 17/SCALE: 1" = 1000
USERF PLT504.PAA
ENDOC
USR = ABC
IOP3 = 1 / -1=USE ALL COMPS; 0=USE PCF LIMITS; 1=USE SPECIFIED BOUNDARY FILE
IOP6 = 0 / 0=FILE #9 ONLY; 1=FILES #9 & #12
IOP11 = 1 / SYMBOL NUMBER FOR COMPOSITE LOCATION (default = 3)
PAR5 = 7.99
PLT = 0.02 0.05 .1 99999 1 2 3 4 / CUTOFFS FOR PEN CHANGE
USR = ABC
IOP3 = -1 / USE ALL COMPOSITES
PAR1 = 10 / SEAM 1 ONLY
PAR2 = .1 1000 / SIZE OF SYMBOL, SCALE
ITM1 = SULF F5.2 -.2 -.07 .1 / PLOT SULFUR
ITM2 = D.H.# A10 -.3 .07 .1 / PLOT DH-IDENT
END
USR = ABC
END
BORDR 1 4000 20000 6000 28000
LINE 1 4000 10000 20000 10000
LINE 1 20000 15000 4000 15000
LINE 1 4000 20000 20000 20000
LINE 1 20000 25000 4000 25000
LINE 1 5000 28000 5000 6000
LINE 1 10000 6000 10000 28000
LINE 1 15000 28000 15000 6000
BORDR 1 3000 21000 5000 29000 0
COORD 1 3400 9950 3400 26000 5000 .1 0 -1 10000 5000
COORD 1 5050 28100 20050 28100 5000 .1 90 -1 5000 5000
USERF PLT504.PAC
ENDOC
2. A common error is to incorrectly specify the output format for label plotting. If you do not specify a large enough format
field, asterisks will be plotted where a number greater than the format is plotted, i.e. plotting a number 1000.0 with a
format of F3.1.
3. Another error is to incorrectly specify the location of the plot label with the xoffset and yoffset. Run a test plot on a
small portion of the data to check the offset values.
PROGRAM FLOW DIAGRAM access Files 11 and 12, pick up the existing codes and
determine a code for the composite using a majority
rule, or
RUN FILE:
input a set of data containing drillhole lengths with
- Item to update appropriate codes and determine the majority code for
- Dh selection
each composite of a drillhole.
- PCF - Add geology or mining Drillholes can be selected using the File 12 data unit numbers.
- File 9 (11, 12 opt.) codes to the Composite
- Optional ASCII input File
The item to be updated is specified on an ITMn line. If File 11
file
is used as the source of the geology codes, then a corresponding
item from File 11 should be specified.
After you run M505V1, list the composite file with M503V1 to
document and check the codes inserted with this program.
A single line of text that will be printed at the top of each page of printer output from the run.
NOTE: 1) Must match value used in M501V1 for PAR4. See the 500 Technical Section.
2) Enter 1.0 in a case of length compositing
Optional item to store minor geologic code within each composite interval (i.e., store two codes for each
interval).
Optional item to store the percentage of the composite interval covered by the major geologic code.
NOTE: DZBOT must be used for 3-D (non-vertical) composites. The label-09 is the same item used as
DEPTH in M501V1. For fixed length composites (horizontal or near-vertical holes) the depth item is
usually -TO-.
NOTE: DZBOT can be used alone with sequential 3-D composites. There is no adequate way to back
calculate drillhole depths from seam composites, hence the drillhole depths must be available from
M501SM in order to run M505V1.
END
End input with a blank line. The Geologic Data Lines must match FMT1.
USR = ABC
I-O = 1
END
SM-008 0.0 1
SM-008 80.0 2
SM-008 620.0 3
SM-008 1460.
SM-009 0.0 1
SM-009 81.0 3
SM-009 1260. 4
SM-009 1600.
USR = ABC
COM
COM DH GEOLOGY FUNCTION
COM
COM -
COM = (DH GEOL INPUT DATA) FMT1 = (10A1, F7.2, 14X, I3)
COM -
END
77-001 0.00 4.10 4.10 1 0
77-001 4.10 92.10 88.00 2 0
77-001 92.10 99.30 7.20 3 0
USR = ABC
COM
COM DH GEOL LOAD FUNCTION
COM
COM
I-O = 2
END
SM-001 0.0 1
SM-001 62.5 2
SM-001 512.5 2
SM-002 0.0 1
SM-002 60.0 2
SM-002 468.8 2
SM-003 0.0 1
SM-003 75.1 2
SM-003 625.1 2
SM-020 0.0 1
SM-020 50.0 3
SM-020 1660.6 3
This error message will appear if a DH-IDENT specified in the input ASCII file does not exist in File 12.
You can only run M505V1 for 3-D block model composites unless the FROM and -TO- depths are stored for each
composite.
A single line of text that will be printed at the top of each page of printer output from the run.
PAR1 = MINIMUM VALUE OF SORT ITEM (default = minimum value for item)
PAR2 = MAXIMUM VALUE OF SORT ITEM (default = YMAX)
NOTE: This default is only useful when sorting on Northing.Set PAR3 so that (PAR2 PAR1) / PAR3
(maximum number of sort subintervals)
END
USR = ABC
END
USR = ABC
IOP1 = 10 / M09N09; FIRST, LAST COMPS TO SORT
ITM2 = ROCKTABLE2122
IOP21 = 24 / TABLE OF ROCK CODES
END
The value you specified in IOP4 is larger than the number of composites currently stored in FILE8. The new composites
will be appended to FILE 8 (equivalent to specifying IOP4=0).
Initially the composites are coarsely sorted into increments of length PAR3 (see requiremetns and considerations). If
PAR3 is not sufficiently small, the number of composites per increment is exceeded. Attempt to lessen the number of
composites per increment by reducing PAR3. Note that PAR3 must be greater than or equal to (PAR2PAR1)/
mamimum number of sort subintervals.
4. ?M506V1 . . . ERROR. The number of words required to store data for 1 composite is xxx.
This value is stressing the limits of the program. Reduce the value of PAR3 and try again.
You will sometimes get this message if your File 8 requires more than 10 words to store data for one composite even
though the number of composites per sort increment is less than the maximum. Determine the number of words needed
to store one data unit (M105TS). If, for example, you need 20 words, then divide your PAR3 in half. If you need 30
words, divide by three.
The next two data fields output normally contain the EAST and
INPUT: CALCS: NORTH coordinates of the composite in whichever order you
- PCF - Select data specify on the GET line. However, if this data will be used in
- File 9 - Create ASCII file the interpolation with M620V1, EAST must precede NORTH.
- File 12 Output fields must be in EAST, NORTH, ELEV. order for use
in M236V2.
NOTE: Composites must have been loaded with the M501 series
of programs. Composites loaded with M500V1 cannot be
extracted.
OUTPUT:
A single line of text that will be printed at the top of each page of printer output from the run.
NOTE: If IOP3=1, the specified boundaries can be in a separate file. Use one of the following command lines:
CMD = INCL LOCPT FILE filename (to include points within the boundary) OR
CMD = EXCL LOCPT FILE filename (to exclude points within the boundary)
1) Using IOP5=2, for rotated models ROW# is meaningless. An artificial grid is used to determine ROW#.
It is not used by M236V1.
2) If IOP5=3, the first three fields are REF#, DH-IDENT, and UNIT#.
3) Using IOP5=4 or IOP5=5, and generating a file for M500V1, then specify the first items as EAST NORTH ELEV.
(3DBM) or EAST NORTH ZTOP ZBOT (GSM) or if generating a file for M502V1, then specify the first item as
ELEV. (3DBM) or SEAM/ZONE/etc (GSM)
NOTES:
1) By default, the first coordinate triplet is in items EAST, NORTH and ELEV. (note the dot), the second in
XTOP, YTOP and ZTOP, and the third in XBOT, YBOT and ZBOT. To override the items use the CMD
statements below.
2) For non-rotated models IOP8 is ignored. For horizontally rotated models only X and Y coordinates are
used to convert to model coordinates (the elevation in model and project coordinates is the same in this
case).
IOP5=0, the order of the items in the output is: REF#, the DH-IDENT and every item on the GET line
Example: FMT1 = (I5, 1X, 10A1, nFw.d)
IOP5=1, the order of the items in the output is: REF# DH-IDENT XCOL YCOL and every item on the GET line
Example: FMT1 = (I5, 1X, 10A1, 2F12.3, nFw.d)
IOP5=2, the order of the items in the output is: REF# , every item on the GET line, ROW# DH-IDENT
Example: FMT1 = (I5, 1X, nFw.d, F7.0, 10A1)
IOP5=3, the order of the items in the output is: REF# DH-IDENT UNIT# and every item on the GET line
Example: FMT1 = (I5,1X,10A1,I7, nFw.d)
IOP5=4, the order of the items in the output is DH-IDENT and every item on the GET line
Example: (10A1,1X, 3F12.3, nFw.d)
IOP5=5, the order of the items in the output is REF# and every item on the GET line
Example: (I5,1X,3F12.3,nFw.d)
NOTE: If IOP5=2, use RANGE to name the bench elevation on which polygons are to be generated, e.g.,
ITM01 = ELEV. RANGE 4300 4300.
END
If the specified boundaries are to be added at the end of the run file rather than a separate file, then use the following format:
where 0 0 ends the input or input ends on end-of-file. Boundary points must close. The maximum number of boundary points
allowed is 2000. Boundary point lines must preceed any drillhole request lines.
DH-IDENT
where DH-IDENT is the drillhole label, one per line. Drillhole request lines must follow any boundary point lines.
USR = ABC
IOP1 = 10 / FILE 12 LIMITS
IOP3 = 0 / DEFAULT TO PCF GEOMETRIC LIMITS
IOP4 = 19 / OUTPUT DATA TO FILE 19
END
NOTE: To access the MSTorque database, use 9=TORQ and include an MSTorque info file (20=mstorq.inf). For information
about the MSTorque info file, please refer to the M500-Technical Section, Access to MSTorque Composite Set Data.
A single line of text that will be printed at the top of each page of printer output from the run.
NOTE: If File 12 is specified on the NAMES LINE, then IOP1 & IOP2 refer to File 12. If IOP5 0 or File 8/9 is
specified on the NAMES LINE, then IOP1 and IOP2 refer to that file.
NOTE: If IOP3=1, the specified boundaries can be in a separate file. Use one of the following command lines:
CMD = INCL LOCPT FILE filename (to include points within the boundary) OR
CMD = EXCL LOCPT FILE filename (to exclude points within the boundary)
IOP6 = 0 OMIT STORING DATA FROM THE USER ROUTINE INTO FILE 9/8 OR MSTORQUE
= 1 STORE DATA FROM THE USER ROUTINE INTO FILE 9/8 OR MSTORQUE
NOTE: You cannot write to both File 12 and 9/8. If both IOP5 and IOP6 are 0 and there are UPD/PUT items,
the items will be calculated but not stored. If File 8 is selected, IOP5 is automatically set to -1.
PUT09 (or 8 or 12)= LABELS OF ITEMS TO BE STORED AFTER M508RP CALCULATIONS (if IOP6=1)
NOTE: If you use File 12 Range Checking you MUST specify the File 12 item with a GET12 command. You
cannot range check on File 9(8) items in the same run. GET12, PUT12, OUT12 and label-12 are not allowed if
9=TORQ.
NOTE: GET/UPD items that are undefined values will be substituted by the UNDEF value before the
RPN expressions are evaluated. UPD/PUT/OUT items that are undefined will be substituted after the RPN
expressions have been evaluated. Undefined values that occur during RPN evaluation cause the result to be
undefined. CMD ALIAS IS is disabled for the TORQ option; there is no need for label substitution because
there are no restrictions for TORQUE label names.
END
If the specified boundaries are to be added at the end of the run file rather than a separate file, then use the following format:
where 0 0 ends the input or input ends on end-of-file. Boundary points MUST close. The maximum number of boundary points
allowed is 2000.
TEMP = TOTCU
You do not need to enclose attributes with blanks and other special characters in quotes, however you may do so for better
readability. You may also use brackets. All three expressions below are valid:
USR = ABC
IOP1 = 0 0 / FIRST AND LAST DATA UNITS FROM FILE 12
END
EQCU = TOTCU + MOLY*5.0
The MEDITM information (I-O=2) should allow you to check the items of information that will be retrieved (GET
items) from File 9/8 or 12. CHECK IF ALL DESIRED ITEMS ARE LISTED.
2. If you make a mistake in entering the variable names in the equations specified in the run file, the calculations will be in
error: IF RESULTS ARE NOT CORRECT, PLEASE CHECK THE EQUATIONS SPECIFIED AND LABELS IN THE
RUN FILE.
3. If there is no PUT item, nothing will be stored in File 9/8 or 12. Set IOP6=1 and make sure you include a PUTnn
statement. If IOP6=1 and no PUTnn was included, the error message
will appear.
4. If M508RP values are stored in File 9 that are either above the PCF maximum or below the minimum, the error message
will appear.
RUN FILE: The user subroutine is called by the main routine (M508V1) a
number of times. A call is always made for initialization. The
- DH selection value of the argument NTRY, determines what action should
- Item selection
be taken in the subroutine:
NTRY=1 - Initialization
Procedure Summary for Use of USR508 IMPORTANT NOTE: If the M508V1 run replaces existing
values in the composite file, there is no way to restore the original
1. Define the desired calculation and draw a flowchart. file. Save a copy of the file (on disk or tape) before you begin
your M508V1 run.
2. Code the subroutine; see the Technical Section
explaining the linkage between M508V1 and USR508. Please refer to the Technical 600 section for more information
Also see the sample subroutine in this section. about how to write MINESIGHT user subroutines.
5. Check output for all drillholes closely with list and maps.
A single line of text information that will be printed at the top of each page of printer output from the run.
IOP2 = N12 LAST SURVEY DATA UNIT TO ACCESS (default = last survey stored)
NOTE: If IOP3=1, the specified boundaries can be in a separate file. Use one of the following command lines:
CMD = INCL LOCPT FILE filename (to include points within the boundary) OR
CMD = EXCL LOCPT FILE filename (to exclude points within the boundary)
IOP4 = 0 CALL THE USER ROUTINE FOR EACH DRILLHOLE IN THE FOLLOWING ORDER:
NOTE: NTRY refers to an argument that is passed to the USR508 subroutine from the M508V1 program (see
example subroutine listing that follows). The USR508 routine must test the NTRY argument and branch to the
appropriate section.
NOTE: If IOP5 = -1, File 12 is still required for the run to determine the location of the composite intervals.
Array B12 will contain valid data, but no calls for NTRY=3 or NTRY=7 are made.
NOTE: If IOP6 = -1, no call for NTRY=2 or NTRY=6 are made. If IOP6=1, specify PUT09 line.
PUT09 (or 12) = LABELS OF ITEMS TO BE STORED IN FILE 9 (If IOP6=1 or IOP5=1)
NOTE: You must specify label-09 on the GET09 line, or label-12 on the GET12 line. If you use File 12 Range
Checking you MUST specify the File 12 item with a GET12 command. You cannot range on File 11 and File 12
items in the same run.
END
If the specified boundaries are to be added at the end of the run file rather than a separate file, then use the following format:
where 0 0 ends the input or input ends on end-of-file. Boundary points must close. The maximum number of boundary points
allowed is 2000.
USR = ABC
END
C
C
C NTRY = 2 : PROCESS COMPOSITE INTERVAL DATA
C
500 CONTINUE
CU = RDAT(K09GET)
IF (CU .LE. 0) CU = 0.0
AM = RDAT(K09GET+1)
IF (AM .LE. 0) AM = 0.0
ECU = CU + FACTR * AM
RDAT(K09PAK) = ECU
C
IF (IBUG .GT. 0) WRITE(LPR,620)IDHN, CU, AM, ECU
620 FORMAT( DH REF# ,I6, CU =,F8.2, MO =,F8.3, EQCU =, F8.3)
RETURN
C
C
C NTRY = 3 : PROCESS END OF SURVEY
C
1000 CONTINUE
RETURN
C
C
C NTRY = 4 : PROCESS OR SUMMARIZE AT END OF DRILLHOLE
C
1500 CONTINUE
NDHS = NDHS + 1
RETURN
C
C
C NTRY = 6 : 2ND PASS FOR DH INTERVAL DATA
C
2500 CONTINUE
RETURN
C
C
C NTRY = 7 : 2ND PASS FOR END OF SURVEY
C
3000 CONTINUE
RETURN
C
C
C NTRY = 8 : 2ND PASS FOR END OF DH
C
3500 CONTINUE
RETURN
C
C
C NTRY = 5 : PROCESS OR SUMMARIZE AT END OF ENTIRE RUN
C
5000 CONTINUE
WRITE(LPR,5010) NDHS
5010 FORMAT(// END RUN # DHS = ,I3/)
RETURN
END
You have selected IOP5=1 but not specified any PUT12 items.
You have selected IOP6=1 but not specified any PUT09 items.
3. If IO=2 is specified in the run file, MineSight routines will print a list of all the items in File 9 followed by the items
to be used, in the form
The MEDITM information (IO = 2) should allow you to check the items of information that will be retrieved (GET
items) from File 9. Check if all desired items are listed.
4. If M508V1 values are stored in File 9 that are either above the PCF maximum or below the minimum, the error message
will appear.
5. A test run should always be performed on a small portion of the drillhole data before actually storing any calculations.
Include an option in your USR508 logic to list enough values to be able to check the calculations.
6. If you want an input or output sequential file, you must specify a 19=filename on the NAMES LINE. IOP7 must be set
appropriately.
The M122V1 run file must specify an area that has X limits
from 0.0 to the section length in user units.
A single line of text that will be printed at the top of each page of printer output from the run.
IOP6 = STANDARD CALCOMP SYMBOL NUMBER (0-13) USED FOR DH COLLARS (default=6)
NOTE: To exclude selected drillholes from the section, set IOP9-1 and put negative drillhole REF#s
on the DH selection lines.
NOTE: Used if PLT is specified or following CMD option.Shade types (0-7) are as follows:
NOTE: Used if PLT is specified or following CMD option.Box types (-1 through 5) are as follows:
-1 =no box
0 = box on right/left side of DH line (See IOP18)
1 = ticks on right side of DH line
2 = ticks on left & right side of DH line
3 = ticks on left side of DH line
4 = box on left & right side of DH line
5 = box on left side of DH line
NOTE: A positive number places the ticks on the right. A negative number places the ticks on the left.
PAR20 and PAR21 define the size and location of the ticks.
IOP29 = PEN# FOR THE GRID LINES EASTING =CONST ON THE COLLAR PLOT (default =1)
IOP30 = PEN# FOR THE GRID LINES NORTHING = CONST ON THE COLLAR PLOT (default =1)
IOP31-34=SHADE TYPES FOR COMPOSITE INTERVAL (or use following CMD option)
NOTE: Cutoffs for shading correspond to PLT commands. You must set IOP13=1.
IOP35-38=BOX TYPES FOR COMPOSITE INTERVAL (or use following CMD option)
NOTE: Cutoffs for box type correspond to PLT command. You must set IOP14=1.
PAR1 = LENGTH OF LINE ON EACH SIDE OF DH PLOT INDICATING INTERSECTION OF DH WITH THE
SECTION (default=1.0")
NOTE: The resultant plot will have the beginning of the section shifted and the end of the section
truncated by the distance specified in PAR2.
PAR11 = SIZE OF COLLAR PLOT SYMBOL (in inches) AND DH ANNOTATION (default=0.1)
PAR12 = HEIGHT (in inches) OF REFERENCE LINE PLOT FOR COLLAR PLOT (default = 5.0")
NOTE: If the reference line plot is not large enough to properly plot a drillhole location, the location
will be plotted on the border. To omit plot, set PAR12=-1.
NOTE: In case of length compositing use command line ITMn = DHTOE. (see below)
PAR15 = OFFSET OF COLLAR PLAN FROM SECTION (in plotter units;default = 0.5)
PAR17 = STARTING GRID EASTING FOR COLLAR GRID LABELING (default = PCF21)
PAR18 = STARTING GRID NORTHING FOR COLLAR GRID LABELING (default = PCF22)
NOTE: PAR17 and PAR18 are for reference. The grid will go both directions from PAR17/PAR18.
where
NOTE: P1, p2, p3 and p4 are optional pen numbers. default is 1, 2, 3, 4. If shading with PEN# 1, use 1 1
1 1. Cutoffs are applied to ITMn item and should be stated in the same precision as the data.
NOTE: Fw.d is the format of the plotted grade label. There must be at least one ITM input.
NOTE: Always list one more pen number than you have cutoffs. The additional pen, box and shade are
for any material greater than the last cutoff.
NOTE: IOP13 and IOP14 no longer need to be set to get the box and shade patterns from the above
commands. The existing PLT command and IOPs 13, 14, and 31-38 work exactly as before. However, if
both the PLT and/or IOP31-38 and the CMD commands are present in a run file, the CMD commands
will override the PLT and IOP31-38 commands.
NOTE: Item labels used in data selection must be specified in the plot ITMn lines above.
NOTE: Label-09 is the depth down the drillhole to the toe of the composite as computed in M501V1
using the command
ITMn = DEPTH label-09 * This option must also be used if you set IOP5= -1 in the
M501V1 run to create a table of bench toes.
To use this option, composites must be continuous. (In M501V1, IOP8 must be set to zero.) This item
must be specified if, for M501V1, the fixed length composites option was used (IOP11=1 in M501V1
run file). Bench composites can also be plotted using this command.
END
X1 Y1 X2 Y2 search-dist. filename
where
where
USR = ABC
END
USR = ABC
END
USERF PLT516.PAA
ENDOC
USR = ABC
2. ?M516V1: ERROR. You MUST enter both X and Z codes in the run file.
3. ?M516V1: ERROR. Exceeded the maximum number of DHs for inclusion: 5000
You can add at most 5000 drillholes using the Optional Drillhole Selection Lines.
Drillhole XXXXX, which was specified using an Optional Drillhole Selection Line, was not found. Check the typing of
the drillhole.
You must enter at least one VALUE item with the command ITM1=label-09 VALUE Fw.d
Therefore, the M122V1 run file must specify an area that has X
limits from 0.0 to the section length in user units. If additional
data with actual section coordinates is to be included on the
same plot, the M122 SHIFT command must be used.
A single line of text that will be printed at the top of each page of printer output from the run.
NOTE: Used if PLT is specified or following CMD option. Shade types (0-7) are as follows:
NOTE: Used if PLT is specified or following CMD option. Box types (-1 through 5) are as follows:
-1 = no box
0 = box on right side of DH line
1 = ticks on right side of DH line
2 = ticks on left & right side of DH line
3 = ticks on left side of DH line
4 = box on left & right side of DH line
5 = box on left side of DH line
The default labels for DHID items are lines DHID1 and DHID2. To override default, use
ITMn lines.
IOP31-34 = SHADE TYPES FOR COMPOSITE INTERVAL (or use the following CMD option)
IOP35-38 = BOX TYPES FOR COMPOSITE INTERVAL (or use the following CMD option)
Cutoffs for box type correspond to PLT command. Must set IOP14=1.
PAR11 = SIZE OF COLLAR PLOT SYMBOL AND ANNOTATION (in inches) (default=0.1)
NOTE: P1, p2, p3 and p4 are optional pen numbers. Default is 1, 2, 3, 4. If shading with PEN# 1, use 1 1
1 1. Cutoffs are applied to ITMn item.
NOTE: There must be a FROM ITM with XYZ coordinate labels for top of interval, a TO ITM with XYZ
coordinate labels for bottom of interval, and at least one VALUE ITM input, any others are optional.
Fw.d is the format of the plotted grade label. If the XYZ labels for FROM and TO items are the same
composites will be plotted as points with grade values plotted at 90.
NOTE: IOP13 and IOP14 no longer need to be set to get the box and shade patterns from the above
commands. The existing PLT command and IOPs 13, 14, and 31-38 work exactly as before. However, if
both the PLT and/or IOP31-38 and the CMD commands are present in a run file, the CMD commands
will override the PLT and IOP31-38 commands.
NOTE: Item labels used in data selection must be specified in the plot ITMn lines above.
Optional labels for DHID items (if different from DHID1 and DHID2):
END
Revised: 19-September-01 Page 516-15
MineSight
M516V2 - CROSS-SECTION OF COMPOSITE DATA (M122) Proprietary Software: Mintec, inc.
X1 Y1 X2 Y2 search-distance filename
where
USR = ABC
END
10500 12750 14500 12750 100 PLT.PAA
USR = ABC
IOP4 = 0 / 0=E-W SECTION, 1=N-S SECTION
IOP7 = 0 / 0=PLOT REF LINE, 1=NO REF LINE
IOP8 = 0 / 0=PLOT TICKS ON BORDER, 1=NO TICKS
IOP18 = 1 / 0=PLOT ON RIGHT SIDE OF LINE, 1=PLOT ON LEFT SIDE OF LINE
END
10500 13200 14500 13200 100 PLT516.PAB
USR = ABC
END
10500 13750 14500 13750 100 PLT516.DAT
will occur if the keyword VALUE does not appear on an ITM line.
PROGRAM FLOW DIAGRAM The program can use VBMs with any orientation (including
non-orthogonal).
RUN FILE:
- VBM planes selection
The program can automatically read VBM orientation from the
- Feature codes/equivalent VBM file (if you have a new VBM file). You can also specify
codes selection VBM orientation with IOP5.
You may use a constant range of influence for each VBM plane
(given by PAR11) setting IOP30 = -1. However, if you are
working with a set of parallel planes the program will
automatically take a volume of influence for each plane half-
way between planes. You also have an option (IOP30=1) to use
INPUT: CALCS: the space between two planes (i.e from one bench toe to the
next). PAR11 in these cases may be used as the maximum
- PCF - Calculate intersections of
- VBM File 25
allowed range of influence.
drillhole projections with
- File 8 or VBM features
- File 9 and optional - Load codes to File 8/9/11 Feature codes and equivalent (geologic) codes are generally
File 12 or entered after the END line in a run-file. However, you may
- Files 11 and 12 choose to consider all 3-digit left and right codes (IOP11=1).
By default (IOP11=0 and IOP12=0) codes are equivalenced to
3-digit right codes. If you choose IOP11=1 equivalent codes
should be entered in a run file. You may equivalence several
OUTPUT: feature codes to the same geologic code.
- Updated File 8 or 9 or 11
You also have an option (IOP12=1) to consider 6-digit codes as
- Printer Output file
codes on the right of the feature. In this case you must enter
equivalent codes in a run file.
A single line of text that will be printed at the top of the printer output from the run.
NOTE: If file 12 is used then IOP1 and IOP2 are the first and last surveys to use. If file 8 is used, or file 9 is
used directly, then IOP1 and IOP2 are the first and the last composites to access.
IOP3 = -1 NO BOUNDARY
= 0 LIMIT DHS TO PCF LIMITS
= 1 READ BOUNDARY FROM THE BOUNDARY FILE
NOTE: Used only if file 12 is used. Checks collars only or checks at beginnings of all surveys, depending on
IOP7. If IOP3=1 boundary file should be entered with one of the following commands:
CMD = INCL LOCPT FILE Filename (to include points within the boundary)
or
CMD = EXCL LOCPT FILE Filename (to exclude points within the boundary).
NOTE: IOP5 must be > 0 if the VBM has no default orientation. The entered value of IOP5>0 will override the
default orientation.
IOP11 = -1 USE ALL FEATURES (do not enter any codes in the run file)
= 0 DO NOT READ EQUIVALENT CODES (use right code if 6-digit codes are entered)
= 1 READ EQUIVALENT CODES FROM THE RUN FILE
IOP12 = 0 EQUIVALENCE 3-DIGIT CODE (right code if 6-digit codes are entered)
= 1 EQUIVALENCE THE WHOLE 6-DIGIT CODE
= 2 TREAT LEFT CODES AS 0
NOTE: Equivalent codes must be an integer from 1 to 999. If IOP12=1 you must enter equivalent codes in the
run file (set IOP11=1).
IOP13 = 0 CHECK BOTH LEFT AND RIGHT CODES WHEN SELECTING FEATURES
= 1 CHECK FULL CODE ONLY
NOTE: If you use IOP15=1 you MUST use IOP30=-1 (constant distance of plane influence).
IOP25 = CODE TO ASSIGN TO UNRESOLVED BLOCKS OF THE FIRST TYPE (default = -1)
IOP26 = CODE TO ASSIGN TO UNRESOLVED BLOCKS OF THE SECOND TYPE (default = -2)
NOTE: Used for files 8 and 9 either if file 12 is used, or if both top and bottom coordinates are entered.
NOTE: If both PAR1 and PAR2 are 0 then all VBM planes will be taken.
PAR3 = INCREMENT FOR ENTERING VBM PLANES (DEFAULT = PLANE NUMBER PRECISION)
NOTE: If PAR3=0 and IOP15=0 all planes in PAR1-PAR2 limits will be taken)
PAR11 = MAX ALLOWED DISTANCE OF INFLUENCE (default = USE ALL SPACE BETWEEN PLANES)
IF IOP30=-1 THEN PAR11 IS USED AS DISTANCE OF INFLUENCE
If file 8 is used or file 9 is used directly enter exactly one of the following two commands
END
ncode
fcode1 ... fcodeN
where
neqcode
eqcode1 ... eqcodeN
where
Free format, not more than 10 per line, lines should have at most 132 characters.Maximum allowed number of DHs is 5000.DH
labels entry ENDS with a BLANK LINE.(You may have blank lines before the start of DH label lines but not between DH label
lines.)
Enter VBM plane numbers, one per line.If you have a non-orthogonal VBM and want to override orientation parameters, enter
themon the same line as a plane number.Order of orientation parameters: X-ref, Y-ref, Z-ref, azimuth, dip.
NOTE: If you are entering plane numbers there should not be anything in the run file except maybe some blank lines after the plane
numbers lines.
PUT11 = ROCK
END
USR = ABC
IOP1 = 0 /FIRST SURVEY
IOP2 = 0 /LAST SURVEY
IOP3 = -1 /BOUNDARY OPTION
IOP4 = 1 /0=USE FILE 9 WITH FILE 12; 1=USE FILE 9 DIRECTLY
IOP5 = 2 /VBM ORIENTATION: 0=USE VBM FILE,
COM 1=W-E, 2=S-N, 3=PLANAR, 4=NON-ORTHOGONAL
IOP6 = 0 /USE SURVEYS IN IOP1-IOP2 LIMITS
IOP7 = 0 /0 = CHECK ALL SERVEYS AGAINST BOUNDARY
COM 1 = CHECK COLLARS ONLY
IOP8 = 0 /1 - CHECK EACH ASSAY AGAINST PAR4-PAR9 LIMITS
IOP11 = 1 /TYPE OF FEATURE CODES ENTRY (-1 - USE ALL)
IOP12 = 0 /1=EQUIVALENCE THE WHOLE 6-DIGIT CODE
COM 2=TREAT ALL LEFT CODES AS 0
IOP13 = 0 /1=CHECK AGAINST FULL CODE ONLY
IOP19 = 0 /1=EQUIVALENCE CODES BEFORE SORTING INTERSECTIONS
IOP30 = 0 /0=USE VOLUNE HALF-WAY BETWEEN PLANES
COM 1= USE VOLUME FROM A PLANE TO THE NEXT ONE
IOP31 = 1 /0=COORDINATES AT MIDPOINT; 1=COORDINATES AT TOE
PUT09 = ROCK
CMD = COORD EAST NORTH ELEV
END
3
501 502 503
3
123
PUT09 = ROCK
END
3
501 502 503
3
123
Note: If 9 = filename is specified on the names line, then composite data will be retrieved from this
file 9. To access the MSTorque database, use 8=TORQ and include an MSTorque info file (e.g.,
19=mstorq.inf). For information about the MSTorque info file, please refer to the M500-Technical
Section, Technical Considerations for Interpolation.
A single line of text that will be printed at the top of each page of printer output from the run.
NOTE: A block will not get interpolated if the number of nearby composites is <IOP7.
NOTE: PAR3 not used for SCAN1-D or SCAN2-D. However, still enter a small non-zero value (e.g. 0.1).
END
Area Selection Lines:
IX1 IX2 IY1 IY2 IZ1 IZ2
Columns Rows Levels
USR = abc /
END
1 100 1 100 1 40
1. The program will list all the composites within the search window of each calculated row if I-O=2 is specified in the run
file. The program will also list all the items in the composite file, the total number of points calculated for each row, and
other debug information.
2. If the ellipsoidal search is used, then the 3-D distance calculated between a composite and the block will be the adjusted
distance based on the anisotropy ratios. This is the distance that is printed on the output file.
3. Do not use composite item labels that conflict with keywords used by the program.
4. If youre using more than one ITM command, number them sequentially, without repetition.
5. Multiple Area Selection Lines are allowed. The program will calculate the points specified in the first line, then will
calculate the points specified in the second line and so on. The only exception is when ASCII input option is used and
the composite data follows the run file.
6. The maximum number of composites that the program can handle for a single row is 64000. If this number is exceeded
the program will print an error message and stop calculating at that row. If this happens, first check if you can decrease
the X, Y or Z search distances used by the program. If you can, then try again starting from the row where the program
stopped and see if this will solve the problem. If it does not, try to specify fewer columns to interpolate by using multiple
Area Selection Lines.
7. The maximum number of composites that the program can handle for a single point is 64000. If this number is
exceeded, the program will give an error message and stop calculating at that point. If this happens, try a similar
approach as explained in Step 6.
8. A separate output file can be obtained for use with M525TS if 30=filename appears in the Names Line. See M525TS
program documentation on how to further evaluate the point kriging results.
The Printer Output File contains the statistics of the original and
INPUT: CALCS: the declustered data values. The histograms for both data sets
are also output. If the declustered data is output to an ASCII file,
- PCF - Find composites within the program M500V1 can be used to load this data to another
- File 9 or 8, or an each cell File 9 or 8 for storing.
MSTorque composite - Compute average cell
set value
- Compute decluster weight The cell size specified in PAR7, PAR8 and PAR9 for X, Y, and
Z directions are very important for the de-clustering results.
The cell size in the Z direction is normally equal to the bench
height unless the user wants to see the effect of increasing the
bench height. The cell sizes in the X and Y directions should
OUTPUT: be equal to the average spacing of the samples. Too small a cell
can result in practically no declustering because each cell may
- File 9 or MSTorque contain no more than one composite. Similarly, too large a cell
database
- Printer Output file
can result in artificial declustering of the composites which, in
with statistics and fact, should not be declustered. Therefore make a few test runs
histograms of original and and study the resulting statistics of the declustered composites
declustered composites before storing the final results.
Note: If 9 = filename is specified on the names line, then composite data will be retrieved from this file 9. To
access the MSTorque database, use 9=TORQ and include an MSTorque info file (e.g., 20=mstorq.inf).
For information about the MSTorque info file, please refer to the M500-Technical Section, Technical
Considerations for Interpolation.
A single line of text information that will be printed at the top of each page of printer output from the run.
END
USR = ABC
I-O = 0 / DEBUG
END
USR = ABC
IOP1 = 10 / FIRST AND LAST DATA UNITS FROM FILE 09
IOP6 = 0 / 1=STORE DECLUSTER WTS TO FILE 9 (PUT09 REQUIRED)
IOP7 = 1 / 1=OUTPUT DECLUSTED COMPOSITES TO AN ASCII FILE
IOP11 = 40 / # OF FREQ. INTERVALS FOR HISTOGRAM (0=NO HIST., MAX=50)
I-O = 0 / DEBUG
END
USR = ABC
I-O = 0 /DEBUG
END
2. If there are less than two composites available, the following message will appear
3. If any of the cell size information (i.e., PAR7, PAR8 and PAR9) is not specified, the program will not execute.
4. If you did not specify the coordinate limits for declustering (i.e., PAR1-PAR6), the program will not default to the PCF
limits and execution ends with the following error message
5. If IOP6=1 to store the declustering weight to File 9 or 8, but you did not specify an item to store these weights, the
following message will appear
6. If IOP7=1 to dump the declustered composites to an ASCII file, but you did not specify this file on the Names Line, the
output will be written to a default File 19 named FOR019.DAT (on most machines).
NOTE: If 9 = filename is specified on the names line, then composite data will be retrieved from this file 9.
To request an ASCII input DO NOT enter 9 = on the names line and specify 8 = CARD. To access
the MSTorque database, use 8=TORQ and include an MSTorque info file (e.g., 19=mstorq.inf). For
information about the MSTorque info file, please refer to the M600-Technical Section, Technical
Considerations for Interpolation.
A single line of text that will be printed at the top of each page of printer output from the run.
NOTE: A block will not get interpolated if the number of nearby composites is <IOP7.
IOP9 = 1 OCTANT
= 2 QUADRANT
NOTE: If IOP13=1-10, one cross reference line is read. If IOP13=11-20, two cross reference lines are read.
IOP17 = THE DESIRED POINT ESTIMATE TO BE STORED (if STORE option is used)
= 1 FOR INVERSE DISTANCE WEIGHTING (IDW) POWER 1
= 2 FOR IDW POWER 1.5
= 3 FOR IDW POWER 2.0
= 4 FOR IDW POWER 2.5
= 5 FOR IDW POWER 3.0
= 6 FOR KRIGE (default)
NOTE: IOP17 can be used to specify the desired point estimate to use. If other than a kriged estimate is to be
stored, then set IOP17 in your run file as above.
NOTE: This option will interpolate the specified block codes (points) without using these codes in the
interpolation.
NOTE: PAR3 not used for SCAN1-D or SCAN2-D. However, still enter a small non-zero value (e.g. 0.1).
There are several customizing options for this program. The following table may be used to reference explanations of the uses for
various IOPs and PARs.
Customizing Page See Sample
Option Description number Run File #
I Primary 3D search parameters defining a box around a block 524-7 1
Optional full 3D (ellipsoidal) defining a trend search within the
II 524-7 2
primary search
III Optional interpolation parameters 524-8 6
VIII Item definition lines for values to be computed and stored 524-12 9
These parameters define a rectangular box around the block to be interpolated. This is illustrated in the 600 - series Technical
Section and Figure 1 below.
Figure 1. Definition of Primary Coordinate Search Parameters (PAR1, PAR2, and PAR3)
All of the composites within this rectangular box are checked against the maximum 3-D search distance (PAR4). The composites
that are less than or equal to this distance from the center of the box are retained for further screening, such as geologic matching.
If there are more than n number of composites for a given block after all the tests, then the composites are sorted by increasing distance.
Only the closest n are retained for calculation of the weights. Based on these weights, the grade of the block is computed and stored.
IMPORTANT NOTE: To interpolate GSM or 2-D surface data, the coordinate selection must be based on X and Y coordinates only.
If this is the case, use one of the following command in your run file:
II. Optional full 3-D (ellipsoidal) search parameters defining a trend search within the primary search
To add a full 3-D search within the primary search explained above include the following command:
ROTN = Rotation from the old north axis to the new north axis (clockwise is a positive rotation)
DIPN = Dip of the new north axis measured horizontally (negative dip is plunging down)
DIPE = Dip easterly from the new north axis (the horizontal component of this vector is 90 degrees clockwise from
the horizontal component of the new axis. The dip is measures as in DIPN)
For a detailed description of rotation parameters see the 600-series Technical Section.
Figures 2, 3, and 4 below and Sample Run File #2 show the relationships of the various parameters for ellipsoidal search.
Figure 2. Definition of ROTN, RY and RX used in Figure 3. Definition of DIPN, RY and RZ used in Figure 4. Definition of DIPE, RX and RZ in 3-D
3-D Ellipsoidal Coordinate Search 3-D Ellipsoidal Coordinate Search Ellipsoidal Coordinate Search
This search forms an ellipsoid around the center of the block based on the search ranges and angles used. The composites within this
ellipsoid are retained for calculations. The 3-D distances (RX, RY, RZ) between the block center and the composites are adjusted
distances based on anisotropy ratios between the ranges. Therefore, if the full 3-D search is used, the value of the maximum 3-D
search distance (PAR4) must reflect the distance in the new coordinates.
NOTE: If the number of composites retained for a block equals IOP7, then PAR8 is tested.
The M524V1 run file must contain the specifications of one variogram model to be used for kriging interpolation. To specify the
nugget and sill values for the variogram, enter the following command lines:
For a detailed description of rotation parameters, see the 600-series Technical Section.
NOTES:
1. Nested variogram models can be entered by adding more Variogram Specification Lines in your run file (max=5). For
example, a nested spherical model could be entered as:
CMD = NUGGET CO
CMD = SPH C1 R1 R2 R3 ...
CMD = SPH C2 R4 R5 R6 ...
where CO=NUGGET, CO+C1+C2=SILL, R1 R2 R3= 3-D ranges of the first structure, R4 R5 R6=3-D ranges of the
second structure, etc.
2. The program allows the interpolation of multiple grades provided that these grades have the same variogram parameters.
3. Horizontal and vertical anisotropy ratios are calculated based on the 3-D ranges entered. For example, if R1=100, R2=50
and R3=40, then horizontal anisotropy ratio = R1/R2 = 2.0 vertical anisotropy ratio = R1/R3 = 2.5
4. If rotation is specified, the 3-D ranges will be in the new coordinate system defined by these angles. For example, by
default R1 is the range in the north-south direction. However, if ROTN=45, then R1 will be the range in the NE-SW
direction. The direction of R2 and R3 will similarly be affected based on these angles. See figures 2, 3 and 4 on page
524-8.
CASE 1 There are multiple items to krige and each item has its own variogram, e.g., kriging gold and silver grades with the same
run file but having different variograms for gold and silver.
Warning: When kriging multiple items, the program assumes that the same number of composites are available for the secondary
items as for the primary item. In other words, inconsistent missing intervals such as a missing value for one item but not the other
ones, will produce erroneous results.
CASE 2 Only one item is to be kriged but there are different variograms for each geologic unit, e.g., kriging gold grades with
a different variogram for each rock type. For this case, the label of the geologic item must be specified with either of the following
commands:
In both cases, parameters for multiple variograms are entered through an external input file. The variogram parameter file name
is entered in the run file on the names line as 31=filename. When there is a file name specified in the run file, the program will
automatically look for this file and try to read it. Any variogram parameter specified in the run file through the CMD function is
ignored. Examples of the variogram file for each case follow.
NOTE: There has to be one set of variogram parameters specified for each item or geologic type to be
kriged even when the variogram parameters are the same for more than one item. As many
as nine items or geologic types can be kriged simultaneously. However, the total number of
model items to be stored cannot exceed ten. option will interpolate the specified block codes
(points) without using these codes in the interpolation.
Line #5, #6, #6A, #6B, #7 Repeat of Lines 2-4 for the variogram parameters to be used for krig item #2.
NOTE: Specify a set of variograms for each item to be kriged. For example, if five items are to be
kriged, there must be five sets of variogram parameters.
If there are more variograms than the number of krige items, the program ignores the extra
variograms. However, if there are more krige items than variograms, the program will stop
with an error message.
Line #2: Nugget value for the first geologic code Integer geologic code
M524V1 can store the estimated grade of a point from kriging and/or the estimation error computed for that point to File 9 or to an
MSTorque database. Include the following lines in your run file to use this option:
1. CALC KRIGE can be used for ordinary kriging. Kriging provides an estimation variance for each block kriged. This
value can be obtained using BLOCK ERROR.
2. ALOG KRIGE can be used in place of CALC KRIGE if you need to use the natural logarithms of grades instead of
actual grades.
3. At least one ITM to CALC KRIGE or to ALOG KRIGE must be included in each run. STORE KRIGE and STORE
ERROR are optional.
4. STORE ERROR can be used only if CALC KRIGE or ALOG KRIGE is used. Normally, ordinary kriging variance
is stored as the error of estimation, but if ALOG KRIGE is used, the stored value is the relative logarithmic kriging
variance.
In order for the program to use the composites with the same geology code in the interpolation, include the following lines in your
run file:
IOP4 = 1
Interpolation can be limited to certain geology types based on specified codes in the run file. Use the following lines in your run file:
BLOCK LIMIT = Keywords for the program indicating the option for geology limiting.
v1, v2 ... vn = Numeric values between 1 and 999999 of the items used for geology limiting.
More than one code line may be entered if needed.
label-09 = Label of composite item with geology codes
LIMIT = Keyword for the program
To assign alias (assumed) geology codes to the composites used in the interpolation of a point, specify the number of cross reference
codes in IOP13 and include a line after the END line in your run file for the cross reference codes. For example, if you have five
geology codes (1 through 5) and you want to interpolate the points with codes 1 or 2. You want to use composites with codes 1 or 3
to interpolate the points with code 1 and you also want to use composites with codes 2, 4 and 5 to interpolate the points with code
2. First include in your run file BLOCK LIMIT and MATCH options. Then set the following option in your run file: IOP13 = n. For
the above example, you would have IOP13=5, then 1 2 1 2 2 after the END line in your run file. Sample run file #6 demonstrates
the use of this option.
The composite assay data input to the interpolation programs is normally from MineSight Data File 8. However, the programs are
also designed to read data from an ASCII file. To use an ASCII input file, follow the steps below:
1. Include REF# (or DHID), X, Y and Z coordinates (if 3-D interpolation), and a grade item on each line of input.
2. Include 8=CARD in the Names Line of your run file. Also include 19=filename if the data is in a separate file.
3. Also include the following line in your run file:
REF#, EAST, NORTH, ELEV. = Keywords for the program (not necessarily in this order)
value1, value2, etc = Grade items to interpolate
geol1, etc. = Geology items to interpolate
The format of the data in the ASCII file is free field format.
4. The following line is optional and can be included to specify a CALC items minimum value. Otherwise the minimum is
assumed to be zero.
where
CALC keyword
MINn indicates which CALC item from the INP line is being set
value the minimum value.
For example:
will set the minimum value of ROCK (the second CALC item) to 7.
If composite data lines are included within the run input file, then only one Area Specification Line can be used. Sample Run File
#5 demonstrates the use of this option.
Use the following line in your run file to use this option:
Up to three lines can be entered. Sample run file #8 demonstrates the use of this option.
Composites used in interpolation can be selected using RANGE, OMIT or TABLE with the ITMn command. See ITMn commands
in the Introduction for an explanation of these options.
USE AIK keywords, factor is the proportion of area of influence. If this proportion is 1.0, a sample will have 100%
influence to the blocks within its area of influence.
END
Example: If IOP13=5 (block codes are 1 to 5) & if 1 1 3 3 3 are input for cross reference, then data with geologic codes 1 & 2 will
be used as the same type, with codes 3, 4 & 5 used together.
See the Technical Section for a complete description of these lines. If Composite Data Lines are included within the run file, then only
one Area Specification Line can be used. If Composite Data Lines are in a separate file, include a 19=filename in the Names Line.
USR = ABC
COM -
COM POINT INTERPOLATION
COM -
END
20 30 23 23 20 20
END
20 30 23 23 20 20
USR = ABC
COM -
COM POINT INTERPOLATION
COM -
END
20 30 23 23 20 20
END
000000
USR = ABC
COM -
COM POINT INTERPOLATION
COM -
END
20 30 24 28 20 20
30 DH-43 12652. 12851. 3600. 0.0580 32
31 DH-47 12720. 12752. 3500. 0.0660 43
32 DH-56 12736. 12430. 3600. 0.0520 32
33 DH-60 12705. 12509. 3600. 0.2510 04
34 DH-62 12719. 12603. 3600. 0.1260 53
35 DH-63 12725. 12550. 3650. 0.4020 92
USR = ABC
COM -
COM POINT INTERPOLATION
COM -
END
12122
20 30 23 23 20 20
USR = ABC
COM -
COM POINT INTERPOLATION
COM -
USR = ABC
COM -
COM POINT INTERPOLATION
COM -
USR = ABC
COM -
COM POINT INTERPOLATION
COM -
USR = ABC
COM -
COM POINT INTERPOLATION
COM -
END
000000
2. If IOP2=n, the program will print the location of composites, the distance to the point, and the grades and weights for n
interpolated points.
3. If the ellipsoidal search is used, then the 3-D distance calculated between a composite and the block will be the adjusted
distance based on the anisotropy ratios. This is the distance that is printed on the output file when the debug option is on.
4. The difference between PAR7 and PAR8 is very subtle. PAR7 is tested at all times. PAR8 is tested only if the number of
composites to interpolate the point is equal to IOP7.
5. Five inverse powers are used by the program. By default these powers are: 1.0, 1.5, 2.0, 2.5 and 3.0. The higher the
inverse power, the more weight is given to the nearest sample.
6. Do not use composite item labels that conflict with keywords used by the program.
7. If youre using more than one ITM command, number them sequentially, without repetition.
8. Multiple Area Selection Lines are allowed. The program will interpolate the points specified in the first line, then will
interpolate the points specified in the second line and so on. The only exception is when ASCII input option is used and
the composite data follows the run file.
9. The maximum number of composites that the program can handle for a single row is 2000. If this number is exceeded
the program will print an error message and stop interpolating at that row. If this happens, first check if you can decrease
the X, Y or Z search distances used by the program. If you can, then try again starting from the row where the program
stopped and see if this will solve the problem. If it does not, try to specify fewer columns to interpolate by using multiple
Area Selection Lines.
10. The maximum number of composites that the program can handle for a single point is 1500. If this number is exceeded,
the program will give an error message and stop interpolating at that point. If this happens, try a similar approach as
explained in Step 9.
11. A separate output file can be obtained for use with M525TS if 30=filename appears in the Names Line. See M525TS
program documentation on how to further evaluate the point kriging results.
will appear if
will appear if you do not specify any ITMn lines that will store data in the model file.
will appear if you specify more than two ITMn = CALC POLYG lines.
will appear if you specify more than two items for geology matching.
17. ** M524V1 WARNING. You can only check for duplicate points on blasthole projects. IOP21 is being reset to 0.
IOP21 is applicable to blasthole projects only. There is no need to check for duplicate composites on non-blasthole
projects. If you try to activate IOP21 for a non-blasthole project it will be reset to zero and the program will continue
normally.
18. There are 15 columns (fields) of information in the output file (30=) for M525TS. These columns are Ref. number,
Northing, Easting, Elev. (or Seam), Number of Points, Actual value, IDW estimates (5 different powers), Kriging
estimate, Kriging variance, the sum of the kriging weights X covariances, and the number from the previous column plus
the Lagrange Multiplier.
RUN FILE: If two items are selected, the program generates the statistics
for these items and their differences. The program will try to fit
- Interactive keyboard entry
a straight line through the data, using least square regression, to
see how these two items are correlated. The correlation results
are then output..
INTERACTIVE OPERATION
1. MEDS M525TS <ENTER> on the PC or m525ts <ENTER> on the Sun.
Enter a filename for printer output or <ENTER> if you want the default print file (FOR003.DAT on most computers).
Enter the name of the input file (normally this file is output from M524V1).
<ENTER> if your input file is from M524V1; otherwise, enter the necessary format to read the data.
If you are using the default format, the program will go the Step 9; otherwise, it will go to Step 7.
Enter the label of the first item to read(max = 10 characters). The label entered is for display only. It should describe the
column of information being read.
The program will repeat Step 8 as many times as the number of items entered at Step7, then it will go to Step 10.
Enter the label of the composite data file item used in point kriging. This label is optional if one item is kriged; otherwise,
it must be specified so that the program will read the correct item. This label is also used for display in output; therefore,
it is recommended. If MSTorque attribute was kriged, you must enter a column width of 20. For file 9 items, column
width (5) is optional.
10. ENTER MIN AND MAX VALUES ACCEPTABLE FOR ITEM# 1 OR <ENTER>:
Enter the minimum and maximum values for the controlling item. Defaults are 0.0 and 99999.0.
NOTE: At this point, the program will display a table of the items used and summary statistics. Use the information in
this table as a guide to select histogram and scatter plot parameters.
Enter the selected item number(s). Each item number is displayed in the summary table on the screen.
First histogram:
# of intervals = 15
starting value = 0.
the increment is computed
If a single item was selected at Step 11, the program will output the histogram and go back to Step 11. Otherwise it will
go the following step.
14. ENTER MIN AND MAX VALUES FOR SCATTER GRAPH OR <ENTER>:
NOTE: If <ENTER> is entered, default values are used or values from previous graph are assigned.
Enter the minimum and maximum values for X and Y axes of the scatter graph. Suggestion: Try to enter (MAXMIN)
divisible by 5 for a reasonable looking plot. This question will appear only if you selected option 1 or 3 at Step 13. If
your option is 1, then the program will go to Step 11. Otherwise, it will go the next step.
Enter the maximum value of the X and Y axes of the scatter plot. This question will appear only if your selection at Step
13 is 2 or 3. The program will go to Step 11 if an <ENTER> or a negative value is entered. Otherwise the program will
go to the next step.
The default plot size is 5.0 X 5.0 for a factor of 1. <ENTER> to use a plot factor of 1. You may increase or reduce the
size of the plot by entering a factor other than 1.0
Enter the name of the file to which scatter plot commands will be written. This file is an ASCII file and may be plotted by
M122V1 by using a USERF <filename> in the run file.
NOTE: If you <ENTER> instead of specifying a file name, the output is written to unit 30. The default filenames differ
depending on platform.
USR = ABC
MAP = 0 10 1 0 10 1
END
PEN #1
USERF PLT525.TST
ENDOC
CALCULATION OF COMPOSITES
Limiting Angle
Two types of fixed length composites can be calculated with M501V1: vertical and horizontal. The limiting angle controls which
method is used for a hole. By default, the limiting angle is:
NOTE: If you choose a non-default limiting angle in M501V1 (PAR4), then you must remember to include this
value in your M505V1 (PAR1) and M516V1 (PAR13) runfiles.
Vertical Composites
If the dip angle of a drillhole is greater than the limiting angle, a vertical composite is calculated. The vertical bench height (DZ) is
used as the control for compositing. The weight used in compositing is the calculated vertical length of each sample interval, not the
slope length of the interval.
Horizontal Composites
When the dip angle of a drillhole is less than the limiting angle, a horizontal composite is calculated. This prevents long vertical
composites being generated for near-horizontal drillholes. The E-W direction block size (DX) is used as the control for compositing.
The weight used in compositing is the calculated horizontal length of each sample interval, not the slope length of the interval.
Once the compositing calculation is begun with either the vertical or horizontal method, it will be used for the entire hole regardless
of subsequent surveys. There may be times when you will want to specify a limiting angle to prevent horizontal compositing of a few
holes that are close to, but less than the default angle. You should check for any holes that have surveys that are close to the limiting
angle and decide if they should be composited horizontally or vertically.
Note about using the DHID option in the composite file (File 9) without the survey file (File 12):
Starting with the 2001 update you can store drillhole labels with composites in files 9 and 8. To use this option, you must reserve two
full-word items in file9. By default, MineSight programs expect those items to be named DHID1 and DHID2, but you can name
them differently. These items are similar to DLAB1 and DLAB2 in the survey file (File 12). To initialize a full-word item enter min
= max = 0., and precision = 1.
NOTE: A 10-character drillhole ID is coded and packed into two words. The actual DH ID can be retrieved only programmatically,
and you must direct MineSight programs to retrieve DHID from file9. The option to use DHID from file9 directly was added to all
the relevant procedures.
1. NTRY is an integer argument specifying entry control as follows: NTRY=1 for subroutine initialization. The call is
made after program initialization but before any drillhole data is accessed.
NTRY=2 for each drillhole interval within the current survey data unit. Interval data is subject to ITM restrictions as are
survey data units. Collar survey data units are also subject to boundary restrictions.
NTRY=4 for user calculations at the end of all survey data units for the current drillhole.
NTRY=5 for user calculations at the end of all drillholes (or end of run).
NTRY=6 same as NTRY=2 but in the optional second pass of the data. It follows an NTRY of 4 for each drillhole.
2. IDHN is an integer argument containing the current drillhole reference number. IDHN is zero for NTRY=l and
NTRY=5.
In addition to other system common blocks, the following common block data is available:
/BUF/B12(1024), B09(1024)
1. B12 is an array containing the File 12 record in which the current survey data unit is located. See the /USR/ common
block for information of addressing data units.
2. B09 is an array containing the File 9 record in which the current interval data unit is located. Most words in B12 are
accessible directly although you may have to equivalence a floating point array to avoid conversion problems. File 9
data is usually packed in array B09 so, unless you understand MineSight packing schemes, you are advised to keep
your paws off this array.
/USR/K09GET, L09GET, K09PAK, L09PAK, LOC09, J09, K09, L09, K12GET, L12GET,
K12PAK, L12PAK, LOC12, J12, K12, L12
1. K09GET is the integer subscript for information in the ITM common block for the first File 9 GET item. Arrays of
interest are ITM(3,99), RDAT(99) and IDAT(99).
4. L09PAK is the subscript for the last File 9 PAK item. If you use UPD09 in your run file, then GET subscripts will fall
within the range of the PUT subscripts.
5. LOC09 is the base subscript of the current interval data unit and is addressed by B09 (LOC09 + 1).
7. K09 is the first data unit number of interval data for the current survey data unit.
8. L09 is the last data unit number for the current survey data unit.
THE REMAINING INTEGERS ARE THE SAME AS ABOVE BUT APPLICABLE TO File 12.
Access to MSTorque
Composite Set Data
A number MineSight programs have access to MSTorque Composites. They include the m620-interpolation programs, m508rp,
m522v1, m523v1, m524v1, m617v1 and a number of programs based on the m508v1 driver (cldiff, gn508crd, gn508grd,
gn508tru, gncmpint, vvc508).
To indicate MSTorque data, on the names line specify 9=TORQ and also add the MSTorque information file, mstorq.inf. For
all programs except interpolation, m522v1 and m524v1, enter 20=mstorq.inf. For m522v1, m524v1 and the m620-interpolation
programs (m620v1/v2, m621v1, m624v1, etc), the torque information file should is entered as file 19: 19=mstorq.inf
Structure of mstorq.inf:
version - 2.0.0.0
connection - This and project are used if you are using Windows authentication. If datasource occurs in this file, then
connection and project are ignored.
datasource - SQL authentication. If you have this keyword in the info file, then connection and project are ignored. Datasource
uses an .xml file - do not edit this file or try to prepare it by hand. The content is created by MSCompass when the MSTorque
database you are using is attached.
type - The type of sample site. To access drillholes only use drillholes for the keyword value, for blastholes use blastholes, to
access other composite sample site types use other, and to get all types of composite sample sites use all. If the keyword type
is not present, all will be used by default.
Example:
version 2.0.0.0
torquepath d:\medexe7.50\mstorque
connection tucn-ws018\sqlexpress2008r2
project osamp_2.5
type drillholes
coverage formation
Attributes that include blanks, apostrophes () and slashes (/) should be enclosed in quotes, e.g.
Keyword Description
East Mid-Span Easting
North Mid-Span Northing
Elev. Mid-Span Elevation
From Depth to the composite start
To or -TO- Depth to the composite end
Length or LNGTH Composite interval Length along the hole
XTop Easting at the top of the interval span
YTop Northing at the top of the interval span
ZTop Elevation at the top of the interval span
XBot Easting at the bottom of the interval span
YBot Northing at the bottom of the interval span
ZBot Elevation at the bottom of the interval span
TotLen Composited length
TotWt Composited weight
Top and bottom coordinates are calculated on the fly when the data is loaded to an MSBasis program.
Note: The names above are treated as keywords by MSBasis programs to access corresponding data and should not be used as
attribute names .
Programs are enabled to retrieve Enumeration Codes for Sample Attributes which have Data Type Enumeration (see MSTorque
documentation for details), modify codes and write them back to MSTorque database.