Point function not working for negative coordinates
Hello yt users, I have been trying to obtain density values at different points using point data object. The problem is that this function does not work for negative coordinates regardless the fact that the computational domain contains those coordinates. This is the error I got: Entering folder ./lr_teq_p_3_nk4 Using rho_0: 6.372606e-24 Using 5.0 snapshots yt : [INFO ] 2018-09-12 17:34:34,371 Parameters: current_time = 0.0 yt : [INFO ] 2018-09-12 17:34:34,371 Parameters: domain_dimensions = [256 256 256] yt : [INFO ] 2018-09-12 17:34:34,372 Parameters: domain_left_edge = [ -1.54283920e+19 -1.54283920e+19 -1.54283920e+19] yt : [INFO ] 2018-09-12 17:34:34,372 Parameters: domain_right_edge = [ 1.54283920e+19 1.54283920e+19 1.54283920e+19] yt : [INFO ] 2018-09-12 17:34:34,372 Parameters: cosmological_simulation = 0.0 domain: [ 10.00000265 10.00000265 10.00000265] pc YTPoint (lr_hdf5_plt_cnt_0000): , p=[ 1. 1. -0.001] cm [] g/cm**3 --------------------------------------------------------------------------- ValueError Traceback (most recent call last) /share/Part1/guido/w_profile2.py in <module>() 84 print(point) 85 print(point['gas','density']) ---> 86 dens_array[snap_order,ipoint]=point['gas','density'] 87 ipoint=ipoint+1 88 print("time :",time) ValueError: setting an array element with a sequence. As you can see, the point p=[ 1. 1. -0.001] cm belongs to the computational domain. When I only use positive coordinates the code runs well. Correct me if I am wrong, but I think this is a bug. Do you guys know another way to get values of any field at a fixed point? Cheers, -- Guido
participants (2)
-
Guido granda muñoz
-
Matthew Turk