#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on 2025
@author: A. Petrenko
"""

#!/usr/bin/env python3
# -*- coding: utf-8 -*-

import numpy as np
import matplotlib.pyplot as plt

# Etude des mesures du transect 1a effectué le 01/01/2011 entre 04:18 et 23:38.

# Extration des données des coefficients d'absorption particulaire des fichiers texte :
data1ap = np.genfromtxt('/Users/apetrenko/Desktop/M2/OPB305_OptiqueMarine/TD_Absorption/OPB305/Tara_ACS_apcp2011_1ap.txt',skip_header=31)

# Extration des données d'incertitudes des absorptions particulaires des fichiers texte :
inc = np.genfromtxt('/Users/apetrenko/Desktop/M2/OPB305_OptiqueMarine/TD_Absorption/OPB305/Tara_ACS_apcp2011_1ap_uncertainty.txt',skip_header=31) 

# Création d'une liste des longueurs d'onde mesurées
LO = [405.5,410.1,414.4,419,423.7,428.6,433.1,437.7,442.5,447,452,456.6,461.1,465.5,470.6,475.7,480.6,485.5,489.9,494.4,499.1,503.8,508.5,513.5,518.4,523.2,527.7,532.3,537.2,541.5,545.9,550.7,555.3,559.8,563.2,567.6,572,576.3,580.5,584.9,589.3,593.9,598.4,603,607.6,612.7,617.4,622,626.5,630.9,635.4,640,644.8,649.3,654.2,658.6,663.4,668.2,672.6,676.8,680.9,685.5,689.7,693.7,697.6,701.7,705.4,709.4,712.8,716.4,720,723.3,726.8,730,733.4,736.7,739.6,742.4,745.6,748.1]

####### PLOT ap ################################################

#   ATTENTION    642 est le nombre de lignes de mesures; à adapter au fichier que vous avez chargé
for i in range(642):

# si les longueurs d'onde correspondent aux colonnes 6 à 86 alors (sinon adapter); verifier que cela correspond aux nombres de valeurs de LO
plt.plot(LO,data1ap[i,6:86],'k',linewidth=0.05)
plt.xlabel('Longueurs d\'onde [nm]')
plt.ylabel('ap [1/m]')