Python - Wetland Mitigation & Evaluation process

This script is created for Wetland Mitigation & Evaluation process. NJ transit assets/property were intersected with wetlands from different organisation.

# Import arcpy module
import arcpy
from arcpy import env
from arcpy.sa import *
import numpy as np

arcpy.env.overwriteOutput = True

arcpy.env.workspace = arcpy.GetParameterAsText(0)
##ws = arcpy.env.workspace = r"W:\Projects\NJ_Transit\WetlandsEvaluation\PhaseI\IdentifiedParcels\Rail\NJCL\NJCL.gdb"

#Local variable
WMA = "WMA"
HUC8 = "HUC8"
HUC14 = "HUC14"

for fc in arcpy.ListFeatureClasses("*Parcels"):
    print fc
    # Process: Intersect
    inFeatures = [fc, "WMA", "HUC8", "HUC14"]
    intersectOutput = "WHH_" + str(fc)
    arcpy.Intersect_analysis(inFeatures, intersectOutput, "ALL", "", "INPUT")

    # Process: Delete Field
    arcpy.DeleteField_management(intersectOutput, "FID_NWI_Wetlands_ClippedToParcels;ACRES;SHAPE_Leng;FID_HUC8;FID_WMA;DEPWMAS_;DEPWMAS_ID")

    # Process: Add Geometry Attributes
    arcpy.AddGeometryAttributes_management(intersectOutput, "AREA;PERIMETER_LENGTH", "FEET_US", "SQUARE_FEET_US", "PROJCS['NAD_1983_UTM_Zone_18N',GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',1640416.666666667],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',-75.0],PARAMETER['Scale_Factor',0.9996],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Foot_US',0.3048006096012192]]")