From 2594a9e2c73c2a3a9acaa78fef8252cf5de3b524 Mon Sep 17 00:00:00 2001 From: Grant Firl Date: Fri, 6 Dec 2024 17:22:11 -0500 Subject: [PATCH] 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