LAPOD 2.0 apps: Media de direcciones paleomagnéticas medias

Calculadora para obtener una media vectorial a partir de dos o más direcciones paleomagnéticas medias.

1. Datos de entrada

Introduce las direcciones medias. Los campos Dec e Inc son obligatorios. Los campos R, n, k y α95 son opcionales.

Si no se proporciona R, el programa puede estimarlo usando esta prioridad:

1. Usar R si fue proporcionado.
2. Si no existe R, estimar R a partir de n y k.
3. Si no existe R ni k, usar R ≈ n.
4. Si solo existen Dec e Inc, usar peso uniforme R = 1.

Estimar automáticamente R cuando no se proporcione

Si esta opción está activada, el programa calculará R usando, en orden: R directo, n y k, R ≈ n, o R = 1. Si está desactivada, cada fila deberá tener un valor de R.

Mostrar stereonet de igual área

Grafica las direcciones introducidas, los círculos α95 cuando estén disponibles y la dirección media resultante.

#	Dec (°)	Inc (°)	R	n	k	α95 (°)	Eliminar

2. Entrada rápida por CSV

Puedes pegar datos en formato CSV. Se aceptan las columnas:

Dec,Inc,R,n,k,alpha95
12.3,45.6,8.75,10,,
18.1,47.2,,10,120,
5.4,44.8,,9,,
350.2,42.1,,,,

3. Método de cálculo

Cada dirección media se convierte a coordenadas cartesianas ponderadas por su longitud resultante:

X = R cos(I) cos(D)
Y = R cos(I) sin(D)
Z = R sin(I)

Posteriormente se suman las componentes:

Xsum = ΣX
Ysum = ΣY
Zsum = ΣZ

La dirección resultante se obtiene mediante:

Rtot = sqrt(Xsum² + Ysum² + Zsum²)
Dec_mean = atan2(Ysum, Xsum)
Inc_mean = atan2(Zsum, sqrt(Xsum² + Ysum²))

Si se dispone de información suficiente sobre n, también se estiman los parámetros Fisherianos:

N = Σn
k = (N - 1) / (N - Rtot)
cos(α95) = 1 - [(N - Rtot) / Rtot] × [0.05^(-1/(N-1)) - 1]

Para la gráfica se usa una proyección igual área tipo Schmidt. Las inclinaciones positivas se muestran como símbolos rellenos y las negativas como símbolos abiertos.