suds:

require(lubridate)
require(dplyr)
require(wasa.ops)
require(hydraulics)
require(sf)


### assuming subbasin 64 as close to reference state.


#input.data <- "~/proj/SUDS_Famalicao/CAPITULO_2/modelo_rede_drenagem/run_model/"
input.data <- "~/proj/suds/input.data/"
step=5*60

### prepare precipitation

p <- read.table(paste0(input.data,"/precipitacao_raw"),sep=",") %>% as_tibble


I0h <- p[,-3:-4] %>%
    mutate(dt=dmy_hm(V1)) %>%
    rename(value=V2) %>%
    select(dt,value) %>%
    mutate(year=year(dt),month=month(dt),day=day(dt),hour=hour(dt)) %>%
    mutate_at(vars(year,month, day, hour),as.integer)

t5m <- seq(I0h$dt[1],I0h$dt[length(I0h$dt)],paste0(step," sec"))

tdf <- data_frame(year=year(t5m),month=month(t5m),day=day(t5m),hour=hour(t5m),minute=minute(t5m),second=second(t5m)) %>% mutate_all(as.integer)

I0 <- left_join(tdf,I0h,by=c("year","month","day","hour")) %>%
    mutate(dt=ymd_hms(paste(year,month,day,hour,minute,second))) %>%
    select(dt,value)


## prepare subbasin

losses.sb <- readRDS(paste0(input.data,"/losses.rds")) %>% rename(name=subbacia) %>% as_tibble %>% mutate(model="referencia")

losses.sb <- losses.sb %>% filter(name==64) %>% select(-name,-bacia) %>% left_join(select(losses.sb,name,bacia,model),.)

#subbasinsf <- readRDS(paste0(input.data,"/parm_sb.rds")) %>%
#    st_as_sf %>%
#    as_tibble %>%
#    mutate(dt=5*60) %>%

subbasin <- readRDS(paste0(input.data,"/parm_sb.rds")) %>%
    .@data %>%
    as_tibble %>%
    mutate(step=step) %>%
    rename(L=longest_path,S=av_slope) %>%
    mutate(volume_lago=ifelse(is.na(volume_lago),0,volume_lago)) %>% ## editing retention structures
    mutate(Qoutmax=ifelse(volume_lago==0,0,0.005)) %>%  ### set Qoutmax for lakes
    rename(name=subbacia) %>%
    left_join(.,select(losses.sb,name,initial.proposta,permanent.proposta),by="name") %>%
    rename(hi.max=initial.proposta,he=permanent.proposta) %>%
    mutate(hi=hi.max) %>%
    mutate(X=0.2) %>%### muskingum
    mutate(affluent=0,effluent=0) %>%
    mutate(c.factor=c_existente) %>%
    mutate(pipe.V=0,pipe.Qin=0,pipe.Qout=0) %>% ## state variables
    mutate(runoff.V=0,runoff=0,runoff.out=0) %>% ## state variables
    mutate(structure.V=0,structure.Qin=0,structure.Qout=0,structure.Qoverflow=0) %>% ## state variables
    mutate(model="referencia",volume_lago=0,length_colector=0) %>%
    mutate(loss=0)


subbasin <- subbasin %>% filter(name==64) %>% select(-name,-bacia,-S,-area,-L) %>% left_join(select(subbasin,name,bacia,S,area,L,model),.)


subbasin.template <- subbasin

##### define network

strahler <- read.table(paste0(input.data,"/strahler.proposta"),header=T,sep="\t")