-
Notifications
You must be signed in to change notification settings - Fork 34
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
fix crash when using clm lake and GFS PBL/sfclay #148
fix crash when using clm lake and GFS PBL/sfclay #148
Conversation
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.
The changes look fine to me, just a small comment on the formatting
@@ -754,6 +756,11 @@ SUBROUTINE clm_lake_run( & | |||
weasd(i) = weasdi(i) | |||
snowd(i) = snodi(c) ! surface_snow_thickness_water_equivalent_over_ice | |||
|
|||
|
|||
if (.not. icy(i)) then | |||
flag_lakefreeze(i)=.true. |
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.
@JiliDong-NOAA Jili, should it be if(icy(i)) then when the lake grid point is frozen?
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.
@JiliDong-NOAA I see now, that this will be true only for new ice points. I agree with this change. Thank you.
This PR is to fix the RRFS ensemble member crash when combining clm lake and GFS PBL/sfclay:
#141
A new flag is created to track new freezing lake ice and let sfc_diff to calculate stability variables over the new icy grids for GFS PBL in the second surface loop, which otherwise will remain as missing values.