Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

Question about aerosol-aware thompson scheme with use_aero_icbc = .false.

yutinghe

New member
Hi all,
I am running WRF model (version 4.2.2) using aerosol-aware thompson scheme (mp_physics=28) with use_aero_icbc = .false. In fact, I want to use this microphysics scheme in the context of a very low aerosol number concentration. Thus, I modify the naCCN0, naCCN1, naIN0 and naIN1 in module_mp_thompson.F as follows:
Code:
! REAL, PARAMETER, PRIVATE:: naIN0 = 1.5E6
! REAL, PARAMETER, PRIVATE:: naIN1 = 0.5E6
! REAL, PARAMETER, PRIVATE:: naCCN0 = 300.0E6
! REAL, PARAMETER, PRIVATE:: naCCN1 = 50.0E6
REAL, PARAMETER, PRIVATE:: naIN0 = 3.0
REAL, PARAMETER, PRIVATE:: naIN1 = 1.0
REAL, PARAMETER, PRIVATE:: naCCN0 = 6.0
REAL, PARAMETER, PRIVATE:: naCCN1 = 1.0
However, I found that the QNWFA variable in wrfout files only at the first output interval with my modified value (i.e. less than 6), and the rests still have large values (about 7 powers of 10). So what should I do to make the aerosol number concentration at each output interval is a very small value?
Besides, I have another question: I found that the QNWFA2D variable in wrfout files only at the first output interval with values, the other intervals are zero, is that correct?

Here I pasted some pieces of my namelist.input:
&physics
mp_physics = 28, 28, 28, 28, 28,
use_aero_icbc = .false.
aer_opt = 3
cu_physics = 1, 0, 0, 0, 0,
ra_lw_physics = 4, 4, 4, 4, 4,
ra_sw_physics = 1, 1, 1, 1, 1,
bl_pbl_physics = 1, 1, 1, 1, 1,
sf_sfclay_physics = 1, 1, 1, 1, 1,
sf_surface_physics = 2, 2, 2, 2, 2,
radt = 10, 10, 10, 10, 10,
bldt = 0, 0, 0, 0, 0,
cudt = 0, 0, 0, 0, 0,
icloud = 1,
num_land_cat = 21,
sf_urban_physics = 1, 1, 1, 1, 1,
sst_update = 1,
tmn_update = 1,
/

I will appreciate any help.
Yuting
 
Last edited:
Top