atmoschem / eixport

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

Error on wrf_meta #59

Closed Schuch666 closed 1 year ago

Schuch666 commented 4 years ago

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 4 years ago

Thanks for sharing this

Em seg., 31 de ago. de 2020 às 19:47, Daniel Schuch < notifications@github.com> escreveu:

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

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/atmoschem/eixport/issues/59, or unsubscribe https://github.com/notifications/unsubscribe-auth/AGRM74DB6ZL5UMLGMMLL5MDSDQR7ZANCNFSM4QRAR5AA .