Found wrong unit in wrfchemi
Closed this issue · 5 comments
what value did you use for the arguments of the wrf_create function?
wrf_create(wrfinput_dir = "WRF/ABRIL/C1/CET/", # ext4
wrfchemi_dir = "WRF/ABRIL/C1/CET", # ext4
domains = 1,
frames_per_auxinput5 = length(sect),
variables = polss,
verbose = T,
separator = "_")
R version 3.6.2 (2019-12-12)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 19.04
Matrix products: default
BLAS: /usr/lib/x86_64-linux-gnu/blas/libblas.so.3.8.0
LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.8.0
locale:
[1] LC_CTYPE=pt_BR.UTF-8 LC_NUMERIC=C LC_TIME=pt_BR.UTF-8
[4] LC_COLLATE=pt_BR.UTF-8 LC_MONETARY=pt_BR.UTF-8 LC_MESSAGES=pt_BR.UTF-8
[7] LC_PAPER=pt_BR.UTF-8 LC_NAME=C LC_ADDRESS=C
[10] LC_TELEPHONE=C LC_MEASUREMENT=pt_BR.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] eixport_0.4.3
loaded via a namespace (and not attached):
[1] Rcpp_1.0.3 rstudioapi_0.10 raster_3.0-12 magrittr_1.5 sfheaders_0.1.0
[6] units_0.6-5 tidyselect_1.0.0 lattice_0.20-38 R6_2.4.1 rlang_0.4.4
[11] cptcity_1.0.4 dplyr_0.8.4 tools_3.6.2 grid_3.6.2 packrat_0.5.0
[16] gibble_0.3.0 KernSmooth_2.23-16 e1071_1.7-3 DBI_1.1.0 class_7.3-15
[21] assertthat_0.2.1 lifecycle_0.1.0 tibble_2.1.3 silicate_0.2.0 crayon_1.3.4
[26] sf_0.8-1 tidyr_1.0.2 purrr_0.3.3 vctrs_0.2.2 codetools_0.2-16
[31] ncdf4_1.17 glue_1.3.1 sp_1.3-2 pillar_1.4.3 compiler_3.6.2
[36] classInt_0.4-2 pkgconfig_2.0.3
You should check the content of your polss variable (passed to variable) and also change n_aero according to the number of aerosols. aerosols must be the last names passed to the variable argument.
This design is based on wrf chem registry convention, do you think I need to change the documentation?
data("emis_opt")
polss <- emis_opt$ecb05_opt1
Do you mean that the units of the aerosols depend on the number of n_aero?
Maybe, more documentation would help. At least, it is not clear for me the meaning and consequences of n_aero :(
So you trying to create a file with:
polss
[1] "E_NO2" "E_XYL" "E_TOL" "E_TERP" "E_SO2" "E_ORA2" "E_OLT" "E_OLI" "E_OL2" "E_NO"
[11] "E_NH3" "E_ISO" "E_HCL" "E_HCHO" "E_ETH" "E_CSL" "E_CO" "E_CH3OH" "E_C2H5OH" "E_ALD"
[21] "E_ALDX" "E_HC3" "E_HC5" "E_HC8" "E_KET" "E_PM25I" "E_PM25J" "E_ECI" "E_ECJ" "E_ORGI"
[31] "E_ORGJ" "E_SO4I" "E_SO4J" "E_NO3I" "E_NO3J" "E_SO4C" "E_NO3C" "E_ORGC" "E_ECC" "E_PM_10"
you need to change n_aero argument from 14 to15