This map was created as part of a preliminary coalbed methane resource assessment conducted by the ISGS.
CODE NAME 1670 TROWBRIDGE 1760 SHELBYVILLE 1780 OPDYKE 1790 LOUDON 1800 BELLE RIVE 1870 FRIENDSVILLE 1910 BRISTOL HILL 2490 DANVILLE 2610 JAMESTOWN 2660 HERRIN 2790 SPRINGFIELD 2840 HOUCHIN CREEK 2940 SURVANT 3020 COLCHESTER 3170 SEELYVILLE 3210 DEKOVEN 3240 WILEY 3250 DAVIS 3272 SEAHORNE COAL 3320 WISE RIDGE 3330 MT. RORAH 3370 MURPHYSBORO 3490 ROCK ISLAND 3500 LITCHFIELD 3510 ASSUMPTION 3590 WILLISProcess work and notes by CC: The procedure was to adapt cpk's amls from various availability projects and calc tonnage for the Keystone writeup to get amls that do what we need.
A couple of false starts involved going to the library to dissolve the various res covers to thickness covers over here in the cbm space. Besides hitting the snag that the res2660 covers wouldn't dissolve or mapjoin, we realized that after dissolving, we could no longer easily get back the areas that should have thickness = 0 (outside the crop, in channels, in areas of insufficient data, and mined-out areas). Thus, we needed to copy the res covers to the cbm space (copy_res.aml), then update them for this purpose.
calc-thick-0.aml reselects for ctype = 3, 5, 6, or 9 and then calcs thick = 0; runcalc0.aml runs calc-thick-0.aml for each specific res cover.
CTYPE96 1 Minable coal 2 Split coal 3 Channel 4 Less than 150' deep (Not used for CTYPE96. This was used for earlier CTYPEs) 5 Insufficient data 6 Mined-out area 7 Outside crop 8 Non-study area 9 Unminable area between mines 10 Coal too thinthen rundis_all.aml invokes dis.aml for each specific res cover that had been copied over; dis.aml dissolves on the item THICK so that only the polygons referring to thickness are left. rundis_all does this for each res cover, so that coverages named tk[fips #]_ [seam] were created.
then join_all.aml was run to mapjoin all thickness coverages created in the last step for each seam. -note, seam 2660 was omitted; due to problems with the input res covers, these need to be treated separately (it makes the process hang; they may not join at all; how 2660 Herrin seam was treated will be described separately)
Then the item THICK needs to be altered in info so that when the different seams are unioned, it will be clear which seam has what thickness. Otherwise, the item THICK will be overwritten with each successive union, and there would be no way to calculate a cumulative thickness. alter.aml alters the item THICK (a 3,3,I item) to seam-specific thick items
** NOTE: this did not work, I did them by hand. first i had an aml that did it by itself, but it would not run due to error in the &data arc info line. so then i copied the format of runcalc0.aml and made an aml that should have worked by invoking an altr.aml for each one, and it did go in & out of info, but it did not make the alterations and i don't know why. check w/ cpk. it was taking longer to make the aml than to do it by hand, so i just did it, but i'd like to know why it did not work.
union all of them
Usage: UNION <in_cover> <union_cover> <out_cover> {fuzzy_tolerance}
{JOIN | NOJOIN}
union tk_1670 tk_1760 tmpa 4 join
union tmpa tk_1780 tmpb 4 join
union tmpb tk_1790 tmpc 4 join
union tmpc tk_1800 tmpd 4 join
union tmpd tk_1870 tmpe 4 join
union tmpe tk_1910 tmpf 4 join
union tmpf tk_2490 tmpg 4 join
union tmpg tk_2610 tmph 4 join
union tmph tk_2790 tmpi 4 join
union tmpi tk_2840 tmpj 4 join
union tmpj tk_2940 tmpk 4 join
union tmpk tk_3020 tmpl 4 join
union tmpl tk_3170 tmpm 4 join
union tmpm tk_3210 tmpn 4 join
union tmpn tk_3240 tmpo 4 join
dropitem TMPO.PAT TMPO.PAT TMPN#
dropitem TMPO.PAT TMPO.PAT TMPN-ID
dropitem TMPO.PAT TMPO.PAT TMPM#
dropitem TMPO.PAT TMPO.PAT TMPM-ID
dropitem TMPO.PAT TMPO.PAT TMPL#
dropitem TMPO.PAT TMPO.PAT TMPL-ID
dropitem TMPO.PAT TMPO.PAT TMPK#
dropitem TMPO.PAT TMPO.PAT TMPK-ID
dropitem TMPO.PAT TMPO.PAT TMPJ#
dropitem TMPO.PAT TMPO.PAT TMPJ-ID
dropitem TMPO.PAT TMPO.PAT TMPI#
dropitem TMPO.PAT TMPO.PAT TMPI-ID
dropitem TMPO.PAT TMPO.PAT TMPH#
dropitem TMPO.PAT TMPO.PAT TMPH-ID
dropitem TMPO.PAT TMPO.PAT TMPG#
dropitem TMPO.PAT TMPO.PAT TMPG-ID
dropitem TMPO.PAT TMPO.PAT TMPF#
dropitem TMPO.PAT TMPO.PAT TMPF-ID
dropitem TMPO.PAT TMPO.PAT TMPE#
dropitem TMPO.PAT TMPO.PAT TMPE-ID
dropitem TMPO.PAT TMPO.PAT TMPD#
dropitem TMPO.PAT TMPO.PAT TMPD-ID
dropitem TMPO.PAT TMPO.PAT TMPC#
dropitem TMPO.PAT TMPO.PAT TMPC-ID
dropitem TMPO.PAT TMPO.PAT TMPB#
dropitem TMPO.PAT TMPO.PAT TMPB-ID
dropitem TMPO.PAT TMPO.PAT TMPA#
dropitem TMPO.PAT TMPO.PAT TMPA-ID
dropitem TMPO.PAT TMPO.PAT TK_1670#
dropitem TMPO.PAT TMPO.PAT TK_1670-ID
dropitem TMPO.PAT TMPO.PAT TK_1760#
dropitem TMPO.PAT TMPO.PAT TK_1760-ID
dropitem TMPO.PAT TMPO.PAT TK_1780#
dropitem TMPO.PAT TMPO.PAT TK_1780-ID
dropitem TMPO.PAT TMPO.PAT TK_1790#
dropitem TMPO.PAT TMPO.PAT TK_1790-ID
dropitem TMPO.PAT TMPO.PAT TK_1800#
dropitem TMPO.PAT TMPO.PAT TK_1800-ID
dropitem TMPO.PAT TMPO.PAT TK_1870#
dropitem TMPO.PAT TMPO.PAT TK_1870-ID
dropitem TMPO.PAT TMPO.PAT TK_1910#
dropitem TMPO.PAT TMPO.PAT TK_1910-ID
dropitem TMPO.PAT TMPO.PAT TK_2490#
dropitem TMPO.PAT TMPO.PAT TK_2490-ID
dropitem TMPO.PAT TMPO.PAT TK_2610#
dropitem TMPO.PAT TMPO.PAT TK_2610-ID
dropitem TMPO.PAT TMPO.PAT TK_2790#
dropitem TMPO.PAT TMPO.PAT TK_2790-ID
dropitem TMPO.PAT TMPO.PAT TK_2840#
dropitem TMPO.PAT TMPO.PAT TK_2840-ID
dropitem TMPO.PAT TMPO.PAT TK_2940#
dropitem TMPO.PAT TMPO.PAT TK_2940-ID
dropitem TMPO.PAT TMPO.PAT TK_3020#
dropitem TMPO.PAT TMPO.PAT TK_3020-ID
dropitem TMPO.PAT TMPO.PAT TK_3170#
dropitem TMPO.PAT TMPO.PAT TK_3170-ID
dropitem TMPO.PAT TMPO.PAT TK_3210#
dropitem TMPO.PAT TMPO.PAT TK_3210-ID
dropitem TMPO.PAT TMPO.PAT TK_3240#
dropitem TMPO.PAT TMPO.PAT TK_3240-ID
union tmpo tk_3250 tmpp 4 join
union tmpp tk_3272 tmpq 4 join
union tmpq tk_3320 tmpr 4 join
union tmpr tk_3330 tmps 4 join
union tmps tk_3370 tmpt 4 join
union tmpt tk_3490 tmpu 4 join
union tmpu tk_3500 tmpv 4 join
union tmpv tk_3510 tmpw 4 join
union tmpw tk_3590 tmpy 4 join
dropitem TMPY.PAT TMPY.PAT TMPW#
dropitem TMPY.PAT TMPY.PAT TMPW-ID
dropitem TMPY.PAT TMPY.PAT TMPV#
dropitem TMPY.PAT TMPY.PAT TMPV-ID
dropitem TMPY.PAT TMPY.PAT TMPU#
dropitem TMPY.PAT TMPY.PAT TMPU-ID
dropitem TMPY.PAT TMPY.PAT TMPT#
dropitem TMPY.PAT TMPY.PAT TMPT-ID
dropitem TMPY.PAT TMPY.PAT TMPS#
dropitem TMPY.PAT TMPY.PAT TMPS-ID
dropitem TMPY.PAT TMPY.PAT TMPR#
dropitem TMPY.PAT TMPY.PAT TMPR-ID
dropitem TMPY.PAT TMPY.PAT TMPQ#
dropitem TMPY.PAT TMPY.PAT TMPQ-ID
dropitem TMPY.PAT TMPY.PAT TMPP#
dropitem TMPY.PAT TMPY.PAT TMPP-ID
dropitem TMPY.PAT TMPY.PAT TMPO#
dropitem TMPY.PAT TMPY.PAT TMPO-ID
dropitem TMPY.PAT TMPY.PAT TK_3210-ID
dropitem TMPY.PAT TMPY.PAT TK_3250#
dropitem TMPY.PAT TMPY.PAT TK_3250-ID
dropitem TMPY.PAT TMPY.PAT TK_3272#
dropitem TMPY.PAT TMPY.PAT TK_3272-ID
dropitem TMPY.PAT TMPY.PAT TK_3320#
dropitem TMPY.PAT TMPY.PAT TK_3320-ID
dropitem TMPY.PAT TMPY.PAT TK_3330#
dropitem TMPY.PAT TMPY.PAT TK_3330-ID
dropitem TMPY.PAT TMPY.PAT TK_3370#
dropitem TMPY.PAT TMPY.PAT TK_3370-ID
dropitem TMPY.PAT TMPY.PAT TK_3490#
dropitem TMPY.PAT TMPY.PAT TK_3490-ID
dropitem TMPY.PAT TMPY.PAT TK_3500#
dropitem TMPY.PAT TMPY.PAT TK_3500-ID
dropitem TMPY.PAT TMPY.PAT TK_3510#
dropitem TMPY.PAT TMPY.PAT TK_3510-ID
dropitem TMPY.PAT TMPY.PAT TK_3590#
dropitem TMPY.PAT TMPY.PAT TK_3590-ID
dropitem TMPY.PAT TMPY.PAT TK_3500#
****** start here, messing w/ #6; PROBLEMS FROM MAPJOINS OF RES COVERS
RECORDED IN join_2660.aml 9it stopped being an aml when it didn't
run and i recorded what worked and what didn't in that file)
/coal14/avail/herrin/all-layers/transfer-to-usgs/thick6
note that herrn_tk in /cbm/thick has codes like 2, 3, 5 for ss channel,
insuff. data, thin coal. these would be treated like inches, so this
thick6 cover is the better one to start with.
need to union with a statewide mines cover
note that I can't reselect the ctype = 6 (Mined-out area) and ctype =
9 (Unminable area between mines) from the res covers because i'd
get the same problem of not being able to put them together as i
did with the whole and dissolved versions of the res covers. I will
have to assume that the unminable area between mines is negligible.
As it turns out, I should be able to get ctype = 9 by using the
type = A in the mine coverage.
then calc the mined areas to thick = 0
(first check that item thick is 3, 3, i), YES!
Arc: copy /coal14/avail/herrin/all-layers/transfer-to-usgs/thick6 thick6
copy /coal11/coalgisdb/herrn_mines mines6
union thick6 mines6 tkm6 4 join
info
SEL TKM6.PAT
RES TYPE CN 'A' OR TYPE CN 'B' OR TYPE CN 'K' OR TYPE CN 'L' OR TYPE CN 'M' OR TYPE CN 'N'
CALC THICK = 0
q stop
dissolve tkm6 tk_2660 THICK poly
union tmpy tk_2660 thickall 4 join
altered item thick to TK2660
dropitem THICKALL.PAT THICKALL.PAT TK_2660#
dropitem THICKALL.PAT THICKALL.PAT TK_2660-ID
dropitem THICKALL.PAT THICKALL.PAT TMPY#
dropitem THICKALL.PAT THICKALL.PAT TMPY-ID
NOTE: 57,203 records
additem
Usage: ADDITEM <in_info_file> <out_info_file> <item_name> <item_width>
<output_width> <item_type> {decimal_places} {start_item}
additem THICKALL.PAT THICKALL.PAT TKA 4 4 I
additem THICKALL.PAT THICKALL.PAT TKB 4 4 I
additem THICKALL.PAT THICKALL.PAT THICKALL 4 4 I
INFO
ARC
SEL THICKALL.PAT
CALC TKA = TK1670 + TK1760 + TK1780 + TK1790 + TK1800 + TK1870 + TK1910 + TK2490 + TK2610 + TK2790 + TK2840 + TK2940 + TK3020
CALC TKB = TK3170 + TK3210 + TK3240 + TK3250 + TK3272 + TK3320 + TK3330 + TK3370 + TK3490 + TK3500 + TK3510 + TK3590 + TK2660
CALC THICKALL = TKA + TKB
coverage thickall is the original gross thickness map.
next, copy this and make a real gross thickness map; there are some
negative thicknesses (range is up to -57) and this may be altering the
cumulative thickness.
copy thickall thickallg
SEL THICKALLG.PAT
RES TK1670 < 0
CALC TK1670 = 0
ASEL
RES TK1760 < 0
CALC TK1760 = 0
ASEL
RES TK1780 < 0
CALC TK1780 = 0
ASEL
RES TK1790 < 0
CALC TK1790 = 0
ASEL
RES TK1800 < 0
CALC TK1800 = 0
ASEL
RES TK1870 < 0
CALC TK1870 = 0
ASEL
RES TK1910 < 0
CALC TK1910 = 0
ASEL
RES TK2490 < 0
CALC TK2490 = 0
ASEL
RES TK2610 < 0 THIS HAD 306
CALC TK2610 = 0
ASEL
RES TK2790 < 0
CALC TK2790 = 0
ASEL
RES TK2840 < 0
CALC TK2840 = 0
ASEL
RES TK2940 < 0
CALC TK2940 = 0
ASEL
RES TK3020 < 0
CALC TK3020 = 0
ASEL
RES TK3170 < 0
CALC TK3170 = 0
ASEL
RES TK3210 < 0
CALC TK3210 = 0
ASEL
RES TK3240 < 0
CALC TK3240 = 0
ASEL
RES TK3250 < 0
CALC TK3250 = 0
ASEL
RES TK3272 < 0
CALC TK3272 = 0
ASEL
RES TK3320 < 0
CALC TK3320 = 0
ASEL
RES TK3330 < 0
CALC TK3330 = 0
ASEL
RES TK3370 < 0
CALC TK3370 = 0
ASEL
RES TK3490 < 0
CALC TK3490 = 0
ASEL
RES TK3500 < 0
CALC TK3500 = 0
ASEL
RES TK3510 < 0
CALC TK3510 = 0
ASEL
RES TK3590 < 0
CALC TK3590 = 0
ASEL
RES TK2660 < 0
CALC TK2660 = 0
ASEL
CALC TKA = TK1670 + TK1760 + TK1780 + TK1790 + TK1800 + TK1870 + TK1910 + TK2490 + TK2610 + TK2790 + TK2840 + TK2940 + TK3020
CALC TKB = TK3170 + TK3210 + TK3240 + TK3250 + TK3272 + TK3320 + TK3330 + TK3370 + TK3490 + TK3500 + TK3510 + TK3590 + TK2660
CALC THICKALL = TKA + TKB
then copy thickall and make a net thickness map, which will drop out all
thicknesses less than or equal to 18" before adding them
copy thickall thickalln
SEL THICKALLN.PAT
RES TK1670 < 18
CALC TK1670 = 0
ASEL
RES TK1760 < 18
CALC TK1760 = 0
ASEL
RES TK1780 < 18
CALC TK1780 = 0
ASEL
RES TK1790 < 18
CALC TK1790 = 0
ASEL
RES TK1800 < 18
CALC TK1800 = 0
ASEL
RES TK1870 < 18
CALC TK1870 = 0
ASEL
RES TK1910 < 18
CALC TK1910 = 0
ASEL
RES TK2490 < 18
CALC TK2490 = 0
ASEL
RES TK2610 < 18
CALC TK2610 = 0
ASEL
RES TK2790 < 18
CALC TK2790 = 0
ASEL
RES TK2840 < 18
CALC TK2840 = 0
ASEL
RES TK2940 < 18
CALC TK2940 = 0
ASEL
RES TK3020 < 18
CALC TK3020 = 0
ASEL
RES TK3170 < 18
CALC TK3170 = 0
ASEL
RES TK3210 < 18
CALC TK3210 = 0
ASEL
RES TK3240 < 18
CALC TK3240 = 0
ASEL
RES TK3250 < 18
CALC TK3250 = 0
ASEL
RES TK3272 < 18
CALC TK3272 = 0
ASEL
RES TK3320 < 18
CALC TK3320 = 0
ASEL
RES TK3330 < 18
CALC TK3330 = 0
ASEL
RES TK3370 < 18
CALC TK3370 = 0
ASEL
RES TK3490 < 18
CALC TK3490 = 0
ASEL
RES TK3500 < 18
CALC TK3500 = 0
ASEL
RES TK3510 < 18
CALC TK3510 = 0
ASEL
RES TK3590 < 18
CALC TK3590 = 0
ASEL
RES TK2660 < 18
CALC TK2660 = 0
ASEL
CALC TKA = TK1670 + TK1760 + TK1780 + TK1790 + TK1800 + TK1870 + TK1910 + TK2490 + TK2610 + TK2790 + TK2840 + TK2940 + TK3020
CALC TKB = TK3170 + TK3210 + TK3240 + TK3250 + TK3272 + TK3320 + TK3330 + TK3370 + TK3490 + TK3500 + TK3510 + TK3590 + TK2660
CALC THICKALL = TKA + TKB
dissolve thickalln netthick THICKALL poly
cleanup:
killer_res.aml kills all the res covers copied over and messed with
killer_tk.aml kills all the individual thickness covers dissolved out of
the res covers
Polygon attributes:
NET_THICKNESS - cumulative (net) coal thickness:
0 - no coal included (e.g no coal present, mined,
cut-out, not mapped, etc.)
60 - 18 inches to 5 feet
120 - 5 to 10 feet
180 - 10 to 15 feet
240 - 15 to 20 feet
300 - 20 to 25 feet
ITEM_LABEL - text label of above thickness ranges
In no event shall the DNR have any liability whatsoever for payment of any consequential, incidental, indirect, special, or tort damages of any kind, including, but not limited to, any loss of profits arising out of use of or reliance on the geographic data or arising out of the delivery, installation, operation, or support by DNR.