diff --git a/src/bmi_lgar.cxx b/src/bmi_lgar.cxx index 9b70990..0708286 100644 --- a/src/bmi_lgar.cxx +++ b/src/bmi_lgar.cxx @@ -808,7 +808,7 @@ GetValuePtr (std::string name) else if (name.compare("smcmin") == 0) return (void*)this->state->lgar_calib_params.theta_r; else if (name.compare("van_genuchten_m") == 0) - return (void*)this->state->lgar_calib_params.vg_m; + return (void*)this->state->lgar_calib_params.vg_n; else if (name.compare("van_genuchten_alpha") == 0) return (void*)this->state->lgar_calib_params.vg_alpha; else if (name.compare("hydraulic_conductivity") == 0) diff --git a/src/lgar.cxx b/src/lgar.cxx index b296134..a877c1d 100755 --- a/src/lgar.cxx +++ b/src/lgar.cxx @@ -86,7 +86,7 @@ extern void lgar_initialize(string config_file, struct model_state *state) // initialize array for holding calibratable parameters state->lgar_calib_params.theta_e = new double[state->lgar_bmi_params.num_layers]; state->lgar_calib_params.theta_r = new double[state->lgar_bmi_params.num_layers]; - state->lgar_calib_params.vg_m = new double[state->lgar_bmi_params.num_layers]; + state->lgar_calib_params.vg_n = new double[state->lgar_bmi_params.num_layers]; state->lgar_calib_params.vg_alpha = new double[state->lgar_bmi_params.num_layers]; state->lgar_calib_params.Ksat = new double[state->lgar_bmi_params.num_layers]; @@ -103,7 +103,7 @@ extern void lgar_initialize(string config_file, struct model_state *state) state->lgar_calib_params.theta_e[i] = state->soil_properties[soil].theta_e; state->lgar_calib_params.theta_r[i] = state->soil_properties[soil].theta_r; - state->lgar_calib_params.vg_m[i] = state->soil_properties[soil].vg_m; + state->lgar_calib_params.vg_n[i] = state->soil_properties[soil].vg_n; state->lgar_calib_params.vg_alpha[i] = state->soil_properties[soil].vg_alpha_per_cm; state->lgar_calib_params.Ksat[i] = state->soil_properties[soil].Ksat_cm_per_h;