supporting¶

Test

aggreg(atmfile, SCOPEspec)¶: Aggregate MODTRAN data over SCOPE bands by averaging (over rectangular band passes)

calczenithangle(Doy, t, Omega_g, Fi_gm, Long, Lat)¶

author: Christiaan van der Tol (c.vandertol@utwente.nl) date: Jan 2003 update: Oct 2008 by Joris Timmermans (j_timmermans@itc.nl):

corrected equation of time

Oct 2012 (CvdT) comment: input time is GMT, not local time!

function [Fi_s,Fi_gs,Fi_g]= calczenithangle(Doy,t,Omega_g,Fi_gm,Long,Lat)

calculates pi/2-the angle of the sun with the slope of the surface.

input: Doy day of the year t time of the day (hours, GMT) Omega_g slope azimuth angle (deg) Fi_gm slope of the surface (deg) Long Longitude (decimal) Lat Latitude (decimal)

output: Fi_s ‘classic’ zenith angle: perpendicular to horizontal plane Fi_gs solar angle perpendicular to surface slope Fi_g projected slope of the surface in the plane through the solar beam and the vertical

calc_rssrbs(SMC, LAI, rbs)¶

count_k(nvars, v, vmax, id)¶: nvars = number of digits v = current vector of digits vmax = maximum values of digits id = starting digit vnew = new vector of digits

fixedp_brent_ari(func, x0, corner, tolFn, verbose)¶

Find a fixed point of func(x) using Brent’s method, as described by Brent 1971

func is a single-argument function, f(x) that returns a value the same size as x:: The goal is to find f(x) = x (or for Brent, f(x) - x = 0).

x0 is the initial guess (or 2 x n matrix if we want to generalize) tol is the tolerance in x (or if two-valued, x, f(x)? ) corner (optional) is a known “edge” in the function that could slow down the algorithm

if specified and the first two points include the corner, the corner will be substituted as a starting point.

Written by: Ari Kornfeld, 2016-10

latin_hypercube_input(tab, n_spectra, outdir)¶

leafangles(a, b)¶

Subroutine FluorSail_dladgen Version 2.3 For more information look to page 128 of “theory of radiative transfer models applied in optical remote sensing of vegetation canopies”

FluorSail for Matlab FluorSail is created by Wout Verhoef, National Aerospace Laboratory (NLR) Present e-mail: w.verhoef@utwente.nl

This code was created by Joris Timmermans, International institute for Geo-Information Science and Earth Observation. (ITC) Email: j.timmermans@utwente.nl

main function

meanleaf(canopy, F, choice, Ps)¶

Monin_Obukhov(constants, meteo, H)¶

Planck(wl, Tb, em)¶

plot_output_one2one()¶: …

satvap(T)¶

function [es,s]= satvap(T) Author: Dr. ir. Christiaan van der Tol Date: 2003

calculates the saturated vapour pressure at temperature T (degrees C) and the derivative of es to temperature s (kPa/C) the output is in mbar or hPa. The approximation formula that is used is: es(T) = es(0)*10^(aT/(b+T)); where es(0) = 6.107 mb, a = 7.5 and b = 237.3 degrees C and s(T) = es(T)*ln(10)*a*b/(b+T)^2

Sint(y, x)¶: Simpson integration x and y must be any vectors (rows, columns), but of the same length x must be a monotonically increasing series

slope_satvap(T)¶

function [es,s]= satvap(T) Author: Dr. ir. Christiaan van der Tol Date: 2003

calculates the saturated vapour pressure at temperature T (degrees C) and the derivative of es to temperature s (kPa/C) the output is in mbar or hPa. The approximation formula that is used is: es(T) = es(0)*10^(aT/(b+T)); where es(0) = 6.107 mb, a = 7.5 and b = 237.3 degrees C and s(T) = es(T)*ln(10)*a*b/(b+T)^2

Soil_Inertia0(cs, rhos, lambdas)¶: soil thermal inertia

Soil_Inertia1(SMC)¶: soil inertia method by Murray and Verhoef (

soil_respiration(Ts)¶

soltir_tp7(filename)¶: soltir_tp7 Reads MODTRAN tp7 file and applies a new MIT algorithm to derive 18 spectral functions for atmospheric correction and simulations at BOA and TOA

tav(alfa, nr)¶

zo_and_d(soil, canopy, constants)¶

function zom_and_d calculates roughness length for momentum and zero plane displacement from vegetation height and LAI

Date: 17 November 2008: 17 April 2013 (structures)
Author: A. Verhoef: implemented into Matlab by C. van der Tol (c.vandertol@utwente.nl)

Source: Verhoef, McNaughton & Jacobs (1997), HESS 1, 81-91

usage:: zo_and_d (soil,canopy)
canopy fields used as inpuyt:: LAI one sided leaf area index hc vegetation height (m)
soil fields used:: Cd Averaged drag coefficient for the vegetation CR Drag coefficient for isolated tree CSSOIL Drag coefficient for soil CD1 Fitting parameter Psicor Roughness layer correction
constants used (as global): kappa Von Karman’s constant
output:: zom roughness lenght for momentum (m) d zero plane displacement (m)