c17q3 changes

fixes that include route geo rebuild issue. update scenario numbering.

Author: nrferguson

apply_future_rail_actions.sas

@@ -78,7 +78,7 @@ filename inrte "&inpath.\Temp\rte_out.txt";
         data a2rpt(keep=descr ln reduc); set a2rpt(where=(layover='99'));
             ln=substr(tr_line,1,3);
             reduc=trv_time*100;
-        proc sort data=a2rpt nodupkey; by ln;
+        proc sort data=a2rpt nodupkey; by descr;
         data __null__; set a2rpt;
             file "&reportpath" mod;
             put "--> " descr ": reduced " ln "travel time along entire line by " reduc "percent";

create_Emme_rail_files_GTFS.sas

@@ -29,7 +29,7 @@ ________________________________________________________________________________
 %let reportpath=&outpath.\&scen.00\rail_changes.txt;
 %let maxzone=1961;                                                        *** highest zn09 POE zone number;
 %let dropnode=48007;                                                      *** Exclude South Bend station until zone system expands to include it;
-%let basescen=200;                                                        *** base year scenario;
+%let basescen=100;                                                        *** base year scenario;
 %let counter=1;
 %let tod=0;
 %let a1=0;%let a2=0;%let a3=0;%let a4=0;%let a5=0;%let a7=0;                        *** initialize action code counts;
@@ -72,7 +72,7 @@ data __null__;
         put "--> removed scenario = " scenario "people mover link anode = " anode "bnode = " bnode "from network";
 
 *---------------------------------------------------------------*;
-    *** READ IN FUTURE CODING DATA FOR SCENARIOS 200 - 600 ***;
+    *** READ IN FUTURE CODING DATA FOR SCENARIOS 200 - 700 ***;
 *---------------------------------------------------------------*;
 %macro future;
   %if &sc>&basescen %then %do;
@@ -223,7 +223,7 @@ run;
 
        * - - - - - - - - - - - - - - - - - - - - - - - - - - *;
              **VERIFY THAT ITINERARIES DO NOT STOP AT JUNCTIONS**;
-        data check; set combine(where=(((itin_b>38999 & itin_b<40000) | (itin_b>48999 & itin_b<50000)) & dw_code=0));
+        data check; set combine(where=(((itin_b>39000 & itin_b<40000) | (itin_b>49000 & itin_b<50000)) & dw_code=0));
             proc print; title "NETWORK &scen.00 ITINERARY SEGMENTS THAT STOP AT A JUNCTION";
        * - - - - - - - - - - - - - - - - - - - - - - - - - - *;
 

create_scenario_files_GTFS.py

@@ -33,14 +33,15 @@
 
 ##======================================================================##
 
-## scenario options: beginning with c10q3 
+## scenario options: beginning with c17q3
 ## --------------------------------------
-##            100 - 2010 network
-##            200 - 2015 network
-##            300 - 2020 network
-##            400 - 2025 network 
-##            500 - 2030 network
+##            100 - 2015 network
+##            200 - 2020 network
+##            300 - 2025 network
+##            400 - 2030 network 
+##            500 - 2035 network
 ##            600 - 2040 network
+##            700 - 2050 network
 
 ##======================================================================##
 
@@ -66,7 +67,7 @@
 f = string.replace(c, "\\", '/') + "/mrn_programs"
 newdir = path + "\\" + scenario
 current = "all_runs"                                   # current routes
-if scenario in ("200","900"):
+if scenario in ("100","900"):
     current = "all_runs_base"                          # base year routes
 future = "future"                                      # future coding
 rail_routes = current

itinerary_node_update_Aug2012.sas

@@ -1,5 +1,9 @@
-*NRF revised 5-30-2017 to use action codes to inform future it_order after split link.;
-
+/*
+NRF 5/30/2017 - corrected future it_order overwrite issue using action codes to inform future it_order after split link.
+NRF 6/2/2017 - corrected error in rebuilding future route geo. was not using node ID of 0 to find coordinates for new node.
+NRF 6/23/2017 - corrected error in rebuiding future route geo. new node coordinates were being used to apply node ids to
+                route geo based on original node coordinates when node is moved. 
+*/
 *================================================================*;
    *** UPDATE ITINERARY TO REFLECT SPLIT LINKS W/ TEMPORARY ANODE-BNODE VALUES (ONLY IF NEW_MILE.DBF EXISTS) ***;
 *================================================================*;
@@ -366,16 +370,18 @@ filename delndf "&dir.\Temp\deleted_node.dbf";
             input line_num tr_line $ x y m;
 
         proc import datafile="&dir.\Temp\new_node.dbf" out=rail_nodes replace;
-        data rail_nodes(keep=node point_x point_y); set rail_nodes;
-        data rail_nodes(drop=point_x point_y); set rail_nodes;
+        data rail_nodes(keep=node point_x point_y point_x0 point_y0); set rail_nodes;
+        data rail_nodes(drop=point_x point_y point_x0 point_y0); set rail_nodes;
             x=point_x;
-	    y=point_y;
+	        y=point_y;
+            x0=point_x0;
+            y0=point_y0;run;
 
 	proc sql noprint;
 	    create table rte_nodes as
-	    select rte_geo.line_num,rte_geo.tr_line,rte_geo.x,rte_geo.y,rte_geo.m,rail_nodes.node
+	    select rte_geo.line_num,rte_geo.tr_line,rail_nodes.x,rail_nodes.y,rte_geo.m,rail_nodes.node
 	    from rte_geo,rail_nodes
-	    where rte_geo.x=rail_nodes.x & rte_geo.y=rail_nodes.y;
+	    where rte_geo.x=rail_nodes.x0 & rte_geo.y=rail_nodes.y0;
 	proc sort data=rte_nodes;
 	    by line_num m;
 
@@ -390,6 +396,8 @@ filename delndf "&dir.\Temp\deleted_node.dbf";
 	        set tmpnode;
 	        itin_a=anode;
 	        itin_b=bnode;
+            if tempa>90000 then itin_a=0;
+            if tempb>90000 then itin_b=0;
 
 	    proc sort data=rte_itin;
 	        by itin_a itin_b;

update_network_edits_Aug2012.py

@@ -27,6 +27,8 @@
 #   08-07-2012: Revised to iterate through & update all route systems.      #
 #   05-30-2017: revised to pass future route action codes to SAS for        #
 #               processing split links.                                     #
+#   06-23-2017: revised to remove point_x0 and point_y0 fields from final   #
+#               node feature class.                                         #
 #                                                                           #
 #############################################################################
 
@@ -185,6 +187,7 @@
 arcpy.AddMessage("---> Updating Node Feature Class")
 arcpy.MakeXYEventLayer_management(new_node_dbf, "point_x", "point_y", temp_node_Layer, "PROJCS['NAD_1927_StatePlane_Illinois_East_FIPS_1201',GEOGCS['GCS_North_American_1927',DATUM['D_North_American_1927',SPHEROID['Clarke_1866',6378206.4,294.9786982]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',-88.33333333333333],PARAMETER['Scale_Factor',0.999975],PARAMETER['Latitude_Of_Origin',36.66666666666666],UNIT['Foot_US',0.3048006096012192]];IsHighPrecision")
 arcpy.FeatureClassToFeatureClass_conversion(temp_node_Layer, mrn_gdb, "temp")
+arcpy.DeleteField_management(temp, ["point_x0", "point_y0"])
 arcpy.SelectLayerByAttribute_management(railnet_node, "CLEAR_SELECTION", "")
 arcpy.DeleteRows_management(railnet_node)
 arcpy.Append_management(temp, railnet_node, "TEST", "", "") 
@@ -239,6 +242,7 @@
 env.workspace = "C:/Master_Rail/mrn.gdb/railnet"             ## point inside feature dataset
 fcs = arcpy.ListFeatureClasses('',"arc")
 fcs.remove("railnet_arc")
+
 i = 0
 for fc in fcs:
     arcpy.AddMessage("---> Updating Geometry for " + fcs[i] + " Route System ...")
@@ -275,7 +279,7 @@
         f -= 1
         arcpy.AddMessage("---> Geometry Written for " + str(f) + " Future Routes")
         outFile.close()
-        
+    
     ## Run SAS to Update Itineraries ##
     # -- finish set up to run SAS
     y2 = c + "$" + orig_itinerary_dbf + "$X$3$X$X$" + orig_future_routes_dbf

update_nodes.sas

@@ -13,7 +13,9 @@
       3. new nodes due to added link - new XY coordinates will be calculated automatically, will need node number  
       4. new nodes due to split links - extra processing required to determine which node(s) are new because arc
                                         table will show duplicate entries of original link with anode & bnode values.
-  */
+
+    6/23/17 NRF: revised to retain original coordinates (point_x0, point_y0) in new_node.dbf for use when rebuilding route geo.
+*/
 
 
 %let tot=0;
@@ -115,13 +117,13 @@ proc sort data=keepnodes; by node;
 
    *** READ IN CURRENT NODE DATA AND ATTACH ATTRIBUTE DATA TO NEW NODE DATASET ***;
 proc import datafile="&dir.\Temp\temp_node.dbf" out=c replace;
-data c(drop=point_x point_y); set c; proc sort; by node;
+data c(rename=(point_x=point_x0 point_y=point_y0)); set c; proc sort; by node;
 data keepnodes(drop=_type_ _freq_); merge keepnodes (in=hit) c; by node; if hit;
   if pspace=. then pspace=0; if pcost=. then pcost=0;
 
 
    *** RE-ORDER VARIABLES TO MATCH FEATURE CLASS ORDER ***;
-data final; retain node label pspace pcost ftr_pspace ftr_pcost point_x point_y; set keepnodes;
+data final; retain node label pspace pcost ftr_pspace ftr_pcost point_x point_y point_x0 point_y0; set keepnodes;
 
 proc export data=final outfile="&dir.\Temp\new_node.dbf" dbms=dbf replace;
 

update_split_itinerary_Aug2012.sas

@@ -2,6 +2,8 @@
     Craig Heither, revised 08/07/2012
 
    Program updates the values of itin_a and itin_b in the itinerary coding for rail routes whose links have been split.
+   
+   NRF 6/2/2017 - added logic to set dwell code and dwell time in itin as no stop if node ID represents junction.
 
   */
 
@@ -38,6 +40,10 @@ data a(drop=c newc); set a;
 data itin(drop=newa newb); merge x (in=hit) a; by itin_a itin_b; if hit;
   if newa>0 then itin_a=newa;
   if newb>0 then itin_b=newb;
+  if not(index(tr_line,'*')) and (39000<itin_b<40000 or 49000<itin_b<50000) then do;
+      dw_code=1;
+      dw_time=0;
+      end;
    proc sort; by tr_line it_order;
 
    *** RE-ORDER VARIABLES ***;