The standard approach to obtaining knowledge about the properties of the surface layer of a comet from observations of gas production consists of two stages. First, various thermophysical models are used to calculate gas production for a few sets of parameters. Second, a comparison of observations and theoretical predictions is performed. This approach is complicated because the values of many model characteristics are known only approximately. Therefore, it is necessary to investigate the sensitivity of the simulated outgassing to variations in the properties of the surface layer. This problem was recently considered by us for aggregates up to tens of microns in size. For millimetre-size aggregates, a qualitative extension of the method us...