gs2crmod 0.11.67 → 0.11.68

data/VERSION

@@ -1 +1 @@
-0.11.67
+0.11.68

data/gs2crmod.gemspec

@@ -5,7 +5,7 @@
 
 Gem::Specification.new do |s|
   s.name = "gs2crmod"
-  s.version = "0.11.67"
+  s.version = "0.11.68"
 
   s.required_rubygems_version = Gem::Requirement.new(">= 0") if s.respond_to? :required_rubygems_version=
   s.authors = ["Edmund Highcock", "Ferdinand van Wyk"]

data/lib/gs2crmod/graphs.rb

@@ -10,53 +10,54 @@ DataKit = GraphKit::DataKit
 
 def auto_axiskits(name, options)
 	hash = cache[:auto_axiskits] ||= {'t' => ['Time', ''],
-                'phi2tot_over_time' => ['Phi^2 Total', ''],
-                'apar2_over_time' => ['Apar^2 Total', ''],
-                'growth_rate_by_ky_over_time' => ['Growth Rate by ky', ''],
-                 'growth_rate_by_kx_over_time' => ['Growth Rate by kx', ''],  
-		 'growth_rate_by_mode_over_time' => ["Growth Rate by mode", ''],
-# <MJL additions 2013-09-19>
-                 'frequency_by_ky_over_time' => ['Real frequency by ky', ''],
-                  'frequency_by_kx_over_time' => ['Real frequency by kx', ''],
-# </MJL>
-                'phi2_by_ky_over_time' => ['Phi^2 by ky', ''],
-                 'phi2_by_kx_over_time' => ['Phi^2 by ky', ''],  
-                'es_heat_by_ky_over_time' => ['Phi^2 by ky', ''],
-                 'es_heat_by_kx_over_time' => ['Phi^2 by kx', ''],  
-		 'phi2_by_mode_over_time' => ["Phi^2 by mode", ''],
-	 'tpar2_by_mode_over_time' => ["(delta T_parallel)^2 by mode", '%'],
- 		'tperp2_by_mode_over_time' => ["(delta T_perp)^2 by mode", '%'],
-                              'hflux_tot' => ['Total Heat Flux', ''],
-                                'es_heat_par' => ['Parallel electrostatic heat flux', ''],
-                                'es_heat_perp' => ['Perpendicular electrostatic heat flux', ''],
-                'ky' => ['ky', "1/rho_#{species_letter}"],
-                'kx' => ['kx', "1/rho_#{species_letter}"],
-	        'kpar' => ['kpar', "2 pi/qR"],
-	        'growth_rate_over_kx' => ['Growth Rate', "v_th#{species_letter}/a", 1],
-	        'growth_rate_over_ky' => ['Growth Rate', "v_th#{species_letter}/a", 1],
-	        'growth_rate_over_kx_slice' => ['Growth Rate', "v_th#{species_letter}/a", 1],
-	        'growth_rate_over_ky_slice' => ['Growth Rate', "v_th#{species_letter}/a", 1],
-	        'growth_rate_over_ky_over_kx' => ["Growth Rate", "v_th#{species_letter}/a", 2],
-          'frequency_over_ky' => ['Frequency', "v_th#{species_letter}/a", 1],
-	        'transient_es_heat_flux_amplification_over_kx' => ['Transient Electrostatic Heat Amplification', "", 1],
-	        'transient_es_heat_flux_amplification_over_ky' => ['Transient Electrostatic Heat Amplification', "", 1],
-	        'transient_amplification_over_kx' => ['Transient Amplification', "", 1],
-	        'transient_amplification_over_ky' => ['Transient Amplification', "", 1],
-	        'spectrum_over_kx' => ["Spectrum at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", '', 1],
-	        'zonal_spectrum' => ["Zonal spectrum at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", '', 1],
-	        'spectrum_over_ky' => ["Spectrum at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", '', 1],
-	        'es_heat_over_kx' => ["Heat Flux at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", 'Q_gB', 1],
-	        'es_heat_over_ky' => ["Heat Flux at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", 'Q_gB', 1],
-	       	'es_heat_flux_over_ky_over_kx' => ["Heat flux at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", '', 2],
-	       	'spectrum_over_kpar' => ["Spectrum at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", '', 1],
-	       	'spectrum_over_ky_over_kx' => ["Spectrum at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", '', 2],
-	       	'spectrum_over_ky_over_kpar' => ["Spectrum at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", '', 2],
-	        #'phi0_over_x_over_y' => ["Phi at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", '', 2],
-	        'phi0_over_x_over_y' => ["Phi at theta = 0", '', 2],
-	        'es_mom_flux_over_time' => ["#{species_type((options[:species_index] or 1)).capitalize} Momentum Flux", '', 1]
-
-                  
-                   }
+    'phi2tot_over_time' => ['Phi^2 Total', ''],
+    'apar2_over_time' => ['Apar^2 Total', ''],
+    'growth_rate_by_ky_over_time' => ['Growth Rate by ky', ''],
+    'growth_rate_by_kx_over_time' => ['Growth Rate by kx', ''],  
+    'growth_rate_by_mode_over_time' => ["Growth Rate by mode", ''],
+    # <MJL additions 2013-09-19>
+    'frequency_by_ky_over_time' => ['Real frequency by ky', ''],
+    'frequency_by_kx_over_time' => ['Real frequency by kx', ''],
+    # </MJL>
+    'phi2_by_ky_over_time' => ['Phi^2 by ky', ''],
+    'phi2_by_kx_over_time' => ['Phi^2 by ky', ''],  
+    'es_heat_by_ky_over_time' => ['Phi^2 by ky', ''],
+    'es_heat_by_kx_over_time' => ['Phi^2 by kx', ''],  
+    'phi2_by_mode_over_time' => ["Phi^2 by mode", ''],
+    'tpar2_by_mode_over_time' => ["(delta T_parallel)^2 by mode", '%'],
+    'tperp2_by_mode_over_time' => ["(delta T_perp)^2 by mode", '%'],
+    'hflux_tot' => ['Total Heat Flux', ''],
+    'es_heat_par' => ['Parallel electrostatic heat flux', ''],
+    'es_heat_perp' => ['Perpendicular electrostatic heat flux', ''],
+    'ky' => ['ky', "1/rho_#{species_letter}"],
+    'kx' => ['kx', "1/rho_#{species_letter}"],
+    'x' => ['x', "rho_#{species_letter}", 1],
+    'kpar' => ['kpar', "2 pi/qR"],
+    'growth_rate_over_kx' => ['Growth Rate', "v_th#{species_letter}/a", 1],
+    'growth_rate_over_ky' => ['Growth Rate', "v_th#{species_letter}/a", 1],
+    'growth_rate_over_kx_slice' => ['Growth Rate', "v_th#{species_letter}/a", 1],
+    'growth_rate_over_ky_slice' => ['Growth Rate', "v_th#{species_letter}/a", 1],
+    'growth_rate_over_ky_over_kx' => ["Growth Rate", "v_th#{species_letter}/a", 2],
+    'frequency_over_ky' => ['Frequency', "v_th#{species_letter}/a", 1],
+    'transient_es_heat_flux_amplification_over_kx' => ['Transient Electrostatic Heat Amplification', "", 1],
+    'transient_es_heat_flux_amplification_over_ky' => ['Transient Electrostatic Heat Amplification', "", 1],
+    'transient_amplification_over_kx' => ['Transient Amplification', "", 1],
+    'transient_amplification_over_ky' => ['Transient Amplification', "", 1],
+    'spectrum_over_kx' => ["Spectrum at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", '', 1],
+    'spectrum_over_ky' => ["Spectrum at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", '', 1],
+    'es_heat_over_kx' => ["Heat Flux at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", 'Q_gB', 1],
+    'es_heat_over_ky' => ["Heat Flux at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", 'Q_gB', 1],
+    'es_heat_flux_over_ky_over_kx' => ["Heat flux at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", '', 2],
+    'spectrum_over_kpar' => ["Spectrum at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", '', 1],
+    'spectrum_over_ky_over_kx' => ["Spectrum at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", '', 2],
+    'spectrum_over_ky_over_kpar' => ["Spectrum at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", '', 2],
+    #'phi0_over_x_over_y' => ["Phi at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", '', 2],
+    'phi0_over_x_over_y' => ["Phi at theta = 0", '', 2],
+    'es_mom_flux_over_time' => ["#{species_type((options[:species_index] or 1)).capitalize} Momentum Flux", '', 1],
+    'zonal_spectrum' => ["Zonal spectrum at t = #{sprintf("%.3f" ,(options[:t] or list(:t)[options[:t_index]] or list(:t).values.max))}", '', 1],
+    'zonal_flow_velocity_over_x' => ['Zonal Flow Velocity', "", 1],
+    'mean_flow_velocity_over_x' => ['Mean Flow Velocity', "", 1]
+  }
 	return hash[name]
 end
 
@@ -165,34 +166,25 @@ GRAPHKIT_OPTIONS_HELP = {
 	no_kpar0: "Don't plot the kpar=0 part (boolean, e.g. no_kpar0: true)",
 	log: "Plot the log of a given quantity (exact meaning varies). boolean",
 	Rmaj: "The major radius in metres. This has no effect on the shape of the graph: it merely multiplies every length",
- n0: " The toroidal mode number of the longest y mode. In effect it is the number of periodic copies of the flux tube that will fit in the torus. Periodicity requires that n0 q  is also an integer. If you specify :n0 where this is not the case, q will automatically be adjusted until it is",
- rho_star: " The ratio of the reference Lamour radius to the GS2 normalising length a. Cannot be specified at the same time as n0. If specified, both n0 and q will be adjusted to ensure periodicity",
- t_index: "The (1-based) time index",
- nakx: "The number of radial wave numbers to include in the plot. In effect, it is a low pass filter which reduces the resolution in the radial direction without changing the shape of the final surface. Minimum value is 4",
- naky: "The number of kys to include in the plot. In effect, it is a low pass filter which reduces the resolution in the y direction without changing the shape of the final surface. Minimum value is 4",
- gs2_coordinate_factor: "When set to 1, plot the graph in GS2 coordinates. When set to  0 plot the graph in real space. Can be set at any value between 0 and 1: the graph will smoothly distort between the two limits",
- xmax: "The (0-based) index of the maximum value of x to include in the plot",
- xmin: "The (0-based) index of the minimum value of x to include in the plot",
- ymax: "The (0-based) index of the maximum value of y to include in the plot",
- ymin: "The (0-based) index of the minimum value of y to include in the plot",
- thetamax: "The (0-based) index of the maximum value of theta to include in the plot",
- thetamin: "The (0-based) index of the minimum value of theta to include in the plot",
- ncopies: " The number of periodic copies of the flux tube to include",
- torphi_values: "An array of two values of the toroidal angle. The graph will be plotted in between those two values with poloidal cross sections at either end",
- magnify: " The magnification factor of the small section. It can take any value greater than or equal to 1",
-
+  n0: " The toroidal mode number of the longest y mode. In effect it is the number of periodic copies of the flux tube that will fit in the torus. Periodicity requires that n0 q  is also an integer. If you specify :n0 where this is not the case, q will automatically be adjusted until it is",
+  rho_star: " The ratio of the reference Lamour radius to the GS2 normalising length a. Cannot be specified at the same time as n0. If specified, both n0 and q will be adjusted to ensure periodicity",
+  t_index: "The (1-based) time index",
+  nakx: "The number of radial wave numbers to include in the plot. In effect, it is a low pass filter which reduces the resolution in the radial direction without changing the shape of the final surface. Minimum value is 4",
+  naky: "The number of kys to include in the plot. In effect, it is a low pass filter which reduces the resolution in the y direction without changing the shape of the final surface. Minimum value is 4",
+  gs2_coordinate_factor: "When set to 1, plot the graph in GS2 coordinates. When set to  0 plot the graph in real space. Can be set at any value between 0 and 1: the graph will smoothly distort between the two limits",
+  xmax: "The (0-based) index of the maximum value of x to include in the plot",
+  xmin: "The (0-based) index of the minimum value of x to include in the plot",
+  ymax: "The (0-based) index of the maximum value of y to include in the plot",
+  ymin: "The (0-based) index of the minimum value of y to include in the plot",
+  thetamax: "The (0-based) index of the maximum value of theta to include in the plot",
+  thetamin: "The (0-based) index of the minimum value of theta to include in the plot",
+	theta_index: "integer, index of theta at which to plot (e.g. theta_index: 20)",
+	kxfac: "float, overrides calculation of kxfac in zonal flow velocity function",
+  ncopies: " The number of periodic copies of the flux tube to include",
+  torphi_values: "An array of two values of the toroidal angle. The graph will be plotted in between those two values with poloidal cross sections at either end",
+  magnify: " The magnification factor of the small section. It can take any value greater than or equal to 1",
 }
 	
-
-# def graphkit(name, options={})
-# 	unless [:Failed, :Complete].include? status
-# 		return get_graphkit(name, options)
-# 	else
-# 		return cache[[:graphkit, name, options]] ||= get_graphkit(name, options)
-# 	end
-# end
-
-
 def graphkit(name, options={})
 	logf :graphkit
 	# If an array of t, kx or ky values is provided, plot one graph for each value and then sum the graphs together
@@ -2375,6 +2367,39 @@ module GraphKits
 		end
 	end
 
+	def zonal_flow_velocity_vs_x_graphkit(options={})
+		case options[:command]
+		when :help
+			return  "zonal_flow_velocity_vs_x: Graph of the zonal flow velocity kxfac*IFT(i k_x phi). kxfac = (qinp/rhoc)*grho(rhoc)."
+		when :options
+			return  [:t, :t_index, :theta_index, :kxfac]
+		else
+      options[:ky_index]=0
+			kit = GraphKit.autocreate({x: axiskit('x', options), y: axiskit("zonal_flow_velocity_over_x", options)})
+			kit.title  = "Zonal Flow Velocity versus x"
+			kit.file_name = options[:graphkit_name] + options[:t_index].to_s
+			kit.data[0].with = 'lp'
+			kit.pointsize = 2.0
+			kit
+		end
+	end
+
+	def mean_flow_velocity_vs_x_graphkit(options={})
+		case options[:command]
+		when :help
+			return  "mean_flow_velocity_vs_x: Graph of the mean flow velocity (x - x(centre))*g_exb"
+		when :options
+			return  [:t, :t_index, :theta_index]
+		else
+      options[:ky_index]=0
+			kit = GraphKit.autocreate({x: axiskit('x', options), y: axiskit("mean_flow_velocity_over_x", options)})
+			kit.title  = "Mean Flow Velocity versus x"
+			kit.file_name = options[:graphkit_name] + options[:t_index].to_s
+			kit.data[0].with = 'lp'
+			kit.pointsize = 2.0
+			kit
+		end
+	end
 
 end
 

data/lib/gs2crmod/gsl_data.rb

@@ -990,18 +990,65 @@ module GSLVectors
 			ysize = ny*2-2+ny%2
 			GSL::Vector.indgen(ysize, 0, ly/ysize)
 		end
+
+    #This function reads in the 'grho' variable from the netcdf file.
+    def grho_gsl_vector(options)
+      grho = GSL::Vector.alloc(netcdf_file.var('grho').get('start' => [0], 'end' => [-1]).to_a)
+      return grho
+    end
+
+    #This function returns the zonal flow velocity as a function of x (the radial coordinate). 
+    #This is v_ZF = kxfac*IFT(i k_x phi_imag), where kxfac = (qinp/rhoc)*grho(rhoc). 
+		def zonal_flow_velocity_over_x_gsl_vector(options)
+			Dir.chdir(@directory) do
+        raise CRFatal.new("Need to specify a theta_index.") unless options[:theta_index]
+        raise CRFatal.new("Need either qinp or pk and epsl specified in order to calculate kxfac.
+                          If using numerical equil use the option :kxfac to override calculation.") unless @qinp or (@pk and @epsl or options[:kxfac])
+				phi = gsl_vector_complex('phi_zonal', options)
+
+        kx = gsl_vector(:kx).to_box_order
+        x = gsl_vector(:x)
+        grho = gsl_vector('grho')[options[:theta_index]]
+        if @qinp
+          kxfac = (@qinp/@rhoc)*grho
+        elsif @pk and @epsl
+          kxfac = (@epsl/@pk)*grho
+        elsif options[:kxfac]
+          kxfac = options[:kxfac]
+        else
+          raise 'Error: Need qinp or pk and epsl to calculate kxfac' 
+        end
+
+        vec_zf_vel = GSL::Vector.alloc(kx.size)
+        #Take imaginary part since i k_x will lead to imaginary part being real
+        vec_zf_vel = kxfac*(phi*kx).backward.imag*kx.size
+        return vec_zf_vel
+			end
+		end
+
+    #This function returns the mean flow velocity as a function of x (the radial coordinate). 
+    #This is v_g_exb = (x - x(centre))*g_exb. The x-x(centre) ensures that the flow is zero 
+    #at the middle of the box.
+		def mean_flow_velocity_over_x_gsl_vector(options)
+			Dir.chdir(@directory) do
+				raise CRFatal.new("Need to have g_exb > 0 to have a mean flow.") unless @g_exb > 0
+        raise CRFatal.new("Need to specify a theta_index.") unless options[:theta_index]
+        x = gsl_vector(:x)
+
+        vec_exb_vel = GSL::Vector.alloc(x.size)
+        #Take imaginary part since i k_x will lead to imaginary part being real
+        vec_exb_vel = (x - x[x.size/2])*@g_exb
+        return vec_exb_vel
+			end
+		end
+
 		def zonal_spectrum_gsl_vector(options)
 			Dir.chdir(@directory) do
-		     gmzf = gsl_matrix('spectrum_over_ky_over_kx',options)
-         veczf = GSL::Vector.alloc(gmzf.shape[1])
-                   # p gmzf.get_row(0).size
-                   # p gmzf.get_row(0)
-		     gmzf.shape[1].times{|i| veczf[i] = gmzf[0,i]}
-		     return veczf
-		#else
-			#raise CRError.new("Unknown gsl_vector requested: #{name}")
+        gmzf = gsl_matrix('spectrum_over_ky_over_kx',options)
+        veczf = GSL::Vector.alloc(gmzf.shape[1])
+        gmzf.shape[1].times{|i| veczf[i] = gmzf[0,i]}
+        return veczf
 			end
-	#  			eputs data; gets
 		end
   
 end	# module GSLVectors	
@@ -1076,6 +1123,32 @@ module GSLVectorComplexes
 	#  			eputs data; gets
 	end
 
+  #This function returns a complex GSL vector of the zonal (ky=0) component of phi/phi_t at a given theta index 
+  #and time index if write_phi_over_time was enabled during the simulation
+	def phi_zonal_gsl_vector_complex(options)
+    Dir.chdir(@directory) do
+      if options[:t_index] or options[:t]
+        #extra option required is t_index
+        raise CRFatal.new("write_phi_over_time is not enabled so this function won't work") unless @write_phi_over_time
+        
+        options.convert_to_index(self, :t)
+        a = netcdf_file.var('phi_t').get({
+          'start' => [0,options[:theta_index],0,0, options[:t_index] - 1], 
+          'end' => [-1,options[:theta_index],-1,0, options[:t_index] - 1]
+        })
+        vector = GSL::Vector::Complex.alloc(GSL::Vector.alloc(a[0,0,0..-1,0,0]), GSL::Vector.alloc(a[1,0,0..-1,0,0]))
+        return vector
+      else
+        a = netcdf_file.var('phi').get({
+          'start' => [0, options[:theta_index], 0, 0], 
+          'end' => [-1, options[:theta_index], -1, 0]
+        })
+        vector = GSL::Vector::Complex.alloc(GSL::Vector.alloc(a[0,0,0..-1,0]), GSL::Vector.alloc(a[1,0,0..-1,0]))
+        return vector
+			end
+		end
+	end
+
 end		
 include GSLVectorComplexes
 

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: gs2crmod
 version: !ruby/object:Gem::Version
-  version: 0.11.67
+  version: 0.11.68
   prerelease: 
 platform: ruby
 authors:
@@ -14,7 +14,7 @@ date: 2014-09-03 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: coderunner
-  requirement: &13799260 !ruby/object:Gem::Requirement
+  requirement: &17929040 !ruby/object:Gem::Requirement
     none: false
     requirements:
     - - ! '>='
@@ -22,10 +22,10 @@ dependencies:
         version: 0.14.15
   type: :runtime
   prerelease: false
-  version_requirements: *13799260
+  version_requirements: *17929040
 - !ruby/object:Gem::Dependency
   name: rubyhacks
-  requirement: &13797960 !ruby/object:Gem::Requirement
+  requirement: &17926880 !ruby/object:Gem::Requirement
     none: false
     requirements:
     - - ! '>='
@@ -33,10 +33,10 @@ dependencies:
         version: 0.1.2
   type: :runtime
   prerelease: false
-  version_requirements: *13797960
+  version_requirements: *17926880
 - !ruby/object:Gem::Dependency
   name: ruby-netcdf-updated
-  requirement: &13812880 !ruby/object:Gem::Requirement
+  requirement: &17925740 !ruby/object:Gem::Requirement
     none: false
     requirements:
     - - ! '>='
@@ -44,10 +44,10 @@ dependencies:
         version: 0.6.6.1
   type: :runtime
   prerelease: false
-  version_requirements: *13812880
+  version_requirements: *17925740
 - !ruby/object:Gem::Dependency
   name: shoulda
-  requirement: &13811080 !ruby/object:Gem::Requirement
+  requirement: &17940940 !ruby/object:Gem::Requirement
     none: false
     requirements:
     - - =
@@ -55,10 +55,10 @@ dependencies:
         version: 3.0.1
   type: :development
   prerelease: false
-  version_requirements: *13811080
+  version_requirements: *17940940
 - !ruby/object:Gem::Dependency
   name: rdoc
-  requirement: &13810040 !ruby/object:Gem::Requirement
+  requirement: &17938700 !ruby/object:Gem::Requirement
     none: false
     requirements:
     - - ~>
@@ -66,10 +66,10 @@ dependencies:
         version: '3.12'
   type: :development
   prerelease: false
-  version_requirements: *13810040
+  version_requirements: *17938700
 - !ruby/object:Gem::Dependency
   name: bundler
-  requirement: &13809160 !ruby/object:Gem::Requirement
+  requirement: &17937200 !ruby/object:Gem::Requirement
     none: false
     requirements:
     - - ! '>'
@@ -77,10 +77,10 @@ dependencies:
         version: 1.0.0
   type: :development
   prerelease: false
-  version_requirements: *13809160
+  version_requirements: *17937200
 - !ruby/object:Gem::Dependency
   name: jeweler
-  requirement: &13808120 !ruby/object:Gem::Requirement
+  requirement: &17936160 !ruby/object:Gem::Requirement
     none: false
     requirements:
     - - ! '>='
@@ -88,10 +88,10 @@ dependencies:
         version: 1.8.4
   type: :development
   prerelease: false
-  version_requirements: *13808120
+  version_requirements: *17936160
 - !ruby/object:Gem::Dependency
   name: minitest
-  requirement: &13807320 !ruby/object:Gem::Requirement
+  requirement: &17935100 !ruby/object:Gem::Requirement
     none: false
     requirements:
     - - ~>
@@ -99,7 +99,7 @@ dependencies:
         version: '4'
   type: :development
   prerelease: false
-  version_requirements: *13807320
+  version_requirements: *17935100
 description: GS2 is a gyrokinetic flux tube initial value turbulence code which can
   be used for fusion or astrophysical plasmas. CodeRunner is a framework for the automated
   running and analysis of large simulations. This module allows GS2 (and its sister