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

make snow tests more comprehensive, fix bug #330

Merged
merged 1 commit into from
Sep 19, 2023
Merged

Conversation

kmdeck
Copy link
Member

@kmdeck kmdeck commented Sep 15, 2023

Purpose

Investigate conservation of energy and water for the bucket model due to issues seen in the Coupler
CliMA/ClimaCoupler.jl#424
CliMA/ClimaCoupler.jl#421

To-do

Content

  • Add to test script for snow bucket conservation: run with nonzero precip, and with nonuniform evaporation field.
  • Fix bug in heaviside function
  • More care needed when tracking energy conservation. The energy in land is \rhoc_soil * \sum T_soil dz - \rho_l LH_f0*SWE. When snow falls on land, the energy changes because SWE changes, and we need to account for that in the test. This is not automatically accounted for because we dont have a snow energy variable (because we track SWE as the prognostic variable and the snow temperature is assumed to be the same as the freezing point). We also do not track the energy in the water on land, so there is no contribution for liquid precip.

Review checklist

I have:

In the Content, I have included

  • relevant unit tests, and integration tests,
  • appropriate docstrings on all functions, structs, and modules, and included relevant documentation.

  • I have read and checked the items on the review checklist.

@kmdeck kmdeck requested a review from LenkaNovak September 15, 2023 18:08
Copy link
Contributor

@LenkaNovak LenkaNovak left a comment

Choose a reason for hiding this comment

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

LGTM, thank you, @kmdeck !

@kmdeck kmdeck force-pushed the investigate_snow_behavior branch from cc2d36f to f6586bd Compare September 19, 2023 00:16
@kmdeck
Copy link
Member Author

kmdeck commented Sep 19, 2023

bors r+

@bors
Copy link
Contributor

bors bot commented Sep 19, 2023

Build succeeded!

The publicly hosted instance of bors-ng is deprecated and will go away soon.

If you want to self-host your own instance, instructions are here.
For more help, visit the forum.

If you want to switch to GitHub's built-in merge queue, visit their help page.

@bors bors bot merged commit 4143c8a into main Sep 19, 2023
@bors bors bot deleted the investigate_snow_behavior branch September 19, 2023 01:31
@juliasloan25 juliasloan25 mentioned this pull request Sep 22, 2023
1 task
mitraA90 pushed a commit that referenced this pull request Dec 22, 2023
330: make snow tests more comprehensive, fix bug r=kmdeck a=kmdeck

## Purpose 
Investigate conservation of energy and water for the bucket model due to issues seen in the Coupler
CliMA/ClimaCoupler.jl#424
CliMA/ClimaCoupler.jl#421



## To-do

## Content
- Add to test script for snow bucket conservation: run with nonzero precip, and with nonuniform evaporation field.
- Fix bug in heaviside function
- More care needed when tracking energy conservation. The energy in land is \rhoc_soil * \sum T_soil dz - \rho_l LH_f0*SWE. When snow falls on land, the energy changes because SWE changes, and we need to account for that in the test. This is not automatically accounted for because we dont have a snow energy variable (because we track SWE as the prognostic variable and the snow temperature is assumed to be the same as the freezing point).  We also do not track the energy in the water on land, so there is no contribution for liquid precip.

Review checklist

I have:
- followed the codebase contribution guide: https://clima.github.io/ClimateMachine.jl/latest/Contributing/
- followed the style guide: https://clima.github.io/ClimateMachine.jl/latest/DevDocs/CodeStyle/
- followed the documentation policy: https://github.com/CliMA/policies/wiki/Documentation-Policy
- checked that this PR does not duplicate an open PR.

In the Content, I have included 
- relevant unit tests, and integration tests, 
- appropriate docstrings on all functions, structs, and modules, and included relevant documentation.


----
- [X] I have read and checked the items on the review checklist.


Co-authored-by: kmdeck <[email protected]>
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.

2 participants