Skip to content

Conversation

trontrytel
Copy link
Member

Hi Folks! This PR adds 2-moment microphysics to diagnostic EDMF.

@sajjadazimi - could I ask you to take a look? I wasn't sure about two things:

  1. What velocity I should be using in the aerosol activation.
  2. How to interact with global aerosol data.

For the second point I temporarily hijacked two parameters: seasalt_calibration_coefficient, reference_seasalt_aerosol_mass_concentration to allow the users to prescribe the assumed mode size and concentration for seasalt aerosol. Once we make the breaking release of CliMA Params I will add proper parameters for it.

We do get some clouds and drizzle and the results differ for high and low concentrations. But in the couple of configurations I tried we have been raining out a lot of the water very early on.

@trontrytel trontrytel force-pushed the aj/diagnostic_edmf_2M branch from be0943e to a25436c Compare August 20, 2025 23:39
@trontrytel trontrytel marked this pull request as ready for review August 20, 2025 23:40
Copy link
Member

@sajjadazimi sajjadazimi left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thank you! LGTM!
The velocity passed to aerosol activation seems correct. And for reading aerosol data we need to add the list of aerosols to yml and then use the function compute_prescribed_aerosol_properties! from microphysics wrappers .

dt,
)
# Aerosol activation
# TODO - right now I'm hijacking the
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

We can read prescribed aerosol data from artifacts (Merra-2 dataset). For example the function microphysics_wrappers/compute_prescribed_aerosol_properties! reads sea salt and sulfate data and computes their number concentrations and the mean radius of seasalt. In order to use the data we need to pass a list of aerosols in the yml file.

max(
0,
w_component.(
Geometry.WVector.(u³ʲ_data_prev_halflevel)
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think this is correct! We need to pass the velocity component along gravity and I think this is it.

Copy link
Member

@oalcabes oalcabes left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

The bits of this that I understand look good to me!

Copy link
Contributor

@ericmpham ericmpham left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I don't have much to add, I think it looks okay. 😄

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

4 participants