Skip to content

add semi-Lagrangian sedimentation of graupel into Thompson MP#770

Merged
climbfuji merged 5 commits into
NCAR:mainfrom
RuiyuSun:semi-lagrangian_sedi_graupel
Nov 23, 2021
Merged

add semi-Lagrangian sedimentation of graupel into Thompson MP#770
climbfuji merged 5 commits into
NCAR:mainfrom
RuiyuSun:semi-lagrangian_sedi_graupel

Conversation

@RuiyuSun
Copy link
Copy Markdown
Contributor

@RuiyuSun RuiyuSun commented Nov 5, 2021

This PR is to add semi-Lagrangian sedimentation of graupel into the Thompson microphysics scheme.

Copy link
Copy Markdown
Collaborator

@climbfuji climbfuji left a comment

Choose a reason for hiding this comment

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

This looks good to me, will approve after the regression testing is done. To confirm, this will not change the answers of the existing regression tests, because the feature is turned off by default?

@RuiyuSun
Copy link
Copy Markdown
Contributor Author

RuiyuSun commented Nov 5, 2021 via email

Copy link
Copy Markdown
Contributor

@gthompsnWRF gthompsnWRF left a comment

Choose a reason for hiding this comment

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

Very minor changes to names requested, otherwise, I'm fine with this.

Comment thread physics/module_mp_thompson.F90
Comment thread physics/module_mp_thompson.F90 Outdated
nr_tmp(:) = nr(:)
call nislfv_rain_ppm(kte,dzq,vtrk,rr,rainsfc,dtcfl,R1)
call nislfv_rain_ppm(kte,dzq,vtnrk,nr,vtr,dtcfl,R2)
call nislfv_precip_ppm(kte,dzq,vtrk,rr,rainsfc,dtcfl,R1)
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

While you changed "rain" to "precip" I really don't understand "nislfv" nor "ppm" and would appreciate something better. These are jumbled letters to me and others. Got anything else?

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

nislfv and ppm are described in the description section of the routine.
nislfv= non-iteration semi-Lagrangain forward advection
ppm = 2nd order interpolation with monotonic piecewise parabolic method

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

Good to know it's documented somehow, but I see it as an improvement to use something like subroutine semi_lagrange_precip (or anything to give a clue by name in the calling routine. It's great to have the more full descript inside the routine.

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

What do you suggest using? Any idea?

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

@Songyou184 What do you think about the names?

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

What do you suggest using? Any idea?

I did give you my suggested name of subroutine in my comment. It is: semi_lagrange_precip or perhaps better is semi_lagrange_sedim

@Songyou184
Copy link
Copy Markdown

Songyou184 commented Nov 6, 2021 via email

@RuiyuSun
Copy link
Copy Markdown
Contributor Author

RuiyuSun commented Nov 6, 2021 via email

@RuiyuSun
Copy link
Copy Markdown
Contributor Author

RuiyuSun commented Nov 6, 2021

The latter, semi_lagrange_sedim, is my preference. I wonder if my previous comments in git environment were delivered correctly. In addition to this naming, in 1st line of the description in the subroutine, ....for cloud is to be replaced with ....for falling hydrometeors. Also, the corresponding reference is to be added. Juang, H.-M., and S.-Y. Hong, 2010: Forward semi-Lagrangian advection with mass conservation and positive definiteness for falling hydrometeors. Mon. Wea. Rev., 138, 1778-1791.

On Sat, Nov 6, 2021 at 9:41 AM gthompsnWRF @.> wrote: @.* commented on this pull request. ------------------------------ In physics/module_mp_thompson.F90 <#770 (comment)>: > @@ -3937,11 +3937,12 @@ subroutine mp_thompson (qv1d, qc1d, qi1d, qr1d, qs1d, qg1d, ni1d, & do n = 1, niter rr_tmp(:) = rr(:) nr_tmp(:) = nr(:) - call nislfv_rain_ppm(kte,dzq,vtrk,rr,rainsfc,dtcfl,R1) - call nislfv_rain_ppm(kte,dzq,vtnrk,nr,vtr,dtcfl,R2) + call nislfv_precip_ppm(kte,dzq,vtrk,rr,rainsfc,dtcfl,R1) What do you suggest using? Any idea? I did give you my suggested name of subroutine in my comment. It is: semi_lagrange_precip or perhaps better is semi_lagrange_sedim — You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub <#770 (comment)>, or unsubscribe https://github.com/notifications/unsubscribe-auth/AVR3U6WNWUTPADXJODCZRNLUKVLE5ANCNFSM5HN4HQXA . 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.

Adjustments were made. Thank you for your reviews and please let me know if you have any other suggestions.

@Songyou184
Copy link
Copy Markdown

Songyou184 commented Nov 6, 2021 via email

Copy link
Copy Markdown
Contributor

@gthompsnWRF gthompsnWRF left a comment

Choose a reason for hiding this comment

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

Thanks Ruiyu for making the small changes suggested!

@RuiyuSun
Copy link
Copy Markdown
Contributor Author

RuiyuSun commented Nov 8, 2021 via email

@RuiyuSun
Copy link
Copy Markdown
Contributor Author

Two additional suggestions were made by Songyou:
SEMI_SEDI_UPDATE and SEMI_SEDI_DECFL options are needed in early stages of developing the explicit semi-lagrangian, and no longer needed as namelist variables. SEDI_SEMI is to be the master namelist variable, I think.
Instead, DECFL = 8.0, is to be moved to a namelist variable, for future evaluation.

Copy link
Copy Markdown
Contributor

@mzhangw mzhangw left a comment

Choose a reason for hiding this comment

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

minor format changes need to be made.

Comment thread physics/module_mp_thompson.F90 Outdated
Comment thread physics/module_mp_thompson.F90 Outdated
Copy link
Copy Markdown
Collaborator

@climbfuji climbfuji left a comment

Choose a reason for hiding this comment

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

I suggest to revert the whitespace changes that I mentioned in a follow-up PR.

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.

5 participants