From 006b54ead6e6c771bad945da2322a5af3a610106 Mon Sep 17 00:00:00 2001
From: Samuel Lazerson <lazersos@gmail.com>
Date: Tue, 5 Mar 2024 10:15:52 +0100
Subject: [PATCH] Added 120kV OPTEMIST beamline

---
 BEAMS3D/w7x_beam_params.m | 31 +++++++++++++++++++++++++++++++
 1 file changed, 31 insertions(+)

diff --git a/BEAMS3D/w7x_beam_params.m b/BEAMS3D/w7x_beam_params.m
index d42140f..fea39d0 100644
--- a/BEAMS3D/w7x_beam_params.m
+++ b/BEAMS3D/w7x_beam_params.m
@@ -32,6 +32,7 @@ function w7x_beam_params( source ,varargin)
 lrudix = 0;
 loptemist = 0;
 loptemist150kev = 0;
+loptemist120kev = 0;
 loptemist300kev = 0;
 lwrite_beams3d = 0;
 grid = 60;
@@ -80,6 +81,11 @@ function w7x_beam_params( source ,varargin)
                     vmec_data.rmax_surf = 6.5; %Maybe should even be after this section, to avoid WALL OUTSIDE GRID DOMAIN! error
                     vmec_data.rmin_surf = 4.25;
                     vmec_data.zmax_surf = 1.25;
+                case {'optemist120kev','OPTEMIST120KEV','Optemist120keV'}
+                    loptemist120kev=1;
+                    vmec_data.rmax_surf = 6.5; %Maybe should even be after this section, to avoid WALL OUTSIDE GRID DOMAIN! error
+                    vmec_data.rmin_surf = 4.25;
+                    vmec_data.zmax_surf = 1.25;
                 case {'optemist150kev','OPTEMIST150KEV','Optemist150keV'}
                     loptemist150kev=1;
                     vmec_data.rmax_surf = 6.5; %Maybe should even be after this section, to avoid WALL OUTSIDE GRID DOMAIN! error
@@ -390,6 +396,31 @@ function w7x_beam_params( source ,varargin)
     end
 end
 
+if loptemist120kev
+    E_OPT = [120] * 1000; % Nominal values
+    PFRAC = [99.8 0.1 0.1]./100;
+    P_OPT = 1E6; %1 MW 
+    for i=1:length(E_OPT)
+        r_beam(1,j) = sqrt(xo_RUDI.^2+yo_RUDI.^2); %RUDIX Geometry
+        r_beam(2,j) = sqrt(xt_RUDI.^2+yt_RUDI.^2);
+        p_beam(1,j) = atan2(yo_RUDI,xo_RUDI);
+        p_beam(2,j) = atan2(yt_RUDI,xt_RUDI);
+        z_beam(1,j) = zo_RUDI;
+        z_beam(2,j) = zt_RUDI;
+        power_beam(j) = P_OPT;
+        energy_beam(j) = E_OPT(i);
+        div_beam(j)    = 0.0125;
+        note{j} = 'OPTEMIST 120keV BEAM';
+        source = [source j]; % Treat as source 9
+        j = j + 1;
+    end
+    if lplots
+        hold(ax, 'on');
+        plot3(ax,[xo_OPT xt_OPT],[yo_OPT yt_OPT],[zo_OPT zt_OPT],'k');
+        hold(ax,'off');
+    end
+end
+
 if loptemist300kev
     E_OPT = [300] * 1000; % Nominal values
     PFRAC = [34 33 33]./100;