Adding Ristic microphysics scheme #1619

davegill · 2021-12-25T18:49:55Z

TYPE: New feature

KEYWORDS: microphysic, forecast, cloud, precipitation

SOURCE: Ivan Ristic (Weather2 - www.weather2.rs)

DESCRIPTION OF CHANGES:
Original PR Adding new Ristic microphysics scheme (#1308) was more than a year out of date, and the effort to "merge everything in" did not go well. This PR started with the develop branch and cherry-picked the two (non-NMM) commits from the original Ristic PR.

In order to improve cloud and precipitation forecast we developed new cloud prediction scheme and we implemented it in WRF model. Fractional cloud cover, cloud liquid water, cloud ice and cloud snow are explicitly predicted by adding three
prognostic equations for fractional cloud cover, cloud mixing ratio and snow per cloud fraction to the model. Sedimentation of ice and snow is also included in parameterization. Precipitation of rain and snow are determined from cloud fields. Clouds predicted like this can be used also in radiation parameterization.

Thermodynamic wet bulb temperature will be used for describing clouds because it is constant during water phase changes. By using this temperature moist static energy of model grid box and cloudy part inside the grid box is the same and principle of energy conservation is satisfied.

A complete description is now found in Ristic I., Kordic I., 2018: Cloud parameterization and cloud prediction scheme in the Eta numerical weather model. NWCM - Serbian Academy of Sciences and Arts which can be found at https://www.sanu.ac.rs/wp-content/uploads/2018/10/11_Cloud%20parameterization%20and%20cloud%20prediction%20scheme%20in%20Eta%20numerical%20weather%20model.pdf .

LIST OF MODIFIED FILES:
modified: Registry/Registry.EM_COMMON
modified: dyn_em/solve_em.F
modified: phys/Makefile
modified: phys/module_microphysics_driver.F
new file: phys/module_mp_ivanr_micro.F

TESTS CONDUCTED:
All the wrf-coop tests passed and the microphysics is being actively used for running a model with NMM.
Integration of the model for test cases indicate that new cloud prediction scheme improved forecast compared to the original model. New fractional cloud cover formula showed good results in practice, since the fractional cloud cover, predicted in this way, was much closer to the real cloud cover values. Significant progress has been made in stratiform precipitation forecast. Positive impact on convection scheme is also noticed.
One test was run on 17. November 2011 with the Wrf model with ECMWF as boundary data. Model run for 72 hours in horizontal resolution of about 22 km and vertical resolution of 38 layers for the Europe domain. Fog was first to be tested.

The second test was run on 24. June 2015 with Wrf model with ECMWF as boundary data. Model was run for 24 hours in
horizontal resolution of about 22 km and a vertical resolution of 38 layers for the Europe domain. The second test situation tested mid-morning precipitation over northern Serbia.

More details can be found in the aforementioned pdf file.

RELEASE NOTE: Ivan Ristic: Custom implementation of microphysics

TYPE: New feature KEYWORDS: microphysic, forecast, cloud, precipitation SOURCE: Ivan Ristic (Weather2 - www.weather2.rs) DESCRIPTION OF CHANGES: In order to improve cloud and precipitation forecast we developed new cloud prediction scheme and we implemented it in WRF model. Fractional cloud cover, cloud liquid water, cloud ice and cloud snow are explicitly predicted by adding three prognostic equations for fractional cloud cover, cloud mixing ratio and snow per cloud fraction to the model. Sedimentation of ice and snow is also included in parameterization. Precipitation of rain and snow are determined from cloud fields. Clouds predicted like this can be used also in radiation parameterization. Thermodynamic wet bulb temperature will be used for describing clouds because it is constant during water phase changes. By using this temperature moist static energy of model grid box and cloudy part inside the grid box is the same and principle of energy conservation is satisfied. A complete description is now found in Ristic I., Kordic I., 2018: Cloud parameterization and cloud prediction scheme in the Eta numerical weather model. NWCM - Serbian Academy of Sciences and Arts which can be found at https://www.sanu.ac.rs/wp-content/uploads/2018/10/11_Cloud%20parameterization%20and%20cloud%20prediction%20scheme%20in%20Eta%20numerical%20weather%20model.pdf . LIST OF MODIFIED FILES: modified: Registry/Registry.EM_COMMON modified: dyn_em/solve_em.F modified: phys/Makefile modified: phys/module_microphysics_driver.F new file: phys/module_mp_ivanr_micro.F TESTS CONDUCTED: All the wrf-coop tests passed and the microphysics is being actively used for running a model with NMM. Integration of the model for test cases indicate that new cloud prediction scheme improved forecast compared to the original model. New fractional cloud cover formula showed good results in practice, since the fractional cloud cover, predicted in this way, was much closer to the real cloud cover values. Significant progress has been made in stratiform precipitation forecast. Positive impact on convection scheme is also noticed. One test was run on 17. November 2011 with the Wrf model with ECMWF as boundary data. Model run for 72 hours in horizontal resolution of about 22 km and vertical resolution of 38 layers for the Europe domain. Fog was first to be tested. ![image](https://user-images.githubusercontent.com/51482208/97575741-7b35d080-19ed-11eb-8310-fc6a6e3c5ac8.png) ![image](https://user-images.githubusercontent.com/51482208/97576039-f303fb00-19ed-11eb-8410-7736be350b6d.png) The second test was run on 24. June 2015 with Wrf model with ECMWF as boundary data. Model was run for 24 hours in horizontal resolution of about 22 km and a vertical resolution of 38 layers for the Europe domain. The second test situation tested mid-morning precipitation over northern Serbia. ![image](https://user-images.githubusercontent.com/51482208/97577542-dff22a80-19ef-11eb-9a86-e4d3b0a05e10.png) More details can be found in the aforementioned pdf file. RELEASE NOTE: Ivan Ristic: Custom implementation of microphysics

davegill · 2021-12-25T19:40:30Z

@Nikolar97
Nikola,
Would you please verify that this PR still is able to work correctly. Since your original PR only had two small commits, they were easier to handle with a "git cherry-pick".

davegill · 2021-12-25T19:41:06Z

jenkins

Please find result of the WRF regression test cases in the attachment. This build is for Commit ID: 49e02f166c778bbbd7c39c757aca13b71697baff, requested by: davegill for PR: https://github.com/scala-computing/WRF/pull/1619. For any query please send e-mail to David Gill.

    Test Type              | Expected  | Received |  Failed
    = = = = = = = = = = = = = = = = = = = = = = = =  = = = =
    Number of Tests        : 23           24
    Number of Builds       : 60           58
    Number of Simulations  : 159           155        0
    Number of Comparisons  : 96           93        0

    Failed Simulations are: 
    None
    Which comparisons are not bit-for-bit: 
    None

davegill · 2021-12-25T23:06:56Z

@Nikolar97 @weiwangncar @dudhia
Is this scheme suitable for use with ARW? We have removed the NMM dynamical core. We will not merge in a MP scheme designed specifically for NMM.

dudhia · 2021-12-25T23:15:27Z

This is not for V4.4. It can be made compatible with V4.3.2 if they need NMM.

davegill · 2021-12-25T23:58:43Z

@Nikolar97 @dudhia @weiwangncar

This is not for V4.4. It can be made compatible with V4.3.2 if they need NMM.

Jimy,
Thanks! I'll close this PR.

Nikolar97 · 2021-12-26T18:44:28Z

We would like to to make it compatible with ARW , is it enough for us to put our CLD variable into the existing CLDFRAC array. Do you have any pointers how we could go along with these changes, it would be greatly appreciated. Thanks in advance.

dudhia · 2021-12-27T08:21:03Z

I have not looked at the code, but use the existing variables as much as possible. CLDFRAC is used by the radiation schemes.

…

On Sun, Dec 26, 2021 at 11:44 AM Nikola Ristic ***@***.***> wrote: We would like to to make it compatible with ARW , is it enough for us to put our CLD variable into the existing CLDFRAC array. Do you have any pointers how we could go along with these changes, it would be greatly appreciated. Thanks in advance. — Reply to this email directly, view it on GitHub <#1619 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AEIZ77G4QDS4JXONTKZPMITUS5PBPANCNFSM5KYB2M7Q> . Triage notifications on the go with GitHub Mobile for iOS <https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675> or Android <https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub>. You are receiving this because you were mentioned.Message ID: ***@***.***>

davegill · 2021-12-29T16:05:23Z

Apparently, we are OPEN for business!

Nikolar97 · 2022-01-29T12:51:42Z

In Ristic microphysics, the cloud fraction is a prognostic quantity. This is a
new approach and we need advice on how to best implement
the cloud fraction in the radiation driver module_radiation_driver.F . We need
brief advice on how to insert a CLD variable from Ristic microphysics into
CLDFRA variable.
The problem is also the ICLOUD variable, whether to assign a new value to be
used by our microphysics or to use one of the existing ones and just add the
part related to the equalization of CLDFRC and CLD variable.
With these changes, Ristic microphysics will be totally compatible with all
radiation parameterizations and the ARW core, which is the goal.
Is there anyone we can talk to about this topic?

davegill · 2022-01-29T19:07:27Z

@weiwangncar @dudhia
Ristic questions

weiwangncar · 2022-01-29T20:02:13Z

@Nikolar97 We will get back to you soon.

Nikolar97 · 2022-04-12T09:25:49Z

Hello, we are currently waiting for an answer in regards to our last question so we can try to continue with our integration of Ristic microphysics with the ARW core. Is there any update on the current situation?

dudhia · 2022-04-12T15:14:25Z

Our plan with this is to have it as a feature branch. We are still finalizing our V4.4 release this month and after that you may be able to upgrade and have an up-to-date branch. Jimy

…

On Tue, Apr 12, 2022 at 3:35 AM Nikola Ristic ***@***.***> wrote: Hello, we are currently waiting for an answer in regards to our last question so we can try to continue with our integration of Ristic microphysics with the ARW core. Is there any update on the current situation? — Reply to this email directly, view it on GitHub <#1619 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AEIZ77BCHZOW7CVW5H2ZKDDVEU7U7ANCNFSM5KYB2M7Q> . You are receiving this because you were mentioned.Message ID: ***@***.***>

Nikolar97 · 2022-04-12T18:45:01Z

@dudhia Thank you for your answer, we will wait for the new release to come out and then start implementing our microphysics.

Nikolar97 · 2022-05-22T18:16:17Z

@davegill @dudhia Hello, we have finished porting our microphysics to the V4.4 release and successfully tested. What are the next steps we should take to have it as a feature branch?

weiwangncar · 2022-05-23T16:14:48Z

@Nikolar97 Thank you for updating your branch to v4.4. Let us discuss some ways to maintain your code outside the main repository, and get back to you, hopefully soon.

Nikolar97 · 2022-08-03T16:19:26Z

@weiwangncar @davegill @dudhia Hello, we would like to know if there are any updates on the current situation. We are eager to continue our discussions on the maintenance of the branch.

dudhia · 2022-08-04T15:27:16Z

We are looking at how best to do this.

…

On Wed, Aug 3, 2022 at 10:19 AM Nikola Ristic ***@***.***> wrote: @weiwangncar <https://github.com/weiwangncar> @davegill <https://github.com/davegill> @dudhia <https://github.com/dudhia> Hello, we would like to know if there are any updates on the current situation. We are eager to continue our discussions on the maintenance of the branch. — Reply to this email directly, view it on GitHub <#1619 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AEIZ77DWOIS3Q7ZUO6G6RE3VXKLZTANCNFSM5KYB2M7Q> . You are receiving this because you were mentioned.Message ID: ***@***.***>

dudhia · 2022-08-04T18:50:09Z

@Nikolar97 do you have your updated code in your own github fork?

Nikolar97 · 2022-08-12T17:46:18Z

@dudhia Sorry for replaying so late, we were on a holiday vacation. We updated the current fork to the newest WRF version and then pushed our changes on it(with the new microphysics integrated with the 4.4 WRF).https://github.com/Nikolar97/WRF

dudhia · 2022-08-12T18:11:44Z

Maybe you can do a new pull request, which will also test the compilation. If you have your own fork, it is preferable to do it from there. I am on vacation now, so my replies may be slow. Jimy

…

On Fri, Aug 12, 2022 at 11:46 AM Nikola Ristic ***@***.***> wrote: @dudhia <https://github.com/dudhia> Sorry for replaying so late, we were on a holiday vacation. We updated the current fork to the newest WRF version and then pushed our changes on it(with the new microphysics integrated with the 4.4 WRF). — Reply to this email directly, view it on GitHub <#1619 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AEIZ77D7454SRPIMFHW3UFTVY2EXJANCNFSM5KYB2M7Q> . You are receiving this because you were mentioned.Message ID: ***@***.***>

Nikolar97 · 2022-08-12T18:24:45Z

@dudhia On which branch should i do my new pull request if it matters?

dudhia · 2022-08-12T21:11:04Z

Since it is just for testing you can use the master branch which is the latest release.

…

On Fri, Aug 12, 2022 at 12:24 PM Nikola Ristic ***@***.***> wrote: @dudhia <https://github.com/dudhia> On which branch should i do my new pull request if it matters? — Reply to this email directly, view it on GitHub <#1619 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AEIZ77CBSRPKCYN6CRNCDZTVY2JHRANCNFSM5KYB2M7Q> . You are receiving this because you were mentioned.Message ID: ***@***.***>

Nikolar97 · 2022-08-13T17:20:51Z

@dudhia Created the pull request for testing at #1757.

dudhia · 2023-01-26T20:45:40Z

Superseded by #1757

davegill requested review from a team as code owners December 25, 2021 18:49

davegill added PHYSICS New Feature Develop Branch Await WPRP labels Dec 25, 2021

Adding descriptions of custom fields

d5a459d

davegill mentioned this pull request Dec 25, 2021

Adding new Ristic microphysics scheme #1308

Closed

davegill closed this Dec 25, 2021

davegill reopened this Dec 29, 2021

davegill removed the Develop Branch label Jan 6, 2022

dudhia closed this Jan 26, 2023

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Adding Ristic microphysics scheme #1619

Adding Ristic microphysics scheme #1619

davegill commented Dec 25, 2021 •

edited

Loading

davegill commented Dec 25, 2021

davegill commented Dec 25, 2021

davegill commented Dec 25, 2021

dudhia commented Dec 25, 2021

davegill commented Dec 25, 2021

Nikolar97 commented Dec 26, 2021

dudhia commented Dec 27, 2021 via email

davegill commented Dec 29, 2021

Nikolar97 commented Jan 29, 2022

davegill commented Jan 29, 2022

weiwangncar commented Jan 29, 2022

Nikolar97 commented Apr 12, 2022

dudhia commented Apr 12, 2022 via email

Nikolar97 commented Apr 12, 2022

Nikolar97 commented May 22, 2022 •

edited

Loading

weiwangncar commented May 23, 2022

Nikolar97 commented Aug 3, 2022

dudhia commented Aug 4, 2022 via email

dudhia commented Aug 4, 2022

Nikolar97 commented Aug 12, 2022 •

edited

Loading

dudhia commented Aug 12, 2022 via email

Nikolar97 commented Aug 12, 2022

dudhia commented Aug 12, 2022 via email

Nikolar97 commented Aug 13, 2022

dudhia commented Jan 26, 2023

Adding Ristic microphysics scheme #1619

Adding Ristic microphysics scheme #1619

Conversation

davegill commented Dec 25, 2021 • edited Loading

davegill commented Dec 25, 2021

davegill commented Dec 25, 2021

davegill commented Dec 25, 2021

dudhia commented Dec 25, 2021

davegill commented Dec 25, 2021

Nikolar97 commented Dec 26, 2021

dudhia commented Dec 27, 2021 via email

davegill commented Dec 29, 2021

Nikolar97 commented Jan 29, 2022

davegill commented Jan 29, 2022

weiwangncar commented Jan 29, 2022

Nikolar97 commented Apr 12, 2022

dudhia commented Apr 12, 2022 via email

Nikolar97 commented Apr 12, 2022

Nikolar97 commented May 22, 2022 • edited Loading

weiwangncar commented May 23, 2022

Nikolar97 commented Aug 3, 2022

dudhia commented Aug 4, 2022 via email

dudhia commented Aug 4, 2022

Nikolar97 commented Aug 12, 2022 • edited Loading

dudhia commented Aug 12, 2022 via email

Nikolar97 commented Aug 12, 2022

dudhia commented Aug 12, 2022 via email

Nikolar97 commented Aug 13, 2022

dudhia commented Jan 26, 2023

davegill commented Dec 25, 2021 •

edited

Loading

Nikolar97 commented May 22, 2022 •

edited

Loading

Nikolar97 commented Aug 12, 2022 •

edited

Loading