From 01c00720eb0dee18a23849f5668be7b339e4819e Mon Sep 17 00:00:00 2001 From: Yessica Dominguez <yessicadomin@gmail.com> Date: Mon, 12 Dec 2022 00:45:44 +0000 Subject: [PATCH] Upload New File --- .../codigos python/DetectorHits.py | 42 +++++++++++++++++++ 1 file changed, 42 insertions(+) create mode 100644 HumedadConNeutrones/codigos python/DetectorHits.py diff --git a/HumedadConNeutrones/codigos python/DetectorHits.py b/HumedadConNeutrones/codigos python/DetectorHits.py new file mode 100644 index 0000000..f5ef192 --- /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() + + + + -- GitLab