Feature/rtodling/aero affected radiance

CHANGELOG.md

@@ -10,6 +10,9 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 ### Added
 
 ### Changed
+06Feb2026:
+- revive knob for aerosol-affected radiances (zero-diff when not used)
+
 29Jan2026:
 - revise ABI QC
 

GEOSaana_GridComp/GEOSgsi_Coupler/CMakeLists.txt

@@ -8,16 +8,17 @@ set (SRCSO
   geos_nstmod.F90
   gsi_fixture_GEOS.F90
   g5_pertmod.F90
+  crtm_aerosol.F90
   cplr_pertmod.F90
   cplr_ensemble.F90
   cplr_nst.F90
   cplr_timermod.F90
+  cplr_Set_CRTM_Aerosol.F90
   m_zeitTimer.F90
 #  stub_GEOSagcmPert_GridCompMod.F90
   )
 
 set (SRCSX
-  crtm_aerosol.F90 cplr_Set_CRTM_Aerosol.F90
   crtm_cloud.F90   cplr_Set_CRTM_Cloud.F90
   )
 

GEOSaana_GridComp/GEOSgsi_Coupler/crtm_aerosol.F90

@@ -14,11 +14,10 @@ module crtm_aerosol
   use CRTM_Aerosol_Define, only: CRTM_Aerosol_type
   use CRTM_Parameters,     only: DUST_AEROSOL, SEASALT_SSAM_AEROSOL, &
                                  SEASALT_SSCM1_AEROSOL, SEASALT_SSCM2_AEROSOL, &
-                                 SEASALT_SSCM3_AEROSOL, SEASALT_SSCM3_AEROSOL, &
+                                 SEASALT_SSCM3_AEROSOL, &
                                  BLACK_CARBON_AEROSOL, ORGANIC_CARBON_AEROSOL, &
                                  SULFATE_AEROSOL
 
-
   use Chem_RegistryMod, only: Chem_Registry, Chem_RegistryCreate, Chem_RegistryDestroy
   use Chem_MieMod,      only: Chem_Mie, Chem_MieCreate, Chem_MieDestroy, &
                               Chem_MieQueryIdx, Chem_MieQuery

GEOSaana_GridComp/GSI_GridComp/CMakeLists.txt

@@ -81,6 +81,7 @@ set (SRCS_UTIL
   gridmod.F90
   gscond_ad.f90
   gsi_4dcouplermod.f90
+  gsi_crtmcoupler_aeromod.f90
   gsi_enscouplermod.f90
   gsi_io.f90
   gsi_nemsio_mod.f90
@@ -139,7 +140,6 @@ set (SRCS_UTIL
   rsearch.F90
   rtlnmc_version3.f90
   satthin.F90
-  set_crtm_aerosolmod.f90
   set_crtm_cloudmod.f90
   sfc_model.f90
   simpin1.f90
@@ -731,6 +731,7 @@ if (EXTENDED_SOURCE)
 endif ()
 
 set (SRCS_STUBS
+  stub_crtm_aerosol.f90
   stub_get_wrf_mass_ensperts.f90
   stub_nstmod.f90
   stub_pertmod.F90

GEOSaana_GridComp/GSI_GridComp/set_crtm_aerosolmod.f90 → GEOSaana_GridComp/GSI_GridComp/cplr_gfs_crtm_aeromod.f90

@@ -1,36 +1,4 @@
-module set_crtm_aerosolmod
-!$$$ module documentation block
-!           .      .    .                                       .
-! module:   set_crtm_aerosolmod
-!  prgmmr: todling          org: gmao                date: 2011-06-01
-!
-! abstract: module providing interface to set-crtm-aerosol procedures
-!
-! program history log:
-!   2011-06-01  todling
-!   2011-09-20  hclin   - separate na and na_crtm for p25 handling
-!   2019-03-21  martin - replaced blank subroutine here with that previously
-!                        in stub_set_crtm_aerosol.f90;
-!                        also moved eff rad for dust to size function
-!
-! subroutines included:
-!   sub Set_CRTM_Aerosol_
-!
-! attributes:
-!   language: f90
-!   machine:
-!
-!$$$ end documentation block
-
-implicit none
-
-private
-
-public Set_CRTM_Aerosol
-
-contains
-
-  subroutine Set_CRTM_Aerosol ( km, na, na_crtm, aero_name, aero_conc, rh, aerosol)
+subroutine Set_CRTM_Aerosol_ ( km, na, na_crtm, aero_name, aero_conc, rh, aerosol)
 
 !$$$  subprogram documentation block
 !                .      .    .                                       .
@@ -213,6 +181,5 @@ function GOCART_Aerosol_size( kk, itype,  & ! Input
     return
     end function GOCART_Aerosol_size
 
-  end subroutine Set_CRTM_Aerosol
+  end subroutine Set_CRTM_Aerosol_
 
-end module set_crtm_aerosolmod

GEOSaana_GridComp/GSI_GridComp/crtm_interface.f90

@@ -81,6 +81,7 @@ module crtm_interface
     n_aerosols_jac,aerosol_names_jac,rad_obs_type,cw_cv,ql_cv
 use control_vectors, only: lcalc_gfdl_cfrac
 use ncepnems_io, only: imp_physics
+use gsi_crtmcoupler_aeromod, only: Set_CRTM_Aerosol 
 
 implicit none
 
@@ -886,7 +887,8 @@ subroutine init_crtm(init_pass,mype_diaghdr,mype,nchanl,nreal,isis,obstype,radmo
  endif
 
 ! Initial GFDL saturation water vapor pressure tables
-  if (n_actual_aerosols_wk>0 .or. n_clouds_fwd_wk>0 .and. imp_physics==11) then
+  if (imp_physics==11) then
+   if (n_actual_aerosols_wk>0 .or. n_clouds_fwd_wk>0) then
 
      if (mype==0) write(6,*)myname_,':initial and load GFDL saturation water vapor pressure tables'
 
@@ -907,6 +909,7 @@ subroutine init_crtm(init_pass,mype_diaghdr,mype,nchanl,nreal,isis,obstype,radmo
      des2 (length) = des2 (length - 1)
      desw (length) = desw (length - 1)
 
+   endif
   endif
 
  return
@@ -1100,7 +1103,6 @@ subroutine call_crtm(obstype,obstime,data_s,nchanl,nreal,ich, &
   use constants, only: zero,half,one,one_tenth,fv,r0_05,r10,r100,r1000,constoz,grav,rad2deg, &
       sqrt_tiny_r_kind,two,three,four,five,t0c,rd,eps,rd_over_cp,rearth
   use constants, only: max_varname_length,pi
-  use set_crtm_aerosolmod, only: set_crtm_aerosol
   use set_crtm_cloudmod, only: set_crtm_cloud
   use crtm_module, only: limit_exp,o3_id,toa_pressure
   use obsmod, only: iadate