-
Notifications
You must be signed in to change notification settings - Fork 19
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Remove 0.5*timestep logic from call to zm #109
Remove 0.5*timestep logic from call to zm #109
Conversation
@@ -68,9 +68,6 @@ subroutine zm_conv_convtran_run(ncol, pver, & | |||
|
|||
real(kind_phys), intent(in) :: dpdry(:,:) ! Delta pressure between interfaces (ncol,pver) | |||
|
|||
real(kind_phys), intent(in) :: dt ! 2 delta t (model time increment) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I presume dt is just not used in convtran?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
That is correct - when I went to see if the dt needed adjusting in this routine, I discovered it was not used, so I removed it.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good to me!
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
one tiny thing!
zhang_mcfarlane/zm_conv_momtran.F90
Outdated
@@ -59,7 +59,7 @@ subroutine zm_conv_momtran_run(ncol, pver, pverp, & | |||
real(kind_phys), intent(in) :: ed(:,:) ! Mass entraining from downdraft (ncol,pver) | |||
real(kind_phys), intent(in) :: dp(:,:) ! Delta pressure between interfaces (ncol,pver) | |||
real(kind_phys), intent(in) :: dsubcld(:) ! Delta pressure from cloud base to sfc (ncol) | |||
real(kind_phys), intent(in) :: dt ! time step in seconds : 2*delta_t | |||
real(kind_phys), intent(in) :: dt |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
can we put the ! time step in seconds
comment back?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
That's fine by me @cacraigucar
ZM has for a very long time passed in .5*timestep and then inside multiplied by 2 in several places (and also done 0.5/input_timestep in several other locations). It also hinders CCPP'ization as we do not have a half timestep standard name. This confusing, unnecessary logic will be removed.
To confirm the results, see this email from @adamrher
Your ZM branch has round-off level differences with cam6_4_018, so it's working as expected (plots attached).
Ignore the top right panel "cam6_4_015" -- that's my clubb+mf branch and something is very wrong with that run.