function [porosity, b, CGsat, p, a, C2ref,C1sat] = soilprop(stype)
%
% This function returns tabulated soil parameters from Noilhan and Planton
% (1989) and porosity from
% https://hess.copernicus.org/articles/22/4935/2018/hess-22-4935-2018-t04.pdf 
% November, 2022
% 
if strcmpi(stype,'Sand')
    porosity=0.437;
    b=4.05;
    CGsat=3.222;
    p=4;
    a=0.387;
    C2ref=3.9;
    C1sat=0.082;
elseif strcmpi(stype,'Loamy sand')
    porosity=0.437;
    b=4.38;
    CGsat=3.057;
    p=4;
    a=0.404;
    C2ref=3.7;
    C1sat=0.098;
elseif strcmpi(stype,'Sandy loam')
    porosity=0.453;
    b=4.9;
    CGsat=3.56;
    p=4;
    a=0.219;
    C2ref=1.8;
    C1sat=0.132;
elseif strcmpi(stype,'Silt loam')
    porosity=0.501;
    b=5.3;
    CGsat=4.418;
    p=6;
    a=0.102;
    C2ref=0.8;
    C1sat=0.153;    
elseif strcmpi(stype,'Loam')
    porosity=0.463;
    b=5.39;
    CGsat=4.111;
    p=6;
    a=0.148;
    C2ref=0.8;
    C1sat=0.191;
elseif strcmpi(stype,'Sandy clay loam')
    porosity=0.398;
    b=7.12;
    CGsat=3.67;
    p=6;
    a=0.135;
    C2ref=0.8;
    C1sat=0.213;
elseif strcmpi(stype,'Silty clay loam')
    porosity=0.471;
    b=7.75;
    CGsat=3.593;
    p=8;
    a=0.127;
    C2ref=0.4;
    C1sat=0.385;    
elseif strcmpi(stype,'Clay loam')
    porosity=0.464;
    b=8.52;
    CGsat=3.995;
    p=10;
    a=0.084;
    C2ref=0.6;
    C1sat=0.227;
elseif strcmpi(stype,'Sandy clay')
    porosity=0.43;
    b=10.4;
    CGsat=3.058;
    p=8;
    a=0.139;
    C2ref=0.3;
    C1sat=0.421;
elseif strcmpi(stype,'Silty clay')
    porosity=0.479;
    b=10.4;
    CGsat=3.729;
    p=10;
    a=0.075;
    C2ref=0.3;
    C1sat=0.375;    
elseif strcmpi(stype,'Clay')
    porosity=0.475;
    b=11.4;
    CGsat=3.6;
    p=12;
    a=0.083;
    C2ref=0.3;
    C1sat=0.342;    
end
end