Procedures

set the following two backflow to set (see Bolker 199x)

vh_js !!

Model initialisation routine for the CABLE offline driver.

Computes the fraction of net photosynthate allocated to leaf, wood and froot using an allocation scheme modified from Pierre Friedlingstein.

revise turnover and NPP and dcplantdt to reflect above adjustments End of adjustments to avoid negative stores Ticket#108

vh_js !! vh_js !! vh_js

Carbon allocation algorithm implemented by Dan Qiu and YP Wang in 2011 based on Adam Wolf 1 and 2

CLN CHARACTER(LEN=99), INTENT(IN) :: fcnpspin vh_js !! CLN901 format(A99) $ WHERE (pop%pop_grid(:)%cmass_sum_old.gt.0.1 .and. pop%pop_grid(:)%cmass_sum.gt.0.1 ) $ casapool%Cplant(Iw,2) = casapool%Cplant(Iw,2)(1.0- min( POP%pop_grid(:)%cat_mortality/(POP%pop_grid(:)%cmass_sum_old),0.99)) $ casapool%Nplant(Iw,2) = casapool%Nplant(Iw,2)(1.0- min( POP%pop_grid(:)%cat_mortality/(POP%pop_grid(:)%cmass_sum_old),0.99)) $ ENDWHERE

Calculates the extinction coefficients for black leaves. It returns:

This SUBROUTINE evaluates the value of the friction velocity $u_*$ as used during the iteration loop of the Monin-Obukhov (MO) similarity theory in define_canopy. $u_*$ is used when quantifying the stability parameter canopy%zetar in update_zetar, and when evaluating the resistance network components rt0 and rt1 in define_canopy. The friction velocity quantifies the magnitude of turbulent mixing that is occurring as well as the flux of momentum from the atmosphere to the land.

Searches for CABLE parameters in the met forcing file, and if it finds any, uses these values to overwrite the values that have already been loaded from the default parameter loading and/or the restart file.

MMY : move from subroutine derived_parameters to here since derived_parameters is now an option but compulsory

This subroutine computes the height of the canopy above ground/snow surface when there is snow present.

Initialize the surface wetness factor of the soil/snow (ssnow_wetfac) array

This subroutine computes the leaf-area index of a canopy when there is snow present. The formulae assume that the leaf area is distributed uniformly in the vertical within the canopy.

Allocates the luc%%var(imland,mvmax) variables

Allocates the luc%var(mp) variables.

Checks the mass balance and writes the output CNP pool sizes for each plant functional type.

Deallocates the luc%var(mland,mvmax) variables

Deallocates the luc%var(mp) variables

Main driver for the land-use change

Maps luc%var_y(mland,mvmax) to lucmp%var(mp)

Assigns lucmp%var(mp) to luc%var_x(mland,mvmax)

Redistribution of the plant functional type transition to ensure that the change in plant functional type fractions from the state data is consistent with the estimates from previous states and transition.

This subroutine applies the land-use changes to the transfer of the C, N and P pools, the biophysical states, the soil texture and the soil order for each patch. Then it will seed the deforested land.

Assigns luc%var_y to luc%var_x

This SUBROUTINE converts the previously evaluated rate of potential evaporation ssnow_potev from the soil into the associated latent heat flux canopy_fes. The evaluation accounts for whether the water flux is in liquid or solid phase (evaporation or sublimation, dew or frost), the direction of the flux (to or from the soil), and any limits due to a lack of water in the surface layer of the soil column. The SUBROUTINE also evaluates the partition of the latent heat flux between soil canopy_fess and puddle canopy_fesp contributions.

This SUBROUTINE checks that input values of leaf area and canopy height lie between set values (defined by PFT) and overwrites if necessary. These limits are needed to ensure that the canopy roughness properties of the land point are well behaved (see ruff_resist).

Defines the priority order of sources of parameter values for CABLE, determines the total number of patches over all grid cells, and writes parameter values to CABLE's parameter arrays.

Evaluates the integrated similarity function for momentum transfer, $\psi_m(\xi)$.

Evaluates the integrated similarity function for turbulent transfer of scalars, $\psi_s(\xi)$.

This SUBROUTINE evaluates the specific humidity (water vapour mixing ratio) at saturation (in kgkg$^{-1}$), at a given air temperature (in $^{\circ}$C) and air pressure (in hPa) for an array of size mp.

This SUBROUTINE evaluates the specific humidity (water vapour mixing ratio) at saturation (in kgkg$^{-1}$), at a given air temperature (in $^{\circ}$C) and air pressure (in hPa).

$

Calculates the roughness parameters and the aerodynamic contribution to the resistances controlling the fluxes between the land and atmosphere for each land point.

Set individual variables to output according to the values of the group options from the namelist entries in output%.

calculate snow processes and thermal soil

CLN CHARACTER(LEN=99), INTENT(IN) :: fcnpspin vh_js !! CLN OPEN(91, file=fcnpspin) CLN read(91,*) myearspin CLN901 format(A99) $ WHERE(xkNlimiting .eq. 0) !Chris Lu 4/June/2012 $ xkNlimiting = 0.001 $ END WHERE $ call totcnppools(1,veg,casamet,casapool,bmcplant,bmnplant,bmpplant,bmclitter,bmnlitter,bmplitter, & $ bmcsoil,bmnsoil,bmpsoil,bmnsoilmin,bmpsoillab,bmpsoilsorb,bmpsoilocc,bmarea)

This SUBROUTINE updates the value of the stability parameter $\xi$ during the iteration loop of the Monin-Obukhov (MO) similarity theory in define_canopy.
$\xi$ quantifies the role that the surface fluxes play in setting the efficiency of turbulent transfer from the land to the atmosphere, and hence the aerodynamic component of the resistance network for those same surface fluxes (an implicit problem which requires iteration to solve).