R/MCTTRP_mod_compartments_creation.R

In LauraDavilaPena/cwmcttrp: Clarke and Wright Multi-Compartment Truck and Trailer Routing Problem (CW-MCTTRP)

boolean_available_compartments_destination_route <- function(input, result, intermediate_solution, inserted_client, 
                                                             destination_route, initial_solution, penalty_max){
  
  #resultado <- intermediate_solution
  
  
  total_truck_hoppers <- length(input$H.camion[1,])
  total_trailer_hoppers <- length(input$H.trailer[1,])
  
  capacity_truck_hoppers <- input$H.camion[1,1]
  capacity_trailer_hoppers <- input$H.trailer[1,1]
  
  # penalizacion
  
  
  n1 <- input$n1
  
  
  avail <- FALSE
  avail_truck <- FALSE
  
  
  client_origin_route_index <- route_of_client(inserted_client, initial_solution)$index
  client_origin_route_info <- initial_solution[[client_origin_route_index]]
  
  
  if(inserted_client <= n1 ){ #vc
    # Habia un fallo. 27/04/2021
    for(ind_vc in 1:length(client_origin_route_info$clients_vc)){
      if(client_origin_route_info$clients_vc[[ind_vc]]$id == inserted_client ){
        client_vc_index <- ind_vc
      }
    }
    # hasta aqui (el client_vc_index no estaba bien, ya que lo calculabamos como en la linea de abajo)
    VC_index <- which(client_origin_route_info$VCs==inserted_client)
    client_vc_info <- client_origin_route_info$clients_vc[[client_vc_index]]
    client_hoppers_trailers <- client_vc_info$hoppers_trailers
    client_hoppers_trucks <- client_vc_info$hoppers_trucks
    client_hoppers_demands <- client_vc_info$demands
    
    
    if(destination_route$type == "PTR"){
      needed_truck_hoppers_per_product <- ceiling(client_hoppers_demands/capacity_truck_hoppers)
      needed_truck_hoppers_total <- sum(needed_truck_hoppers_per_product)
      # Fallo 27/04/2021. La penalty_max estaba mal colocada en el siguiente condicional
      if(needed_truck_hoppers_total <=  total_truck_hoppers - destination_route$used_hoppers_truck + floor(penalty_max/capacity_truck_hoppers)){
        
        avail <- TRUE
        
        
      }
      
    }else{ #if destination_route is either CVR or PVR
      
      # Fallo 27/04/2021. Puse <= en vez de ==
      if( total_trailer_hoppers - destination_route$used_hoppers_trailer <= 0){ # no hay tolvas de trailer suficientes: chequeo trucks
        needed_truck_hoppers_per_product <- ceiling(client_hoppers_demands/capacity_truck_hoppers)
        needed_truck_hoppers_total <- sum(needed_truck_hoppers_per_product)
        #  Fallo 27/04/2021. La penalty_max estaba mal colocada en el siguiente condicional
        if(needed_truck_hoppers_total <=  total_truck_hoppers - destination_route$used_hoppers_truck + floor(penalty_max/capacity_truck_hoppers)){ # me caben
          
          avail <- TRUE 
          avail_truck <- TRUE
          
        }
      }else{ # en caso de que si haya tolvas de trailer disponibles (en destination_route) hay que ver si me llegan (junto con las que haya disponibles de truck)
        
        # 28/04/2021. Para ver si toda su demanda cabria en el camion
        needed_truck_hoppers_per_product <- ceiling(client_hoppers_demands/capacity_truck_hoppers)
        needed_truck_hoppers_total <- sum(needed_truck_hoppers_per_product)
        if(needed_truck_hoppers_total <=  total_truck_hoppers - destination_route$used_hoppers_truck + floor(penalty_max/capacity_truck_hoppers)){
          avail_truck <- TRUE
        }
        # hasta aqui lo añadido
        
        needed_trailer_hoppers_per_product <- ceiling(client_hoppers_demands/capacity_trailer_hoppers)
        needed_trailer_hoppers_total <- sum(needed_trailer_hoppers_per_product)
        
       
        #  Fallo 27/04/2021. La penalty_max estaba mal colocada en el siguiente condicional
        if(needed_trailer_hoppers_total <=  total_trailer_hoppers - destination_route$used_hoppers_trailer + floor(penalty_max/capacity_trailer_hoppers)){ # me caben
          
          avail <- TRUE
          
          
        }else{
          # Fallo 27/04/2021. Añadi penalty_max como input de la siguiente funcion
            checking_compartments_available <- is_there_available_compartments(total_truck_hoppers, total_trailer_hoppers, destination_route, 
                                                                               client_hoppers_demands, needed_trailer_hoppers_per_product, 
                                                                               capacity_truck_hoppers, capacity_trailer_hoppers, penalty_max)
            
            available_compartments <- checking_compartments_available$available_compartments
            
            if(available_compartments){
              needed_truck_hoppers_per_product <- checking_compartments_available$needed_truck_hoppers_per_product
              needed_truck_hoppers_total <- sum(needed_truck_hoppers_per_product)
              
              
              avail <- TRUE 
              
            }
            
           
        }
        
        
      }
    }
    
    
  }else if(inserted_client > n1){ #tc
    # Fallo 27/04/2021. El indice estaba mal calculado
    for(ind_tc in 1:length(client_origin_route_info$clients_tc)){
      if(client_origin_route_info$clients_tc[[ind_tc]]$id == inserted_client ){
        client_tc_index <- ind_tc
      }
    }
    TC_index <- which(client_origin_route_info$TCs==inserted_client)
    client_tc_info <- client_origin_route_info$clients_tc[[client_tc_index]]
    client_hoppers_trucks <- client_tc_info$hoppers_trucks
    client_hoppers_demands <- client_tc_info$demands
    
    
    #Como el cliente es TC, la "destination_route" si o si sera una PTR o una subruta
    # (en ambos casos toda su demanda tiene que ir en truck)
    needed_hoppers_per_product <- ceiling(client_hoppers_demands/capacity_truck_hoppers)
    needed_hoppers_total <- sum(needed_hoppers_per_product)
    
    #  Fallo 27/04/2021. La penalty_max estaba mal colocada en el siguiente condicional
    if(needed_hoppers_total <=  total_truck_hoppers - destination_route$used_hoppers_truck + floor(penalty_max/capacity_truck_hoppers)){
      
      avail <- TRUE
      avail_truck <- TRUE
      
    }
    
    
    
    
  }
  
  
  
  
  return(list(avail=avail, avail_truck=avail_truck))
}




create_only_one_vehicle_type_hoppers_for_vc <- function(inserted_client, destination_route, client_hoppers_the_other_vehicles_type, 
                                                        client_hoppers_vehicles_type, capacity_vehicle_type_hoppers,
                                                        client_vc_info, client_hoppers_demands, 
                                                        needed_vehicle_type_hoppers_per_product, used_hoppers_vehicle_type,
                                                        hoppers_the_other_vehicles_type, hoppers_vehicles_type){
  
  new_clients_vc <- list()
  new_clients_vc <- destination_route$clients_vc
  if( length(client_hoppers_the_other_vehicles_type)==0 ){ # no tiene tolvas de truck; todas son de trailer
    
    destination_route[[used_hoppers_vehicle_type]] <- destination_route[[used_hoppers_vehicle_type]] + length(client_hoppers_vehicles_type)
    destination_route$total_load <- destination_route$total_load + sum(client_hoppers_demands)
    new_clients_vc[[length(destination_route$clients_vc)+1]] <- client_vc_info
    destination_route$clients_vc <- new_clients_vc
    # Fallo 27/04/2021. Faltaba por añadir lo siguiente
    destination_route$VCs <- c(destination_route$VCs, inserted_client)
    
    
  }else{ 
    new_hoppers_vehicles_type <- list()
    t <- 1
    client_hoppers_demands_res <- client_hoppers_demands
    for (f in 1:length(client_hoppers_demands)){
      if(needed_vehicle_type_hoppers_per_product[f]!=0){
        while(client_hoppers_demands_res[f]>0){
          quantity <- min(client_hoppers_demands_res[f], capacity_vehicle_type_hoppers)
          new_hoppers_vehicles_type[[t]] <- data.frame(f, quantity, quantity/capacity_vehicle_type_hoppers)
          colnames(new_hoppers_vehicles_type[[t]]) <- c("Pienso", "Cantidad", "Proporcion")
          client_hoppers_demands_res[f] <- client_hoppers_demands_res[f] - quantity
          t <- t + 1
        }
      }
    }
    
    destination_route$total_load <- destination_route$total_load  + sum(client_hoppers_demands)
    destination_route[[used_hoppers_vehicle_type]] <- destination_route[[used_hoppers_vehicle_type]] + length(new_hoppers_vehicles_type)
    new_clients_vc[[length(destination_route$clients_vc)+1]] <- client_vc_info
    new_clients_vc[[length(destination_route$clients_vc)+1]][[hoppers_the_other_vehicles_type]] <- list()
    new_clients_vc[[length(destination_route$clients_vc)+1]][[hoppers_vehicles_type]] <- new_hoppers_vehicles_type
    destination_route$clients_vc <- new_clients_vc
    destination_route$VCs <- c(destination_route$VCs, inserted_client)
    
  }
  
  
  return(destination_route)
}



# Fallo 27/04/2021. Añadi penalty_max como input en la siguiente funcion
is_there_available_compartments <- function(total_truck_hoppers, total_trailer_hoppers, destination_route, 
                                            client_hoppers_demands, needed_trailer_hoppers_per_product, 
                                            capacity_truck_hoppers, capacity_trailer_hoppers, penalty_max){
  
  available_compartments <- FALSE
  
  # Fallo 27/04/2021. Añadi el penalty_max aqui
  available_trailer_hoppers <- total_trailer_hoppers - destination_route$used_hoppers_trailer + floor(penalty_max/capacity_trailer_hoppers)
  
  needed_trailer_hoppers_per_product_res <- needed_trailer_hoppers_per_product
  needed_trailer_hoppers_per_product_res2 <- needed_trailer_hoppers_per_product_res
  load_res <- client_hoppers_demands
  
  
  for (i in 1:length(needed_trailer_hoppers_per_product)){
    while(available_trailer_hoppers>0 & needed_trailer_hoppers_per_product_res2[i]>0){
      needed_trailer_hoppers_per_product_res2[i] <- max(needed_trailer_hoppers_per_product_res[i]-available_trailer_hoppers,0 )
      load_res[i] <- max(0, load_res[i] - (needed_trailer_hoppers_per_product_res[i]-needed_trailer_hoppers_per_product_res2[i])*capacity_trailer_hoppers)
      available_trailer_hoppers <- available_trailer_hoppers - (needed_trailer_hoppers_per_product_res[i] - needed_trailer_hoppers_per_product_res2[i])
      needed_trailer_hoppers_per_product_res <- needed_trailer_hoppers_per_product_res2
    }
  }
  
  
  # Ahora he "consumido" todas las tolvas de trailer, y todavia no he acoplado load_res de mercancia. Tengo que ver si esto me cabe en las tolvas de truck que hay libres
  
  
  # Fallo 27/04/2021. La penalty_max estaba mal en la siguiente linea
  available_truck_hoppers <- total_truck_hoppers - destination_route$used_hoppers_truck + floor(penalty_max/capacity_truck_hoppers)
  needed_truck_hoppers_per_product <- ceiling(load_res/capacity_truck_hoppers)
  needed_truck_hoppers_total <- sum(needed_truck_hoppers_per_product)
  
  # Fallo 27/04/2012. Cambie el condicional poniendo available...
  if(needed_truck_hoppers_total <=  available_truck_hoppers){
    available_compartments <- TRUE
  }
  
  
  
  return(list(available_compartments=available_compartments, needed_truck_hoppers_per_product=needed_truck_hoppers_per_product)) 
}




create_truck_and_trailer_hoppers_for_vc <- function(inserted_client, needed_trailer_hoppers_total, total_truck_hoppers, total_trailer_hoppers, destination_route, client_hoppers_trailers, client_hoppers_trucks,
                                                    client_vc_info, client_hoppers_demands, needed_trailer_hoppers_per_product, capacity_truck_hoppers, capacity_trailer_hoppers, needed_truck_hoppers_per_product){
  
  
  available_trailer_hoppers <- total_trailer_hoppers - destination_route$used_hoppers_trailer
  client_hoppers_demands_res <- client_hoppers_demands
  
  
  new_hoppers_trailers <- list()
  t <- 1
  
  for (f in 1:length(client_hoppers_demands)){
    if(needed_trailer_hoppers_per_product[f]!=0){
      while(client_hoppers_demands_res[f]>0 & available_trailer_hoppers > 0){
        quantity <- min(client_hoppers_demands_res[f], capacity_trailer_hoppers)
        new_hoppers_trailers[[t]] <- data.frame(f, quantity, quantity/capacity_trailer_hoppers)
        colnames(new_hoppers_trailers[[t]]) <- c("Pienso", "Cantidad", "Proporcion")
        client_hoppers_demands_res[f] <- client_hoppers_demands_res[f] - quantity
        available_trailer_hoppers <- available_trailer_hoppers - 1
        t <- t + 1
      }
    }
  }
  
  
  # Ya he llenado las tolvas de trailer. Ahora vamos con las de truck
  available_truck_hoppers <- total_truck_hoppers - destination_route$used_hoppers_truck
  new_hoppers_trucks <- list()
  tt <- 1
  
  for (f in 1:length(client_hoppers_demands)){
    if(needed_truck_hoppers_per_product[f]!=0){
      while(client_hoppers_demands_res[f]>0 & available_truck_hoppers > 0){
        quantity <- min(client_hoppers_demands_res[f], capacity_truck_hoppers)
        new_hoppers_trucks[[tt]] <- data.frame(f, quantity, quantity/capacity_truck_hoppers)
        colnames(new_hoppers_trucks[[tt]]) <- c("Pienso", "Cantidad", "Proporcion")
        client_hoppers_demands_res[f] <- client_hoppers_demands_res[f] - quantity
        available_truck_hoppers <- available_truck_hoppers - 1
        tt <- tt + 1
      }
    }
  }
  
  
  
  new_clients_vc <- list()
  new_clients_vc <- destination_route$clients_vc
  
  
  destination_route$total_load <- destination_route$total_load  + sum(client_hoppers_demands)
  destination_route$used_hoppers_trailer <- destination_route$used_hoppers_trailer + length(new_hoppers_trailers)
  destination_route$used_hoppers_truck <- destination_route$used_hoppers_truck + length(new_hoppers_trucks)
  new_clients_vc[[length(destination_route$clients_vc)+1]] <- client_vc_info
  new_clients_vc[[length(destination_route$clients_vc)+1]]$hoppers_trailers <- new_hoppers_trailers
  new_clients_vc[[length(destination_route$clients_vc)+1]]$hoppers_trucks <- new_hoppers_trucks
  destination_route$clients_vc <- new_clients_vc
  destination_route$VCs <- c(destination_route$VCs, inserted_client)
  
  
  
  
  return(destination_route)
}



# Fallo 27/04/2021. Añadi el penalty_max como input de esta funcion
check_available_compartments <- function(input, result, intermediate_solution, inserted_client, 
                                         destination_route, initial_solution, penalty_max, client_in_subtour){
  
  #resultado <- intermediate_solution
  
  total_truck_hoppers <- length(input$H.camion[1,])
  total_trailer_hoppers <- length(input$H.trailer[1,])
  
  capacity_truck_hoppers <- input$H.camion[1,1]
  capacity_trailer_hoppers <- input$H.trailer[1,1]
  
  
  n1 <- input$n1
  avail <- FALSE
  
  client_origin_route_index <- route_of_client(inserted_client, initial_solution)$index
  client_origin_route_info <- initial_solution[[client_origin_route_index]]
  
  
  if(inserted_client <= n1 ){ #vc
    # Fallo 27/04/2021. Estaba mal calculado el indice
    for(ind_vc in 1:length(client_origin_route_info$clients_vc)){
      if(client_origin_route_info$clients_vc[[ind_vc]]$id == inserted_client ){
        client_vc_index <- ind_vc
      }
    }
    VC_index <- which(client_origin_route_info$VCs==inserted_client)
    client_vc_info <- client_origin_route_info$clients_vc[[client_vc_index]]
    client_hoppers_trailers <- client_vc_info$hoppers_trailers
    client_hoppers_trucks <- client_vc_info$hoppers_trucks
    client_hoppers_demands <- client_vc_info$demands
    
    
    if(destination_route$type == "PTR"){
      needed_truck_hoppers_per_product <- ceiling(client_hoppers_demands/capacity_truck_hoppers)
      needed_truck_hoppers_total <- sum(needed_truck_hoppers_per_product)
      
      # Fallo 27/04/2021. La penalty_max estaba mal en el siguiente condicional
      if(needed_truck_hoppers_total <=  total_truck_hoppers - destination_route$used_hoppers_truck + floor(penalty_max/capacity_truck_hoppers)){
        #destination_route <- create_only_truck_hoppers_for_vc(needed_hoppers_total, total_truck_hoppers, destination_route, client_hoppers_trailers, client_hoppers_trucks,
        #                                                       client_vc_info, client_hoppers_demands, needed_hoppers_per_product, capacity_truck_hoppers)
        
        destination_route <- create_only_one_vehicle_type_hoppers_for_vc(inserted_client, destination_route, client_hoppers_trailers, client_hoppers_trucks, capacity_truck_hoppers,
                                                                         client_vc_info, client_hoppers_demands, needed_truck_hoppers_per_product, "used_hoppers_truck",
                                                                         "hoppers_trailers", "hoppers_trucks")
          
        avail <- TRUE
        #client_origin_route_info$clients_vc <- client_origin_route_info$clients_vc[-client_vc_index]  
        #client_origin_route_info$clients_vc <- client_origin_route_info$clients_vc[-client_vc_index]  
        
      }
      
    }else{ #if destination_route is either CVR or PVR
      
      # 29/04/2021
      # añadir aqui alguna condicion para cuando el vc lo queremos en subtour y, 
      # en consecuencia, toda su mercancia debe ir en los trucks
      # AÑADIR. 
      if(client_in_subtour){
        
        needed_truck_hoppers_per_product <- ceiling(client_hoppers_demands/capacity_truck_hoppers)
        needed_truck_hoppers_total <- sum(needed_truck_hoppers_per_product)
        if(needed_truck_hoppers_total <=  total_truck_hoppers - destination_route$used_hoppers_truck + floor(penalty_max/capacity_truck_hoppers)){ # me caben
          
          #destination_route <- create_only_truck_hoppers_for_vc(needed_truck_hoppers_total, total_truck_hoppers, destination_route, client_hoppers_trailers, client_hoppers_trucks,
          #                                                 client_vc_info, client_hoppers_demands, needed_truck_hoppers_per_product, capacity_truck_hoppers)
          destination_route <- create_only_one_vehicle_type_hoppers_for_vc(inserted_client, destination_route, client_hoppers_trailers, client_hoppers_trucks, capacity_truck_hoppers,
                                                                           client_vc_info, client_hoppers_demands, needed_truck_hoppers_per_product, "used_hoppers_truck",
                                                                           "hoppers_trailers", "hoppers_trucks")
          
          
          avail <- TRUE 
          
        }
      }else{
      
      # Fallo 27/04/2021. Puse <= en vez de ==
      if( total_trailer_hoppers - destination_route$used_hoppers_trailer <= 0){ # no hay tolvas de trailer suficientes: chequeo trucks
        needed_truck_hoppers_per_product <- ceiling(client_hoppers_demands/capacity_truck_hoppers)
        needed_truck_hoppers_total <- sum(needed_truck_hoppers_per_product)
        
        # Fallo 27/04/2021. La penalty_max estaba mal en el siguiente condicional
        if(needed_truck_hoppers_total <=  total_truck_hoppers - destination_route$used_hoppers_truck + floor(penalty_max/capacity_truck_hoppers)){ # me caben
          
          #destination_route <- create_only_truck_hoppers_for_vc(needed_truck_hoppers_total, total_truck_hoppers, destination_route, client_hoppers_trailers, client_hoppers_trucks,
          #                                                 client_vc_info, client_hoppers_demands, needed_truck_hoppers_per_product, capacity_truck_hoppers)
          destination_route <- create_only_one_vehicle_type_hoppers_for_vc(inserted_client, destination_route, client_hoppers_trailers, client_hoppers_trucks, capacity_truck_hoppers,
                                                                           client_vc_info, client_hoppers_demands, needed_truck_hoppers_per_product, "used_hoppers_truck",
                                                                           "hoppers_trailers", "hoppers_trucks")
          
          
          avail <- TRUE 
          
        }
      }else{ # en caso de que si haya tolvas de trailer disponibles (en destination_route) hay que ver si me llegan (junto con las que haya disponibles de truck)

        needed_trailer_hoppers_per_product <- ceiling(client_hoppers_demands/capacity_trailer_hoppers)
        needed_trailer_hoppers_total <- sum(needed_trailer_hoppers_per_product)
        
        # Fallo 27/04/2021. La penalty_max estaba mal colocada en el siguiente condicional
        if(needed_trailer_hoppers_total <=  total_trailer_hoppers - destination_route$used_hoppers_trailer + floor(penalty_max/capacity_trailer_hoppers)){ # me caben
          #destination_route <- create_only_trailer_hoppers_for_vc(needed_trailer_hoppers_total, total_trailer_hoppers, destination_route, client_hoppers_trailers, client_hoppers_trucks,
          #                                                                     client_vc_info, client_hoppers_demands, needed_trailer_hoppers_per_product, capacity_truck_hoppers)
          
          
          destination_route <- create_only_one_vehicle_type_hoppers_for_vc(inserted_client, destination_route, client_hoppers_trucks, client_hoppers_trailers, capacity_trailer_hoppers,
                                                                           client_vc_info, client_hoppers_demands, needed_trailer_hoppers_per_product, "used_hoppers_trailer",
                                                                           "hoppers_trucks", "hoppers_trailers")
          
          
          avail <- TRUE
        }else{
          
          # Fallo 27/04/2021. Tuve que añadir el penalty_max como input en la siguiente funcion
          checking_compartments_available <- is_there_available_compartments(total_truck_hoppers, total_trailer_hoppers, destination_route, 
                                                                             client_hoppers_demands, needed_trailer_hoppers_per_product, 
                                                                             capacity_truck_hoppers, capacity_trailer_hoppers, penalty_max)
          
          available_compartments <- checking_compartments_available$available_compartments
          
          if(available_compartments){
            needed_truck_hoppers_per_product <- checking_compartments_available$needed_truck_hoppers_per_product
            needed_truck_hoppers_total <- sum(needed_truck_hoppers_per_product)
            
            destination_route <- create_truck_and_trailer_hoppers_for_vc(inserted_client, needed_trailer_hoppers_total, total_truck_hoppers, total_trailer_hoppers, destination_route, client_hoppers_trailers, client_hoppers_trucks,
                                                                         client_vc_info, client_hoppers_demands, needed_trailer_hoppers_per_product, capacity_truck_hoppers, capacity_trailer_hoppers, needed_truck_hoppers_per_product)
            
            
            avail <- TRUE 
            
          }
          
          
        }
        
        
      }
      }
    }
    
    
  }else if(inserted_client > n1){ #tc
    # Fallo 27/04/2021. El indice estaba mal calculado
    for(ind_tc in 1:length(client_origin_route_info$clients_tc)){
      if(client_origin_route_info$clients_tc[[ind_tc]]$id == inserted_client ){
        client_tc_index <- ind_tc
      }
    }
    TC_index <- which(client_origin_route_info$TCs==inserted_client)
    client_tc_info <- client_origin_route_info$clients_tc[[client_tc_index]]
    client_hoppers_trucks <- client_tc_info$hoppers_trucks
    client_hoppers_demands <- client_tc_info$demands
    
    
    #Como el cliente es TC, la "destination_route" si o si sera una PTR o el subtour de una CVR,
    # pero la modificacion del subtour se realiza en la funcion perturbacion
    needed_hoppers_per_product <- ceiling(client_hoppers_demands/capacity_truck_hoppers)
    needed_hoppers_total <- sum(needed_hoppers_per_product)
    
    # Fallo 27/04/2021. La penatly_max estaba mal colocada en el siguiente condicional
    if(needed_hoppers_total <=  total_truck_hoppers - destination_route$used_hoppers_truck + floor(penalty_max/capacity_truck_hoppers)){
      new_clients_tc <- list()
      new_clients_tc <- destination_route$clients_tc
      
      destination_route$used_hoppers_truck <- destination_route$used_hoppers_truck + length(client_hoppers_trucks)
      destination_route$total_load <- destination_route$total_load + sum(client_hoppers_demands)
      new_clients_tc[[length(destination_route$clients_tc)+1]] <- client_tc_info
      destination_route$clients_tc <- new_clients_tc
      destination_route$TCs <- c(destination_route$TCs, inserted_client)
      
      avail <- TRUE
      
    }
    
    
    
    
  }
  
  
  
  
  
  
  
  return(list(avail = avail, destination_route = destination_route))
}