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

Using rad_diag with cam5-mam3 in cesm1.1.1

S

s_albani@cornell_edu

New Member
Hi!
I have a question on using the diagnostic radiation calculations enabled in cam5.2
The namelist variables rad_diag_* allow specifying a list of species.
I'm using FC5 -trop_mam3, radiation_scheme = 'rrtmg', and the default namelist is defining the 3 modes through the variable "mode_defs".
Category "M" variable are also usead in the default rad_climate definition.
So my question is how can I still use rad_diag_* where I selectively remove some species insetad of a mode?
I tried to specify rad_climate in terms of each species insetad of modes through user_nl_cam.

  rad_climate            = 'A:Q:H2O', 'N:O2:O2', 'N:CO2:CO2', 'N:ozone:O3', 'N:N2O:N2O', 'N:CH4:CH4', 'N:CFC11:CFC11',
  'N:CFC12:CFC12',
  'A:so4_a1:/glade/p/cesm/cseg//inputdata/atm/cam/physprops/sulfate_rrtmg_c080918.nc',
  'A:pom_a1:/glade/p/cesm/cseg//inputdata/atm/cam/physprops/ocpho_rrtmg_c101112.nc',
  'A:soa_a1:/glade/p/cesm/cseg//inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc',
  'A:bc_a1:/glade/p/cesm/cseg//inputdata/atm/cam/physprops/bcpho_rrtmg_c100508.nc',
  'A:dst_a1:/glade/p/cesm/cseg//inputdata/atm/cam/physprops/dust4_rrtmg_c090521.nc',
  'A:ncl_a1:/glade/p/cesm/cseg//inputdata/atm/cam/physprops/ssam_rrtmg_c100508.nc',
  'A:so4_a2:/glade/p/cesm/cseg//inputdata/atm/cam/physprops/sulfate_rrtmg_c080918.nc',
  'A:soa_a2:/glade/p/cesm/cseg//inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc',
  'A:ncl_a2:/glade/p/cesm/cseg//inputdata/atm/cam/physprops/ssam_rrtmg_c100508.nc',
  'A:dst_a3:/glade/p/cesm/cseg//inputdata/atm/cam/physprops/dust4_rrtmg_c090521.nc',
  'A:ncl_a3:/glade/p/cesm/cseg//inputdata/atm/cam/physprops/ssam_rrtmg_c100508.nc',
  'A:so4_a3:/glade/p/cesm/cseg//inputdata/atm/cam/physprops/sulfate_rrtmg_c080918.nc'

The model builds but crashes with the message:

 rad_cnst_get_mode_props: mode_idx=            1   nmodes=            0
  ENDRUN:rad_cnst_get_mode_props: mode list index out of range

Any suggestion is welcome, thanks!
Sam
 
E

eaton

CSEG and Liaisons
To do a diagnostic calculation with species removed from the modes you need to supply new mode definitions.  For example, to remove black carbon which only occurs in mode 1 you can add a definition for a modified mode 1 to the list of modes defined in the mode_defs variable.  Suppose the new mode is called mam3_mode1_noBC.  Then append the following to the list of strings in mode_defs:mode_defs = all the modes that impact climate,
'mam3_mode1_noBC:accum:=',
'A:num_a1:N:num_c1:num_mr:+',
'A:so4_a1:N:so4_c1:sulfate:/inputdata/atm/cam/physprops/sulfate_rrtmg_c080918.nc:+',
'A:pom_a1:N:pom_c1:p-organic:/inputdata/atm/cam/physprops/ocpho_rrtmg_c101112.nc:+',
'A:soa_a1:N:soa_c1:s-organic:/inputdata/atm/cam/physprops/ocphi_rrtmg_c100508.nc:+',
'A:dst_a1:N:dst_c1:dust:/inputdata/atm/cam/physprops/dust_aeronet_rrtmg.nc:+',
'A:ncl_a1:N:ncl_c1:seasalt:/inputdata/atm/cam/physprops/ssam_rrtmg_c100508.nc'

Then assuming the diagnostic calculation for BC forcing is defined in
rad_diag_1:

rad_diag_1 = 'A:Q:H2O', 'N:O2:O2', 'N:CO2:CO2', 'N:ozone:O3',
'N:N2O:N2O', 'N:CH4:CH4', 'N:CFC11:CFC11', 'N:CFC12:CFC12',
'M:mam3_mode1_noBC:/inputdata/atm/cam/physprops/mam3_mode1_rrtmg_c110318.nc',
'M:mam3_mode2:/inputdata/atm/cam/physprops/mam3_mode2_rrtmg_c110318.nc',
'M:mam3_mode3:/inputdata/atm/cam/physprops/mam3_mode3_rrtmg_c110318.nc'

 
S

s_albani@cornell_edu

New Member
Hi!Thank you very much for the clarification.I still have one doubt: from the example you made in defining  rad_diag_1: 'M:mam3_mode1_noBC:/inputdata/atm/cam/physprops/mam3_mode1_rrtmg_c110318.nc',Shouldn't I re-define mam3_mode1_rrtmg_c110318.nc at this point?Thanks,Samuel
 
E

eaton

CSEG and Liaisons
I asked Steve Ghan this same question.  His response:On Thu, Mar 28, 2013 at 11:19:22AM -0700, Ghan, Steven J wrote:
> I don't think a new optics file is needed, because the width of the modes
> is not changed, and removing constituents doesn't increase the range of
> refractive index that is spanned by the constituents in the mode. It makes
> no difference if the range is reduced.

 
C

chen63@illinois_edu

Member
To continue with the discussion about radiation diagnostics. I am running CESM 1.0.4, which I believe is running CAM 5.1.1. The rad_climate variable is specified with species, not modes, as shown below:  rad_climate            = 'P_Q:H2O', 'D_O2:O2', 'D_CO2:CO2', 'D_ozone:O3', 'D_N2O:N2O', 'D_CH4:CH4', 'D_CFC11:CFC11',  'D_CFC12:CFC12', 'P_so4_a1:$DIN_LOC_ROOT/atm/cam/physprops/sulfate_rrtmg_c080918.nc',  'P_pom_a1:$DIN_LOC_ROOT/atm/cam/physprops/ocpho_rrtmg_c101112.nc',  'P_soa_a1:$DIN_LOC_ROOT/atm/cam/physprops/ocphi_rrtmg_c100508.nc',  'P_bc_a1:$DIN_LOC_ROOT/atm/cam/physprops/bcpho_rrtmg_c100508.nc',  'P_dst_a1:$DIN_LOC_ROOT/atm/cam/physprops/dust4_rrtmg_c090521.nc',  'P_ncl_a1:$DIN_LOC_ROOT/atm/cam/physprops/ssam_rrtmg_c100508.nc',  'P_so4_a2:$DIN_LOC_ROOT/atm/cam/physprops/sulfate_rrtmg_c080918.nc',  'P_soa_a2:$DIN_LOC_ROOT/atm/cam/physprops/ocphi_rrtmg_c100508.nc',  'P_ncl_a2:$DIN_LOC_ROOT/atm/cam/physprops/ssam_rrtmg_c100508.nc',  'P_dst_a3:$DIN_LOC_ROOT/atm/cam/physprops/dust4_rrtmg_c090521.nc',  'P_ncl_a3:$DIN_LOC_ROOT/atm/cam/physprops/ssam_rrtmg_c100508.nc',  'P_so4_a3:$DIN_LOC_ROOT/atm/cam/physprops/sulfate_rrtmg_c080918.nc' In this case, if I want to run radiation diagnostics, can I just remove the species I am interested in? For example, if I want to run the radiation diagnostics for black carbon, I just remove bc from the above list, and set rad_diag_1 in the user_nl_cam as follows: rad_diag_1            = 'P_Q:H2O', 'D_O2:O2', 'D_CO2:CO2', 'D_ozone:O3', 'D_N2O:N2O', 'D_CH4:CH4', 'D_CFC11:CFC11',  'D_CFC12:CFC12', 'P_so4_a1:$DIN_LOC_ROOT/atm/cam/physprops/sulfate_rrtmg_c080918.nc',  'P_pom_a1:$DIN_LOC_ROOT/atm/cam/physprops/ocpho_rrtmg_c101112.nc',  'P_soa_a1:$DIN_LOC_ROOT/atm/cam/physprops/ocphi_rrtmg_c100508.nc',  'P_dst_a1:$DIN_LOC_ROOT/atm/cam/physprops/dust4_rrtmg_c090521.nc',  'P_ncl_a1:$DIN_LOC_ROOT/atm/cam/physprops/ssam_rrtmg_c100508.nc',  'P_so4_a2:$DIN_LOC_ROOT/atm/cam/physprops/sulfate_rrtmg_c080918.nc',  'P_soa_a2:$DIN_LOC_ROOT/atm/cam/physprops/ocphi_rrtmg_c100508.nc',  'P_ncl_a2:$DIN_LOC_ROOT/atm/cam/physprops/ssam_rrtmg_c100508.nc',  'P_dst_a3:$DIN_LOC_ROOT/atm/cam/physprops/dust4_rrtmg_c090521.nc',  'P_ncl_a3:$DIN_LOC_ROOT/atm/cam/physprops/ssam_rrtmg_c100508.nc',  'P_so4_a3:$DIN_LOC_ROOT/atm/cam/physprops/sulfate_rrtmg_c080918.nc' Is this the correct way to do it? Is there any other changes I have to make? I am not running higher version than CESM1.0.4 (CAM5.1.1) and don't know what the rad_climate in the later versions look like.  Thanks! Yanju
 
E

eaton

CSEG and Liaisons
cam-5.1.1 (cesm1_0_4) only supports a special case of diagnostic radiation calculations when modal aerosols are used, and that is the case of no aerosols.  This capability can be used to look at the total aerosol forcing, but is not able to look at forcings of individual species in the modes.  For that you should use the latest cesm1_2_0 release in which the radiation diagnostic calculations for modal aerosols have been improved over the version in cesm1_1_1 by accounting for changes in the water uptake and wet radius of the modes when species are removed.
 
You must log in or register to reply here.
Top