Procedures

Allocate and initialise arrays in the climate_data type

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

Associate the climate type pointers to the climate_data arrays.

vh_js !!

Model initialisation routine for the CABLE offline driver.

Model initialisation routine (CRU specific). TRENDY experiment using CRU-NCEP.

Model initialisation routine (default met specific).

Model initialisation routine (GSWP specific).

Model initialisation routine (PLUME specific). PLUME experiment setup using WATCH.

Model initialisation routine (site met specific). Site experiment, e.g. AmazonFace (spinup or transient run type).

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.

Compare reference values for quasi-bitwise reproducibility checks and write reference value on failed reproducibility.

Deallocate arrays in the climate_data type

canopy%fes_cor = 0.0

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.

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.

Reading LU input data, zeroing biomass in empty secondary forest tiles and tranferring LUC-based age weights for secondary forest to POP structure

Calculates the energy and water balances at each time step as well as the cumulative balance over the simulation.

Check if an MPI return code signaled an error. If so, print the corresponding message and abort the execution.

Class method to abort execution of an MPI group.

Contructor for mpi_grp_t class.

Overload the default construct for mpi_grp_t

Finalise MPI.

Initialise MPI and set default communicator.

Stub for when MPI is not available

Stub for when MPI is not available

Nullify the climate type pointers.

Select the correct files given the year for filenames following the gswp format

Select the correct files given the year for filenames following the princeton format

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).

$

Replace the year in the filename with the value of ncciy for the gswp file format.

Replace the year in the filename with the value of ncciy for the princeton format

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

Offline serial driver.

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).