Found wrong unit in wrfchemi

ibarraespinosa commented 4 years ago

It should be ug/m2/s

Schuch666 commented 4 years ago

what value did you use for the arguments of the wrf_create function?

ibarraespinosa commented 4 years ago

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

Schuch666 commented 4 years ago

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?

ibarraespinosa commented 4 years ago

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 :(

Schuch666 commented 4 years ago

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

atmoschem / eixport

Found wrong unit in wrfchemi #53