Skip to content

Claude/rh definition mixactivate ajs0 x #2266

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

Closed
liranpeng wants to merge 8 commits into from
+2,436 −150
Closed

Claude/rh definition mixactivate ajs0 x #2266

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
1e6d9da
Liran add extra output for parcel model
liranpeng Jan 3, 2025
4a9fa0d
add aerosol volume and lower/upper size limit
liranpeng Jan 6, 2025
be57f13
replace eq
liranpeng Apr 4, 2025
5a90cc4
Update module_mixactivate.F
liranpeng Oct 23, 2025
b7bd48f
switch mode
liranpeng Nov 10, 2025
78ce2ec
Update module_mixactivate.F
liranpeng Nov 17, 2025
aae37c5
Add namelist flag to control ML vs traditional aerosol activation
claude Jan 3, 2026
6944eaa
Fix missing config_flags in remaining mixactivate wrapper calls
claude Jan 5, 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 Registry/Registry.EM
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -20,6 +20,7 @@ state real SST_INPUT ij misc 1 - rh "SST_IN

# Chem Scalars
rconfig integer chem_opt namelist,physics max_domains 0 rh "chem_opt" "" ""
rconfig integer use_ml_activate namelist,physics max_domains 1 rh "use_ml_activate" "Use ML activation (1) or traditional activation (0)" ""
state real - ikjftb chem 1 - - -
#
# Tracer Scalars: If you want full resolved and non-resolved dispersion, compile WRF-Chem
Expand Down
132 changes: 126 additions & 6 deletions Registry/registry.chem
View file
Open in desktop

Large diffs are not rendered by default.

20 changes: 17 additions & 3 deletions chem/chem_driver.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1021,7 +1021,7 @@ end SUBROUTINE sum_pm_driver
end if
#endif
endif

call wrf_debug(15,'Liran start vertical mixing with dry deposition')
!
! do vertical mixing with dry deposition
!
Expand All @@ -1048,8 +1048,22 @@ end SUBROUTINE sum_pm_driver
call dry_dep_driver(grid%id,curr_secs,ktau,grid%dt,config_flags, &
grid%gmt,ijulian,t_phy,moist,scalar,p8w,t8w,vvel, &
rri,p_phy,chem,tracer,rho,dz8w,rh,grid%exch_h,grid%hfx,grid%dx, &
grid%cldfra, grid%cldfra_old,grid%raincv_b,seasin,dustin, &
grid%ccn1, grid%ccn2, grid%ccn3, grid%ccn4, grid%ccn5, grid%ccn6, &
grid%cldfra,grid%cldfra_old,grid%raincv_b,seasin,dustin, &
grid%fn11,grid%fn12,grid%fn13,grid%fn14,grid%fn21,grid%fn22,grid%fn23,grid%fn24,&
grid%fn31,grid%fn32,grid%fn33,grid%fn34,grid%fn41,grid%fn42,grid%fn43,grid%fn44,&
grid%na11,grid%na12,grid%na13,grid%na14,grid%na21,grid%na22,grid%na23,grid%na24,&
grid%na31,grid%na32,grid%na33,grid%na34,grid%na41,grid%na42,grid%na43,grid%na44,&
grid%vo11,grid%vo12,grid%vo13,grid%vo14,grid%vo21,grid%vo22,grid%vo23,grid%vo24,&
grid%vo31,grid%vo32,grid%vo33,grid%vo34,grid%vo41,grid%vo42,grid%vo43,grid%vo44,&
grid%dl11,grid%dl12,grid%dl13,grid%dl14,grid%dl21,grid%dl22,grid%dl23,grid%dl24,&
grid%dl31,grid%dl32,grid%dl33,grid%dl34,grid%dl41,grid%dl42,grid%dl43,grid%dl44,&
grid%dh11,grid%dh12,grid%dh13,grid%dh14,grid%dh21,grid%dh22,grid%dh23,grid%dh24,&
grid%dh31,grid%dh32,grid%dh33,grid%dh34,grid%dh41,grid%dh42,grid%dh43,grid%dh44,&
grid%sg11,grid%sg12,grid%sg13,grid%sg14,grid%sg21,grid%sg22,grid%sg23,grid%sg24,&
grid%sg31,grid%sg32,grid%sg33,grid%sg34,grid%sg41,grid%sg42,grid%sg43,grid%sg44,&
grid%hg11,grid%hg12,grid%hg13,grid%hg14,grid%hg21,grid%hg22,grid%hg23,grid%hg24,&
grid%hg31,grid%hg32,grid%hg33,grid%hg34,grid%hg41,grid%hg42,grid%hg43,grid%hg44,&
grid%ccn1, grid%ccn2, grid%ccn3, grid%ccn4, grid%ccn5, grid%ccn6, &
grid%qndropsource,grid%ivgtyp,grid%tsk,grid%gsw,grid%vegfra,pbl_h, &
grid%rmol,grid%ust,grid%znt,grid%xlat,grid%xlong, &
zmid,z_at_w,grid%xland,grid%ash_fall, &
Expand Down
98 changes: 90 additions & 8 deletions chem/dry_dep_driver.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -14,7 +14,21 @@ MODULE module_dry_dep_driver
subroutine dry_dep_driver(id,curr_secs,ktau,dtstep,config_flags, &
gmt,julday,t_phy,moist,scalar,p8w,t8w,w,alt, &
p_phy,chem,tracer,rho_phy,dz8w,rh,exch_h,hfx,dx, &
cldfra, cldfra_old,raincv,seasin,dustin, &
cldfra,cldfra_old,raincv,seasin,dustin,&
fn11,fn12,fn13,fn14,fn21,fn22,fn23,fn24,&
fn31,fn32,fn33,fn34,fn41,fn42,fn43,fn44,&
na11,na12,na13,na14,na21,na22,na23,na24,&
na31,na32,na33,na34,na41,na42,na43,na44,&
vo11,vo12,vo13,vo14,vo21,vo22,vo23,vo24,&
vo31,vo32,vo33,vo34,vo41,vo42,vo43,vo44,&
dl11,dl12,dl13,dl14,dl21,dl22,dl23,dl24,&
dl31,dl32,dl33,dl34,dl41,dl42,dl43,dl44,&
dh11,dh12,dh13,dh14,dh21,dh22,dh23,dh24,&
dh31,dh32,dh33,dh34,dh41,dh42,dh43,dh44,&
sg11,sg12,sg13,sg14,sg21,sg22,sg23,sg24,&
sg31,sg32,sg33,sg34,sg41,sg42,sg43,sg44,&
hg11,hg12,hg13,hg14,hg21,hg22,hg23,hg24,&
hg31,hg32,hg33,hg34,hg41,hg42,hg43,hg44,&
ccn1, ccn2, ccn3, ccn4, ccn5, ccn6, nsource, &
ivgtyp,tsk,gsw,vegfra,pbl,rmol,ust,znt,xlat,xlong,z,z_at_w,&
xland,ash_fall,h2oaj,h2oai,nu3,ac3,cor3,asulf,ahno3, &
Expand Down Expand Up @@ -144,7 +158,20 @@ subroutine dry_dep_driver(id,curr_secs,ktau,dtstep,config_flags, &

REAL, INTENT(OUT), dimension(ims:ime,kms:kme,jms:jme) :: nsource, &
ccn1,ccn2,ccn3,ccn4,ccn5,ccn6 ! number conc of aerosols activated at supersat

REAL, INTENT(OUT), dimension(ims:ime,kms:kme,jms:jme) :: fn11,fn12,fn13,fn14,fn21,fn22,fn23,fn24,&
fn31,fn32,fn33,fn34,fn41,fn42,fn43,fn44,&
na11,na12,na13,na14,na21,na22,na23,na24,&
na31,na32,na33,na34,na41,na42,na43,na44,&
vo11,vo12,vo13,vo14,vo21,vo22,vo23,vo24,&
vo31,vo32,vo33,vo34,vo41,vo42,vo43,vo44,&
dl11,dl12,dl13,dl14,dl21,dl22,dl23,dl24,&
dl31,dl32,dl33,dl34,dl41,dl42,dl43,dl44,&
dh11,dh12,dh13,dh14,dh21,dh22,dh23,dh24,&
dh31,dh32,dh33,dh34,dh41,dh42,dh43,dh44,&
sg11,sg12,sg13,sg14,sg21,sg22,sg23,sg24,&
sg31,sg32,sg33,sg34,sg41,sg42,sg43,sg44,&
hg11,hg12,hg13,hg14,hg21,hg22,hg23,hg24,&
hg31,hg32,hg33,hg34,hg41,hg42,hg43,hg44
REAL, INTENT(IN ) :: &
dtstep,gmt,dx

Expand Down Expand Up @@ -1068,18 +1095,31 @@ subroutine dry_dep_driver(id,curr_secs,ktau,dtstep,config_flags, &
end where

qsrflx(:,:,:) = 0.0

CALL wrf_debug(15,'Liran test aerosol')
mixactivate_select: SELECT CASE(config_flags%chem_opt)

CASE (RADM2SORG_AQ, RADM2SORG_AQCHEM, RACMSORG_AQ, RACMSORG_AQCHEM_KPP, RACM_ESRLSORG_AQCHEM_KPP, CBMZSORG_AQ, &
CB05_SORG_AQ_KPP)
CALL wrf_debug(15,'call mixactivate for sorgam aerosol')
CALL wrf_debug(15,'call mixactivate for LR sorgam aerosol')
call sorgam_mixactivate ( &
id, ktau, dtstep, config_flags, idrydep_onoff, &
dryrho_phy, t_phy, w, cldfra, cldfra_old, &
ddvel, z, dz8w, p8w, t8w, exch_h, &
moist(ims,kms,jms,P_QV), moist(ims,kms,jms,P_QC), moist(ims,kms,jms,P_QI), &
scalar(ims,kms,jms,P_QNDROP), f_qc, f_qi, chem, &
scalar(ims,kms,jms,P_QNDROP), f_qc, f_qi, chem,&
fn11,fn12,fn13,fn14,fn21,fn22,fn23,fn24,&
fn31,fn32,fn33,fn34,fn41,fn42,fn43,fn44,&
na11,na12,na13,na14,na21,na22,na23,na24,&
na31,na32,na33,na34,na41,na42,na43,na44,&
vo11,vo12,vo13,vo14,vo21,vo22,vo23,vo24,&
vo31,vo32,vo33,vo34,vo41,vo42,vo43,vo44,&
dl11,dl12,dl13,dl14,dl21,dl22,dl23,dl24,&
dl31,dl32,dl33,dl34,dl41,dl42,dl43,dl44,&
dh11,dh12,dh13,dh14,dh21,dh22,dh23,dh24,&
dh31,dh32,dh33,dh34,dh41,dh42,dh43,dh44,&
sg11,sg12,sg13,sg14,sg21,sg22,sg23,sg24,&
sg31,sg32,sg33,sg34,sg41,sg42,sg43,sg44,&
hg11,hg12,hg13,hg14,hg21,hg22,hg23,hg24,&
hg31,hg32,hg33,hg34,hg41,hg42,hg43,hg44,&
ccn1, ccn2, ccn3, ccn4, ccn5, ccn6, nsource, &
ids,ide, jds,jde, kds,kde, &
ims,ime, jms,jme, kms,kme, &
Expand All @@ -1091,7 +1131,21 @@ subroutine dry_dep_driver(id,curr_secs,ktau,dtstep,config_flags, &
dryrho_phy, t_phy, w, cldfra, cldfra_old, &
ddvel, z, dz8w, p8w, t8w, exch_h, &
moist(ims,kms,jms,P_QV), moist(ims,kms,jms,P_QC),moist(ims,kms,jms,P_QI), &
scalar(ims,kms,jms,P_QNDROP), f_qc, f_qi, chem, &
scalar(ims,kms,jms,P_QNDROP), f_qc, f_qi, chem,&
fn11,fn12,fn13,fn14,fn21,fn22,fn23,fn24,&
fn31,fn32,fn33,fn34,fn41,fn42,fn43,fn44,&
na11,na12,na13,na14,na21,na22,na23,na24,&
na31,na32,na33,na34,na41,na42,na43,na44,&
vo11,vo12,vo13,vo14,vo21,vo22,vo23,vo24,&
vo31,vo32,vo33,vo34,vo41,vo42,vo43,vo44,&
dl11,dl12,dl13,dl14,dl21,dl22,dl23,dl24,&
dl31,dl32,dl33,dl34,dl41,dl42,dl43,dl44,&
dh11,dh12,dh13,dh14,dh21,dh22,dh23,dh24,&
dh31,dh32,dh33,dh34,dh41,dh42,dh43,dh44,&
sg11,sg12,sg13,sg14,sg21,sg22,sg23,sg24,&
sg31,sg32,sg33,sg34,sg41,sg42,sg43,sg44,&
hg11,hg12,hg13,hg14,hg21,hg22,hg23,hg24,&
hg31,hg32,hg33,hg34,hg41,hg42,hg43,hg44,&
ccn1, ccn2, ccn3, ccn4, ccn5, ccn6, nsource, &
ids,ide, jds,jde, kds,kde, &
ims,ime, jms,jme, kms,kme, &
Expand All @@ -1104,6 +1158,20 @@ subroutine dry_dep_driver(id,curr_secs,ktau,dtstep,config_flags, &
ddvel, z, dz8w, p8w, t8w, exch_h, &
moist(ims,kms,jms,P_QV), moist(ims,kms,jms,P_QC), moist(ims,kms,jms,P_QI), &
scalar(ims,kms,jms,P_QNDROP), f_qc, f_qi, chem, &
fn11,fn12,fn13,fn14,fn21,fn22,fn23,fn24,&
fn31,fn32,fn33,fn34,fn41,fn42,fn43,fn44,&
na11,na12,na13,na14,na21,na22,na23,na24,&
na31,na32,na33,na34,na41,na42,na43,na44,&
vo11,vo12,vo13,vo14,vo21,vo22,vo23,vo24,&
vo31,vo32,vo33,vo34,vo41,vo42,vo43,vo44,&
dl11,dl12,dl13,dl14,dl21,dl22,dl23,dl24,&
dl31,dl32,dl33,dl34,dl41,dl42,dl43,dl44,&
dh11,dh12,dh13,dh14,dh21,dh22,dh23,dh24,&
dh31,dh32,dh33,dh34,dh41,dh42,dh43,dh44,&
sg11,sg12,sg13,sg14,sg21,sg22,sg23,sg24,&
sg31,sg32,sg33,sg34,sg41,sg42,sg43,sg44,&
hg11,hg12,hg13,hg14,hg21,hg22,hg23,hg24,&
hg31,hg32,hg33,hg34,hg41,hg42,hg43,hg44,&
ccn1, ccn2, ccn3, ccn4, ccn5, ccn6, nsource, &
ids,ide, jds,jde, kds,kde, &
ims,ime, jms,jme, kms,kme, &
Expand All @@ -1118,7 +1186,21 @@ subroutine dry_dep_driver(id,curr_secs,ktau,dtstep,config_flags, &
dryrho_phy, t_phy, w, cldfra, cldfra_old, &
ddvel, z, dz8w, p8w, t8w, exch_h, &
moist(ims,kms,jms,P_QV), moist(ims,kms,jms,P_QC), moist(ims,kms,jms,P_QI), &
scalar(ims,kms,jms,P_QNDROP), f_qc, f_qi, chem, &
scalar(ims,kms,jms,P_QNDROP), f_qc, f_qi, chem,&
fn11,fn12,fn13,fn14,fn21,fn22,fn23,fn24,&
fn31,fn32,fn33,fn34,fn41,fn42,fn43,fn44,&
na11,na12,na13,na14,na21,na22,na23,na24,&
na31,na32,na33,na34,na41,na42,na43,na44,&
vo11,vo12,vo13,vo14,vo21,vo22,vo23,vo24,&
vo31,vo32,vo33,vo34,vo41,vo42,vo43,vo44,&
dl11,dl12,dl13,dl14,dl21,dl22,dl23,dl24,&
dl31,dl32,dl33,dl34,dl41,dl42,dl43,dl44,&
dh11,dh12,dh13,dh14,dh21,dh22,dh23,dh24,&
dh31,dh32,dh33,dh34,dh41,dh42,dh43,dh44,&
sg11,sg12,sg13,sg14,sg21,sg22,sg23,sg24,&
sg31,sg32,sg33,sg34,sg41,sg42,sg43,sg44,&
hg11,hg12,hg13,hg14,hg21,hg22,hg23,hg24,&
hg31,hg32,hg33,hg34,hg41,hg42,hg43,hg44,&
ccn1, ccn2, ccn3, ccn4, ccn5, ccn6, nsource, &
qsrflx, &
ids,ide, jds,jde, kds,kde, &
Expand Down
Loading