RSGInc · americalexander · Jan 16, 2026 · Jan 20, 2026 · Jan 26, 2026 · Jan 26, 2026
diff --git a/resident/configs/accessibility.yaml b/resident/configs/accessibility.yaml
@@ -1,2 +1,3 @@
 # columns from land_use table to add to df
 land_use_columns: ['EMP_RET', 'EMP_TOTAL', 'TOTHHS', 'walk_dist_local_bus']
+explicit_chunk: 0.2
diff --git a/resident/configs/annotate_disaggregate_accessibility.csv b/resident/configs/annotate_disaggregate_accessibility.csv
@@ -1,4 +1,4 @@
-# annotating the proto_disaggregate_accessibilty table
+# annotating the proto_disaggregate_accessibilty table,,
 Description,Target,Expression
 workplace location for zero auto only,workplace_location_accessibility_0,"np.where(df.auto_ownership == 0, df.workplace_location_accessibility, np.nan)"
 workplace location for auto deficient,workplace_location_accessibility_1,"np.where(df.auto_ownership == 1, df.workplace_location_accessibility, np.nan)"

diff --git a/resident/configs/annotate_households.csv b/resident/configs/annotate_households.csv
@@ -31,4 +31,5 @@ home_is_rural,home_is_rural,"reindex(land_use.pseudomsa, households.home_zone_id
 ,num_hh_in_zone,"reindex(households.groupby('home_zone_id').size(), households.home_zone_id)"
 #,multiple_auto_hh_in_zone,"reindex_i(df.groupby('home_zone_id').multiple_auto_hh.sum(), df.index)"
 #,multiple_auto_hh_in_zone_share,multiple_auto_hh_in_zone/num_hh_in_zone
-hh owns an ebike,ebike_owner,"rng.random_for_df(households)[:,0] < ebikeownership"
+hh owns an ebike,ebike_owner,"rng.random_for_df(households)[:,0] < ebikeownership"
+Retired Adults Only Households,retired_adults_only_hh,(households.hhsize > 0) & (households.hhsize == num_retired_adults)
diff --git a/resident/configs/annotate_landuse.csv b/resident/configs/annotate_landuse.csv
@@ -4,10 +4,10 @@ household_density,household_density,land_use.TOTHHS / land_use.ACRES.clip(lower=
 population_density,population_density,land_use['TOTPOP'] / land_use.ACRES.clip(lower=1)
 employment_density,employment_density,(land_use['EMP_TOTAL'] / land_use.ACRES.clip(lower=1)).clip(upper=500)
 density_index,density_index,(household_density * employment_density) / (household_density + employment_density).clip(lower=1)
+effective active acres,ACTIVE_ACRES,land_use.ACRES * land_use.PARKATTRACT
 # FIXME pseudomsa is used in a few places but there is no equivalent in Metro landuse,,
 ,pseudomsa,0
 # FIXME missing in landuse but used in SANDAG -- setting to 0 for now,,
-,ACTIVE_ACRES,0
 ,micro_dist_local_bus,0
 ,microtransit,0
 ,nev,0

diff --git a/resident/configs/annotate_persons.csv b/resident/configs/annotate_persons.csv
@@ -11,14 +11,14 @@ employment status type,pemploy,"pd.Series(np.zeros(len(persons)), index=persons.
 ,pemploy,"np.where((persons.age >= 16) & ((persons.ESR == 3) | (persons.ESR == 6)), PEMPLOY_NOT, pemploy)"
 ,pemploy,"np.where((persons.age>=16) & (persons.ESR != 3) & (persons.ESR != 6) & (persons.WKHP >= 35) & (persons.WKW >= 1) & (persons.WKW <= 4), PEMPLOY_FULL, pemploy)"
 student category,pstudent,"pd.Series(np.zeros(len(persons)), index=persons.index)"
-,pstudent,"np.where((pemploy == 1) & (persons.age >= 16), PSTUDENT_NOT, pstudent)"
-,pstudent,"np.where((pemploy == 1) & (persons.age < 16), PSTUDENT_GRADE_OR_HIGH, pstudent)"
+,pstudent,"np.where((pemploy == PEMPLOY_FULL) & (persons.age >= 16), PSTUDENT_NOT, pstudent)"
+,pstudent,"np.where((pemploy == PEMPLOY_FULL) & (persons.age < 16), PSTUDENT_GRADE_OR_HIGH, pstudent)"
 ,pstudent,"np.where((persons.SCHG < 1) & (persons.age >= 16), PSTUDENT_NOT, pstudent)"
 ,pstudent,"np.where((persons.SCHG < 1) & (persons.age < 16), PSTUDENT_GRADE_OR_HIGH, pstudent)"
-,pstudent,"np.where((persons.SCHG >= 15) & (persons.age >= 16) & (pemploy != 1), PSTUDENT_UNIVERSITY, pstudent)"
-,pstudent,"np.where((persons.SCHG >= 15) & (persons.age < 16) & (pemploy != 1), PSTUDENT_GRADE_OR_HIGH, pstudent)"
-,pstudent,"np.where((persons.age <= 19) & (pemploy != 1) & (persons.SCHG >=1) & (persons.SCHG<=14), PSTUDENT_GRADE_OR_HIGH, pstudent)"
-,pstudent,"np.where((persons.age > 19) & (pemploy != 1) & (persons.SCHG >=1) & (persons.SCHG<=14),  PSTUDENT_UNIVERSITY, pstudent)"
+,pstudent,"np.where((persons.SCHG >= 15) & (persons.age >= 16), PSTUDENT_UNIVERSITY, pstudent)"
+,pstudent,"np.where((persons.SCHG >= 15) & (persons.age < 16), PSTUDENT_GRADE_OR_HIGH, pstudent)"
+,pstudent,"np.where((persons.age <= 19) & (persons.SCHG >=1) & (persons.SCHG<=14), PSTUDENT_GRADE_OR_HIGH, pstudent)"
+,pstudent,"np.where((persons.age > 19) & (persons.SCHG >=1) & (persons.SCHG<=14),  PSTUDENT_UNIVERSITY, pstudent)"
 ,pstudent,"np.where(pstudent == 0, 3, pstudent)"
 person type,ptype,"pd.Series(PTYPE_NONWORK, index=persons.index)"
 ,ptype,"np.where(persons.age >= 65, PTYPE_RETIRED, ptype)"
@@ -62,4 +62,4 @@ travel time sensitivity factor for non-work travel,time_factor_nonwork,"np.minim
 ,naics_code,_naics_code.astype('int') if all(_naics_code!=0) else 0
 #,,
 # FIXME does filling with -1 make sense here? depends on what the missing values mean,,
-occupation categories mapped to workplace segments,occupation,persons.OCCP.fillna(-1).map(occupation_xwalk)
+occupation categories mapped to workplace segments,occupation,persons.OCCP.fillna(999).map(occupation_xwalk)
diff --git a/resident/configs/auto_ownership.csv b/resident/configs/auto_ownership.csv
@@ -14,7 +14,7 @@ util_hh_income_verylow,household income <15k,@(df.income<15000),coef_household_i
 util_hh_income_low,household income 15k-30k,@((df.income>=15000) & (df.income<30000)),coef_household_income_1530k_0_CARS,,coef_household_income_1530k_2_CARS,coef_household_income_1530k_3_CARS,coef_household_income_1530k_4_CARS
 util_hh_income_mid,household income 30-60k,@((df.income>=30000) & (df.income<60000)),coef_household_income_3060k_0_CARS,,coef_household_income_3060k_2_CARS,coef_household_income_3060k_3_CARS,coef_household_income_3060k_4_CARS
 util_hh_income_veryhigh,household income 100k+,@(df.income>=100000),coef_household_income_100k_0_CARS,,coef_household_income_100k_2_CARS,coef_household_income_100k_3_CARS,coef_household_income_100k_4_CARS
-# FIXME no bldgsz equivalent in the Metro data,,,,,,
+# FIXME no bldgsz equivalent in the Metro data,,,,,,,
 # util_attached,attached dwelling,"@np.where(df.bldgsz == 2, 0, 1)",coef_attached_0_CARS,,coef_attached_2_CARS,coef_attached_3_CARS,coef_attached_4_CARS
 # 0.785 sq miles in 1/2 mile radius * 640 acres per square mile = 502,,,,,,,
 util_intersection_density,Intersection count in 1/2 mile radius of household,totint/502,coef_intersection_density_0_CARS,,coef_intersection_density_2_CARS,coef_intersection_density_3_CARS,coef_intersection_density_4_CARS

diff --git a/resident/configs/auto_ownership_preprocessor.csv b/resident/configs/auto_ownership_preprocessor.csv
@@ -1,7 +1,7 @@
 Description,Target,Expression
 ,_PERSON_COUNT,"lambda query, persons, households: persons.query(query).groupby('household_id').size().reindex(households.index).fillna(0).astype(np.int8)"
 ,num_old_retirees,"_PERSON_COUNT('80 <= age', persons, households)"
-# averaging disaggregate accessibilities across all household members
+# averaging disaggregate accessibilities across all household members,,
 ,shopping_accessibility_0,persons_merged.groupby('household_id').shopping_accessibility_0.mean().reindex(df.index)
 ,shopping_accessibility_1,persons_merged.groupby('household_id').shopping_accessibility_1.mean().reindex(df.index)
 ,shopping_accessibility_2,persons_merged.groupby('household_id').shopping_accessibility_2.mean().reindex(df.index)

diff --git a/resident/configs/constants.yaml b/resident/configs/constants.yaml
@@ -99,7 +99,6 @@ occupation_xwalk:
   5: nat_res_cnstr_maint
   6: prod_trans_move
   999: missing # no occupation 
-  -1: missing # no occupation
 
 # RIDEHAIL Settings - all costs in dollars, times in minutes
 Taxi_baseFare: 3.00
@@ -355,4 +354,4 @@ OpVTCr: 0.1442
 
 # model region constants
 KNR_available: True  # SKATS does not have KNR, so set to False in configs_skats/constants.yaml
-model_region_name: "Metro"
+model_region_name: "Metro"
diff --git a/resident/configs/disaggregate_accessibility.yaml b/resident/configs/disaggregate_accessibility.yaml
@@ -94,7 +94,7 @@ CREATE_TABLES:
                 # FIXME this should be updated when occupation_xwalk is updated
                 OCCP: # occupation codes
                     1: 1   # Management Occupations
-                    2: -1    # Non-worker
+                    2: 999    # Non-worker
 
     PROTO_TOURS:
         index_col: proto_tour_id

diff --git a/resident/configs/external_non_mandatory_identification.csv b/resident/configs/external_non_mandatory_identification.csv
@@ -1,6 +1,6 @@
 Label,Description,Expression,external_tour,internal_tour
 util_dist_to_nearest_ext_station,Distance to nearest external station,dist_to_external_zone,coef_dist_to_nearest_ext_station,
-# FIXME need actual external counts at external station - fix in destination choice size terms too
+# FIXME need actual external counts at external station - fix in destination choice size terms too,,,,
 # util_log_size,Log size of nearest station,"@np.log1p(reindex(land_use.external_nonwork, df.closest_external_zone))",coef_log_size_of_nearest_ext_station,
 util_log_size,Log size of nearest station,"@np.log1p(reindex(land_use.EXTERNAL, df.closest_external_zone))",coef_log_size_of_nearest_ext_station,
 util_escort,escort tour ASC,"@np.where(df.tour_type == 'escort', 1, 0)",coef_escort,

diff --git a/resident/configs/external_non_mandatory_identification.yaml b/resident/configs/external_non_mandatory_identification.yaml
@@ -9,3 +9,5 @@ LOGIT_TYPE: MNL
 EXTERNAL_COL_NAME: is_external_tour
 # set to True if not external but CHOOSER_FILTER_COLUMN_NAME is True
 INTERNAL_COL_NAME: is_internal_tour
+
+EXTERNAL_TOUR_ALT: 1
diff --git a/resident/configs/external_worker_identification.csv b/resident/configs/external_worker_identification.csv
@@ -1,16 +1,16 @@
-Label,Description,Expression,work_external,work_internal
-util_dist_to_nearest_ext_station,Distance to nearest external station,dist_to_external_zone,coef_dist_to_nearest_ext_station,
+Label,Description,Expression,work_internal,work_external
+util_dist_to_nearest_ext_station,Distance to nearest external station,dist_to_external_zone,,coef_dist_to_nearest_ext_station
 # FIXME just have a single external indicator not counts at the station here and in destination choice size terms,,,,
-#util_size_of_nearest_ext_station,Size of nearest external station,"@np.log1p(reindex(land_use.external_work,df.closest_external_zone))",coef_size_of_nearest_ext_station,
-util_size_of_nearest_ext_station,Size of nearest external station,"@np.log1p(reindex(land_use.EXTERNAL,df.closest_external_zone))",coef_size_of_nearest_ext_station,
-util_dist_lt_2p5,Distance less than 2.5 miles,"@np.where(df.dist_to_external_zone<2.5,1,0)",coef_dist_lt_2p5,
-util_part_time,Part time worker,"@np.where(df.pemploy == 2,1,0)",coef_part_time,
-util_inc_lt15,Household Income less than $15k,@(df.income<15000),coef_inc_lt15,
-util_inc_15_25,Household income $15k-$25k,@(df.income>=15000) & (df.income<25000) ,coef_inc_15_25,
-util_inc_25_50,Household Income $25k-$50k,@(df.income>=25000) & (df.income<50000) ,coef_inc_25_50,
-util_inc_150_250,Household Income $150k-$50k,@(df.income>=150000) * (df.income<250000) ,coef_inc_150_250,
-util_inc_250plus,Household Income $250k+,@(df.income>=250000),coef_inc_250plus,
-util_asc,Alternative-specific constant for external worker,1,asc_external_worker,
-util_global_switch,Global switch to have everything internal,@NO_EXTERNAL,-999,
-util_2016,Constant for pre-COVID conditions,@PRE_COVID,asc_external_2016,
-util_calib,Constant for calibration,1,0.35,
+#util_size_of_nearest_ext_station,Size of nearest external station,"@np.log1p(reindex(land_use.external_work,df.closest_external_zone))",,coef_size_of_nearest_ext_station
+util_size_of_nearest_ext_station,Size of nearest external station,"@np.log1p(reindex(land_use.EXTERNAL,df.closest_external_zone))",,coef_size_of_nearest_ext_station
+util_dist_lt_2p5,Distance less than 2.5 miles,"@np.where(df.dist_to_external_zone<2.5,1,0)",,coef_dist_lt_2p5
+util_part_time,Part time worker,"@np.where(df.pemploy == 2,1,0)",,coef_part_time
+util_inc_lt15,Household Income less than $15k,@(df.income<15000),,coef_inc_lt15
+util_inc_15_25,Household income $15k-$25k,@(df.income>=15000) & (df.income<25000) ,,coef_inc_15_25
+util_inc_25_50,Household Income $25k-$50k,@(df.income>=25000) & (df.income<50000) ,,coef_inc_25_50
+util_inc_150_250,Household Income $150k-$50k,@(df.income>=150000) * (df.income<250000) ,,coef_inc_150_250
+util_inc_250plus,Household Income $250k+,@(df.income>=250000),,coef_inc_250plus
+util_asc,Alternative-specific constant for external worker,1,,asc_external_worker
+util_global_switch,Global switch to have everything internal,@NO_EXTERNAL,,-999
+util_2016,Constant for pre-COVID conditions,@PRE_COVID,,asc_external_2016
+util_calib,Constant for calibration,1,,0.35
diff --git a/resident/configs/external_worker_identification.yaml b/resident/configs/external_worker_identification.yaml
@@ -4,6 +4,8 @@ COEFFICIENTS: external_worker_identification_coeffs.csv
 
 LOGIT_TYPE: MNL
 
+EXTERNAL_WORKER_ALT: 1
+
 # boolean column to filter choosers (True means keep)
 # will only expose these people to the model
 CHOOSER_FILTER_COLUMN_NAME: is_out_of_home_worker

diff --git a/resident/configs/joint_tour_frequency_composition.csv b/resident/configs/joint_tour_frequency_composition.csv
@@ -21,7 +21,7 @@ util_constant_for_2_visiting_tour,Constant for 2 visiting tour,@df.social==2,coe
 util_constant_for_1_visiting_tour,Constant for 1 visiting tour,@df.social==2,coef_constant_for_1_visiting_tour
 util_1_visiting_and_1_discretionary_tour,1 Visiting and 1 Discretionary Tour,@((df.social==1) & (df.othdiscr==1)),coef_1_visiting_and_1_discretionary_tour
 util_constant_for_2_discretionary_tour,Constant for 2 discretionary tour,othdiscr==2,coef_constant_for_2_discretionary_tour
-util_constant_for_2_discretionary_tour,Constant for 2 discretionary tour,othdiscr==1,coef_constant_for_1_discretionary_tour
+util_constant_for_1_discretionary_tour,Constant for 1 discretionary tour,othdiscr==1,coef_constant_for_1_discretionary_tour
 util_number_of_active_full_time_workers_shopping,Number of Active Full time workers /Shopping,num_travel_active_full_time_workers * shopping,coef_number_of_active_full_time_workers_shopping
 util_number_of_active_nonworkers_shopping,Number of Active Non-workers /Shopping,num_travel_active_non_workers * shopping,coef_number_of_active_nonworkers_shopping
 util_number_of_active_pre_driving_age_school_children_shopping,Number of Active Pre- Driving Age School Children /Shopping,num_travel_active_pre_driving_age_school_kids * shopping,coef_number_of_active_pre_driving_age_school_children_shopping

diff --git a/resident/configs/license_holding_status.csv b/resident/configs/license_holding_status.csv
@@ -1,5 +1,5 @@
 Label,Description,Expression,has_license,does_not_have_license
-# dummy specification for license holding status model,,,
+# dummy specification for license holding status model,,,,
 util_constant,alternative specific constant,1,coef_has_license_ASC,0
 util_unavailable_for_under_16,unavailable for persons under 16 years old,@df.age < 16,coef_unavailable,0
 util_age,Age of person in years,@df.age,coef_age,0

diff --git a/resident/configs/logging.yaml b/resident/configs/logging.yaml
@@ -46,12 +46,12 @@ logging:
       class: logging.Formatter
       # format: '%(levelname)s - %(name)s - %(message)s'
       format: '%(levelname)s - %(message)s'
-      datefmt: '%d/%m/%Y %H:%M:%S'
+      datefmt: '%Y-%m-%dT%H:%M:%S%z'
 
     fileFormatter:
       class: logging.Formatter
       format: '%(asctime)s - %(levelname)s - %(name)s - %(message)s'
-      datefmt: '%d/%m/%Y %H:%M:%S'
+      datefmt: '%Y-%m-%dT%H:%M:%S%z'
 
     elapsedFormatter:
       (): activitysim.core.tracing.ElapsedTimeFormatter

diff --git a/resident/configs/mandatory_tour_scheduling_annotate_alts_preprocessor.csv b/resident/configs/mandatory_tour_scheduling_annotate_alts_preprocessor.csv
@@ -6,7 +6,7 @@ duration_reference_bin,durationRefBin,20
 departure_shift,departureLinearShift1,"(9-df.start)*(df.start<=9) + (df.start-9)*(df.start>9)"
 arrival_shift,arrivalLinearShift1,"(29-df.end)*(df.end<=29) + (df.end-29)*(df.end>29)"
 duration_shift,durationShift,"(20-df.duration)*(df.duration<=20) + (df.duration-20)*(df.duration>20)"
-# school specific bins
+# school specific bins,,
 departure_reference_bin_school,departureRefBin_school,9
 arrival_reference_bin_school,arrivalRefBin_school,29
 duration_reference_bin_school,durationRefBin_school,16

diff --git a/resident/configs/non_mandatory_tour_destination.csv b/resident/configs/non_mandatory_tour_destination.csv
@@ -71,4 +71,4 @@ util_joint_0_2_ASC,joint tours in 0_2bin calibration constant,"@(_DIST.between(0
 util_joint_2_5_ASC,joint tours in 2_5bin calibration constant,"@(_DIST.between(2,5)) * (df.get('tour_category', default=False) == 'joint')",,coef_abm3_dist_2_5jointmaint_asc,coef_abm3_dist_2_5jointdisc_asc,coef_abm3_dist_2_5jointmaint_asc,coef_abm3_dist_2_5jointdisc_asc,coef_abm3_dist_2_5jointdisc_asc
 util_joint_5_10_ASC,joint tours in 5_10bin calibration constant,"@(_DIST.between(5,10)) * (df.get('tour_category', default=False) == 'joint')",,coef_abm3_dist_5_10jointmaint_asc,coef_abm3_dist_5_10jointdisc_asc,coef_abm3_dist_5_10jointmaint_asc,coef_abm3_dist_5_10jointdisc_asc,coef_abm3_dist_5_10jointdisc_asc
 util_joint_10_30_ASC,joint tours in 10_30bin calibration constant,"@(_DIST.between(10,30)) * (df.get('tour_category', default=False) == 'joint')",,coef_abm3_dist_10_30jointmaint_asc,coef_abm3_dist_10_30jointdisc_asc,coef_abm3_dist_10_30jointmaint_asc,coef_abm3_dist_10_30jointdisc_asc,coef_abm3_dist_10_30jointdisc_asc
-util_indiv_dist_ASC,indiv tours in distance calibration constant,"@(_DIST) * (df.get('tour_category', default=False) == 'joint')",,0.10,,0.10,,
+util_joint_dist_ASC,indiv tours in distance calibration constant,"@(_DIST) * (df.get('tour_category', default=False) == 'joint')",,0.10,,0.10,,
diff --git a/resident/configs/non_mandatory_tour_destination.yaml b/resident/configs/non_mandatory_tour_destination.yaml
@@ -3,6 +3,7 @@ SPEC: non_mandatory_tour_destination.csv
 COEFFICIENTS: non_mandatory_tour_destination_coefficients.csv
 
 SAMPLE_SIZE: 30
+ESTIMATION_SAMPLE_SIZE: 30
 
 SIZE_TERM_SELECTOR: non_mandatory
 

diff --git a/resident/configs/non_mandatory_tour_frequency_alternatives.csv b/resident/configs/non_mandatory_tour_frequency_alternatives.csv
@@ -201,4 +201,5 @@ escort,shopping,othmaint,eatout,social,othdiscr
 2,1,1,0,1,1
 2,1,2,1,0,0
 2,0,2,1,1,0
+2,0,2,0,1,1
 2,1,1,1,1,0
diff --git a/resident/configs/non_mandatory_tour_scheduling.yaml b/resident/configs/non_mandatory_tour_scheduling.yaml
@@ -30,6 +30,7 @@ SIMULATE_CHOOSER_COLUMNS:
   - num_non_escort_tours
   - num_add_shop_maint_tours
   - num_add_soc_discr_tours
+  - retired_adults_only_hh
 
 LOGSUM_SETTINGS: tour_mode_choice.yaml
 

diff --git a/resident/configs/school_escorting.yaml b/resident/configs/school_escorting.yaml
@@ -12,6 +12,8 @@ ALTS: school_escorting_alts.csv
 
 LOGIT_TYPE: MNL
 
+ESCORTEE_AGE_CUTOFF: 19
+
 NUM_ESCORTEES: 3
 NUM_CHAPERONES: 2
 

diff --git a/resident/configs/school_escorting_coefficients_inbound.csv b/resident/configs/school_escorting_coefficients_inbound.csv
@@ -6,10 +6,10 @@ coef_child_age_u5_noes,-1.36718,F
 coef_ln_dist_from_school_noes,-0.01787,F
 coef_ln_dist_from_school_u6_noes,-0.23304,F
 coef_ln_dist_from_school_6to9_noes,-0.07286,F
-coef_child_age_16p_rs,1.97184
-coef_child_age_10to15_rs,1.73544
-coef_child_age_6to9_rs,1.73544
-coef_child_age_u6_rs,1.34996
+coef_child_age_16p_rs,1.97184,
+coef_child_age_10to15_rs,1.73544,
+coef_child_age_6to9_rs,1.73544,
+coef_child_age_u6_rs,1.34996,
 coef_hh_inc_u25k_noes,0.0,F
 coef_hh_inc_25to50k_noes,0.0,F
 coef_zero_auto_hh_noes,0.13165,F

diff --git a/resident/configs/school_escorting_preprocessor_outbound.csv b/resident/configs/school_escorting_preprocessor_outbound.csv
@@ -12,13 +12,13 @@ Daily activity pattern - chauffer 2,cdap_chauf2,"reindex(persons.cdap_activity,
 Age - child 1,age_child1,"reindex(persons.age, df.child_id1)"
 Age - child 2,age_child2,"reindex(persons.age, df.child_id2)"
 Age - child 3,age_child3,"reindex(persons.age, df.child_id3)"
-# Departure times to school and work
+# Departure times to school and work,,
 Preferred departure time to school - child 1,pref_depart_time_school1,"reindex(tours[(tours.tour_type == 'school') & (tours.tour_num == 1)].set_index('person_id').start, df.child_id1)"
 Preferred departure time to school - child 2,pref_depart_time_school2,"reindex(tours[(tours.tour_type == 'school') & (tours.tour_num == 1)].set_index('person_id').start, df.child_id2)"
 Preferred departure time to school - child 3,pref_depart_time_school3,"reindex(tours[(tours.tour_type == 'school') & (tours.tour_num == 1)].set_index('person_id').start, df.child_id3)"
 Preferred departure time to work / univ - chauffer 1,pref_depart_time_chauf1,"reindex(tours[(tours.tour_category == 'mandatory') & (tours.tour_num == 1)].set_index('person_id').start, df.chauf_id1)"
 Preferred departure time to work / univ - chauffer 2,pref_depart_time_chauf2,"reindex(tours[(tours.tour_category == 'mandatory') & (tours.tour_num == 1)].set_index('person_id').start, df.chauf_id2)"
-# Distances and times to school and work
+# Distances and times to school and work,,
 School location - child 1,school_location_child1,"reindex(persons.school_zone_id, df.child_id1)"
 School location - child 2,school_location_child2,"reindex(persons.school_zone_id, df.child_id2)"
 School location - child 3,school_location_child3,"reindex(persons.school_zone_id, df.child_id3)"