September 2016 Review

Reading info

Site records

site_prcp_sept_tbl <- read_table(here::here("data","csv","SEPTIEMBRE.txt"), 
    skip = 1)

site_time_tbl<-paste0(site_prcp_sept_tbl$Date," ",site_prcp_sept_tbl$Time,site_prcp_sept_tbl[,3] %>% pull()) %>% 
        str_replace_all(.,"a"," AM") %>% 
        str_replace_all(.,"p"," PM")

site_time_tbl2<-fast_strptime(x = site_time_tbl,format = "%d/%m/%y %H:%M:%S %p",tz = "UTC") %>% 
        as.POSIXct()

site_prcp_z<-zoo(site_prcp_sept_tbl$Rate,site_time_tbl2)

PERSIANN - CCS Source

site_persian_ccs_tbl <- read_csv(here::here("data","csv","CCS_2023090100_2023091623.csv"))


site_persian_ccs_time_tbl1<-site_persian_ccs_tbl$Time %>% 
        str_remove_all(.,"CCS_1h") %>% 
        ymd_h() 

site_persian_ccs_z<-zoo(site_persian_ccs_tbl$`Rain(mm)`,site_persian_ccs_time_tbl1)

PERSIANN Source

site_persian_tbl <- read_csv(here::here("data","csv","PERSIANN_2023090100_2023091623.csv"))

site_persiann_time_tbl<-site_persian_tbl$Time %>% 
        str_remove_all(.,"PERSIANN_1h") %>% 
        ymd_h() 

site_persiann_z<-zoo(site_persian_tbl$`Rain(mm)`,site_persiann_time_tbl)

Results

The results were compared at hourly and daily scale

Hourly

merge(SITE=site_prcp_z,
        PERSIANN_CCS=site_persian_ccs_z,
      PERSIANN=site_persiann_z) %>% dygraph(group=1) %>% 
        dyBarChart()

Daily

merge(SITE=site_prcp_z,
        PERSIANN_CCS=site_persian_ccs_z,
      PERSIANN=site_persiann_z) %>% 
        subdaily2daily(FUN=sum) %>% dygraph(ylab = "mm/d",group=1) %>% 
        dyBarChart()

merge(SITE=site_prcp_z,
        PERSIANN_CCS=site_persian_ccs_z,
      PERSIANN=site_persiann_z) %>% 
        subdaily2daily(FUN=sum) %>%
        fortify.zoo() %>% 
        pander::pandoc.table(style="simple",caption="daily precipitation based on different sources [mm/d]")

daily precipitation based on different sources [mm/d]

	Index	SITE	PERSIANN_CCS	PERSIANN
1	2023-09-01	0	0	0
49	2023-09-02	0	0	0
97	2023-09-03	0	0	0
145	2023-09-04	0	0	0
193	2023-09-05	0	0	0
241	2023-09-06	0	0	0
289	2023-09-07	0.5	0	0
337	2023-09-08	0	0	0
385	2023-09-09	0	0	0
433	2023-09-10	0	0	0
481	2023-09-11	0	0	0
529	2023-09-12	0	0	0
577	2023-09-13	0	0	0
625	2023-09-14	0.25	0	0
673	2023-09-15	4.77	0	1.13
721	2023-09-16	3.51	14	1.78
769	2023-09-17	0	0	0
817	2023-09-18	0	0	0