A question regarding the aerosol emission file in CESM

[yuanwang@jpl_nasa_gov]

I'm trying to run CESM with cyclical aerosol forcing in year 1970. I found there was a problem with the aerosol surface emission files. My setting details are listed below. The model always crashes with the error message "ENDRUN:open_trc_datafile: cycle year not found" after it opened the file 'aerocom_mam3_dms_surf_1849-2006_c090804.nc'. Does anyone happen to know what is the problem here?

Thanks,
Yuan

srf_emis_type = 'CYCLICAL'
srf_emis_cycle_yr = 1970
srf_emis_specifier = 'DMS -> /CSMDATA/inputdata/atm/cam/chem/trop_mozart_aero/emis/aerocom_mam3_dms_surf_1849-2006_c090804.nc',
'SO2 -> /CSMDATA/inputdata/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_so2_surf_1850-2005_c090804.nc',
'SOAG -> /CSMDATA/inputdata/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_soag_1.5_surf_1850-2005_c100429.nc',
'bc_a1 -> /CSMDATA/inputdata/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_bc_surf_1850-2005_c090804.nc',
'num_a1 -> /CSMDATA/inputdata/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_num_a1_surf_1850-2005_c090804.nc',
'num_a2 -> /CSMDATA/inputdata/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_num_a2_surf_1850-2005_c090804.nc',
'pom_a1 -> /CSMDATA/inputdata/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_oc_surf_1850-2005_c090804.nc',
'so4_a1 -> /CSMDATA/inputdata/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_so4_a1_surf_1850-2005_c090804.nc',
'so4_a2 -> /CSMDATA/inputdata/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_so4_a2_surf_1850-2005_c090804.nc'

 

[eaton]

In 'CYCLICAL' mode the tracer_data module assumes that the cycle year is in the dataset.  The file for dms only contains data for years 1849 and 2006.  The easiest way to fix this problem is to create a dms dataset that contains the 1970 data.  Just linearly interpolate for each month between the years in the existing file. 

 

[yuanwang@jpl_nasa_gov]

Thanks for your reply, Eaton,
About the new DMS dataset, should I just insert one year 1970, or I should do interpolation every ten year like the other emission data files? Do I need to modify the tracer_data module code to let model know that?Yuan

 

[eaton]

The cyclical mode only needs the 12 monthly averages for 1970.  They could be in a file by themselves, or could be inserted into a copy of the existing file.  There are no modifications needed in tracer_data. 

 

[yuanwang@jpl_nasa_gov]

Create a new dms data file by interpolating the original dms data and setting the new date as 1970****. Now problem has been fixed. Thanks Eaton.Btw, just found another DMS emission data file available on the server 'DMS_emissions_1.9x2.5_1850-2000_surface_c070320', which has temperal series of DMS. Can we use that file?

 

[tilmes]

you could also use this file:inputdata/atm/cam/chem/emis/accmip_mz4_1850_2000/IPCC_emissions_DMS_surface_1850-2000_1.9x2.5.nc Simone

 

[loulij@gmail_com]

Hi, I wonder which compset did you use? I am also trying to do a transient run with cyclical emissions but didn't find the proper compset. Thank you so much! Jing

 

[yuanwang@jpl_nasa_gov]

Hi,
Since I wanted to get an equilibrium climate state, I used F_2000_CAM5. If I understand you correctly, you want a transient run with fixed aerosol forcing but time-varying all other forcings? If that is the case, you can build a case using F_AMIP_CAM5 and then specify the aerosol focing types as 'cyclical' in user_nl_cam.  Yuan

 

[loulij@gmail_com]

Hi Yuan, Thanks a lot for your reply! I am trying to do a transient run with climate driving aerosols with emission fixed to a certain year. For CAM5, it seems aerosols are not interactive with climate?  Thanks!Jing

 

[yuanwang@jpl_nasa_gov]

Hi Jing,
In CAM5, aerosols are fully interactive with the grid-scale clouds, radation and climate. You can find more details in the technical note of CAM5. Yuan

 

[loulij@gmail_com]

I see. I 'll look into the document. Many thanks for your help!