Giter VIP home page Giter VIP logo

Comments (6)

gitpeterwind avatar gitpeterwind commented on September 14, 2024

Hi Kevin,
We do not normally distribute our FINN forest fire netcdf files. If you have got a file that has been converted by us, which emission year was it, and what is the creation date? (global attributes created_date and created_hour). We have had some erroneous versions, and maybe your copy is one of them.

from emep-ctm.

kevinbasoa avatar kevinbasoa commented on September 14, 2024

Hi Peter,

I used a inventory that a made myself, with the database from NCAR (https://www2.acom.ucar.edu/modeling/finn-fire-inventory-ncar), I don't have access to any "official" forest fires inventory. Here is a dump from the inventory that I made.

I think that emissions from aerosol in molecules are or could be wrong calculated. This is a inventory for last days of 2018 and January of 2019.

dimensions:
lat = 880 ;
lon = 880 ;
time = 43 ;
variables:
float lat(lat) ;
lat:units = "degrees_north" ;
lat:long_name = "latitude" ;
float lon(lon) ;
lon:units = "degrees_east" ;
lon:long_name = "longitude" ;
float time(time) ;
time:units = "days since 1990-01-01 00:00:00" ;
time:long_name = "time" ;
time:calendar = "Gregorian" ;
float BC(time, lat, lon) ;
BC:units = "mol/day" ;
float CO(time, lat, lon) ;
CO:units = "mol/day" ;
float ACET(time, lat, lon) ;
ACET:units = "mol/day" ;
float ALK4(time, lat, lon) ;
ALK4:units = "mol/day" ;
float BENZ(time, lat, lon) ;
BENZ:units = "mol/day" ;
float C2H4(time, lat, lon) ;
C2H4:units = "mol/day" ;
float C2H6(time, lat, lon) ;
C2H6:units = "mol/day" ;
float C3H8(time, lat, lon) ;
C3H8:units = "mol/day" ;
float CH20(time, lat, lon) ;
CH20:units = "mol/day" ;
float GLYC(time, lat, lon) ;
GLYC:units = "mol/day" ;
float HAC(time, lat, lon) ;
HAC:units = "mol/day" ;
float MEK(time, lat, lon) ;
MEK:units = "mol/day" ;
float MGLY(time, lat, lon) ;
MGLY:units = "mol/day" ;
float PRPE(time, lat, lon) ;
PRPE:units = "mol/day" ;
float TOLU(time, lat, lon) ;
TOLU:units = "mol/day" ;
float XYLE(time, lat, lon) ;
XYLE:units = "mol/day" ;
float CH3CHO(time, lat, lon) ;
CH3CHO:units = "mol/day" ;
float NH3(time, lat, lon) ;
NH3:units = "mol/day" ;
float NO(time, lat, lon) ;
NO:units = "mol/day" ;
float NO2(time, lat, lon) ;
NO2:units = "mol/day" ;
float OC(time, lat, lon) ;
OC:units = "mol/day" ;
float PM25(time, lat, lon) ;
PM25:units = "mol/day" ;
float SO2(time, lat, lon) ;
SO2:units = "mol/day" ;

// global attributes:
:title = "FINN_v2.5_2018" ;

Since this inventory I made it myself, I think is builded wrong.

Could you guide me if there something wrong with this inventory?

Kevin

from emep-ctm.

gitpeterwind avatar gitpeterwind commented on September 14, 2024

and you use the correct unit conversion? they specify
"NOTE on units: In the text files with emissions for each fire, the aerosols (BC, OC, PM2.5, PM10), as well as the gas-phase lumped species 'NOXasNO' and NMOC, are in units of kg/day, and the gas-phase species are in moles/day. In these files 'AREA' is the area burned in m^2; 'BMASS' is the biomass burned per area burned in kg/m^2. 'GENVEG' specifies the vegetation type assumed for the fire: 1 = grasslands and savanna; 2 = woody savanna/shrublands; 3 = tropical forest; 4 = temperate forest; 5 = boreal forest; 6 = temperate evergreen forest; 9 = croplands; 0 = no vegetation.
In the 0.1x0.1 degree gridded files, all compounds are converted to molecules/cm^2/sec, where aerosols have been converted from kg to molecules using molecular weight of 12 g/mole.
"

from emep-ctm.

kevinbasoa avatar kevinbasoa commented on September 14, 2024

I read it and I think the second apply to my case, I only use the netCDF files, not the text files. I convert these to molescules/day, so I did not change molecules to mass. I do not know if EMEP convert "molecules of PM" as they do (NCAR).

Maybe this is making my results bad?

from emep-ctm.

gitpeterwind avatar gitpeterwind commented on September 14, 2024

Yes, the emep model assumes kg/day for PM25, given you are using the "CMX_BiomassBurning_FINNv2p5.txt" conversion file

from emep-ctm.

kevinbasoa avatar kevinbasoa commented on September 14, 2024

Thanks for the clarification.

I will change the units and come back with you.

Best regards

from emep-ctm.

Related Issues (20)

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❤️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.