diff --git a/examples/domains/standardLayers/domainsStandardLayersScript.m b/examples/domains/standardLayers/domainsStandardLayersScript.m
index 27aacf398..aefdd2a7e 100644
--- a/examples/domains/standardLayers/domainsStandardLayersScript.m
+++ b/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'});
 
diff --git a/targetFunctions/+domainsTF/domainsReflectivity.m b/targetFunctions/+domainsTF/domainsReflectivity.m
index 507a1aa86..ee75899de 100644
--- a/targetFunctions/+domainsTF/domainsReflectivity.m
+++ b/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