qmDoc = sys.argv[7]
    currDir = os.path.dirname(geoDB)

    arcpy.AddMessage("COC Name: " + COCName)
    
    # Local variables...
    desc = arcpy.Describe(COCLayer)
    scriptPath = sys.path[0]
    xmlTemp = scriptPath + "\\filtered_metadata_template.xml"
    COCLayerN = COCLayer.split(os.sep)[-1]
    if desc.DataType == "FeatureLayer":
        COCLayerBase = COCLayerN
    else:
        COCLayerBase = COCLayerN.split(".")[0]
    #COCLayerBase = ARD_HEA_Tools.sanitize(COCLayerBase)
    COCLayerOut = geoDB + "\\" + ARD_HEA_Tools.sanitize(COCLayerBase) + "_Layer"
    COCStats = geoDB + "\\" + ARD_HEA_Tools.sanitize(COCName) + "_Stats"
    COCFieldStat = COCField + " " + STATType
    COCStatsLyr = ARD_HEA_Tools.sanitize(COCName) + "_Stats_Layer"
    COCFiltered = geoDB + "\\" + ARD_HEA_Tools.sanitize(COCName) + "_filtered"
    COCFilteredLyr = ARD_HEA_Tools.sanitize(COCName) + "_filtered"
    SpatRef = arcpy.Describe(COCLayer).SpatialReference
    SAField = STATType + "_" + COCField
    InventTable = geoDB + "\\COC_INVENTORY"
    fltrString = arcpy.AddFieldDelimiters(COCLayer, COCField) + " > -999.99"
    chkString1 = fltrString + " AND " + arcpy.AddFieldDelimiters(COCLayer, COCField) + " <= 0"
    chkString2 = arcpy.AddFieldDelimiters(COCLayer, COCField) + " > 0 AND " + arcpy.AddFieldDelimiters(COCLayer, COCField) + " < 1"

    # Check for projected coordinate system...
    project = desc.SpatialReference.Name
    units = desc.SpatialReference.LinearUnitName.upper()
# Load required toolboxes...
sub_folder = "ArcToolbox/Toolboxes/"
install_dir = arcpy.GetInstallInfo("desktop")['InstallDir'].replace("\\","/")
tbx_home = os.path.join(install_dir, sub_folder)
arcpy.AddToolbox(tbx_home+"Data Management Tools.tbx")

try:
    # Report version...
    ver = ARD_HEA_Tools.version()
    arcpy.AddMessage("ARD HEA Tools Version: " + ver)

    # Script arguments...
    projDir = sys.argv[1]
    projNameIn = sys.argv[2]
    analystName = sys.argv[3]
    projName = ARD_HEA_Tools.sanitize(projNameIn)

    # Local variables...
    geoDBname = projName + "_GIS.mdb"
    geoDBfolder = projDir + "\\" + projName
    geoDB = geoDBfolder + "\\" + geoDBname
    prjAttr = geoDB + "\\PROJECT_ATTRIBUTES"
    COCData = geoDB + "\\COC_DATA"
    COCInvent = geoDB + "\\COC_INVENTORY"
    SiteAttr = geoDB + "\\SITE_ATTRIBUTES"
    COCAnalysis = geoDB + "\\ANALYSIS_TABLE"
    Footprint = geoDB + "\\FOOTPRINTS"

    # Create analysis folder
    arcpy.CreateFolder_management(projDir, projName)
    arcpy.AddMessage("ARD HEA Tools Version: " + ver)

    # Script arguments...
    geoDB = sys.argv[1]
    inLayer = sys.argv[2]
    habType = sys.argv[3]
    conType = sys.argv[4]
    remStat = sys.argv[5]
    subSite = sys.argv[6]
    siteDoc = sys.argv[7]
    depth = "-not applicable-"

    # Local variables...
    desc = arcpy.Describe(inLayer)
    if desc.dataType == "FeatureLayer":
        inBase = ARD_HEA_Tools.sanitize(inLayer.split(os.sep)[-1])
    else:
        inBase = desc.BaseName
    AnalysisGrid = geoDB + "\\ANALYSIS_GRID"
    AnalysisPnts = geoDB + "\\ANALYSIS_PNTS"

    inJoin = geoDB + "\\" + inBase + "_ident"
    inJoinLyr = inBase + "_ident_lyr"
    SiteAttr = geoDB + "\\SITE_ATTRIBUTES"
    tmpJoin = geoDB + "\\SITE_JOIN_VIEW"
    prjAttr = geoDB + "\\PROJECT_ATTRIBUTES"

    # Set the geoprocessing environment...
    env.overwriteOutput = 1
    env.XYTolerance = "0.00000000000000000001"