Skip to content
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

[production/RRFS.v1] physics updates for RRFS.v1 code freeze #176

Conversation

haiqinli
Copy link
Collaborator

@haiqinli haiqinli commented Feb 23, 2024

This PR includes the following updates:
1). MYNN updates;
2). RUC LSM updates (#163);
3). GF updates for humidity DA; restore scale-awareness in the 1st hour; suppress weak radar reflectivity over water;
4). smoke plume rise updates for stability on WCOSS2;
5). GWD updates to reduce the high 10m wind speed bias.

…nd it will be consistent with convective precipitation"
@MatthewPyle-NOAA
Copy link
Collaborator

@grantfirl Does more need to be done before this PR will be reviewed/merged? Still trying to understand the process with the production release branch. Thanks!

@grantfirl
Copy link
Collaborator

@MatthewPyle-NOAA We just need to review it and @haiqinli needs to upload the test_changes.list file in ufs-community/ufs-weather-model#2147, assuming that we don't request changes. I'm also working on getting #147 and #159 in

@MatthewPyle-NOAA
Copy link
Collaborator

Thanks for the explanation, @grantfirl

@@ -21,11 +21,11 @@ module GFS_MP_generic_post
subroutine GFS_MP_generic_post_run( &
im, levs, kdt, nrcm, nncl, ntcw, ntrac, imp_physics, imp_physics_gfdl, imp_physics_thompson, imp_physics_nssl, &
imp_physics_mg, imp_physics_fer_hires, cal_pre, cplflx, cplchm, cpllnd, progsigma, con_g, rhowater, rainmin, dtf, &
frain, rainc, rain1, rann, xlat, xlon, gt0, gq0, prsl, prsi, phii, tsfc, ice, phil, htop, refl_10cm, &
frain, rainc, rain1, rann, xlat, xlon, gt0, gq0, prsl, prsi, phii, tsfc, ice, phil, htop, refl_10cm, &
Copy link
Collaborator

Choose a reason for hiding this comment

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

These are only whitespace changes. Is that intended? If not, please revert this change.

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

@grantfirl I will revert this. Thanks.

@@ -123,7 +123,7 @@ subroutine rrfs_smoke_wrapper_run(im, kte, kme, ktau, dt, garea, land, jdate,
ebu_smoke,fhist,min_fplume, &
max_fplume, hwp, hwp_ave, wetness, ndvel, ddvel_inout, &
peak_hr_out,lu_nofire_out,lu_qfire_out, &
fire_heat_flux_out, frac_grid_burned_out, kpbl,oro, &
fire_heat_flux_out, frac_grid_burned_out, kpbl,oro, totprcp, &
Copy link
Collaborator

Choose a reason for hiding this comment

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

Why are you adding totprcp when you're not using it?

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

@grantfirl I will remove this unused variable. Thanks.

@@ -916,6 +916,14 @@
type = real
kind = kind_phys
intent = in
[totprcp]
Copy link
Collaborator

Choose a reason for hiding this comment

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

Same question. Please remove this change if you're not using totprcp.

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

@grantfirl I will remove this unused variable. Thanks.

@grantfirl
Copy link
Collaborator

@haiqinli Looks fine with the exception of some minor typo-like things. Once you address those, I'll approve. Please remember to commit the test_changes.list file that is generated now in the ufs-weather-model/tests directory for that PR.

@grantfirl
Copy link
Collaborator

This contains the changes from #163.

@@ -1363,7 +1363,8 @@ subroutine drag_suite_run( &
DO k=kts,km
wsp=SQRT(uwnd1(i,k)**2 + vwnd1(i,k)**2)
! alpha*beta*Cmd*Ccorr*2.109 = 12.*1.*0.005*0.6*2.109 = 0.0759
var_temp = 0.0759*EXP(-(zl(i,k)/H_efold)**1.5)*a2* &
! Change alpha to 35 -- 0.0759 becomes 0.2214
var_temp = 0.2214*EXP(-(zl(i,k)/H_efold)**1.5)*a2* &
Copy link
Collaborator

Choose a reason for hiding this comment

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

@haiqinli @mdtoyNOAA Will this change affect global and HAFS applications as well ? Have any tests been done to check the impact on applicaitons other than RRFS ?

Copy link
Collaborator

Choose a reason for hiding this comment

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

We tested this change in the HR3 C768 FV3GFS and detected improvements to the 500hPa ACC scores.

Copy link
Collaborator

Choose a reason for hiding this comment

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

Mike, that is good know. I'd assume this has major impact on surface winds and turbulent exchange coefficients as well. Did you also see positive impact on grid-to-obs stats ?

Copy link
Collaborator

Choose a reason for hiding this comment

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

Mike, since this change is hardwired and will be include in HR4. Please record this change in our tracking spreadsheet.

Copy link
Collaborator

Choose a reason for hiding this comment

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

@yangfanglin Keep in mind that this PR is just targeting the RRFS.v1 branch.

Copy link
Collaborator

Choose a reason for hiding this comment

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

@MatthewPyle-NOAA I see. So the change is made locally and will not be committed back to the main repo.

Copy link
Collaborator

Choose a reason for hiding this comment

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

@yangfanglin It may be included in a future PR targeting the main repo, but that will be a separate, future step.

Copy link
Collaborator

Choose a reason for hiding this comment

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

Mike, that is good know. I'd assume this has major impact on surface winds and turbulent exchange coefficients as well. Did you also see positive impact on grid-to-obs stats ?

In the 10m winds we see significant improvements in RMSE in both west and east CONUS. However, a slight degradation of windspeed bias in west (little change in bias in east).

Copy link
Collaborator

Choose a reason for hiding this comment

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

For VSDB verification, go to:
https://www.emc.ncep.noaa.gov/gmb/jhan/vsdbw/hr3gwd01/
The modifications to TOFD are the experiment "HR3GWD01". "HR2D13" is HR2.

@grantfirl grantfirl changed the title physics updates for RRFS.v1 code freeze [production/RRFS.v1] physics updates for RRFS.v1 code freeze Feb 28, 2024
@jkbk2004
Copy link

jkbk2004 commented Mar 1, 2024

tests are done at ufs-community/ufs-weather-model#2147. @grantfirl can you merge this pr?

@jkbk2004
Copy link

jkbk2004 commented Mar 1, 2024

@MatthewPyle-NOAA You may go ahead to merge this pr if you have permission.

@MatthewPyle-NOAA
Copy link
Collaborator

@jkbk2004 I don't have that permission, at least not yet.

@climbfuji
Copy link

I'll merge it for you

@climbfuji climbfuji merged commit 8a4b102 into ufs-community:production/RRFS.v1 Mar 1, 2024
3 checks passed
@jkbk2004
Copy link

jkbk2004 commented Mar 1, 2024

I'll merge it for you

Thanks! @climbfuji

@grantfirl
Copy link
Collaborator

@climbfuji Thanks. I got a late start this morning.

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.

7 participants