Changeset 3398 for TOOLS

Timestamp:

10/02/17 12:07:53 (7 years ago)

Author:

oboucher

Message:

Hopefully last change regarding N>Ox units.
Now emissions are provided in kg NO for INCA.

Location:

TOOLS/CMIP6_FORCING/AER_TROP_EMISSIONS

Files:

• README (modified) (2 diffs)
• REGRID/regrid.sh (modified) (3 diffs)

TOOLS/CMIP6_FORCING/AER_TROP_EMISSIONS/README

@@ -12 +12 @@
 grouped together, and FF and BB emissions are also grouped together by species. 
 FF emissions of NOx are reported as NO2, while BB emissions of NOx are reported as NO. 
-For AER configuration of INCA, we report NOx as NO2 (this may be different in other versions of INCA)
+For AER configuration of INCA, we report NOx as NO (this may be different in other versions of INCA)
 It is unclear as of now if BB emissions should be used as such (with their exact year) or time-filtered.
 
@@ -28 +28 @@
 /homedata/obolmd/AEROSOL
 
-Update: 26/09/2017
-report NOx emissions as NO2
+Update: 02/10/2017
+report NOx emissions as NO for both FF and BB emissions
 
 **********************************************************

TOOLS/CMIP6_FORCING/AER_TROP_EMISSIONS/REGRID/regrid.sh

@@ -5 +5 @@
 #--corrected on 22/6/2017 for BB: undef values zeroed before remapping
 #--corrected on 26/09/2017 for NOx units: PNNL dataset is kg NO2, VUA is kg NO. 
-#  INCA expects NO2 in its AER version, this may have to be changed for the full INCA version
+#  INCA expects NO in its AER version.
 #
 #INCA  conc_dms flx_nox flx_bc flx_pom flx_bbbc flx_bbpom flx_so2 flx_so4 flx_h2s flx_nh3
@@ -71 +71 @@
 fi
 
-#--input file PNNL
+#--input file PNNL (fossil fuel emissions)
 filename=${dirinPNNL}${species}_em_anthro/gn/v20170519/${species}-em-anthro_input4MIPs_emissions_CMIP_CEDS-2017-05-18_gn_${year1}01-${year2}12.nc
 
@@ -470 +470 @@
 rm -f ${dirout}sflx_lmdz_cmip6_${year}.nc
 #--merging all files into a single one
-#--PNNL NOx is NO2 so no change in unit 
-#--VUA NOx is NO so 46/30 ratio to convert to NO2
-cdo merge -expr,'flx_bc=FLX_BC' ${dirout}flux_vector_bc_${year}.nc -expr,'flx_bbbc=FLX_BBBC' ${dirout}flux_vector_bcbb_${year}.nc -expr,'flx_pom=FLX_POM' ${dirout}flux_vector_pom_${year}.nc -expr,'flx_bbpom=FLX_BBPOM' ${dirout}flux_vector_pombb_${year}.nc -expr,'flx_nox=FLX_NOX' ${dirout}flux_vector_nox_${year}.nc -expr,'flx_bbnox=46.*FLX_BBNOX/30.;' ${dirout}flux_vector_noxbb_${year}.nc -expr,'flx_so2=FLX_SO2' ${dirout}flux_vector_so2_${year}.nc -expr,'flx_bbso2=FLX_BBSO2' ${dirout}flux_vector_so2bb_${year}.nc -expr,'flx_nh3=FLX_NH3' ${dirout}flux_vector_nh3_${year}.nc -expr,'flx_bbnh3=FLX_BBNH3' ${dirout}flux_vector_nh3bb_${year}.nc  ${dirout}flux_vector_h2s_${year}.nc ${dirout}flux_vector_so4_${year}.nc -selname,conc_dms ${fileINCAex} ${dirout}sflx_lmdz_cmip6_${year}.nc
+#--PNNL NOx is NO2 so 30/46 scaling factor to change to NO 
+#--VUA NOx is NO so no change in unit
+cdo merge -expr,'flx_bc=FLX_BC' ${dirout}flux_vector_bc_${year}.nc -expr,'flx_bbbc=FLX_BBBC' ${dirout}flux_vector_bcbb_${year}.nc -expr,'flx_pom=FLX_POM' ${dirout}flux_vector_pom_${year}.nc -expr,'flx_bbpom=FLX_BBPOM' ${dirout}flux_vector_pombb_${year}.nc -expr,'flx_nox=30.*FLX_NOX/46.' ${dirout}flux_vector_nox_${year}.nc -expr,'flx_bbnox=FLX_BBNOX;' ${dirout}flux_vector_noxbb_${year}.nc -expr,'flx_so2=FLX_SO2' ${dirout}flux_vector_so2_${year}.nc -expr,'flx_bbso2=FLX_BBSO2' ${dirout}flux_vector_so2bb_${year}.nc -expr,'flx_nh3=FLX_NH3' ${dirout}flux_vector_nh3_${year}.nc -expr,'flx_bbnh3=FLX_BBNH3' ${dirout}flux_vector_nh3bb_${year}.nc  ${dirout}flux_vector_h2s_${year}.nc ${dirout}flux_vector_so4_${year}.nc -selname,conc_dms ${fileINCAex} ${dirout}sflx_lmdz_cmip6_${year}.nc
 
 ncrename -d VECTOR,vector -v VECTOR,vector ${dirout}sflx_lmdz_cmip6_${year}.nc