Abstract aerosol interface; bulk aerosols (BAM); partial MAM; radiative_aerosol and physical props infra; optics core code; pfUnit tests

cime_config/buildlib

@@ -103,6 +103,7 @@ def _build_cam():
                               case.get_value("COMP_INTERFACE")),
                  os.path.join(atm_root, "src", "dynamics", "utils"),
                  os.path.join(atm_root, "src", "physics", "utils"),
+                 os.path.join(atm_root, "src", "aerosol"),
                  os.path.join(atm_root, "src", "history"),
                  os.path.join(atm_root, "src", "history", "buffers", "src"),
                  os.path.join(atm_root, "src", "history", "buffers", "src", "hash"),

cime_config/namelist_definition_cam.xml

@@ -403,6 +403,182 @@
     </values>
   </entry>
 
+<!-- Radiatively active aerosol -->
+  <entry id="mode_defs">
+    <type>char*256(120)</type>
+    <category>radiation</category>
+    <group>rad_aer_nl</group>
+    <units>none</units>
+    <desc>
+      Definitions for the aerosol modes that may be used in the rad_climate and
+      rad_aer_diag_* variables.
+    </desc>
+    <values>
+      <value>UNSET</value>
+    </values>
+  </entry>
+
+  <entry id="bin_defs">
+    <type>char*256(640)</type>
+    <category>radiation</category>
+    <group>rad_aer_nl</group>
+    <units>none</units>
+    <desc>
+      Definitions for the aerosol bins that may be used in the rad_climate and
+      rad_aer_diag_* variables.
+    </desc>
+    <values>
+      <value>UNSET</value>
+    </values>
+  </entry>
+
+  <entry id="rad_aer_climate">
+    <type>char*256(80)</type>
+    <category>radiation</category>
+    <group>rad_aer_nl</group>
+    <units>none</units>
+    <desc>
+      A list of the radiatively active species, i.e., species that affect the
+      climate simulation via the radiative heating rate calculation.
+    </desc>
+    <values>
+      <value>UNSET</value>
+    </values>
+  </entry>
+
+  <entry id="rad_aer_diag_1">
+    <type>char*256(80)</type>
+    <category>radiation</category>
+    <group>rad_aer_nl</group>
+    <units>none</units>
+    <desc>
+      A list of species to be used in the first diagnostic radiative heating rate
+      calculation. These species are not the ones affecting the climate
+      simulation. This is a hook for performing radiative forcing calculations.
+    </desc>
+    <values>
+      <value>UNSET</value>
+    </values>
+  </entry>
+
+  <entry id="rad_aer_diag_2">
+    <type>char*256(80)</type>
+    <category>radiation</category>
+    <group>rad_aer_nl</group>
+    <units>none</units>
+    <desc>
+      Analogous to rad_aer_diag_1, but for the 2nd diagnostic calculation.
+    </desc>
+    <values>
+      <value>UNSET</value>
+    </values>
+  </entry>
+
+  <entry id="rad_aer_diag_3">
+    <type>char*256(80)</type>
+    <category>radiation</category>
+    <group>rad_aer_nl</group>
+    <units>none</units>
+    <desc>
+      Analogous to rad_aer_diag_1, but for the 3rd diagnostic calculation.
+    </desc>
+    <values>
+      <value>UNSET</value>
+    </values>
+  </entry>
+
+  <entry id="rad_aer_diag_4">
+    <type>char*256(80)</type>
+    <category>radiation</category>
+    <group>rad_aer_nl</group>
+    <units>none</units>
+    <desc>
+      Analogous to rad_aer_diag_1, but for the 4th diagnostic calculation.
+    </desc>
+    <values>
+      <value>UNSET</value>
+    </values>
+  </entry>
+
+  <entry id="rad_aer_diag_5">
+    <type>char*256(80)</type>
+    <category>radiation</category>
+    <group>rad_aer_nl</group>
+    <units>none</units>
+    <desc>
+      Analogous to rad_aer_diag_1, but for the 5th diagnostic calculation.
+    </desc>
+    <values>
+      <value>UNSET</value>
+    </values>
+  </entry>
+
+  <entry id="rad_aer_diag_6">
+    <type>char*256(80)</type>
+    <category>radiation</category>
+    <group>rad_aer_nl</group>
+    <units>none</units>
+    <desc>
+      Analogous to rad_aer_diag_1, but for the 6th diagnostic calculation.
+    </desc>
+    <values>
+      <value>UNSET</value>
+    </values>
+  </entry>
+
+  <entry id="rad_aer_diag_7">
+    <type>char*256(80)</type>
+    <category>radiation</category>
+    <group>rad_aer_nl</group>
+    <units>none</units>
+    <desc>
+      Analogous to rad_aer_diag_1, but for the 7th diagnostic calculation.
+    </desc>
+    <values>
+      <value>UNSET</value>
+    </values>
+  </entry>
+
+  <entry id="rad_aer_diag_8">
+    <type>char*256(80)</type>
+    <category>radiation</category>
+    <group>rad_aer_nl</group>
+    <units>none</units>
+    <desc>
+      Analogous to rad_aer_diag_1, but for the 8th diagnostic calculation.
+    </desc>
+    <values>
+      <value>UNSET</value>
+    </values>
+  </entry>
+
+  <entry id="rad_aer_diag_9">
+    <type>char*256(80)</type>
+    <category>radiation</category>
+    <group>rad_aer_nl</group>
+    <units>none</units>
+    <desc>
+      Analogous to rad_aer_diag_1, but for the 9th diagnostic calculation.
+    </desc>
+    <values>
+      <value>UNSET</value>
+    </values>
+  </entry>
+
+  <entry id="rad_aer_diag_10">
+    <type>char*256(80)</type>
+    <category>radiation</category>
+    <group>rad_aer_nl</group>
+    <units>none</units>
+    <desc>
+      Analogous to rad_aer_diag_1, but for the 10th diagnostic calculation.
+    </desc>
+    <values>
+      <value>UNSET</value>
+    </values>
+  </entry>
+
+
 <!-- Vertical Coordinate -->
 
   <entry id="pver">

cime_config/testdefs/testlist_cam.xml

@@ -229,6 +229,17 @@
     </options>
   </test>
 
+  <!-- microp_aero BAM test -->
+  <test compset="FPHYStest" grid="ne3pg3_ne3pg3_mg37" name="SMS_Ln2" testmods="cam/outfrq_microp_aero_bam_derecho">
+    <machines>
+      <machine name="derecho" compiler="gnu" category="aux_sima"/>
+    </machines>
+    <options>
+      <option name="wallclock">00:10:00</option>
+      <option name="comment">Test for microp_aero with BAM bulk aerosols, including ndrop_bam, nucleate_ice, and inline logic in microp_aero</option>
+    </options>
+  </test>
+
   <!-- CAM7 physics tests -->
   <test compset="FCAM7" grid="ne3pg3_ne3pg3_mg37" name="SMS_Ln9" testmods="cam/outfrq_se_cslam_analy_ic">
     <machines>

cime_config/testdefs/testmods_dirs/cam/outfrq_analy_ic_cam4/user_nl_cam

@@ -27,3 +27,9 @@ hist_add_inst_fields;h1:slflx_PBL,qtflx_PBL,uflx_PBL,vflx_PBL,slflx_cg_PBL,qtflx
 hist_add_inst_fields;h1:slten_PBL,qtten_PBL,uten_PBL,vten_PBL,qvten_PBL,qlten_PBL,qiten_PBL,tten_PBL,rhten_PBL
 hist_add_inst_fields;h1:qt_pre_PBL,sl_pre_PBL,slv_pre_PBL,u_pre_PBL,v_pre_PBL,qv_pre_PBL,ql_pre_PBL,qi_pre_PBL,t_pre_PBL,rh_pre_PBL
 hist_add_inst_fields;h1:DTV,DUV,DVV,DTVKE
+
+! QPC4 is an aquaplanet case - disable prescribed aerosol data
+prescribed_aero_file = ''
+prescribed_aero_specifier = ''
+prescribed_volcaero_file = ''
+aerodep_flx_file = ''

cime_config/testdefs/testmods_dirs/cam/outfrq_microp_aero_bam_derecho/shell_commands

@@ -0,0 +1 @@
+ ./xmlchange CAM_CONFIG_OPTS="--dyn none --physics-suites microp_aero_bam"

cime_config/testdefs/testmods_dirs/cam/outfrq_microp_aero_bam_derecho/user_nl_cam

@@ -0,0 +1,68 @@
+! this snapshot was constructed by:
+! (1) set -chem none in CAM_CONFIG_OPTS in a FHIST_LTso (or any other CAM7 case)
+!     n.b. only CAM7 is supported as its PUMAS accepts tendencies from microp_aero/nucleate_ice
+!          and CAM7 PUMAS is the target for CAM-SIMA CCPPization, not pumas-frozen.
+! (2) isolating the microp_aero section in physpkg.F90
+! as the default 'microp_section' snapshot point will also include the actual microphysics (e.g., PUMAS)
+! (3) cld_macmic_num_steps = 1 also has to be set in CAM user_nl_cam when creating the snapshot.
+
+! this is a CAM7 LT nochem snapshot
+pver = 58
+
+debug_output = 1
+
+! tolerances for testing
+ncdata_check_err = .true.
+min_difference = 2e-15
+
+ncdata = '/glade/campaign/cesm/community/amwg/sima_baselines/cam_sima_test_snapshots/cam_ne3pg3_fhistc_ltso_microp_aero_derecho_gnu_before_c20260518.nc'
+ncdata_check = '/glade/campaign/cesm/community/amwg/sima_baselines/cam_sima_test_snapshots/cam_ne3pg3_fhistc_ltso_microp_aero_derecho_gnu_after_c20260518.nc'
+
+use_hetfrz_classnuc            = .false.
+
+! nucleate_ice:
+nucleate_ice_incloud           = .false.
+nucleate_ice_strat             = 1.0D0
+nucleate_ice_subgrid           = 1.2D0
+nucleate_ice_subgrid_strat             = 1.2D0
+nucleate_ice_use_troplev               = .true.
+use_preexisting_ice            = .true.
+
+! microp_aero:
+microp_aero_wsub_min           =   0.1D0
+microp_aero_wsub_scale         =   0.9D0
+microp_aero_wsubi_min          =   0.001D0
+microp_aero_wsubi_scale                =   1.5D0
+
+! BAM aerosols:
+rad_aer_climate = 'N:sulf:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/sulfate_rrtmg_c080918.nc', 'N:dust1:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/dust1_rrtmg_c080918.nc',
+         'N:dust2:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/dust2_rrtmg_c080918.nc', 'N:dust3:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/dust3_rrtmg_c080918.nc',
+         'N:dust4:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/dust4_rrtmg_c080918.nc', 'N:bcar1:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/bcpho_rrtmg_c080918.nc',
+         'N:bcar2:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/bcphi_rrtmg_c080918.nc', 'N:ocar1:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/ocpho_rrtmg_c080918.nc',
+         'N:ocar2:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/ocphi_rrtmg_c080918.nc', 'N:sslt1:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/seasalt1_rrtmg_c080918.nc',
+         'N:sslt2:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/seasalt2_rrtmg_c080918.nc', 'N:sslt3:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/seasalt3_rrtmg_c080918.nc',
+         'N:sslt4:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/seasalt4_rrtmg_c080918.nc'
+
+! diagnostics for microp_aero BAM suite
+hist_output_frequency;h1: 1*nsteps
+hist_write_nstep0;h1: .true.
+hist_precision;h1: REAL64
+hist_max_frames;h1: 1
+
+! scale_subgrid_vertical_velocity diagnostics
+hist_add_inst_fields;h1: WSUB,WSUBI
+
+! nucleate_ice diagnostics
+hist_add_inst_fields;h1: NIHFTEN,NIIMMTEN,NIDEPTEN,NIMEYTEN
+hist_add_inst_fields;h1: NIREGM,NISUBGRID,NITROP_PD
+hist_add_inst_fields;h1: FHOM,WICE,WEFF
+hist_add_inst_fields;h1: INnso4TEN,INnbcTEN,INndustTEN,INondustTEN,INhetTEN,INhomTEN
+hist_add_inst_fields;h1: INFrehom,INFreIN
+
+! ndrop_bam diagnostics: CCN at fixed supersaturation levels
+hist_add_inst_fields;h1: CCN1,CCN2,CCN3,CCN4,CCN5,CCN6
+
+! ndrop_bam diagnostics: aerosol number concentration (names from physprops files)
+hist_add_inst_fields;h1: SULFATE_m3,DUST1_m3,DUST2_m3,DUST3_m3,DUST4_m3
+hist_add_inst_fields;h1: BCPHO_m3,BCPHI_m3,OCPHO_m3,OCPHI_m3
+hist_add_inst_fields;h1: SEASALT1_m3,SEASALT2_m3,SEASALT3_m3,SEASALT4_m3