Changeset 3398


Ignore:
Timestamp:
10/02/17 12:07:53 (4 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:
2 edited

Legend:

Unmodified
Added
Removed
  • TOOLS/CMIP6_FORCING/AER_TROP_EMISSIONS/README

    r3394 r3398  
    1212grouped together, and FF and BB emissions are also grouped together by species.  
    1313FF emissions of NOx are reported as NO2, while BB emissions of NOx are reported as NO.  
    14 For AER configuration of INCA, we report NOx as NO2 (this may be different in other versions of INCA) 
     14For AER configuration of INCA, we report NOx as NO (this may be different in other versions of INCA) 
    1515It is unclear as of now if BB emissions should be used as such (with their exact year) or time-filtered. 
    1616 
     
    2828/homedata/obolmd/AEROSOL 
    2929 
    30 Update: 26/09/2017 
    31 report NOx emissions as NO2 
     30Update: 02/10/2017 
     31report NOx emissions as NO for both FF and BB emissions 
    3232 
    3333********************************************************** 
  • TOOLS/CMIP6_FORCING/AER_TROP_EMISSIONS/REGRID/regrid.sh

    r3397 r3398  
    55#--corrected on 22/6/2017 for BB: undef values zeroed before remapping 
    66#--corrected on 26/09/2017 for NOx units: PNNL dataset is kg NO2, VUA is kg NO.  
    7 #  INCA expects NO2 in its AER version, this may have to be changed for the full INCA version 
     7#  INCA expects NO in its AER version. 
    88# 
    99#INCA  conc_dms flx_nox flx_bc flx_pom flx_bbbc flx_bbpom flx_so2 flx_so4 flx_h2s flx_nh3 
     
    7171fi 
    7272 
    73 #--input file PNNL 
     73#--input file PNNL (fossil fuel emissions) 
    7474filename=${dirinPNNL}${species}_em_anthro/gn/v20170519/${species}-em-anthro_input4MIPs_emissions_CMIP_CEDS-2017-05-18_gn_${year1}01-${year2}12.nc 
    7575 
     
    470470rm -f ${dirout}sflx_lmdz_cmip6_${year}.nc 
    471471#--merging all files into a single one 
    472 #--PNNL NOx is NO2 so no change in unit  
    473 #--VUA NOx is NO so 46/30 ratio to convert to NO2 
    474 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 
     472#--PNNL NOx is NO2 so 30/46 scaling factor to change to NO  
     473#--VUA NOx is NO so no change in unit 
     474cdo 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 
    475475 
    476476ncrename -d VECTOR,vector -v VECTOR,vector ${dirout}sflx_lmdz_cmip6_${year}.nc 
Note: See TracChangeset for help on using the changeset viewer.