Merge pull request #27 from nrferguson/master

batchin creation updates

Author: nrferguson

Master Rail Network Management.tbx

@@ -13,12 +13,14 @@
                     types if desired.
         05-09-2012: Logic revised to process updated future route coding procedures.
         06-19-2013: Added logic to choose people mover links.
-		08-07-2014: Added check for stops at junctions before itinerary batchin creation
+        08-07-2014: Added check for stops at junctions before itinerary batchin creation
         01-08-2016: Corrected itinerary batchin formatting
         02-01-2016: Metra headways coded as '99' are replaced by the TOD period duration.
         12-20-2016: Updated 'basescen' variable to '200' after base year moved to 2015.
         NRF 8/23/2017: Added support for updated future coding of additional service (action = 1)
                        to include appropriate additional service in all TOD periods.
+        12/9/2019 NRF: hardcoded OHare Express layover in TOD itinerary file.
+        12/10/2019 NRF: updated itinerary file formatting for case where there is only one itinerary segment.
 
 -------------                                                             -------------
 __________________________________________________________________________________________________________________________  */
@@ -29,7 +31,7 @@ ________________________________________________________________________________
 %let ct_ramp=%scan(&sysparm,4,$);
 %let scen=%eval(&sc/100);
 %let reportpath=&outpath.\&scen.00\rail_changes.txt;
-%let maxzone=1961;                                                        *** highest zn09 POE zone number;
+%let maxzone=3649;                                                        *** highest zone17 POE zone number;
 %let dropnode=48007;                                                      *** Exclude South Bend station until zone system expands to include it;
 %let basescen=100;                                                        *** base year scenario;
 %let counter=1;
@@ -50,7 +52,7 @@ data __null__;
        *** READ IN NODE INFORMATION ***;
   proc import datafile="&inpath.\Temp\temp_rlnode_zone.dbf" out=nodes replace; 
   data nodes(keep=node label pspace pcost ftr_pspace ftr_pcost point_x point_y zone); set nodes; 
-    rename zone09=zone;
+    rename zone17=zone;
 
       *** READ IN LINK INFORMATION ***;
   proc import datafile="&inpath.\Temp\temp_arc.dbf" out=network replace;
@@ -233,14 +235,15 @@ run;
       set itins(where=(itin_a not in (&dropnode) & itin_b not in (&dropnode)));
       by tr_line it_order;
       retain o 1;
-	  it_order=o;
-	  o=o+1;
-	  if first.tr_line then do;
+      it_order=o;
+      o=o+1;
+      if first.tr_line then do;
           o=1;
-		  it_order=o;
-		  o=o+1;
-		  end;
+          it_order=o;
+          o=o+1;
+          end;
       if last.tr_line then layover='3';
+      if (last.tr_line & substr(tr_line, 1, 3) = 'oes') then layover='0';    *** hardcoded layover for OHare Express - needs better solution ***;
 
 
   data combine; set routes itins; proc sort; by tr_line it_order;
@@ -265,15 +268,17 @@ run;
       else if dw_code=5 then d=compress('dwt=*'||dw_time);
       else d=compress('dwt='||dw_time);
       lag_d=lag(d);
+      
 
      file out1;
      if _n_=1 then put "c RAIL TRANSIT BATCHIN FILE FOR SCENARIO &scen.00 TOD &tod" / "c &sysdate" / "c us1 holds segment travel time, us2 holds zone fare" / "t lines init";
      if first.tr_line then put 'a' +1 name +2 mode +2 veh_type +2 headway +2 speed +2 desc / +2 'path=no';
-     else if last.tr_line then put +4 lag_d +(10-length(left(trim(lag_d)))) itin_a +2 'ttf=1' +3 'us1=' +0 trv_time +(6-length(left(trim(trv_time)))) 'us2=' +0 zn_fare / +4 'dwt=0.01' +3 itin_b +2 'lay=' +0 layover;
+     else if (it_order=1 & last.tr_line) then put +4 'dwt=0.01' +3 itin_a +2 'ttf=1' +3 'us1=' +0 trv_time +(6-length(left(trim(trv_time)))) 'us2=' +0 zn_fare / +4 'dwt=0.01' +3 itin_b +2 'lay=' +0 layover;
      else if it_order=1 then put +4 'dwt=0.01' +3 itin_a +2 'ttf=1' +3 'us1=' +0 trv_time +(6-length(left(trim(trv_time)))) 'us2=' +0 zn_fare;
+     else if last.tr_line then put +4 lag_d +(10-length(left(trim(lag_d)))) itin_a +2 'ttf=1' +3 'us1=' +0 trv_time +(6-length(left(trim(trv_time)))) 'us2=' +0 zn_fare / +4 'dwt=0.01' +3 itin_b +2 'lay=' +0 layover;
      else if (layov ne '0' and layov ne '') then put +4 lag_d +(10-length(left(trim(lag_d)))) itin_a +2 'ttf=1' +3 'us1=' +0 trv_time +(6-length(left(trim(trv_time)))) 'us2=' +0 zn_fare +2 'lay=' +0 layov;
      else put +4 lag_d +(10-length(left(trim(lag_d)))) itin_a +2 'ttf=1' +3 'us1=' +0 trv_time +(6-length(left(trim(trv_time)))) 'us2=' +0 zn_fare;
-
+    
      %if &scen=9 %then %do;
          proc export data=combine outfile="&outpath.\&scen.00\rail_itinerary_&tod..dbf" dbms=dbf replace;
      %end; 

create_Emme_rail_files_GTFS.sas

@@ -84,7 +84,7 @@
 people_mover = d + "\\mrn.gdb\\" + "people_mover"
 scen_itin_dbf = d + "\\Temp\\scen_itin.dbf"
 ftr_itin_dbf = d + "\\Temp\\ftr_itin.dbf"
-zones = d + "\\zones_2009\\polygon"
+zones = d + "\\zones17\\zones17.shp"
 temp_rlnode_zone_shp = Temp + "\\temp_rlnode_zone.shp"
 outFl = Temp + "\\rte_out.txt"
 link_list = newdir + '\\rail_links_all.csv'