Skip to content

Add Anthropogenic Emissions #215

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
haiqinli wants to merge 2 commits into
base: gsl/develop
Choose a base branch
from
Open

Add Anthropogenic Emissions #215

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
8c891bc
"add the anthropogenic emission"
haiqinli Feb 17, 2026
241727a
"update for enhanced mixing of anthropogenic emission"
haiqinli Feb 20, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
1 change: 1 addition & 0 deletions src/core_atmosphere/CMakeLists.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -136,6 +136,7 @@ set(ATMOSPHERE_CORE_PHYSICS_SMOKE_SOURCES
seas_data_mod.F90
seas_ngac_mod.F90
ssalt_mod.F90
module_anthro_emissions.F90
)
list(TRANSFORM ATMOSPHERE_CORE_PHYSICS_SMOKE_SOURCES PREPEND physics/physics_noaa/SMOKE/)

Expand Down
59 changes: 45 additions & 14 deletions src/core_atmosphere/physics/mpas_atmphys_driver_pbl.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -131,6 +131,7 @@ subroutine allocate_pbl(configs)
if(.not.allocated(tend_chem_settle_p)) allocate(tend_chem_settle_p(ims:ime,kms:kme,jms:jme,1:num_chem))
if(.not.allocated(chem_p)) allocate(chem_p(ims:ime,kms:kme,jms:jme,1:num_chem))
if(.not.allocated(frp_out_p)) allocate(frp_out_p(ims:ime,jms:jme))
if(.not.allocated(aero_emis_for_enhmix_p)) allocate(aero_emis_for_enhmix_p(ims:ime,jms:jme))
endif

pbl_select: select case (trim(pbl_scheme))
Expand Down Expand Up @@ -311,10 +312,11 @@ subroutine deallocate_pbl(configs)

!deallocation of chemistry arrays
if (do_chemistry) then
if(allocated(ddvel_p) ) deallocate(ddvel_p )
if(allocated(chem_p) ) deallocate(chem_p )
if(allocated(tend_chem_settle_p)) deallocate(tend_chem_settle_p)
if(allocated(frp_out_p) ) deallocate(frp_out_p )
if(allocated(ddvel_p) ) deallocate(ddvel_p )
if(allocated(chem_p) ) deallocate(chem_p )
if(allocated(tend_chem_settle_p) ) deallocate(tend_chem_settle_p )
if(allocated(frp_out_p) ) deallocate(frp_out_p )
if(allocated(aero_emis_for_enhmix_p)) deallocate(aero_emis_for_enhmix_p )
endif

pbl_select: select case (trim(pbl_scheme))
Expand Down Expand Up @@ -499,8 +501,12 @@ subroutine pbl_from_MPAS(configs,state,mesh,sfc_input,diag_physics,tend_physics,
real(kind=RKIND),dimension(:,: ),pointer :: ddvel
real(kind=RKIND),dimension(:,:,:),pointer :: chem
real(kind=RKIND),dimension(:),pointer :: frp_out
real(kind=RKIND),dimension(:),pointer :: aero_emis_for_enhmix
integer,pointer:: drydep_opt
logical,pointer:: enh_mix
character(len=StrKIND),pointer :: config_smoke_scheme, &
config_anthro_scheme, &
config_rwc_scheme

!local pointers for MYJ scheme:
real(kind=RKIND),dimension(:),pointer :: chlowq,thz0,qz0,uz0,vz0,ct,akhs,akms,lh,mixht
Expand Down Expand Up @@ -555,17 +561,42 @@ subroutine pbl_from_MPAS(configs,state,mesh,sfc_input,diag_physics,tend_physics,
call mpas_pool_get_dimension(state, 'chemistry_end', chemistry_end)
chem => scalars(chemistry_start:chemistry_end,:,:)
call mpas_pool_get_config(configs,'config_mynn_enh_mix',enh_mix)

if (enh_mix) then
call mpas_pool_get_array(diag_physics,'frp_out',frp_out)
end if
call mpas_pool_get_config(configs,'config_smoke_scheme',config_smoke_scheme)
if ( config_smoke_scheme .ne. 'off' ) then
call mpas_pool_get_array(diag_physics,'frp_out',frp_out)
do j = jts,jte
do i = its,ite
frp_out_p(i,j) = frp_out(i)
enddo
enddo
else
do j = jts,jte
do i = its,ite
frp_out_p(i,j) = 0._RKIND
enddo
enddo
endif

call mpas_pool_get_config(configs,'config_anthro_scheme',config_anthro_scheme)
call mpas_pool_get_config(configs,'config_rwc_scheme',config_rwc_scheme)
if ( config_anthro_scheme .ne. 'off' .or. config_rwc_scheme .ne. 'off') then
call mpas_pool_get_array(diag_physics,'aero_emis_for_enhmix',aero_emis_for_enhmix)
do j = jts,jte
do i = its,ite
aero_emis_for_enhmix_p(i,j) = aero_emis_for_enhmix(i)
enddo
enddo
else
do j = jts,jte
do i = its,ite
aero_emis_for_enhmix_p(i,j) = 0._RKIND
enddo
enddo
endif
endif ! If enh

do j = jts,jte
do i = its,ite
if ( enh_mix ) then
frp_out_p(i,j) = frp_out(i)
endif
enddo
enddo

do n = 1,num_chem
do j = jts,jte
Expand Down Expand Up @@ -1622,7 +1653,7 @@ subroutine driver_pbl(itimestep,configs,state,mesh,sfc_input,diag_physics,tend_p
settle3d = tend_chem_settle_p , &
enh_mix = enh_mix , &
frp_mean = frp_out_p , &
emis_ant_no = xland_p , &
emis_ant_no = aero_emis_for_enhmix_p , &
ids = ids , ide = ide , jds = jds , jde = jde , kds = kds , kde = kde , &
ims = ims , ime = ime , jms = jms , jme = jme , kms = kms , kme = kme , &
its = its , ite = ite , jts = jts , jte = jte , kts = kts , kte = kte , &
Expand Down
Loading