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

Bug fix for COARE3.5 (not used by default) and move some parameters t… #143

Merged
merged 2 commits into from
May 18, 2022

Conversation

joeolson42
Copy link
Collaborator

Two changes:

  1. Bug fix for COARE3.5 - this option is not used by default so it will not cause a change in results.
  2. Moved 4 internal parameters to namelist options:
  • isftcflx (default = 0)
  • iz0tlnd (default = 0)
  • sfclay_compute_flux (default = .false.)
  • sfclay_compute_diag (default = .false.)

@@ -2418,7 +2419,7 @@ SUBROUTINE edson_etal_2013(Z_0,ustar,wsp10,visc,zu)
REAL, INTENT(IN) :: ustar, visc, wsp10, zu
REAL, INTENT(OUT) :: Z_0
REAL, PARAMETER :: G=9.81
REAL, PARAMETER :: m=0.017, b=-0.005
REAL, PARAMETER :: m=0.0017, b=-0.005
Copy link
Collaborator

Choose a reason for hiding this comment

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

Is this change the reason the regression tests failed?

Copy link
Collaborator

Choose a reason for hiding this comment

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

I don't see how the code could reach this point since COARE_OPT is hard-coded to 3.0, so this subroutine is never called.

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

Yes, it's currently hard-coded to the COARE 3.0 algorithm, so this isn't the reason for the failure. This change does make COARE 3.5 competitive with COARE 3.0. Beforehand, COARE 3.5 has a large negative 10-m wind speed bias:

Screen Shot 2022-03-29 at 5 26 11 PM

@SamuelTrahanNOAA
Copy link
Collaborator

These changes are in NCAR#911 and will be merged back to gsl/develop after the community repositories are updated.

@SamuelTrahanNOAA
Copy link
Collaborator

The latest batch of changes (5c6243d) are also in NCAR#911

@SamuelTrahanNOAA
Copy link
Collaborator

#143 will merge these changes

@SamuelTrahanNOAA SamuelTrahanNOAA merged commit cf9e003 into NOAA-GSL:gsl/develop May 18, 2022
TH1D(I)=T1D(I)*THCON(I) !(Theta, K)
TC1D(I)=T1D(I)-273.15 !(T, Celsius)
TH1D(I)=T1D(I)**(100000./P1D(I))**ROVCP !(Theta, K)
TC1D(I)=T1D(I)-273.15 !(T, Celsius)
Copy link
Collaborator

Choose a reason for hiding this comment

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

@joeolson42 Joe, this PR got the "buggy" MYNN surface layer scheme in. It looks like it did not get into the RRFS_dev branch

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.

4 participants