Jarvis Stewart#

class abcmodel.land.jarvis_stewart.JarvisStewartState(alpha, wg, temp_soil, temp2, surf_temp, wl, wq, wtheta, rs=<factory>, rssoil=<factory>, esat=<factory>, qsat=<factory>, dqsatdT=<factory>, e=<factory>, qsatsurf=<factory>, wCO2=<factory>, cliq=<factory>, temp_soil_tend=<factory>, wgtend=<factory>, wltend=<factory>, le_veg=<factory>, le_liq=<factory>, le_soil=<factory>, le=<factory>, hf=<factory>, gf=<factory>, le_pot=<factory>, le_ref=<factory>, vpd=<factory>)[source]#

Bases: StandardLandState

Jarvis-Stewart model state.

class abcmodel.land.jarvis_stewart.JarvisStewartModel(*args: Any, **kwargs: Any)[source]#

Bases: AbstractStandardLandModel

Jarvis-Stewart land surface model with empirical surface resistance.

Parameters:

**kwargs – additional keyword arguments to pass to the base class.

init_state(alpha=0.25, wg=0.21, temp_soil=285.0, temp2=286.0, surf_temp=290.0, wl=0.0, wq=0.0001, wtheta=0.1, rs=1000000.0, rssoil=1000000.0)[source]#

Initialize the model state.

Parameters:
  • alpha (float) – albedo [-]. Default is 0.25.

  • wg (float) – Volumetric soil moisture [m3 m-3]. Default is 0.21.

  • temp_soil (float) – Soil temperature [K]. Default is 285.0.

  • temp2 (float) – Deep soil temperature [K]. Default is 286.0.

  • surf_temp (float) – Surface temperature [K]. Default is 290.0.

  • wl (float) – Canopy water content [m]. Default is 0.0000.

  • wq (float) – Kinematic moisture flux [kg/kg m/s]. Default is 1e-4.

  • wtheta (float) – Kinematic heat flux [K m/s]. Default is 0.1.

  • rs (float) – Surface resistance [s m-1]. Default is 1.0e6.

  • rssoil (float) – Soil resistance [s m-1]. Default is 1.0e6.

Returns:

The initial land state.

update_surface_resistance(state)[source]#

Update the surface resistance rs in the state using the Jarvis-Stewart model.

Parameters:

state (AbstractCoupledState) – CoupledState.

Returns:

CoupledState (with updated land component).

compute_f1(in_srad)[source]#

Compute rad factor f1.

compute_f2(wg)[source]#

Compute soil moisture factor f2.

compute_f3(esat, e)[source]#

Compute VPD factor f3.

compute_f4(theta)[source]#

Compute temperature factor f4.

update_co2_flux(state)[source]#

No CO2 flux is computed using this model. See AgsModel.