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

Land surface upgrades for HR4 #219

Merged
merged 10 commits into from
Jul 23, 2024

Conversation

HelinWei-NOAA
Copy link
Collaborator

To address some concerns from HR3 especially the surface cold biases we have finalized some land surface upgrades after carrying out extensive tests and evaluation.

  1. remove vegetation type-dependent for the diagnostic option (opt_diag) used in HR3
  2. use radiative temperature as skin temperature
  3. use single loop alternative approach to address some flux inconsistence below and above the canopy
  4. fix the thermal roughness bug (the one for bared portion in the composition is not the same as the one used in the bared title)
  5. adjust some momentum roughness length
  6. use the interpolated temp between the skin temp and the first soil temp to get the initial snow temp based on snow depth
  7. use constant snow thermal conductivity to reduce the snow insulation effects
  8. fix radiation parameter bugs for grassland/cropland
  9. adjust albedo parameter for each soil color category

@yangfanglin
Copy link
Collaborator

@barlage @cenlinhe @dustinswales @grantfirl We'd appreciate it if you can review this PR. We'd like to get this PR committed to start GFS HR4 experiments at the earliest time possible.

do is = isnow,0
tsnoxy(ix,is) = tgxy(ix)
! tsnoxy(ix,is) = tgxy(ix)
! tsnoxy(ix,is) = tgxy(ix) + (( sum(dzsno(isnow:is)) -0.5*dzsno(is) )/snd)*(tgxy(ix)-stc(ix,1))
Copy link
Collaborator

Choose a reason for hiding this comment

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

probably best to clean this up by removing all the comments

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

Done

@@ -1989,7 +1989,7 @@ subroutine energy (parameters,ice ,vegtyp ,ist ,nsnow ,nsoil , & !in

real (kind=kind_phys), parameter :: mpe = 1.e-6
real (kind=kind_phys), parameter :: psiwlt = -150. !metric potential for wilting point (m)
real (kind=kind_phys), parameter :: z0 = 0.002 ! bare-soil roughness length (m) (i.e., under the canopy)
real (kind=kind_phys), parameter :: z0 = 0.015 ! bare-soil roughness length (m) (i.e., under the canopy)
Copy link
Collaborator

Choose a reason for hiding this comment

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

where does "0.015" come from? Is it tuned? I did not find this value in the Noah-MP.

Copy link
Collaborator

Choose a reason for hiding this comment

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

it is tuned based on wind speed bias over bare soil and vegetated areas with low vegetation fraction

Copy link
Collaborator

Choose a reason for hiding this comment

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

got it. thanks!

Copy link
Collaborator

@cenlinhe cenlinhe left a comment

Choose a reason for hiding this comment

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

The changes look good to me.

@jkbk2004
Copy link

@grantfirl @dustinswales any chance for final reviews?

@FernandoAndrade-NOAA
Copy link

FernandoAndrade-NOAA commented Jul 23, 2024

Testing for #2354 has completed successfully, please continue with the merge process. @grantfirl

@grantfirl grantfirl merged commit 2a50ccc into ufs-community:ufs/dev Jul 23, 2024
3 checks passed
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