trying to implement surface emission file

[jzweifel]

Hello,

I am a new user to CESM2 and am trying for my first time to adjust my model. Specifically, I am trying to implement a change in sulfate emissions. I tracked down a netcdf file (srf_emis_specifier for num_a2) and manipulated it locally to create a large crude plume of sulfae aersols on the east coast of the U.S. I then reuploaded it into my case directory under the name "aerosol_blob".

Like I said I am trying to implement this change now so I followed the format I found here: Change emissions input - CAM-Chem - wiki.ucar.edu

My user_nl_cam now looks like this:

! Users should add all user specific namelist changes below in the form of
! namelist_var = new_namelist_value

&chem_inparm
srf_emis_specifier = 'num_a2 /glade/work/jzweifel/cases/control_F2000/aerosol_blob.nc'
/


The errors that I am receiving back via email look like this:

PBS Job Id: 1773677.desched1
Job Name: st_archive.control_F2000
Aborted by PBS Server
Job deleted as result of dependency on job 1773676.desched1

I have tried to investigate this a few times after some failed runs and am finding messages that looks like this in my cesm logs:

dec0789.hsn.de.hpc.ucar.edu 1: Opened existing file /glade/work/jzweifel/cases/control_F2000/aerosol_blob.nc
65536
dec0789.hsn.de.hpc.ucar.edu 1: NetCDF: Variable not found
pio_support::pio_die:: myrank= -1 : ERROR: nf_mod.F90: 730 :
NetCDF: Variable not found
dec0790.hsn.de.hpc.ucar.edu 129: NetCDF: Variable not found

Could anyone help me figure out what to do next and how to correctly implement this new emission file I have?

Thanks!

Jack

 

[erik]

This is an issue with CAM, so I'm moving to that forum.
Also alerting @peverley for advice here.

 

[peverley]

Hi Jack,

I took a look at your case directory, and the most useful error appears in your atm.log* file:

srf_emis_inti: spc_name is not included in the simulation
ERROR: srf_emis_inti: invalid surface emission specification

I looked at your namelist changes and saw you have:

&chem_inparm
srf_emis_specifier = 'num_a2 /glade/work/jzweifel/cases/control_F2000/aerosol_blob.nc'
/

The correct syntax for srf_emis_specifier (based on CESM2 CAM6.3 Namelist Definitions) is:

&chem_inparm
srf_emis_specifier = 'num_a2 -> /glade/work/jzweifel/cases/control_F2000/aerosol_blob.nc'
/

I'd suggest making this change and trying again. Let me know if you're still getting an error.

Courtney

 

[jzweifel]

Hi and thank you Courtney,

I actually tried to run this model first with the correct syntax (including the little arrowhead) and was still getting some errors. I thought maybe the arrowhead was actually not supposed to be included so I took it out and got some errors as well.

I tried putting it back in and am still having some issues submitting my case.

Let me know how to best assist you in assisting me!

thanks again,

Jack

 

[peverley]

Hi Jack,

Gotcha. The error you're getting in the cesm.log file is pointing to tracer_data.F90, line 1911. If you look at that line (in /glade/work/jzweifel/cases/derecho_cesm_sandbox/components/cam/src/chemistry/utils/tracer_data.F90) you'll see that it's trying to read in a 'date' variable from your aerosol_blob.nc file and isn't finding it (because it's not there). Can you add a 'date' variable to your file?

Courtney

 

[jzweifel]

I can try to look into this, I'm not exactly sure how this should be done but I will do some investigating! (any tips would be helpful)

Thanks!

Jack

 

[peverley]

I'm certainly not an expert with these files, but you could try doing your same manipulations to a different file (assuming those steps are reproducible) that does contain the date. Here's an example: /glade/campaign/cesm/cesmdata/inputdata/atm/cam/chem/emis/CMIP6_emissions_2000climo/emissions-cmip6_num_so4_a2_anthro-res_surface_2000climo_0.9x1.25_c20170616.nc

Courtney

 

[jzweifel]

I've implemented the date correctly and now everything is running smoothly again. Thanks for the help!