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

Calibrate global bucket LHF and SHF #835

Draft
wants to merge 3 commits into
base: main
Choose a base branch
from

Conversation

AlexisRenchon
Copy link
Member

@AlexisRenchon AlexisRenchon commented Oct 9, 2024

closes #871

This PR goal is to build a global calibration pipeline (minimal working example - MWE) to calibrate a land model globally.
In this commit, we calibrate a bucket model (the framework for full land model would be the same).
Two different approaches are added to this commit: the same calibration is done with just EnsembleKalmanProcesses.jl - EKP, or with ClimaCalibrate.jl (ClimaCalibrate.jl uses EKP, but has different backend, optimised for different HPCs).
In our MWE, we use ERA5 latent heat - LH and sensible heat - SH flux as target observations for the calibration. We stack monthly observation at n (e.g., 50) locations on land.
We calibrate the following parameters with following priors:

prior_κ_soil = EKP.constrained_gaussian("κ_soil", 2, 1, 0, Inf);
prior_ρc_soil = EKP.constrained_gaussian("ρc_soil", 4e6, 2e6, 0, Inf);
prior_f_bucket = EKP.constrained_gaussian("f_bucket", 0.5, 0.3, 0, 1);
prior_W_f = EKP.constrained_gaussian("W_f", 0.4, 0.4, 0, Inf);
prior_p = EKP.constrained_gaussian("p", 2, 1, 1, Inf);
prior_z_0m = EKP.constrained_gaussian("z_0m", 0.01, 0.1, 0, Inf);

and we use:

ensemble_size = 10
n_iterations = 5

Future work will use this pipeline to calibrate the land model globally.

TEMPORARY NOTE - delete before merge:

  • EnsembleKalmanProcesses.jl - done!
  • ClimaCalibrate.jl Script testing WIP... error: NetCDF error: HDF error. From Gabriele: "I see this error when NC files are already open somewhere, for example, if it is still being written. Make sure that writer are closed when the workers are done"

@AlexisRenchon AlexisRenchon self-assigned this Oct 22, 2024
@AlexisRenchon AlexisRenchon force-pushed the kd/ar/calibrate_global_bucket branch from f5876ad to 2079dff Compare October 29, 2024 17:07
@AlexisRenchon AlexisRenchon force-pushed the kd/ar/calibrate_global_bucket branch 16 times, most recently from ff3839c to 07e76c5 Compare November 11, 2024 23:05
@AlexisRenchon AlexisRenchon force-pushed the kd/ar/calibrate_global_bucket branch 3 times, most recently from 60612f3 to 175f5fe Compare November 12, 2024 21:28
src/Artifacts.jl Outdated Show resolved Hide resolved
@AlexisRenchon AlexisRenchon force-pushed the kd/ar/calibrate_global_bucket branch 3 times, most recently from 662deff to 1c233c5 Compare December 2, 2024 18:46
@AlexisRenchon AlexisRenchon force-pushed the kd/ar/calibrate_global_bucket branch 2 times, most recently from 74bb016 to 51f14b4 Compare December 5, 2024 17:57
@AlexisRenchon AlexisRenchon force-pushed the kd/ar/calibrate_global_bucket branch 8 times, most recently from dea2ee2 to bbfe99b Compare January 15, 2025 18:54
@AlexisRenchon AlexisRenchon force-pushed the kd/ar/calibrate_global_bucket branch from bbfe99b to 26f1d76 Compare January 15, 2025 21:57
@AlexisRenchon AlexisRenchon force-pushed the kd/ar/calibrate_global_bucket branch 11 times, most recently from cfae5f5 to be1f539 Compare January 31, 2025 17:10
@AlexisRenchon AlexisRenchon force-pushed the kd/ar/calibrate_global_bucket branch from be1f539 to e2d9569 Compare January 31, 2025 19:25
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

Global cabration MWE bucket model LHF
3 participants