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.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • 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.