Scheduled Downtime
On Tuesday 24 October 2023 @ 5pm MT the forums will be in read only mode in preparation for the downtime. On Wednesday 25 October 2023 @ 5am MT, this website will be down for maintenance and expected to return online later in the morning.
Normal Operations
The forums are back online with normal operations. If you notice any issues or errors related to the forums, please reach out to help@ucar.edu

Turn off default aerosols and use my own

nuvolet

nuvolet

Toni Viudez
Member
Greetings,

Is there anyway to tell CESM/CAM to not use the aerosols by default, and instead use my own files?

Thanks in advance.
 
nuvolet

nuvolet

Toni Viudez
Member
Hi Rebecca,

I was wondering if this emissions could be pointed to a daily files, instead of monthly or climatological (CMIP6 for cam emissions) as the example shows?
Another question is that looking at the atm_in file in ~/CaseDocs, I can see how those files are used, if I do not specify some of the files but I only use my own, the ones I do not specify are also used?.
Are those new files need to be in the same spatial resolution as the original ones?.
I also want update those files to more recent ones because all the files for CMIP6 end in 2015, and I would like to make runs beyond that year.

Thanks in advance.

Cheers.

Toni
 
R

rrbuchholz

Rebecca Buchholz
CSEG and Liaisons
Staff member
Hi Toni,
  • Yes you can use emissions files on daily time resolution instead of monthly.
  • You must specify all of the emissions files to be used. You can double check what will be used in your simulation in the atm_in file in the scratch run space once the model is built.
  • It is recommended to have the emissions on the same resolution as the model grid.
  • We do already have some emissions available past 2015 - check out what is available on our Wiki https://wiki.ucar.edu/display/camchem/Emission+Inventories
Hope that helps!
Best,
Rebecca
 
nuvolet

nuvolet

Toni Viudez
Member
Hi Rebecca,

I followed your suggestions from above.
I found some emissions updates from the Quick Fire Emissions Dataset.
The emissions from those files contain, daily data for these species:

"acet","ald2","alk4","bc","c2h6","c3h6","c3h8","ch2o","ch4","co","co2","mek","nh3","no","oc","pm25","so2"

Then, I re-gridded them (0.9x1.25) to be used in CAM6.
I listed those in the nameless files and when I ran a test case I got this error after do the .cas.submit:

Code: 
Creating component namelists
  2023-06-12 11:24:56 atm
   Calling ~CESM2.2/components/cam//cime_config/buildnml
     ...calling cam buildcpp to set build time options
ERROR: Command ~CESM2.2/components/cam/bld/build-namelist -ntasks 12 -csmdata ~CESM2.2/data/input -infile ~CESM2.2/cime/scripts/test-nudged-new-emissions-2/Buildconf/camconf/namelist -start_ymd 20180701 -ignore_ic_date -use_case hist_trop_strat_nudged_cam6 -inputdata ~CESM2.2/cime/scripts/test-nudged-new-emissions-2/Buildconf/cam.input_data_list -namelist " &atmexp /"  failed rc=255
out=
err=Smartmatch is experimental at ~CESM2.2/components/cam/bld/perl5lib/Build/ChemNamelist.pm line 274.
ERROR(Build::Namelist::_parse_next): expect a F90 constant for a namelist instead got:  ALD2 -> ~/QFED/QFED_reg/qfed2.emis_ald2_L288x192_201807.nc4',

It looks that there is problem in another intermediate script (ChemNamelist.pm) at line 274, but I do not know what I should do to avoid this error.

Thanks in advance

Toni
 
nuvolet

nuvolet

Toni Viudez
Member
rrbuchholz said:
Hi Toni,

ALD2 is not a recognized name for a species in CAM-chem. A list of the species can be found in the supplement of the Emmons et al., 2020 description paper: https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019MS001882

Hope that helps!
Rebecca
Click to expand...
Hi Rebecca,

I looked at the supplement you suggested and I revised the list compared to the emissions species I have from QFED.
After replacing and/or deleting those species that are not in that list I ran CESM and still having same error.
For instance this time for black carbon:

Code: 
|------------------> ./case.submit

Creating component namelists
  2023-06-13 11:08:26 atm
   Calling ~/CESM2.2/components/cam//cime_config/buildnml
     ...calling cam buildcpp to set build time options
ERROR: Command ~/CESM2.2/components/cam/bld/build-namelist -ntasks 12 -csmdata ~/CESM2.2/data/input -infile ~/CESM2.2/cime/scripts/test-nudged-new-emissions/Buildconf/camconf/namelist -start_ymd 20180701 -ignore_ic_date -use_case hist_trop_strat_nudged_cam6 -inputdata ~/CESM2.2/cime/scripts/test-nudged-new-emissions/Buildconf/cam.input_data_list -namelist " &atmexp /"  failed rc=255
out=
err=Smartmatch is experimental at ~/CESM2.2/components/cam/bld/perl5lib/Build/ChemNamelist.pm line 274.
ERROR(Build::Namelist::_parse_next): expect a F90 constant for a namelist instead got:  bc_a1 -> ~/QFED/QFED_reg/qfed2.emis_bc_L288x192_201807.nc4',

On the other hand, I do not what source are those BC from QFED, and I assumed they're bc_a1 (accumulation mode)

Thanks

Toni
 
nuvolet

nuvolet

Toni Viudez
Member
rrbuchholz said:
Hi Toni,

ALD2 is not a recognized name for a species in CAM-chem. A list of the species can be found in the supplement of the Emmons et al., 2020 description paper: https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019MS001882

Hope that helps!
Rebecca
Click to expand...
Hi Rebecca,

I think that now I am able to modify the emissions, but now I got an error from a files from WACCM, because one of the files is not updated:

Code: 
flbc_inti: time out of bounds for dataset = ~/data/input/atm/waccm/lb/LBC_1750-2014_CMIP6_0p5degLat_c170126.nc
 ERROR: Unknown error submitted to shr_abort_abort.

Therefore I found these new files should be used to replace that:

Code: 
LBC_20140116-25001216_CMIP6_SSP119_0p5degLat_c190514.nc*
LBC_20140116-25001216_CMIP6_SSP126_0p5degLat_c180905.nc*
LBC_20140116-25001216_CMIP6_SSP245_0p5degLat_c180905.nc*
LBC_20140116-25001216_CMIP6_SSP370_0p5degLat_c180905.nc*
LBC_20140116-25001216_CMIP6_SSP434_0p5degLat_c190514.nc*
LBC_20140116-25001216_CMIP6_SSP460_0p5degLat_c180905.nc*
LBC_20140116-25001216_CMIP6_SSP534os_0p5degLat_c180905.nc*
LBC_20140116-25001216_CMIP6_SSP585_0p5degLat_c180905.nc*

I replaced the file LBC_1750-2014_CMIP6_0p5degLat_c170126.nc by the list above at the same path. But when I run again CESM is still trying to look at the file replaced.

Should I also do the change in the nameless file?.
How should I do that like in https://wiki.ucar.edu/display/camchem/Change+emissions+input?
How do update WACCM files for future runs beyond 2014?.
Thanks in advance
 
nuvolet

nuvolet

Toni Viudez
Member
rrbuchholz said:
Hi Toni,

Great progress! I think the page on changing dates in CAM-chem might help:
https://wiki.ucar.edu/display/camchem/Changing+Dates+of+a+Run

Best,
Rebecca
Click to expand...
Rebecca,

I used one of the LBC files indicated before:

Code: 
&chem_inparm

srf_emis_specifier        = 'bc_a1 -> ~QFED/QFED_reg/qfed2.emis_bc_L288x192_201807.nc4',
    C2H6 -> ~QFED/QFED_reg/qfed2.emis_c2h6_L288x192_201807.nc4',
    C3H6 -> ~QFED/QFED_reg/qfed2.emis_c3h6_L288x192_201807.nc4',
    C3H8 -> ~QFED/QFED_reg/qfed2.emis_c3h8_L288x192_201807.nc4',
    CH20-> ~QFED/QFED_reg/qfed2.emis_ch2o_L288x192_201807.nc4',
    CH4 -> ~QFED/QFED_reg/qfed2.emis_ch4_L288x192_201807.nc4',
    CO -> ~QFED/QFED_reg/qfed2.emis_co_L288x192_201807.nc4',
    CO2 -> ~QFED/QFED_reg/qfed2.emis_co2_L288x192_201807.nc4',
    MEK -> ~QFED/QFED_reg/qfed2.emis_mek_L288x192_201807.nc4',
    NH3 -> ~QFED/QFED_reg/qfed2.emis_nh3_L288x192_201807.nc4',
    NO -> ~QFED/QFED_reg/qfed2.emis_no_L288x192_201807.nc4',
    SO2 -> ~QFED/QFED_reg/qfed2.emis_so2_L288x192_201807.nc4',
&chem_surfvals_nl
flbc_file = '~CESM2.2/data/input/atm/waccm/lb/LBC_20140116-25001216_CMIP6_SSP585_0p5degLat_c180905.nc'

Then when I run this case it gets stuck at the building step:

Code: 
Generating component namelists as part of build
Creating component namelists
  2023-06-28 12:25:09 atm
   Calling ~/CESM2.2/components/cam//cime_config/buildnml
     ...calling cam buildcpp to set build time options

Is there any way to know why and how to solve it?.

Thanks a lot for you help and advice.
 
nuvolet

nuvolet

Toni Viudez
Member
rrbuchholz said:
Hi Toni,

Did you take a look at the troubleshooting page on the wiki? This error occurs sometimes when there is an error in the user_nl_cam file:
https://wiki.ucar.edu/display/camchem/Troubleshooting#Troubleshooting-(2)Builderror-freezing

Best,
Rebecca
Click to expand...
Hi Rebecca,

I look at that link and I had some errors in the user_nl_cam.
I made a csh script to create all the cases plus modifying the user_nl_cam to indicate the modifications I want to do in the nudging module plus aerosols and reanalysis files.
With these files I founder the emissions (QFED) to update runs after 2015 and after re-grideding them to the specific spatial resolution I have that:

Code: 
    echo " srf_emis_specifier        = 'bc_a1 -> ~QFED/QFED_reg/qfed2.emis_bc_L288x192_201807.nc4'," >> user_nl_cam   
    echo "         ' C2H6 -> ~QFED/QFED_reg/qfed2.emis_c2h6_L288x192_201807.nc4'," >> user_nl_cam
    echo "         ' C3H6 -> ~QFED/QFED_reg/qfed2.emis_c3h6_L288x192_201807.nc4'," >> user_nl_cam
    echo "         ' C3H8 -> ~QFED/QFED_reg/qfed2.emis_c3h8_L288x192_201807.nc4'," >> user_nl_cam
    echo "         ' CH20-> ~QFED/QFED_reg/qfed2.emis_ch2o_L288x192_201807.nc4'," >> user_nl_cam
    echo "         ' CH4 -> ~QFED/QFED_reg/qfed2.emis_ch4_L288x192_201807.nc4'," >> user_nl_cam
    echo "         ' CO -> ~QFED/QFED_reg/qfed2.emis_co_L288x192_201807.nc4'," >> user_nl_cam
    echo "         ' CO2 -> ~QFED/QFED_reg/qfed2.emis_co2_L288x192_201807.nc4'," >> user_nl_cam
    echo "         ' MEK -> ~QFED/QFED_reg/qfed2.emis_mek_L288x192_201807.nc4'," >> user_nl_cam
    echo "         ' NH3 -> ~QFED/QFED_reg/qfed2.emis_nh3_L288x192_201807.nc4'," >> user_nl_cam
    echo "         ' NO -> ~QFED/QFED_reg/qfed2.emis_no_L288x192_201807.nc4'," >> user_nl_cam
    echo "         ' SO2 -> ~QFED/QFED_reg/qfed2.emis_so2_L288x192_201807.nc4'" >> user_nl_cam
    echo " " >> user_nl_cam               
    echo "&chem_surfvals_nl" >> user_nl_cam
    echo " flbc_file =        '~CESM2.2/data/input/atm/waccm/lb/LBC_20140116-25001216_CMIP6_SSP585_0p5degLat_c180905.nc'" >> user_nl_cam

Then CESM shows this error in the cesm.log file:

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

Commenting line where is CH2O to avoid that, then new error in cesm.log:

Code: 
ERROR: srf_emis_inti: ERROR -- cannot specify both fixed LBC and emissions for the same species: CH4

Now, if I comment the line where is specified CH4, new error:

Code: 
ERROR: srf_emis_inti: ERROR -- cannot specify both fixed LBC and emissions for the same species: CO2

I do not want to comment more lines to see when the case runs but I want to know how to avoid those errors.
Initially I am running only one month (July 2018).
Do you have any suggestion?.

Thanks in advance
 
You must log in or register to reply here.
Top