atmoschem / eixport

Export Emissions to Atmospheric Models
https://atmoschem.github.io/eixport/
Other
27 stars 8 forks source link

ncview not reading E_CO from wrf_emissions example #1

Closed ibarraespinosa closed 6 years ago

ibarraespinosa commented 6 years ago

wrfchemis not being read by ncview despite color scale are read:

dir.create("EMISS")
wrf_create(wrfinput_dir = system.file("extdata", package = "eixport"),
           wrfchemi_dir = "EMISS",
           frames_per_auxinput5 = 24)
# get the name of created file
files <- list.files(path = "EMISS",pattern = "wrfchemi",full.names = T)

# open, put some numbers and write
CO <- wrf_get(file = files[1],name = "E_CO")
CO <- CO[,,1]
CO[] = rnorm(length(CO))

perfil <- c(0.010760058, 0.005280596, 0.002883553, 0.002666932,
            0.005781312, 0.018412838, 0.051900411, 0.077834636,
            0.067919758, 0.060831614, 0.055852868, 0.052468599,
            0.050938043, 0.051921718, 0.052756244, 0.052820165,
            0.058388406, 0.072855890, 0.075267137, 0.063246412,
            0.042713523, 0.029108975, 0.022091855, 0.015298458)

wrf_emission(CO,files[1],profile = perfil,names = "E_CO")
}

And now reading wrchemis with ncview

ncview wrfchemi_d01_2011-08-01_00\:00\:00 Error