ALLOCATE(soil_param%VolHeatCapacity(domain%nland))
    ALLOCATE(soil_param%ThermalDiffusivity(domain%nland))
    WHERE(FLOOR(fao) == 0)
       soil_param%VolHeatCapacity = 2.25E+06_dp
       soil_param%ThermalDiffusivity = 7.4E-7_dp
    ELSEWHERE(FLOOR(fao) == 1)
       soil_param%VolHeatCapacity = 1.93E+06_dp
       soil_param%ThermalDiffusivity = 8.7E-7_dp
    ELSEWHERE(FLOOR(fao) == 2)
       soil_param%VolHeatCapacity = 2.10E+06_dp
       soil_param%ThermalDiffusivity = 8.0E-7_dp
    ELSEWHERE(FLOOR(fao) == 3)
       soil_param%VolHeatCapacity = 2.25E+06_dp
       soil_param%ThermalDiffusivity = 7.4E-7_dp
    ELSEWHERE(FLOOR(fao) == 4)
       soil_param%VolHeatCapacity = 2.36E+06_dp
       soil_param%ThermalDiffusivity = 7.1E-7_dp
    ELSEWHERE(CEILING(fao) == 5)
       soil_param%VolHeatCapacity = 2.48E+06_dp
       soil_param%ThermalDiffusivity = 6.7E-7_dp
    ELSEWHERE
       soil_param%VolHeatCapacity = -1.0_dp
       soil_param%ThermalDiffusivity = -1.0_dp
    END WHERE
    IF (ANY(soil_param%VolHeatCapacity < 0.0_dp)) &
         CALL finish('init_soil','FAO soil types must be between 0 and 5')
    DEALLOCATE(fao)