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

Black carbon radiative forcing

Yuan Sun

Yuan Sun
Active Member
What version of the code are you using?
CESM2.2.2

Have you made any changes to files in the source tree?
No

Describe your problem or question:
Hi all,

I am running FCnudged and interested in black carbon radiative forcing.

In the user guideline, 8. Physics modifications via the namelist — camdoc documentation, there are instructions for mam4_mode1_noBC and mam4_mode4_noBC.

I see the note: "The current version of the modal aerosol code does not support doing diagnostic radiation calculations with aerosol modes when the model is run with modal_strat_sulfate set to true. This option is not used withcam6 physics, but it is used with WACCM."

Does this mean I should not calculate the BC direct radiative effect using FCnudged with cam6 physics? Does it support now?

Look forward to any comments. Thank you.

Best,
Yuan
 
Hi Yuan-
ACOM people talked about this. No it still does not support. People are not so sure, but Francis Vitt remembered that modal_strat_sulfate is a switch for allowing stratospheric sulfate aerosols in the modal model (MAM). It looks like modal_strat_sulfate is set to true for most FC compsets. If you don't care about stratospheric aerosols, maybe try setting modal_strat_sulfate = .false. (and make sure tropospheric aerosols look pretty much the same to you) and getting your BC radiative forcing from there. See if this works...
Danny
 
Ilaria Quaglia just showed me an FW case where modal_strat_sulfate = .true. and radiative forcings were outputted the same way as the user guideline instructed. So, maybe you can try to see if you can output radiative forcing even with modal_strat_sulfate = .true.
Danny
 

Yuan Sun

Yuan Sun
Active Member
Ilaria Quaglia just showed me an FW case where modal_strat_sulfate = .true. and radiative forcings were outputted the same way as the user guideline instructed. So, maybe you can try to see if you can output radiative forcing even with modal_strat_sulfate = .true.
Danny
Dear Danny,

Thanks for your response. I ran FCnudged in CAM6.3, adding diagnostic lines below in the namelist.

The error was: ERROR: modal_aero_wateruptake_dr cannot be called with optional arguments and having modal_strat_sulfate set to true

After I deleted the diagnostic lines, the simulations ran successfully. Thus, it seemed that the diagnostic lines called modal_aero_wateruptake module, but conflicted with modal_strat_sulfate.

Thanks for your solutions. I would try them.

Best,
Yuan

# calculate BC direct radiative effect:
echo "mode_defs = 'mam4_mode1:accum:=',
'A:num_a1:N:num_c1:num_mr:+',
'A:so4_a1:N:so4_c1:sulfate:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/sulfate_rrtmg_c080918.nc:+',
'A:pom_a1:N:pom_c1:p-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocpho_rrtmg_c130709.nc:+',
'A:soa1_a1:N:soa1_c1:s-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc:+',
'A:soa2_a1:N:soa2_c1:s-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc:+',
'A:soa3_a1:N:soa3_c1:s-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc:+',
'A:soa4_a1:N:soa4_c1:s-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc:+',
'A:soa5_a1:N:soa5_c1:s-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc:+',
'A:bc_a1:N:bc_c1:black-c:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/bcpho_rrtmg_c100508.nc:+',
'A:dst_a1:N:dst_c1:dust:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/dust_aeronet_rrtmg_c141106.nc:+',
'A:ncl_a1:N:ncl_c1:seasalt:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ssam_rrtmg_c100508.nc',
'mam4_mode2:aitken:=',
'A:num_a2:N:num_c2:num_mr:+',
'A:so4_a2:N:so4_c2:sulfate:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/sulfate_rrtmg_c080918.nc:+',
'A:soa1_a2:N:soa1_c2:s-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc:+',
'A:soa2_a2:N:soa2_c2:s-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc:+',
'A:soa3_a2:N:soa3_c2:s-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc:+',
'A:soa4_a2:N:soa4_c2:s-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc:+',
'A:soa5_a2:N:soa5_c2:s-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc:+',
'A:ncl_a2:N:ncl_c2:seasalt:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ssam_rrtmg_c100508.nc:+',
'A:dst_a2:N:dst_c2:dust:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/dust_aeronet_rrtmg_c141106.nc',
'mam4_mode3:coarse:=',
'A:num_a3:N:num_c3:num_mr:+',
'A:dst_a3:N:dst_c3:dust:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/dust_aeronet_rrtmg_c141106.nc:+',
'A:ncl_a3:N:ncl_c3:seasalt:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ssam_rrtmg_c100508.nc:+',
'A:so4_a3:N:so4_c3:sulfate:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/sulfate_rrtmg_c080918.nc',
'mam4_mode4:primary_carbon:=',
'A:num_a4:N:num_c4:num_mr:+',
'A:pom_a4:N:pom_c4:p-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocpho_rrtmg_c130709.nc:+',
'A:bc_a4:N:bc_c4:black-c:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/bcpho_rrtmg_c100508.nc',
'mam4_mode1_nobc:accum:=',
'A:num_a1:N:num_c1:num_mr:+',
'A:so4_a1:N:so4_c1:sulfate:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/sulfate_rrtmg_c080918.nc:+',
'A:pom_a1:N:pom_c1:p-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocpho_rrtmg_c130709.nc:+',
'A:soa1_a1:N:soa1_c1:s-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc:+',
'A:soa2_a1:N:soa2_c1:s-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc:+',
'A:soa3_a1:N:soa3_c1:s-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc:+',
'A:soa4_a1:N:soa4_c1:s-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc:+',
'A:soa5_a1:N:soa5_c1:s-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc:+',
'A:dst_a1:N:dst_c1:dust:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/dust_aeronet_rrtmg_c141106.nc:+',
'A:ncl_a1:N:ncl_c1:seasalt:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ssam_rrtmg_c100508.nc',
'mam4_mode4_nobc:primary_carbon:=',
'A:num_a4:N:num_c4:num_mr:+',
'A:pom_a4:N:pom_c4:p-organic:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/ocpho_rrtmg_c130709.nc:+'" >> user_nl_cam
echo "rad_diag_1 = 'A:Q:H2O', 'N:O2:O2', 'N:CO2:CO2',
'A:O3:O3', 'N:N2O:N2O', 'N:CH4:CH4',
'N:CFC11:CFC11', 'N:CFC12:CFC12',
'M:mam4_mode1_nobc:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/mam4_mode1_rrtmg_aeronetdust_sig1.6_dgnh.48_c140304.nc',
'M:mam4_mode2:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/mam4_mode2_rrtmg_aitkendust_c141106.nc',
'M:mam4_mode3:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/mam4_mode3_rrtmg_aeronetdust_sig1.2_dgnl.40_c150219.nc',
'M:mam4_mode4_nobc:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/mam4_mode4_rrtmg_c130628.nc',
'N:VOLC_MMR1:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/volc_camRRTMG_byradius_sigma1.6_mode1_c170214.nc',
'N:VOLC_MMR2:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/volc_camRRTMG_byradius_sigma1.6_mode2_c170214.nc',
'N:VOLC_MMR3:/work/n02/n02/yuansun/cesm/cesm_inputdata/atm/cam/physprops/volc_camRRTMG_byradius_sigma1.2_mode3_c170214.nc'" >> user_nl_cam
 
Top