An emission data problem about CAM

[fkc000@163_com]

Hello there!    I tried to change the emission of cesm1.2.2, CAM5. And one of the emission file is“cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_so2_surf_2000_c090726.nc”. Here's the problem: there’re only five sections(“agricultural burning emission”, “domestic emission”, “transportation emission”, “waste treatment emission” and “ship emission”) in it, missing “energy emission” and “industry emission”, which I thought are the main sources of SO2 . Is it a mistake when writing emission data or because of absence of abundant data? Because there’s no such problem in BC emission data.    Looking forward to your reply. Thank you!

 

[tilmes]

The engery and industry emissions are placed in the elevated emission files for example see: ar5_mam3_so2_surf_2000_c120315.ncSimone 

 

[fkc000@163_com]

Here's my atm_in file: clim_soilw_file = '/disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart/dvel/clim_soilw.nc'
depvel_file = '/disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart/dvel/depvel_monthly.nc'
depvel_lnd_file = '/disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart/dvel/regrid_vegetation.nc'
exo_coldens_file = '/disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart/phot/exo_coldens.nc'
ext_frc_specifier =
'SO2 ->/disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_so2_elev_2000_c090726.nc',
'bc_a1 ->/disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_bc_elev_2000_c090726.nc',
'num_a1 -> /disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_num_a1_elev_2000_c090726.nc',
'num_a2 -> /disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_num_a2_elev_2000_c090726.nc',
'pom_a1 -> /disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_oc_elev_2000_c090726.nc',
'so4_a1 -> /disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_so4_a1_elev_2000_c090726.nc',
'so4_a2 -> /disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_so4_a2_elev_2000_c090726.nc'
ext_frc_type = 'CYCLICAL'
fstrat_list = ' '
rsf_file = '/disk6/fkc/cesm1_2_2/input_data/atm/waccm/phot/RSF_GT200nm_v3.0_c080416.nc'
season_wes_file = '/disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart/dvel/season_wes.nc'
sol_facti_cloud_borne = 1.0D0
srf_emis_specifier = 'DMS -> /disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/aerocom_mam3_dms_surf_2000_c090129.nc',
'SO2 -> /disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_so2_surf_2000_c090726.nc',
'SOAG -> /disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_soag_1.5_surf_2000_c100217.nc',
'bc_a1 -> /disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_bc_surf_2000_c090726.nc',
'num_a1 -> /disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_num_a1_surf_2000_c090726.nc',
'num_a2 -> /disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_num_a2_surf_2000_c090726.nc',
'pom_a1 -> /disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_oc_surf_2000_c090726.nc',
'so4_a1 -> /disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_so4_a1_surf_2000_c090726.nc',
'so4_a2 -> /disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_so4_a2_surf_2000_c090726.nc'

The only file "ar5_mam3_so2_surf_2000_c090726.nc"containing so2 emissions doesn't have energy and industry emissions.

And my compset is FC5. Any suggestions?

 

[tilmes]

Those emissions are in the elevated emissions (called external forcings in this model):

<span>'SO2   ->/disk6/fkc/cesm1_2_2/input_data/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_so2_elev_2000_c090726.nc'<br /><br />Sorry I pasted in the wrong file in my earlier response,<br /><br />Simone<br /></span>

 

 

[fkc000@163_com]

Problem solved.Thank you very much! 

 

[thibaut_lurton@cnrs-orleans_fr]

Hi,Regarding this surface/elevation emission files, I had a question:What happens if, for example, I set up the SO2 surface emissions along with the elevation emissions (to take every source into consideration), and on top of that I want to perform an extra external forcing which is not included in neither of my two surface/elevation emission files (e.g. SO2 of planes)? It seems that one both srf_emis_specifier and ext_frc_specifier fields are specified for the same species, no other input can be added.Would it work to have two file calls in the namelist, for the same species? Like:ext_frc_specifier = 'SO2 -> /first_file_SO2.nc'
                         'SO2 -> /second_file_SO2.nc'and they would be summed up? (I doubt it... My guess is that one would override the other...)Or would I need to merge my external forcing data into a single netcdf file?Thanks,Thibaut.

 

[thibaut_lurton@cnrs-orleans_fr]

Hello to all,An additional question related to the previous one: in any of the surface/elevation emission files, whenever an emission matrix is split into several variables (for example, "emiss_wst" + "emiss_tra", etc.), are those dummy names or do they actually correspond to something in the interpreting code?Put it in an other way: is the code sensitive to the actual names and number of the input variables? (would an only variable "total" have the same effect for example?) My point is to determine whether this whole group of variables is summed up automatically (without respect to their names). In which case adding a new variable with an extra source of emission should not be an issue.Thanks!Thibaut.

 

[tilmes]

Hi Thibaut,all the emission types will be automatically added in the model. You can add additional types of emission sources in your file and that should work, Simone