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

Default parameters values include Infs + missing parameters/catke_basic_conv_adj_parameters.jld2 #2

Open
navidcy opened this issue Nov 6, 2022 · 3 comments

Comments

@navidcy
Copy link
Collaborator

navidcy commented Nov 6, 2022

I'm trying to run the global simulations but seems to be missing parameters/catke_basic_conv_adj_parameters.jld2.

Also, some of the parameters have default Inf values. As far as I remember this creates issues? (CliMA/ParameterEstimocean.jl#295 (comment))

neutral_default_mixing_length_parameters = (Cᵇu = Inf, Cᵇc = Inf, Cᵇe = Inf,
Cˢu = Inf, Cˢc = Inf, Cˢe = Inf,
CᴷRiᶜ = Inf, CᴷRiʷ = 0.0)
neutral_default_tke_parameters = (CᵂwΔ = 0.0, Cᵂu★ = 0.0,
Cᴰ⁻ = 0.0, Cᴰ⁺ = 0.0, CᴰRiᶜ = Inf, CᴰRiʷ = 0.0)
mixing_length = MixingLength(; neutral_default_mixing_length_parameters...)
turbulent_kinetic_energy_equation = TurbulentKineticEnergyEquation(; neutral_default_tke_parameters...)
neutral_catke = CATKEVerticalDiffusivity(; mixing_length, turbulent_kinetic_energy_equation)
catke_parameters_name = "catke_basic_conv_adj_parameters.jld2"
catke_parameters_filename = joinpath("..", "parameters", catke_parameters_name)
catke_parameters = load(catke_parameters_filename, "mean") # load ensemble mean parameters from EKI calibration

cc @glwagner, @simone-silvestri

@navidcy navidcy changed the title Int in default parameters + missing parameters/catke_basic_conv_adj_parameters.jld2 Default parameters values include Infs + missing parameters/catke_basic_conv_adj_parameters.jld2 Nov 6, 2022
@glwagner
Copy link
Member

glwagner commented Nov 7, 2022

CATKE doesn't work yet, so try running without it.

@navidcy
Copy link
Collaborator Author

navidcy commented Nov 7, 2022

I did, and that worked fine!

@glwagner
Copy link
Member

glwagner commented Nov 8, 2022

Also should mention that Inf is the "neutral" value for some parameters (ie if C=Inf then it does not affect the outcome). That's because they are combined with other parameters in a construction that's something like a = min(C1, C2, C3); therefore if C2=Inf it does not affect model outcome.

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

No branches or pull requests

2 participants