Hi,
I am currently working on a WRF-Chem experiment to evaluate the impact of individual dust size bins on cloud microphysical properties. I intend to use chem_opt = 202, and I plan to provide IC/BC data using the mozbc utility.
However, I am a bit confused about how to properly represent the dust size bins in the spc_map file. Specifically, I would like to know if the following approach is valid for distributing CESM/MAM4 dust (dst_a3) into GOCART bins:
spc_map =
'BC1 -> 1.0*bc_a4;1.e9',
'BC2 -> 1.0*bc_a1;1.e9',
'OC1 -> 1.0*pom_a4;1.e9',
'OC2 -> 1.0*pom_a1;1.e9',
'SEAS_1 -> 1.0*ncl_a1+1.0*ncl_a2;1.e9',
'SEAS_2 -> 0.5*ncl_a3;1.e9',
'SEAS_3 -> 0.5*ncl_a3;1.e9',
'SEAS_4 -> 0.0*ncl_a3;1.e9',
'DUST_1 -> 0.30*dst_a3;1.e9',
'DUST_2 -> 0.40*dst_a3;1.e9',
'DUST_3 -> 0.20*dst_a3;1.e9',
'DUST_4 -> 0.07*dst_a3;1.e9',
'DUST_5 -> 0.03*dst_a3;1.e9'
/
Additionally, for sensitivity experiments, I am considering isolating each bin by modifying the mapping, e.g., setting DUST_1 -> 1.0*dst_a3;1.e9 and all others to zero. Would this approach be scientifically sound for assessing the contribution of each size bin?
I would greatly appreciate any guidance on this or any recommendations you might have to better implement such experiments.
I am currently working on a WRF-Chem experiment to evaluate the impact of individual dust size bins on cloud microphysical properties. I intend to use chem_opt = 202, and I plan to provide IC/BC data using the mozbc utility.
However, I am a bit confused about how to properly represent the dust size bins in the spc_map file. Specifically, I would like to know if the following approach is valid for distributing CESM/MAM4 dust (dst_a3) into GOCART bins:
spc_map =
'BC1 -> 1.0*bc_a4;1.e9',
'BC2 -> 1.0*bc_a1;1.e9',
'OC1 -> 1.0*pom_a4;1.e9',
'OC2 -> 1.0*pom_a1;1.e9',
'SEAS_1 -> 1.0*ncl_a1+1.0*ncl_a2;1.e9',
'SEAS_2 -> 0.5*ncl_a3;1.e9',
'SEAS_3 -> 0.5*ncl_a3;1.e9',
'SEAS_4 -> 0.0*ncl_a3;1.e9',
'DUST_1 -> 0.30*dst_a3;1.e9',
'DUST_2 -> 0.40*dst_a3;1.e9',
'DUST_3 -> 0.20*dst_a3;1.e9',
'DUST_4 -> 0.07*dst_a3;1.e9',
'DUST_5 -> 0.03*dst_a3;1.e9'
/
Additionally, for sensitivity experiments, I am considering isolating each bin by modifying the mapping, e.g., setting DUST_1 -> 1.0*dst_a3;1.e9 and all others to zero. Would this approach be scientifically sound for assessing the contribution of each size bin?
I would greatly appreciate any guidance on this or any recommendations you might have to better implement such experiments.