From 5d1a4aa9075dc7c03fb9c98eadcc15cc4d0bf743 Mon Sep 17 00:00:00 2001
From: "jordan.schnell" <jordan.schnell@noaa.gov>
Date: Thu, 28 Mar 2024 16:03:36 +0000
Subject: [PATCH 1/2] Fix improperly assigned fire emissions for ebb_dcycle==1
 for retrospectives (NOT operational!)

---
 physics/smoke_dust/module_add_emiss_burn.F90 | 9 ---------
 1 file changed, 9 deletions(-)

diff --git a/physics/smoke_dust/module_add_emiss_burn.F90 b/physics/smoke_dust/module_add_emiss_burn.F90
index 80d91bb0e..036ea8525 100755
--- a/physics/smoke_dust/module_add_emiss_burn.F90
+++ b/physics/smoke_dust/module_add_emiss_burn.F90
@@ -76,15 +76,6 @@ subroutine add_emis_burn(dtstep,dz8w,rho_phy,pi,                  &
 !                   the eastern US,   max_ti= 20.041288* 3600.
 !                   else   max_ti= 18.041288* 3600.
 !                   endif
-! RAR: for option #1 ebb and frp are ingested for 24 hours. No modification is
-! applied! 
-    if (ebb_dcycle==1) then
-      do k=kts,kte
-        do i=its,ite
-         ebu(i,k,1)=ebu_in(i,1)   ! RAR:          
-        enddo
-      enddo
-     endif
 
     if (ebb_dcycle==2) then
 

From 2594a9e2c73c2a3a9acaa78fef8252cf5de3b524 Mon Sep 17 00:00:00 2001
From: Grant Firl <grant.firl@noaa.gov>
Date: Fri, 6 Dec 2024 17:22:11 -0500
Subject: [PATCH 2/2] add back Anning's change to evaporation in Thompson when
 using MERRA2 aerosols

---
 physics/MP/Thompson/module_mp_thompson.F90   |  2 +-
 physics/smoke_dust/module_add_emiss_burn.F90 | 17 -----------------
 2 files changed, 1 insertion(+), 18 deletions(-)

diff --git a/physics/MP/Thompson/module_mp_thompson.F90 b/physics/MP/Thompson/module_mp_thompson.F90
index ac5e69a4b..3114859b5 100644
--- a/physics/MP/Thompson/module_mp_thompson.F90
+++ b/physics/MP/Thompson/module_mp_thompson.F90
@@ -3615,7 +3615,7 @@ subroutine mp_thompson (qv1d, qc1d, qi1d, qr1d, qs1d, qg1d, ni1d,     &
 
    !+---+-----------------------------------------------------------------+ !  EVAPORATION
                elseif (clap .lt. -eps .AND. ssatw(k).lt.-1.E-6 .AND.     &
-                        (is_aerosol_aware .or. merra2_aerosol_aware)) then  
+                        is_aerosol_aware) then  
                   tempc = temp(k) - 273.15
                   otemp = 1./temp(k)
                   rvs = rho(k)*qvs(k)
diff --git a/physics/smoke_dust/module_add_emiss_burn.F90 b/physics/smoke_dust/module_add_emiss_burn.F90
index 831f11139..ff1fb09b6 100755
--- a/physics/smoke_dust/module_add_emiss_burn.F90
+++ b/physics/smoke_dust/module_add_emiss_burn.F90
@@ -66,23 +66,6 @@ subroutine add_emis_burn(dtstep,dz8w,rho_phy,pi,ebb_min,          &
     timeq= mod(timeq,timeq_max)
     coef_con=1._kind_phys/((2._kind_phys*pi)**0.5)
 
-
-! RAR: Grasslands (29% of ther western HRRR CONUS domain) probably also need to
-! be added below, check this later
-! RAR: In the HRRR CONUS domain (western part) crop 11%, 2% cropland/natural
-! vegetation and 0.4% urban of pixels
-!.OR. lu_index(i,j)==14)  then ! Croplands(12), Urban and Built-Up(13),
-!cropland/natural vegetation (14) mosaic in MODI-RUC vegetation classes
-! Peak hours for the fire activity depending on the latitude
-!                   if (xlong(i,j)<-130.) then  max_ti= 24.041288* 3600.    !
-!                   peak at 24 UTC, fires in  Alaska
-!                   elseif (xlong(i,j)<-100.) then   max_ti= 22.041288* 3600.
-!                   ! peak at 22 UTC, fires in the western US
-!                   elseif (xlong(i,j)<-70.) then   ! peak at 20 UTC, fires in
-!                   the eastern US,   max_ti= 20.041288* 3600.
-!                   else   max_ti= 18.041288* 3600.
-!                   endif
-
     if (ebb_dcycle==2) then
     
      do j=jts,jte