create.geotop.inpts.keyword: Creates an 'geotop.inpts' files the keyword and their values...

In ecor/geotopbricks: An R Plug-in for the Distributed Hydrological Model GEOtop

Creates an 'geotop.inpts' files the keyword and their values of a date.frame like the one returned by declared.geotop.inpts.keywords

Description

Creates an 'geotop.inpts' files the keyword and their values of a date.frame like the one returned by declared.geotop.inpts.keywords

Usage

create.geotop.inpts.keyword(
  df,
  wpath = NULL,
  comment.lines = "default",
  header = "default",
  like.meteo.keywords = inherits(df, "sf"),
  coords_keywords = c("MeteoStationCoordinateX", "MeteoStationCoordinateY"),
  num_keyword = "NumberOfMeteoStations",
  ...
)

Arguments

df	data frame returend by declared.geotop.inpts.keywords
wpath	complere path to file (optional). Default is NULL.
comment.lines	string or vector of strings to add as comments for each keyword. If it is NULL the comment lines are omitted.
header	string or vector of strings to add as a header. If it is NULL the header is omitted.
like.meteo.keywords	logical. It is used in case df is a sf object (geospatial), it often used for meteorological stations
coords_keywords	keywords for station coordinates in case like.meteo.keywords==TRUE.
num_keyword	keyword for number of stations in case like.meteo.keywords==TRUE.
...	further arguments for writeLines

Details

In case comment.lines and header are set equal to "default", they are suitably modified within the function code. See the example output.

See Also

writeLines,declared.geotop.inpts.keywords

Examples

library(geotopbricks)
#Simulation working path
wpath <- 
'https://raw.githubusercontent.com/ecor/geotopbricks_doc/master/simulations/panola13_run2xC_test3'
## URL path (RAW VERSION) of 
## https://github.com/ecor/geotopbricks_doc/tree/master/simulations/panola13_run2xC_test3
df <- declared.geotop.inpts.keywords(wpath=wpath)
create.geotop.inpts.keyword(df=df)


df <- list(
## ! START METEO 

NumberOfMeteoStations=18,
MeteoStationCode = c("DWD_01735","DWD_01831","DWD_01832","DWD_01833","DWD_04354","DWD_04755","DWD_05306","DWD_05800","DWD_05802","DWD_06191","CHMI_C1BLAD01","CHMI_C1CHUR01","CHMI_C1HUSI01","CHMI_C1JELE01","CHMI_C1KHOR01","CHMI_C1STRZ01","CHMI_C1VOLR01","CHMI_L1HOJS01"),
MeteoStationElevation = c(628,1449,1436,1307,457,847,940,615,615,684,893,1117.8,492,810,728,811,749,866),
MeteoStationLatitude = c(48.7894,49.1126,49.1129,49.0851,48.7832,48.6819,48.9293,49.028,49.0007,48.9023,48.9911111111111,49.0683333333333,49.0525,48.8030555555556,49.1455555555556,48.9088888888889,48.9088888888889,49.2125),
MeteoStationLongitude = c(13.629,13.1353,13.1338,13.2801,13.3146,13.7351,13.4641,13.2385,13.2137,13.1442,13.6683333333333,13.6152777777778,13.99,13.8897222222222,13.55,13.7202777777778,13.8866666666667,13.1972222222222),
MeteoStationSource = c("DWD","DWD","DWD","DWD","DWD","DWD","DWD","DWD","DWD","DWD","CHMI","CHMI","CHMI","CHMI","CHMI","CHMI","CHMI","CHMI"),
MeteoStationWindVelocitySensorHeight = c(10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10),
MeteoStationTemperatureSensorHeight = c(2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2),
MeteoStationCoordinateX = c(399300.751483554,363922.260272129,363813.623704389,374419.289552317,376193.768289692,406895.244369364,387503.305318229,371234.500519548,369350.262272179,364000.012870334,402582.414145364,398857.930427203,426204.259500581,418472.448448941,394254.345037428,406228.32523286,418420.252345735,368705.300573528),
MeteoStationCoordinateY = c(5404951.24399018,5441647.41925558,5441683.4621225,5438340.55274848,5404725.536756,5392867.0291128,5420733.31978365,5432063.12861282,5429070.90371837,5418255.08734482,5427322.00086173,5435975.56352357,5433783.28344618,5406157.18835262,5444648.8862922,5418116.78091222,5417924.86335169,5452642.98443844)
###! END METEO
)


out1 <- create.geotop.inpts.keyword(df=df,comment.lines="")

###
library(terra)

meteo_sts <- terra::vect(system.file("ex_data/meteo_sts3.rds",package="geotopbricks")) |>
   st_as_sf()
   
   
out2 <- create.geotop.inpts.keyword(df=meteo_sts,comment.lines="")

ecor/geotopbricks documentation built on Jan. 26, 2025, 8:50 p.m.