Index: sansmodels/src/include/refl.h
===================================================================
--- sansmodels/src/include/refl.h	(revision 101065a9bd217be61773f8eba38d0c461efee3fd)
+++ sansmodels/src/include/refl.h	(revision 101065a9bd217be61773f8eba38d0c461efee3fd)
@@ -0,0 +1,154 @@
+// The original code, of which work was not DANSE funded,
+// was provided by J. Cho.
+#if !defined(o_h)
+#define refl_h
+#include "parameters.hh"
+
+/**
+ * Structure definition for sphere parameters
+ */
+ //[PYTHONCLASS] = ReflModel
+ //[DISP_PARAMS] = thick_inter0
+ //[DESCRIPTION] =<text>Calculate neutron reflectivity using the Parratt iterative formula
+ //				Parameters:
+ //				background:background
+ //				scale: scale factor
+ //				sld_bottom0: the SLD of the substrate
+ //				sld_medium: the SLD of the incident medium
+ //					or superstrate
+ //				sld_flatN: the SLD of the flat region of
+ //					the N'th layer
+ //				thick_flatN: the thickness of the flat
+ //					region of the N'th layer
+ //				func_interN: the function used to describe
+ //					the interface of the N'th layer
+ //				thick_interN: the thickness of the interface
+ //					of the N'th layer
+ //				Note: the layer number starts to increase
+ //					from the bottom (substrate) to the top.
+ //		</text>
+ //[FIXED]=  <text></text>
+ //[NON_FITTABLE_PARAMS]= <text>n_layers;func_inter0;func_inter1;func_inter2;func_inter3;func_inter4;func_inter5;func_inter5;func_inter7;func_inter8;func_inter9;func_inter10 </text>
+ //[ORIENTATION_PARAMS]= <text> </text>
+
+class ReflModel{
+public:
+  // Model parameters
+  /// number of layers
+  //  [DEFAULT]=n_layers=1
+  Parameter n_layers;
+    /// Scale factor
+    //  [DEFAULT]=scale= 1.0
+  Parameter scale;
+    /// thick_inter0 [A]
+    //  [DEFAULT]=thick_inter0=1.0 [A]
+  Parameter thick_inter0;
+  /// func_inter0
+  //  [DEFAULT]=func_inter0= 0
+  Parameter func_inter0;
+  /// sld_bottom0 [1/A^(2)]
+  //  [DEFAULT]=sld_bottom0= 2.07e-6 [1/A^(2)]
+  Parameter sld_bottom0;
+  /// sld_medium [1/A^(2)]
+  //  [DEFAULT]=sld_medium= 1.0e-6 [1/A^(2)]
+  Parameter sld_medium;
+  /// Background
+  //  [DEFAULT]=background=0
+  Parameter background;
+
+    //  [DEFAULT]=sld_flat1=4.0e-06 [1/A^(2)]
+    Parameter sld_flat1;
+    //  [DEFAULT]=sld_flat2=3.5e-06 [1/A^(2)]
+    Parameter sld_flat2;
+    //  [DEFAULT]=sld_flat3=4.0e-06 [1/A^(2)]
+    Parameter sld_flat3;
+    //  [DEFAULT]=sld_flat4=3.5e-06 [1/A^(2)]
+    Parameter sld_flat4;
+    //  [DEFAULT]=sld_flat5=4.0e-06 [1/A^(2)]
+    Parameter sld_flat5;
+    //  [DEFAULT]=sld_flat6=3.5e-06 [1/A^(2)]
+    Parameter sld_flat6;
+    //  [DEFAULT]=sld_flat7=4.0e-06 [1/A^(2)]
+    Parameter sld_flat7;
+    //  [DEFAULT]=sld_flat8=3.5e-06 [1/A^(2)]
+    Parameter sld_flat8;
+    //  [DEFAULT]=sld_flat9=4.0e-06 [1/A^(2)]
+    Parameter sld_flat9;
+    //  [DEFAULT]=sld_flat10=3.5e-06 [1/A^(2)]
+    Parameter sld_flat10;
+
+    //  [DEFAULT]=thick_inter1=1 [A]
+    Parameter thick_inter1;
+    //  [DEFAULT]=thick_inter2=1 [A]
+    Parameter thick_inter2;
+    //  [DEFAULT]=thick_inter3=1 [A]
+    Parameter thick_inter3;
+    //  [DEFAULT]=thick_inter4=1 [A]
+    Parameter thick_inter4;
+    //  [DEFAULT]=thick_inter5=1 [A]
+    Parameter thick_inter5;
+    //  [DEFAULT]=thick_inter6=1 [A]
+    Parameter thick_inter6;
+    //  [DEFAULT]=thick_inter7=1 [A]
+    Parameter thick_inter7;
+    //  [DEFAULT]=thick_inter8=1 [A]
+    Parameter thick_inter8;
+    //  [DEFAULT]=thick_inter9=1 [A]
+    Parameter thick_inter9;
+    //  [DEFAULT]=thick_inter10=1 [A]
+    Parameter thick_inter10;
+
+    //  [DEFAULT]=thick_flat1=10 [A]
+    Parameter thick_flat1;
+    //  [DEFAULT]=thick_flat2=100 [A]
+    Parameter thick_flat2;
+    //  [DEFAULT]=thick_flat3=100 [A]
+    Parameter thick_flat3;
+    //  [DEFAULT]=thick_flat4=100 [A]
+    Parameter thick_flat4;
+    //  [DEFAULT]=thick_flat5=100 [A]
+    Parameter thick_flat5;
+    //  [DEFAULT]=thick_flat6=100 [A]
+    Parameter thick_flat6;
+    //  [DEFAULT]=thick_flat7=100 [A]
+    Parameter thick_flat7;
+    //  [DEFAULT]=thick_flat8=100 [A]
+    Parameter thick_flat8;
+    //  [DEFAULT]=thick_flat9=100 [A]
+    Parameter thick_flat9;
+    //  [DEFAULT]=thick_flat10=100 [A]
+    Parameter thick_flat10;
+
+    //  [DEFAULT]=func_inter1=0
+    Parameter func_inter1;
+    //  [DEFAULT]=func_inter2=0
+    Parameter func_inter2;
+    //  [DEFAULT]=func_inter3=0
+    Parameter func_inter3;
+    //  [DEFAULT]=func_inter4=0
+    Parameter func_inter4;
+    //  [DEFAULT]=func_inter5=0
+    Parameter func_inter5;
+    //  [DEFAULT]=func_inter6=0
+    Parameter func_inter6;
+    //  [DEFAULT]=func_inter7=0
+    Parameter func_inter7;
+    //  [DEFAULT]=func_inter8=0
+    Parameter func_inter8;
+    //  [DEFAULT]=func_inter9=0
+    Parameter func_inter9;
+    //  [DEFAULT]=func_inter10=0
+    Parameter func_inter10;
+
+
+  // Constructor
+  ReflModel();
+
+  // Operators to get I(Q)
+  double operator()(double q);
+  double operator()(double qx, double qy);
+  double calculate_ER();
+  double evaluate_rphi(double q, double phi);
+};
+
+#endif
