function pbeDef = paracetamolGrowth_2_volume()

dataApprox = load('./paracetamolGrowth_2_volume_relmax_r1000.matx', '-mat');
dataApprox = dataApprox.dataApprox;
pbeDef = paracetamolGrowth;

pbeDef.dec = dataApprox.dec;

pbeDef.seeding.positionMatrix = ...
    pbeDef.dec.properData.mappingProperToNew*pbeDef.seeding.positionMatrix;

pbeDef.growthRate = @(H, P, y, t)(...
    pbeDef.dec.properData.mappingProperToNew*...
    pbeDef.growthRate(pbeDef.dec.properData.mappingNewToProper*H, ...
        P, y, t));
  
pbeDef.pivotProperty = @(H, up)([up; ...
  computeMeasureDec(pbeDef.dec.properData, 'hc', ...
  {up, pbeDef.dec.properData.mappingNewToProper*H}, 'volume', [], 0)]);
