atmoschem/eixport

Error on wrf_meta

Schuch666 opened this issue · 1 comments

Schuch666 commented

Hi,

The function wrf_meta() is a very useful idea, but I think the output is not correct. For example the code bellow and output:

>library(eixport)
>wrf_meta()
      vars     description                memory_order field_type stagger coordinates
4   E_ACET mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
5   E_ALD2 mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
6   E_ALDX mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
7   E_ALK3 mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
8   E_ALK4 mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
9   E_ALK5 mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
10  E_BALD mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
11 E_CCOOH mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
12    E_CO mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
13  E_CRES    ug m^-2 s^-1                   EMISSIONS        XYZ     104          NA
14   E_ECC    ug m^-2 s^-1                   EMISSIONS        XYZ     104          NA
15   E_ECI    ug m^-2 s^-1                   EMISSIONS        XYZ     104          NA
16   E_ECJ mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
17   E_ETH mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
18  E_ETHA mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
19  E_ETOH mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
20  E_FORM mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
21   E_GLY mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
22   E_HCL mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
23 E_HCOOH mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
24  E_IOLE mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
25 E_IPROD mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
26  E_ISOP mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
27  E_MACR mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
28   E_MEK mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
29  E_MEO2 mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
30  E_MEOH mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
31  E_MGLY mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
32   E_NH3 mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
33    E_NO mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
34   E_NO2    ug m^-2 s^-1                   EMISSIONS        XYZ     104          NA
35  E_NO3C    ug m^-2 s^-1                   EMISSIONS        XYZ     104          NA
36  E_NO3I    ug m^-2 s^-1                   EMISSIONS        XYZ     104          NA
37  E_NO3J mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
38   E_OLE    ug m^-2 s^-1                   EMISSIONS        XYZ     104          NA
39  E_ORGC    ug m^-2 s^-1                   EMISSIONS        XYZ     104          NA
40  E_ORGI    ug m^-2 s^-1                   EMISSIONS        XYZ     104          NA
41  E_ORGJ mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
42   E_PAR mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
43  E_PHEN    ug m^-2 s^-1                   EMISSIONS        XYZ     104          NA
44  E_PM10    ug m^-2 s^-1                   EMISSIONS        XYZ     104          NA
45 E_PM25I    ug m^-2 s^-1                   EMISSIONS        XYZ     104          NA
46 E_PM25J mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
47 E_PROD2 mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
48 E_PSULF mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
49   E_SO2    ug m^-2 s^-1                   EMISSIONS        XYZ     104          NA
50  E_SO4C    ug m^-2 s^-1                   EMISSIONS        XYZ     104          NA
51  E_SO4I    ug m^-2 s^-1                   EMISSIONS        XYZ     104          NA
52  E_SO4J mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
53  E_TERP mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
54   E_TOL mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
55   E_XYL    degree north LATITUDE, SOUTH IS NEGATIVE         XY     104          NA
1    Times            <NA>                        <NA>       <NA>      NA          NA
2     XLAT     degree east LONGITUDE, WEST IS NEGATIVE         XY     104          NA
3    XLONG mol km^-2 hr^-1                   EMISSIONS        XYZ     104          NA
> library(hackWRF)
> meta(variable = 'E_NO2',attname = 'units')
choose a file:
... 
E_NO2 attribute units:
mol km^-2 hr^-1
> meta(variable = 'E_NO2', attname = 'description')
choose a file:
...
E_NO2 attribute description:
EMISSIONS
> meta(variable = 'E_NO2', attname = 'MemoryOrder')
choose a file:
...
E_NO2 attribute MemoryOrder:
XYZ
> meta(variable = 'E_NO2', attname = 'FieldType')
choose a file:
...
E_NO2 attribute FieldType:
104
> meta(variable = 'E_NO2', attname = 'stagger')
choose a file:
...
E_NO2 attribute stagger:
Z

I think the column names are wrong (should be vars, units, description, memory order, FieldType and ??). Also the units on the line of the E_NO2 variable, for example, don't match the (correct) units for this variable in the file.

If you want help to fix this let me know.

Another note is that the metadata is not always return the correct/updated information about the data, it’s good to use this kind of information with caution.

Thanks

ibarraespinosa commented