I created an Arcpy script to generate point shapefiles for different watersheds.
I used split by features to create single watersheds. Each watershed has unique 3 letter code and I used that code to split features. I was not able to create single watershed shapefiles but I was able to create feature class in a geodatabase.
Then, I used these 3 letter codes feature class to create point shapefiles as boundary for each watershed in a loop. I realized when it comes to 'CON' loop stopped.
<SPAN class="" style="color: #101094; border: 0px; font-weight: inherit; font-size: 13px;">import</SPAN><SPAN class="" style="color: #303336; border: 0px; font-weight: inherit; font-size: 13px;"> arcpy</SPAN><SPAN class="" style="color: #101094; border: 0px; font-weight: inherit; font-size: 13px;">from</SPAN><SPAN class="" style="color: #303336; border: 0px; font-weight: inherit; font-size: 13px;"> collections </SPAN><SPAN class="" style="color: #101094; border: 0px; font-weight: inherit; font-size: 13px;">import</SPAN><SPAN class="" style="color: #303336; border: 0px; font-weight: inherit; font-size: 13px;"> defaultdict
geo_db </SPAN><SPAN class="" style="color: #303336; border: 0px; font-weight: inherit; font-size: 13px;">=</SPAN><SPAN class="" style="color: #303336; border: 0px; font-weight: inherit; font-size: 13px;"> </SPAN><SPAN class="" style="color: #7d2727; border: 0px; font-weight: inherit; font-size: 13px;">'C:\\CHIA_Routine\\CHIA_Watersheds.gdb'
</SPAN><SPAN class="" style="color: #303336; border: 0px; font-weight: inherit; font-size: 13px;">smis_wq_resultsXY </SPAN><SPAN
