# Generar achivo de observaciones para el modelo CRHM

require(lubridate)
require(readxl)

dir_main  <-  'C:/Users/maria/Documents/PEGA/U/HIDRO DE NIEVES/AUX 8'


# Cargar valores

data  <-  read_xlsx(paste(dir_main,'Data_ejemplo.xlsx', sep='/'))

fecha  <-  as.POSIXct(data$`Fecha-Hora de Medicion`)
fecha_init  <-  as.character(as.Date(fecha[1]))
fecha_fin   <-  as.character(as.Date(fecha[length(fecha)]))

URH_1  <-  data
URH_2  <-  data
URH_2$`Temp.del Aire (ºC)`  <-  URH_2$`Temp.del Aire (ºC)`+2

obs_df  <-  createObsDataframe(fecha_init, fecha_fin, reps=2, timezone='UTC')

obs_df$t.1  <-  URH_1$`Temp.del Aire (ºC)`
obs_df$t.2  <-  URH_2$`Temp.del Aire (ºC)`
obs_df$rh.1  <-  URH_1$`Humedad (%)`
obs_df$rh.2  <-  URH_2$`Humedad (%)`
obs_df$u.1  <-  URH_1$`Velocidad del Vto. (m/s)`
obs_df$u.2  <-  URH_2$`Velocidad del Vto. (m/s)`
obs_df$p.1  <-  URH_1$`Pptación Instantánea (mm)`
obs_df$p.2  <-  URH_2$`Pptación Instantánea (mm)`

writeObsFile(obs_df, obsfile = paste(dir_main,"ejemplo_forzantes.obs", sep='/'))