Skip to content

P3 microphysics scheme logic#139

Merged
isilber merged 25 commits intoARM-DOE:masterfrom
isilber:p3_dev
May 15, 2025
Merged

P3 microphysics scheme logic#139
isilber merged 25 commits intoARM-DOE:masterfrom
isilber:p3_dev

Conversation

@isilber
Copy link
Copy Markdown
Collaborator

@isilber isilber commented May 15, 2025

This PR adds the P3 microphysics scheme logic to both the microphysics and radiation scheme logic approaches of EMC².
The implementation includes the addition of new LUTs consistent in essence with those used by E3SMv3.
Integrating P3 ice treatment generally reduces the microphysics-radiation gap and enables robust application of EMC² to models utilizing this microphysics approach.

This PR also includes:

  1. Potential bug fix in cases of truncated-value mixing ratios, potentially resulting in bogus subcolumn values.
  2. Some metadata corrections and fixes.

MNT: refactoring of terminal velocity parameterization selection
For a given instrument, all LUTs are stored in single NetCDF file
Also, set `mie_for_ice` to `False` by default if using `P3` since ice shape
is baked into the microphysics, which impacts PSD parameters and
hence the radiation scheme treatment.
…microphys to liq processing

Note: we did not implement ice treatment yet, conditioning on the `mcphys_scheme` `Model` class attribute, which will impact both scattering efficiencies and terminal velocity calculations.
That module was redundant for EMC^2 practices (following the production of the LUTs)
This will cut calculation time later.
This commit implements the P3 physics logic (V_t parameterization,
rime properties, etc.) in radar moment calculations.
MNT: some code organization including generalizing PSD parameter
calculations under a single method.
The MH2005 terminal velocity implementation cannot resolve diameters < 23 um.
This fix fills those invalid values with the minimum resolved velocity (~1.7 cm/s)
This commit implements the P3 physics logic (rime properties, etc.)
in lidar calculations.
Conflicts:
	emc2/simulator/lidar_moments.py
	emc2/simulator/psd.py
	emc2/simulator/radar_moments.py
In detail:
ENH: P3 in lidar calculations (microphysics)
ENH: "hybrid" approach in P3 calculations (using equivalent volume diameter and surface area for scattering calculations)
ENH: enable Mie calculations in "P3" mode (while still using P3 parameterization for Vt calculations)
FIX: Ni normalization to 1 m-3 instead of 1 kg-1
DOC: update metadata and correct some print messages to user
This commit includes updates to the P3-related LUTs to enable Mie option in bulk calculations
As part of the potential bug fix in cases of truncated-value mixing ratios,
reference to the input hydrometeor type (conv/strat) was removed and had to be
restored for cases with deep-convection model output
@isilber isilber merged commit 6143afc into ARM-DOE:master May 15, 2025
7 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant