Enables empty contrast models for domains calculations (#442)

Author: DrPaulSharp

examples/domains/standardLayers/domainsStandardLayersScript.m

@@ -86,7 +86,7 @@
     'resample',     false,...
     'BulkIn',       'SLD Air',...        
     'BulkOut',      'SLD D2O',...        
-    'domainRatio',  'Domain Ratio 1',....
+    'domainRatio',  'Domain Ratio 1',...
     'data',         'Simulation',...
     'model',        {'Domain1', 'Domain2'});
 

targetFunctions/+domainsTF/domainsReflectivity.m

@@ -100,10 +100,18 @@
                 % layers list are required for this contrast, and put them
                 % in the correct order according to geometry. We run it
                 % twice, once for each domain
+                if isempty(inputContrastLayers{i})
+                    inputDomainContrastLayers1 = [];
+                    inputDomainContrastLayers2 = [];
+                else
+                    inputDomainContrastLayers1 = domainContrastLayers{inputContrastLayers{i}(1)};
+                    inputDomainContrastLayers2 = domainContrastLayers{inputContrastLayers{i}(2)};
+                end
+
                 contrastLayers1{i} = allocateLayersForContrast(...
-                    domainContrastLayers{inputContrastLayers{i}(1)}, layerValues);
+                    inputDomainContrastLayers1, layerValues);
                 contrastLayers2{i} = allocateLayersForContrast(...
-                    domainContrastLayers{inputContrastLayers{i}(2)}, layerValues);
+                    inputDomainContrastLayers2, layerValues);
 
             end