diff --git a/HumedadConNeutrones/codigos python/DetectorHits.py b/HumedadConNeutrones/codigos python/DetectorHits.py new file mode 100644 index 0000000000000000000000000000000000000000..f5ef192420ab79efc5ae0a26df382b1e560fda9b --- /dev/null +++ b/HumedadConNeutrones/codigos python/DetectorHits.py @@ -0,0 +1,42 @@ +# -*- coding: utf-8 -*- +""" +Spyder Editor + +This is a temporary script file. +""" +import numpy as np +import matplotlib.pyplot as plt +import pandas as pd + +#Dectector hits. +#Los datos corresponden a 10% de humedad en el suelo con porosidad del 50% + +dlayer= pd.read_table('/work/SALIDASdetectorNeutronHitData.dat') +datareal=dlayer[(dlayer['maximum_Depth_[m]']>0)] +data1=datareal.assign(hola=(datareal['maximum_Depth_[m]']*(-1))) + +#Graficamos vista aerea de el suelo + +datareal.plot(x='x_at_interface_[m]', y='y_at_Interface_[m]', kind="scatter", s = 4, fontsize=15, alpha=0.3) +plt.xlabel('x(m)', size=18) +plt.ylabel('y(m)', size=18) +plt.scatter(0,0, c='red') +plt.xlim(-300,300) +plt.ylim(-300,300) +plt.grid() +plt.show() +########### +#Graficamos vista transversal del suelo +######### +data1.plot(y='hola', x='x_at_interface_[m]', kind="scatter", s = 3, fontsize=15, alpha=0.3) +plt.xlim(-300,300) +plt.ylim(-1,0) +plt.scatter(0,0, c='red') +plt.xlabel('x(m)', size=18) +plt.ylabel('z(m)', size=18) +plt.grid() +plt.show() + + + +