diff --git a/hermesv3_gr/modules/emision_inventories/point_gfas_emission_inventory.py b/hermesv3_gr/modules/emision_inventories/point_gfas_emission_inventory.py index b874c5f27f500eab89c6019408be083e519ef1c3..e34c578a63f6cdc1b0a41b80ae5d5ebefaaf7ca7 100755 --- a/hermesv3_gr/modules/emision_inventories/point_gfas_emission_inventory.py +++ b/hermesv3_gr/modules/emision_inventories/point_gfas_emission_inventory.py @@ -346,9 +346,8 @@ class PointGfasEmissionInventory(EmissionInventory): pollutant['name'], num + 1, len(self.pollutant_dicts)), level=3) print(('\t\tPollutant {0} ({1}/{2})'.format(pollutant['name'], num + 1, len(self.pollutant_dicts)))) aux = netcdf.variables[pollutant['name']][:].flatten()[gdf['src_index']] - gdf[pollutant['name']] = aux * gdf['src_area'] - # print 'masa {0}: {1} '.format(pollutant['name'], gdf[pollutant['name']].sum()) - gdf[pollutant['name']] = (aux / gdf['dst_area']) * netcdf.variables['cell_area'][:].flatten()[ + + gdf[pollutant['name']] = (aux / gdf['dst_area'].values) * netcdf.variables['cell_area'][:].flatten()[ gdf['src_index']] # print netcdf.variables['bc'][:].sum()