<div dir="ltr"><div><div><div><div><div>Hi Lei,<br><br></div>I had spent sometime on reading the GFEDS v4 data using NCL version 6.2.1 and have written the attached NCL script that reads the data for each month, and emission source for 18 years and a particular aerosol type. I do not if you need those calculations as I did inside the loops to obtain final aerosol emissions which is summed over all 6 sources. To check the accuracy of calculation it prints "glob_mn" which can be compared with table values provided on the data site.<br></div><br></div>I hope that may work for you, at least should provide some sense how deal with grouped hdf5 data.<br></div><div>Let us how it works.<br></div><div><br></div>Good Luck!<br></div>Rashed<br><br><br>Date: Tue, 9 Jun 2015 13:35:18 -0400<br>
From: dream916 <<a href="mailto:dream916@gmail.com">dream916@gmail.com</a>><br>
Subject: [ncl-talk] Read HDF5 file<br>
To: ncl-talk forum <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>
<<a href="mailto:CABOT6NTdtm4HYqW5HV%2BV6whq9jzD8xyKVOYxiAQi%2BZgQP%2B%2BuJQ@mail.gmail.com">CABOT6NTdtm4HYqW5HV+V6whq9jzD8xyKVOYxiAQi+ZgQP++uJQ@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br><br>Hi,<br>
<br>
I was trying to read global fire emission dataset v4 in hdf4 format. I<br>
searched online and couldn't find any examples how to read groups within<br>
group. My code is as follows:<br>
<br>
; read hdf5 fire dataset<br>
begin<br>
setfileoption("h5", "FileStructure", "Advanced")<br>
f = addfile("GFED4.0s_1997.hdf5","r")<br>
grps = getfilegroups(f,"/",0)<br>
print(grps)<br>
grps2 = getfilegroups(f,"/emissions/",0)<br>
print(grps2)<br>
lat = f->lat<br>
lon = f->lon<br>
printVarSummary(lat)<br>
printVarSummary(lon)<br>
g = f=>ancill<br>
gridarea = g->grid_cell_area<br>
printVarSummary(gridarea)<br>
emissions = f=>emissions<br>
<br>
in the dataset, emissions is a group which contains two subgroups (month<br>
and partitioning) (see below).<br>
<br>
group: emissions {<br>
<br>
<br>
group: \01 {<br>
<br>
dimensions:<br>
<br>
phony_dim_52 = 720 ;<br>
<br>
phony_dim_53 = 1440 ;<br>
<br>
variables:<br>
<br>
float C(phony_dim_52, phony_dim_53) ;<br>
<br>
string C:long_name = "GFED4s biomass burning carbon emissions" ;<br>
<br>
string C:units = "g C / m^2 / month" ;<br>
<br>
float DM(phony_dim_52, phony_dim_53) ;<br>
<br>
string DM:long_name = "GFED4s biomass burning dry matter emissions" ;<br>
<br>
string DM:units = "kg DM / m^2 / month" ;<br>
<br>
float small_fire_fraction(phony_dim_52, phony_dim_53) ;<br>
<br>
string small_fire_fraction:long_name = "Fraction of total GFED4s<br>
biomass burning emissions originating from \"small fire burned area\"" ;<br>
<br>
string small_fire_fraction:units = "Unitless (fraction)" ;<br>
<br>
<br>
group: partitioning {<br>
<br>
dimensions:<br>
<br>
phony_dim_50 = 720 ;<br>
<br>
phony_dim_51 = 1440 ;<br>
<br>
variables:<br>
<br>
float C_AGRI(phony_dim_50, phony_dim_51) ;<br>
<br>
string C_AGRI:long_name = "Contribution of agricultural waste<br>
burning to total monthly biomass burning carbon emissions" ;<br>
<br>
string C_AGRI:units = "Unitless" ;<br>
<br>
float C_BORF(phony_dim_50, phony_dim_51) ;<br>
<br>
string C_BORF:long_name = "Contribution of boreal forest fires to<br>
total monthly biomass burning carbon emissions" ;<br>
<br>
string C_BORF:units = "Unitless" ;<br>
<br>
So my question is how to read the monthly emissions and partitioning<br>
associated with each month into an array.<br>
Thanks,<br>
Lei Meng<br></div>