# devtools::install_github("mccreigh/rwrfhydro")

library(raster)
library(mapview)
library(htmlwidgets)
library(rwrfhydro)

# HGT.nc is wrfout netcdf file



demdf=data.frame(x=getValues(brick("HGT.nc",var="XLONG")$layer),
                                    y=getValues(brick("HGT.nc",var="XLAT")$layer),
                                    hgt=getValues(brick("HGT.nc",var="HGT")[[1]]))

# from copygb's outcome in grib format

heigth=brick("HGT.grb")
topo_3km=heigth[[91]]
mapview(topo_3km)
saveRDS(topo_3km,"topo_ARW_3km.rds")

hgt_nc=raster("HGT.nc",varname="HGT")
rwrfhydro::ExportGeogrid("HGT.nc","HGT", "geogrid_hgt.tif")
topowrfout=projectRaster(raster("geogrid_hgt.tif"),CRS("+init=epsg:4326"))
mapview(topowrfout)
saveRDS(topowrfout,"topofromwrfout.rds")