Repository navigation
Casim hail diagnostic - #27
James Bruten (james-bruten-mo) merged 9 commits into
Conversation
…y and Mason melting
|
added trac.log now - all passed |
|
Hello James - thanks for picking this up. |
| da=(Kw0*t*dz/Lf/rhoi/v)/ & | ||
| ((a0-a)*a/a0+(Kw0/(C*(Ka+Lv*Dv*beta)))*a**2/a0) | ||
|
|
||
| IF ( Re .gt. 0.0 ) THEN |
There was a problem hiding this comment.
not sure why Re could be negative - is this just sign convention for particles terminal velocity vs updraughts? Every other parameter must be postive?
| da=(Kw0*t*dz/Lf/rhoi/v)/ & | ||
| ((a0-a)*a/a0+(Kw0/(C*(Ka+Lv*Dv*beta)))*a**2/a0) | ||
|
|
||
| IF ( Re .gt. 0.0 ) THEN |
There was a problem hiding this comment.
why would Reynolds number be negative - is this just fixing some numerical errors, or is it when "fall"speed is upwards?
There was a problem hiding this comment.
could just have an a .gt. 0. check before calling the mason_melt routine?
There was a problem hiding this comment.
yep - done that now in the elf4299 change. Moved the if out of mason_melting and into the loop that calls mason_melting.
The ensemble test still passes.
There was a problem hiding this comment.
updated the rose-stem trac.log with new branch. All passes still.
There was a problem hiding this comment.
great, looks good to me
Paul Barrett (paul-barrett)
left a comment
There was a problem hiding this comment.
looks good to me, ready for CR James Bruten (@james-bruten-mo)
| v=a_v*(2.0*a*100.0)**b_v * sqrt(rho0/rhoa(k)) !diam in cm | ||
| da=mason_melt(a,a0,dz_in(k),t(k)-T0,v,rhoi,rhoa(k)) | ||
| a=a-da | ||
| if ( a .gt. 0.0 ) then |
There was a problem hiding this comment.
maintains style of the hail check (a .gt. 0.0), happy with that
James Bruten (james-bruten-mo)
left a comment
There was a problem hiding this comment.
The changes look sensible to me. I'll rerun the testing suites then get it on
deab1a0
into
MetOffice:main
PR Summary
** Fix a bug in the hail diagnostic **
Sci/Tech Reviewer: Paul Barrett (@paul-barrett)
Code Reviewer: James Bruten (@james-bruten-mo)
Added 2 limiters to the hail diagnostic. One to set the change of particle size to 0 if the size goes negative. The second is in the melting equation which just limits the delta temperature to positive values in the profile. This fixed the failure from PS48 that i was given.
Code Quality Checklist
(Some checks are automatically carried out via the CI pipeline)
Testing
stem suites
acceptable (eg. kgo changes)
tests, unit tests, etc.)
trac.log
Test Suite Results - um - casim_hail_diagnostic/run4
Suite Information
Approvals
Code Owners
Config Owners
No UM Config Owners Required
Task Information
✅ succeeded tasks - 979
Security Considerations
Performance Impact
AI Assistance and Attribution
Documentation
Sci/Tech Review
Please alert the code reviewer via a tag when you have approved the SR
Code Review