diff --git a/input/beam.data b/input/beam.data
new file mode 100755
index 0000000..7cc9799
--- /dev/null
+++ b/input/beam.data
@@ -0,0 +1,3 @@
+170               : f (GHz)
+950 0 62 :x00  y00 z00 (cm)  : mirror position
+2.1 2.1 162. 162.  0 :w0xt w0yt d0xt d0yt phiwr  (cm)  [d0>0 towards plasma]
diff --git a/input/gray_params.data b/input/gray_params.data
new file mode 100755
index 0000000..a696a9e
--- /dev/null
+++ b/input/gray_params.data
@@ -0,0 +1,28 @@
+8 12 1           : nray ktx,rhomx
+1   1   0        : igrad (0  rtr - 1 gauss beam) ipass [=1,2 passes in the plasma] ipol
+.2 5000 0        : dst (cm), nsteps < = 8000, , ist (0,1,2 integration in s,ct,Sr)
+2  4  -20        : iwarm= 0-3 ECRH&CD computation,  ilarm  order of larmor expans, imx
+11               : ieccd 0/1 NO/YES  ECCD calculation
+
+5 20             x alfac, betac poloidal and toroidal launching angles 
+1                x P0 (MW)   injected power
+1   0. 0.        x IOX = 1  :  O-MODE  IOX = 2  :  X-MODE, psipol, chipol
+
+0                x ibeam=0/1/2 beam file with single beam shape (0) or with 1D/2D (1/2) table
+beam.data        x filename for beam data
+
+2                x iequil= 0 vacuum, 1 analytical equilibrium, 2 EQDSK
+equiln.txt       x filename for equilibrium
+11 0  1 0        x icocos [=0..8,10..18], ipsinorm =0/1 [psi normalized (1) or not (0) in EQDSK], idesc =0/1, ifreefmt= 0/1
+0  0.005  0.01   : indXpoint, sspl spline coeff psi, ssplf spline coeff fpol
+-1 -1  1         x signum B_phi I_phi +1=counterclockwise , -1=clockwise from above, factb
+
+4.               x rwall [inner wall radius, m]
+
+1  2             x iprof= 0/1 analytical/numerical profiles, irho=0/1/2 radial coordinate rhot/rhop/psi
+profilesn.txt    x filename for profiles
+1.05     0.005   : psi plasma boundary, sspld spline coeff density
+1 1 1            x factT, factn, iscal=0/1/2 nustar=const/ngreenw=const/no_scal.
+
+1  501           x ipec=±1/2 profiles rhop/rhot; nnd=no. radial points (no. of radial intervals +1)
+5 5              : istprj,istpl