This makes perfect sense to me Melinda Morang. I am still getting an error for some reason.
'''python
import arcpy
from arcpy import env
import os, datetime
import numpy as np
import math
arcpy.CheckOutExtension("network")
path = r"C:\afclark82\ChildCare_Final"
gdb = r"C:\afclark82\ChildCare_Final\ChildCare_Final.gdb"
nds = gdb + r"\Network\NetworkDataSet"
facilities = gdb + r"\Locations\ChildCareLocations"
incidents = gdb + r"\Locations\ChildHomeLocations"
nd_close_fac = "ShortRoute"
nd_service_area = "ServiceArea"
Scratch = r"C:\afclark82\ChildCare_Final\Scratch.gdb"
shortPath = Scratch + r"\ShortestPaths"
serviceAreas = Scratch + r"\ServiceAreas"
travel_mode = "Driving"
if arcpy.Exists(shortPath)==False:
arcpy.CreateFeatureDataset_management(path+ r"\Scratch.gdb", "ShortestPaths")
if arcpy.Exists(Scratch + r"\ServiceAreas")==False:
arcpy.CreateFeatureDataset_management(Scratch,"ServiceAreas")
########################################
### CALCULATE SERVICE AREAS FOR MAPS ###
########################################
if len(arcpy.ListFields(facilities,"cutoff"))==0:
arcpy.AddField_management(facilities, "cutoff", "TEXT")
arcpy.management.CalculateField(facilities, "Per50", "np.around(!Per50!,2)", "PYTHON3", "import numpy as np")
arcpy.management.CalculateField(facilities, "Per62", "np.around(!Per62!,2)", "PYTHON3", "import numpy as np")
arcpy.management.CalculateField(facilities, "Per75", "np.around(!Per75!,2)", "PYTHON3", "import numpy as np")
arcpy.management.CalculateField(facilities, "cutoff", '"{} {} {}".format(!Per50!,!Per62!,!Per75!)', "PYTHON3", None)
print("Create a new service area")
#############################
## INITIALIZE THE ANALYSIS ###
##############################
arcpy.nax.MakeNetworkDatasetLayer(nds, nd_service_area)
ServArea = arcpy.nax.ServiceArea(nd_service_area)
#####################
## SET PROPERTIES ###
#####################
ServArea.distanceUnits = arcpy.nax.DistanceUnits.Kilometers
ServArea.travelMode = travel_mode
ServArea.outputType = arcpy.nax.ServiceAreaOutputType.Polygons
ServArea.geometryAtOverlap = arcpy.nax.ServiceAreaOverlapGeometry.Overlap
######################
## LOAD THE INPUTS ###
######################
# MAP FIELDS FOR FACILITIES
fm_sa = ServArea.fieldMappings(arcpy.nax.ServiceAreaInputDataType.Facilities, True)
fm_sa["Name"].mappedFieldName = "CentreID"
fm_sa["Breaks"].mappedFieldName = "cutoff"
whereclause = "cutoff <> 'None None None'"
fac = arcpy.SelectLayerByAttribute_management(facilities,"NEW_SELECTION", whereclause)
arcpy.CopyFeatures_management(fac, Scratch + r"\fac_p")
ServArea.load(arcpy.nax.ServiceAreaInputDataType.Facilities, Scratch + r"\fac_p", fm_sa)
#########################
## SOLVE THE ANALYSIS ###
#########################
result_sa = ServArea.solve()
result_sa.export(arcpy.nax.ServiceAreaOutputDataType.Facilities,gdb + "\FinalAnalysis\SrvArea_Facil")
# Export the results to a feature class
if result_sa.solveSucceeded:
result_sa.export(arcpy.nax.ServiceAreaOutputDataType.Polygons,gdb + "\FinalAnalysis\SrvArea_Final")
else:
arcpy.AddMessage("Solved failed")
arcpy.AddMessage(result_sa.solverMessages(arcpy.nax.MessageSeverity.All))
'''