Added 120kV OPTEMIST beamline

lazersos · Mar 5, 2024 · 006b54e · 006b54e
1 parent 61bf430
commit 006b54e
Showing 1 changed file with 31 additions and 0 deletions.
diff --git 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;