This function creates a CONTROL file for HYSPLIT model. It uses inputs from a data.frame with the receptor information.
Examples
{
# Do not run
obs <- system.file("data-raw", package = "rtorf")
index <- obs_summary(obs)
dt <- obs_read(index)
df <- dt[1]
control_file <- tempfile()
obs_hysplit_control(df, control = control_file)
ff <- readLines(control_file)
cat(ff, sep = "\n")
obs_hysplit_control_read(control_file)
}
#> Number of files of index: 1
#> sector N
#> <char> <int>
#> 1: flask 1
#> 2: Total sectors 1
#> Detected 0 files with agl
#> Detected 1 files without agl
#> Searching flask...
#> 1: ch4_aoa_aircraft-flask_19_allvalid.txt
#> 11 02 16 03
#> 1
#> 34.4290 141.0370 5800.0
#> -240
#> 5
#> 10000.0
#> 3 11
#> /work/noaa/lpdm/metfiles/hrrr/2011/
#> 20110216_hrrr
#> /work/noaa/lpdm/metfiles/hrrr/2011/
#> 20110215_hrrr
#> /work/noaa/lpdm/metfiles/hrrr/2011/
#> 20110214_hrrr
#> /work/noaa/lpdm/metfiles/hrrr/2011/
#> 20110213_hrrr
#> /work/noaa/lpdm/metfiles/hrrr/2011/
#> 20110212_hrrr
#> /work/noaa/lpdm/metfiles/hrrr/2011/
#> 20110211_hrrr
#> /work/noaa/lpdm/metfiles/hrrr/2011/
#> 20110210_hrrr
#> /work/noaa/lpdm/metfiles/hrrr/2011/
#> 20110209_hrrr
#> /work/noaa/lpdm/metfiles/hrrr/2011/
#> 20110208_hrrr
#> /work/noaa/lpdm/metfiles/hrrr/2011/
#> 20110207_hrrr
#> /work/noaa/lpdm/metfiles/hrrr/2011/
#> 20110206_hrrr
#> /work/noaa/lpdm/metfiles/nams/2011/
#> 20110216_hysplit.t00z.namsa
#> /work/noaa/lpdm/metfiles/nams/2011/
#> 20110215_hysplit.t00z.namsa
#> /work/noaa/lpdm/metfiles/nams/2011/
#> 20110214_hysplit.t00z.namsa
#> /work/noaa/lpdm/metfiles/nams/2011/
#> 20110213_hysplit.t00z.namsa
#> /work/noaa/lpdm/metfiles/nams/2011/
#> 20110212_hysplit.t00z.namsa
#> /work/noaa/lpdm/metfiles/nams/2011/
#> 20110211_hysplit.t00z.namsa
#> /work/noaa/lpdm/metfiles/nams/2011/
#> 20110210_hysplit.t00z.namsa
#> /work/noaa/lpdm/metfiles/nams/2011/
#> 20110209_hysplit.t00z.namsa
#> /work/noaa/lpdm/metfiles/nams/2011/
#> 20110208_hysplit.t00z.namsa
#> /work/noaa/lpdm/metfiles/nams/2011/
#> 20110207_hysplit.t00z.namsa
#> /work/noaa/lpdm/metfiles/nams/2011/
#> 20110206_hysplit.t00z.namsa
#> /work/noaa/lpdm/metfiles/gfs0p25/2011/
#> 20110216_gfs0p25
#> /work/noaa/lpdm/metfiles/gfs0p25/2011/
#> 20110215_gfs0p25
#> /work/noaa/lpdm/metfiles/gfs0p25/2011/
#> 20110214_gfs0p25
#> /work/noaa/lpdm/metfiles/gfs0p25/2011/
#> 20110213_gfs0p25
#> /work/noaa/lpdm/metfiles/gfs0p25/2011/
#> 20110212_gfs0p25
#> /work/noaa/lpdm/metfiles/gfs0p25/2011/
#> 20110211_gfs0p25
#> /work/noaa/lpdm/metfiles/gfs0p25/2011/
#> 20110210_gfs0p25
#> /work/noaa/lpdm/metfiles/gfs0p25/2011/
#> 20110209_gfs0p25
#> /work/noaa/lpdm/metfiles/gfs0p25/2011/
#> 20110208_gfs0p25
#> /work/noaa/lpdm/metfiles/gfs0p25/2011/
#> 20110207_gfs0p25
#> /work/noaa/lpdm/metfiles/gfs0p25/2011/
#> 20110206_gfs0p25
#> 1
#> Foot
#> 0
#> 0.01
#> 11 02 16 02 17
#> 1
#> 0.0 0.0
#> 0.10 0.10
#> 30.0 30.0
#> ./
#> cdump
#> 1
#> 50
#> 00 00 00 00 00
#> 00 00 00 00 00
#> 00 240 00
#> 1
#> 0.0 0.0 0.0
#> 0.0 0.0 0.0 0.0 0.0
#> 0.0 0.0 0.0
#> 0.0
#> 0.0
#> $receptor
#> control time year month day hour minute second
#> <char> <POSc> <int> <char> <char> <char> <num> <num>
#> 1: 11 02 16 03 2011-02-16 03:00:00 2011 02 16 03 0 0
#> track_hours recep1 lat lon alt nlocations
#> <char> <list> <char> <char> <char> <char>
#> 1: -240 2011, 2, 16, 3 34.4290 141.0370 5800.0 1
#> vertical_motion top_model id
#> <char> <char> <char>
#> 1: 5 10000.0 2011x02x16x03x00x34.4290Nx141.0370Ex05800
#>
#> $met
#> $met$models
#> [1] 3
#>
#> $met$files
#> [1] "/work/noaa/lpdm/metfiles/hrrr/2011/"
#> [2] "20110216_hrrr"
#> [3] "/work/noaa/lpdm/metfiles/hrrr/2011/"
#> [4] "20110215_hrrr"
#> [5] "/work/noaa/lpdm/metfiles/hrrr/2011/"
#> [6] "20110214_hrrr"
#> [7] "/work/noaa/lpdm/metfiles/hrrr/2011/"
#> [8] "20110213_hrrr"
#> [9] "/work/noaa/lpdm/metfiles/hrrr/2011/"
#> [10] "20110212_hrrr"
#> [11] "/work/noaa/lpdm/metfiles/hrrr/2011/"
#> [12] "20110211_hrrr"
#> [13] "/work/noaa/lpdm/metfiles/hrrr/2011/"
#> [14] "20110210_hrrr"
#> [15] "/work/noaa/lpdm/metfiles/hrrr/2011/"
#> [16] "20110209_hrrr"
#> [17] "/work/noaa/lpdm/metfiles/hrrr/2011/"
#> [18] "20110208_hrrr"
#> [19] "/work/noaa/lpdm/metfiles/hrrr/2011/"
#> [20] "20110207_hrrr"
#> [21] "/work/noaa/lpdm/metfiles/hrrr/2011/"
#> [22] "20110206_hrrr"
#> [23] "/work/noaa/lpdm/metfiles/nams/2011/"
#> [24] "20110216_hysplit.t00z.namsa"
#> [25] "/work/noaa/lpdm/metfiles/nams/2011/"
#> [26] "20110215_hysplit.t00z.namsa"
#> [27] "/work/noaa/lpdm/metfiles/nams/2011/"
#> [28] "20110214_hysplit.t00z.namsa"
#> [29] "/work/noaa/lpdm/metfiles/nams/2011/"
#> [30] "20110213_hysplit.t00z.namsa"
#> [31] "/work/noaa/lpdm/metfiles/nams/2011/"
#> [32] "20110212_hysplit.t00z.namsa"
#> [33] "/work/noaa/lpdm/metfiles/nams/2011/"
#> [34] "20110211_hysplit.t00z.namsa"
#> [35] "/work/noaa/lpdm/metfiles/nams/2011/"
#> [36] "20110210_hysplit.t00z.namsa"
#> [37] "/work/noaa/lpdm/metfiles/nams/2011/"
#> [38] "20110209_hysplit.t00z.namsa"
#> [39] "/work/noaa/lpdm/metfiles/nams/2011/"
#> [40] "20110208_hysplit.t00z.namsa"
#> [41] "/work/noaa/lpdm/metfiles/nams/2011/"
#> [42] "20110207_hysplit.t00z.namsa"
#> [43] "/work/noaa/lpdm/metfiles/nams/2011/"
#> [44] "20110206_hysplit.t00z.namsa"
#> [45] "/work/noaa/lpdm/metfiles/gfs0p25/2011/"
#> [46] "20110216_gfs0p25"
#> [47] "/work/noaa/lpdm/metfiles/gfs0p25/2011/"
#> [48] "20110215_gfs0p25"
#> [49] "/work/noaa/lpdm/metfiles/gfs0p25/2011/"
#> [50] "20110214_gfs0p25"
#> [51] "/work/noaa/lpdm/metfiles/gfs0p25/2011/"
#> [52] "20110213_gfs0p25"
#> [53] "/work/noaa/lpdm/metfiles/gfs0p25/2011/"
#> [54] "20110212_gfs0p25"
#> [55] "/work/noaa/lpdm/metfiles/gfs0p25/2011/"
#> [56] "20110211_gfs0p25"
#> [57] "/work/noaa/lpdm/metfiles/gfs0p25/2011/"
#> [58] "20110210_gfs0p25"
#> [59] "/work/noaa/lpdm/metfiles/gfs0p25/2011/"
#> [60] "20110209_gfs0p25"
#> [61] "/work/noaa/lpdm/metfiles/gfs0p25/2011/"
#> [62] "20110208_gfs0p25"
#> [63] "/work/noaa/lpdm/metfiles/gfs0p25/2011/"
#> [64] "20110207_gfs0p25"
#> [65] "/work/noaa/lpdm/metfiles/gfs0p25/2011/"
#> [66] "20110206_gfs0p25"
#>
#>
#> $emi
#> ngases gas emissions_rate hours_emissions control
#> <char> <char> <char> <char> <char>
#> 1: 1 Foot 0 0.01 11 02 16 02 17
#> time year month day hour min sec
#> <POSc> <int> <char> <char> <char> <char> <num>
#> 1: 2011-02-16 02:17:00 2011 02 16 02 17 0
#> hyemit.start ident
#> <list> <char>
#> 1: 2011, 2, 16, 2 2011x02x16x02x17x34.4290Nx141.0370Ex05800
#>
#> $grid
#> ngrids center_grid_lat center_grid_lon spacing_y spacing_x span_x span_y
#> <char> <char> <char> <char> <char> <char> <char>
#> 1: 1 0.0 0.0 0.10 0.10 30.0 30.0
#> dir name_conc num_vertical_vels height_vertical_levels
#> <char> <char> <char> <char>
#> 1: ./ cdump 1 50
#> sampling_start_time sampling_end_time sampling_interval
#> <char> <char> <char>
#> 1: 00 00 00 00 00 00 00 00 00 00 00 240 00
#>
#> $deposition
#> particle_diamter_um density_gcc shape dep_velocity molecular_weight
#> <char> <char> <char> <char> <char>
#> 1: 0.0 0.0 0.0 0.0 0.0
#> surface_reactivity_ratio diffusitivity_ratio henry_constant
#> <char> <char> <char>
#> 1: 0.0 0.0 0.0
#>
#> $wet_removal
#> actual_henry in_cloud below_cloud radioactive_decay_half_life
#> <char> <char> <char> <char>
#> 1: 0.0 0.0 0.0 0.0
#> pollutant_resuspension
#> <char>
#> 1: 0.0
#>