From 5d1a4aa9075dc7c03fb9c98eadcc15cc4d0bf743 Mon Sep 17 00:00:00 2001 From: "jordan.schnell" Date: Thu, 28 Mar 2024 16:03:36 +0000 Subject: [PATCH 1/3] 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 66a3fea2a590dc8432f418068a275a722426eff9 Mon Sep 17 00:00:00 2001 From: "Haiqin.Li" Date: Wed, 26 Jun 2024 04:23:18 +0000 Subject: [PATCH 2/3] "update deep soil temperature of RUC LSM" --- physics/SFC_Models/Land/RUC/module_sf_ruclsm.F90 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/physics/SFC_Models/Land/RUC/module_sf_ruclsm.F90 b/physics/SFC_Models/Land/RUC/module_sf_ruclsm.F90 index f1647ef81..9677e7bf1 100644 --- a/physics/SFC_Models/Land/RUC/module_sf_ruclsm.F90 +++ b/physics/SFC_Models/Land/RUC/module_sf_ruclsm.F90 @@ -1077,8 +1077,6 @@ SUBROUTINE LSMRUC(xlat,xlon, & tso(i,k,j) = tso1d(k) enddo - tso(i,nzs,j) = tbot(i,j) - do k=1,nzs smfr3d(i,k,j) = smfrkeep(k) keepfr3dflag(i,k,j) = keepfr (k) @@ -1104,8 +1102,10 @@ SUBROUTINE LSMRUC(xlat,xlon, & if(snow(i,j)==zero) EMISSL(i,j) = EMISBCK(i,j) EMISS (I,J) = EMISSL(I,J) ! SNOW is in [mm], SNWE is in [m]; CANWAT is in mm, CANWATR is in m - SNOW (i,j) = SNWE*1000._kind_phys - SNOWH (I,J) = SNHEI + !-- 17 may 2024 - cap snow for points at high elevations where all year round skin temperatures are close to 0 C + !-- Snow density for these points will be 3000/7.5=400 [kg/m^3] + SNOW (i,j) = min(3._kind_phys,SNWE)*1000._kind_phys ! cap to be < 3 m + SNOWH (I,J) = min(7.5_kind_phys,SNHEI) ! cap to be < 7.5 m CANWAT (I,J) = CANWATR*1000._kind_phys if (debug_print) then From 2594a9e2c73c2a3a9acaa78fef8252cf5de3b524 Mon Sep 17 00:00:00 2001 From: Grant Firl Date: Fri, 6 Dec 2024 17:22:11 -0500 Subject: [PATCH 3/3] 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