981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
1001
1002
1003
1004
1005
1006
1007
1008
1009
1010
1011
1012
1013
1014
1015
1016
1017
1018
1019
1020
1021
1022
1023
1024
1025
1026
1027
1028
1029
1030
1031
1032
1033
1034
1035
1036
1037
|
# File 'lib/gs2crmod/gsl_data.rb', line 981
def phi_along_field_line_gsl_vector_complex(options)
Dir.chdir(@directory) do
options.convert_to_index(self, :ky)
if options[:t_index] or options[:t]
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)
case @grid_option
when "single"
temp = GSL::Vector.alloc(netcdf_file.var('phi_t').get({'start' => [0,0,0,0, options[:t_index] - 1], 'end' => [-1,-1,0,0, options[:t_index] - 1]}).to_a[0][0][0].flatten)
when "range"
a = netcdf_file.var('phi_t').get({'start' => [0, 0, options[:kx_index]-1, options[:ky_index] - 1, options[:t_index] - 1], 'end' => [-1, -1, options[:kx_index]-1, options[:ky_index] - 1, options[:t_index]-1]})
temp = GSL::Vector.alloc(a.to_a[0][0].flatten)
when "box"
options.convert_to_index(self, :ky, :kx)
kx_elements = gsl_vector('linked_kx_elements', options).to_a
a = netcdf_file.var('phi_t').get({
'start' => [0,0,0,options[:ky_index] - 1, options[:t_index] - 1],
'end' => [-1,-1,-1, options[:ky_index] - 1, options[:t_index] - 1]
}).to_a[0][0].values_at(*kx_elements).flatten
temp = GSL::Vector.alloc(a)
end
return GSL::Vector::Complex.alloc(temp.subvector_with_stride(0, 2), temp.subvector_with_stride(1, 2))
else
case @grid_option
when "single"
temp = GSL::Vector.alloc(netcdf_file.var('phi').get({'start' => [0,0, 0, 0], 'end' => [-1,-1,0,0]}).to_a.flatten)
when "range"
a = netcdf_file.var('phi').get({'start' => [0, 0, 0, options[:ky_index] - 1], 'end' => [-1, -1, -1, options[:ky_index] - 1]})
temp = GSL::Vector.alloc(a.to_a[0][0].flatten)
when "box"
ep 'kx_elements', kx_elements = gsl_vector('linked_kx_elements', options).to_a
a = netcdf_file.var('phi').get({'start' => [0, 0, 0, options[:ky_index] - 1], 'end' => [-1, -1, -1, options[:ky_index] - 1]})
temp = GSL::Vector.alloc(a.to_a[0].values_at(*kx_elements).flatten)
else
raise "invalid grid option"
end
vector = GSL::Vector::Complex.alloc(temp.subvector_with_stride(0, 2), temp.subvector_with_stride(1, 2))
return vector
end
end
end
|