updating map and newton cotes test case. Updating ROM and analytic aLIGO tilt cavity scan test

kat_test/physics/aligo_BH_tilt.kat

+%--------------------------------------------------------------------------
+% This is a cut down version fo the
+% Finesse model for the full dual recycled inteferometer, with Fabry-Perot
+% arm cavities and OMC.
+%
+% This tests 2 ROM files. Each has a 1nrad in xbeta and -0.5nrad in ybeta
+% on both ITM and ETM, there is also a square aperture included.
+%
+% analytic counterpart to the test file aligo_rom_tilt.kat
+%--------------------------------------------------------------------------
+
+l l1 2500 0 nITMYs2
+
+cav cav_yarm ITMYHR nITMY2 ETMYHR nETMY1
+
+m1 ITMYHR 0.014 $Mloss $phi_ITMY nITMYs2 nITMY2
+attr ITMYHR Rc -1934
+
+# Y arm length
+s LYarm 3994.5 nITMY2 nETMY1
+
+# Y arm end mirror
+m1 ETMYHR 5u $Mloss $phi_ETMY nETMY1 dump
+attr ETMYHR Rc 2245
+
+const Mloss 37.5u
+const phi_ETMY 0
+const phi_ITMY 0
+
+maxtem 6
+
+# map ITMYHR ITM.rom
+# map ETMYHR ETM.rom
+
+attr ETMYHR xbeta 1n
+attr ETMYHR ybeta -0.5n
+conf ETMYHR square_aperture 1
+conf ETMYHR knm_flags 0
+attr ETMYHR r_ap 0.16
+attr ITMYHR xbeta 1n
+attr ITMYHR ybeta -0.5n
+conf ITMYHR square_aperture 1
+conf ITMYHR knm_flags 0
+attr ITMYHR r_ap 0.16
+
+
+ad a00 0 0 0 nETMY1
+ad a10 1 0 0 nETMY1
+ad a01 0 1 0 nETMY1
+ad a20 2 0 0 nETMY1
+ad a02 0 2 0 nETMY1
+ad a40 4 0 0 nETMY1
+ad a04 0 4 0 nETMY1
+ad a11 1 1 0 nETMY1
+ad a22 2 2 0 nETMY1
+
+xaxis ETMYHR phi lin 0 180 1000
+
+yaxis log abs:deg
+gnuterm no
+time

kat_test/physics/aligo_rom_tilt.kat

@@ -26,7 +26,6 @@ m1 ETMYHR 5u $Mloss $phi_ETMY nETMY1 dump
 attr ETMYHR Rc 2245
 
 const Mloss 37.5u
-const nsilica 1.44963098985906
 const phi_ETMY 0
 const phi_ITMY 0
 
@@ -59,8 +58,6 @@ ad a22 2 2 0 nETMY1
 
 xaxis ETMYHR phi lin 0 180 1000
 
-yaxis log abs
-
+yaxis log abs:deg
+gnuterm no
 time
-
-trace 8

kat_test/physics/tilt_newton_cotes.kat

@@ -18,15 +18,17 @@ conf m1 integration_method 4
 conf m1 integration_NC_order 1
 
 attr m1 xbeta -10u
-attr m1 ybeta -10u
+attr m1 ybeta -5u
 
-beam b1 0 n0
-beam b2 0 n3
+beam b1x 0 n0
+beam b2x 0 n3
+beam b1y 0 n0
+beam b2y 0 n3
 
-xaxis b1 x lin -3 3 100
-put b2 x $x1
 
-x2axis b1 y lin -3 3 100
-put b2 y $x2
+xaxis b1x x lin -3 3 100
+put b2x x $x1
+put b1y y $x1
+put b2y y $x1
 
 gnuterm no