Read Collection of Tellier DATfiles.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "90e695c2",
   "metadata": {},
   "source": [
    "### Read cryo DAT files and plots X,Y,Z components"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "07671363",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "from decimal import Decimal\n",
    "from tkinter import *\n",
    "from tkinter import filedialog"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "20486b64",
   "metadata": {},
   "outputs": [],
   "source": [
    "'''\n",
    "read SAMPLE INFORMATION AND DEMAG DATA from binary DAT file from magnetometer (format of 2G-Barcelona)\n",
    "'''\n",
    "def readDATfile(filename):                    \n",
    "    with open(filename, 'r') as f:\n",
    "        data = f.read()\n",
    "        sample={}\n",
    "        sample[\"name\"]= data[15:24].rstrip('\\x00')                #sample name\n",
    "        #sample[\"vol\"]= \"{:5.2f}\".format(float(data[24:29].strip().strip('\\x00')))   #sample volume\n",
    "        sample[\"vol\"]= float(data[24:29].strip().strip('\\x00'))   #sample volume\n",
    "        sample[\"comment\"] = data[30:57].strip().strip('\\x00')\n",
    "        sample[\"CoreAz\"] = float(data[111:114].strip().strip('\\x00'))       # in- CoreAzimut\n",
    "        sample[\"CoreInc\"] = float(data[114:118].strip().strip('\\x00'))      #az - CoreInc\n",
    "        sample[\"BedAz\"] = float(data[118:123].strip('\\x00'))                #str - Bed dip directiom\n",
    "        sample[\"BedInc\"] = float(data[123:127].strip('\\x00'))               #dip - Bed inclination\n",
    "        sample[\"mag_Dec\"] = float(data[142:147].strip('\\x00'))\n",
    "        sample[\"com\"] = 50                                      #50 microTesla\n",
    "\n",
    "        sample_demag=[]\n",
    "        measurement = {}\n",
    "        for n in range(len(data)):\n",
    "            if data[n]==\"Í\":\n",
    "                line = n+1\n",
    "                step = data[line:line+4].strip().strip('\\x00')\n",
    "                treatment = data[line+4:line+5].strip().strip('\\x00')\n",
    "                X = float(data[line+65:line+76].strip().strip('\\x00'))/1E+3     # Am2  \n",
    "                Y = float(data[line+97:line+108].strip().strip('\\x00'))/1E+3    # Am2\n",
    "                Z = float(data[line+129:line+140].strip().strip('\\x00'))/1E+3   # Am2\n",
    "                Intensity = np.sqrt(X**2 + Y**2 + Z**2)\n",
    "                measurement ={\"step\":step,\"treat\":treatment,\n",
    "                          \"X\":X,\"Y\":Y,\"Z\":Z, \"Int (Am2)\":Intensity}\n",
    "                sample_demag.append(measurement)\n",
    "    return(sample,sample_demag)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "8defb126",
   "metadata": {},
   "outputs": [],
   "source": [
    "'''\n",
    "WRITES (ADD) SAMPLE INFORMATION TO TXT FILE\n",
    "'''\n",
    "def WriteTXT_heading_Tellier(sample,output):\n",
    "    with open(output, 'a') as f:        # a= to append (not overwritting)\n",
    "        heading_line = \"Id:{0}in:{1:05.1f} az:{2:04.1f}  dip:{3:04.1f}  str:{4:05.1f} m:{5:05.2f} com:{6}\".format(\n",
    "                                        sample[\"name\"].ljust(13),sample[\"CoreAz\"],sample[\"CoreInc\"],\n",
    "                                        sample[\"BedInc\"],sample[\"BedAz\"],sample[\"vol\"],sample[\"com\"])\n",
    "        print(heading_line)\n",
    "        print()\n",
    "        f.write(heading_line)\n",
    "        f.write('\\n')\n",
    "        f.write('\\n')    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "9b8b7a57",
   "metadata": {},
   "outputs": [],
   "source": [
    "'''\n",
    "WRITES DEMAG DATA TO TXT FILE\n",
    "'''    \n",
    "def WriteTXT_steps_Tellier(sample_demag, step,output):\n",
    "    global step_index\n",
    "    step_index = step\n",
    "    with open(output, 'a') as f:\n",
    "        for n in range(len(sample_demag)):\n",
    "            if sample_demag[n]['step'] == 'NRM':\n",
    "                step = \"0N0\"\n",
    "            if sample_demag[n]['step'][-1:] == '0':\n",
    "                step_index += 1\n",
    "                step = sample_demag[n]['step']+'R'+chr(64+step_index)  #chr(65) is 'A'\n",
    "            if sample_demag[n]['step'][-1:] == '1':\n",
    "                step = str(int(sample_demag[n]['step'])-1)+'V'+chr(64+step_index)\n",
    "            if sample_demag[n]['step'][-1:] == '2':\n",
    "                step = str(int(sample_demag[n]['step'])-2)+'P'+chr(64+step_index)\n",
    "            if sample_demag[n]['step'][-1:] == '3':\n",
    "                step = str(int(sample_demag[n]['step'])-3)+'X+'\n",
    "            if sample_demag[n]['step'][-1:] == '4':\n",
    "                step = str(int(sample_demag[n]['step'])-4)+'X-'\n",
    "            if sample_demag[n]['step'][-1:] == '5':\n",
    "                step = str(int(sample_demag[n]['step'])-5)+'Y+'\n",
    "            if sample_demag[n]['step'][-1:] == '6':\n",
    "                step = str(int(sample_demag[n]['step'])-6)+'Y-'\n",
    "            if sample_demag[n]['step'][-1:] == '7':\n",
    "                step = str(int(sample_demag[n]['step'])-7)+'Z+'\n",
    "\n",
    "            step_line = \"{0}{1:11}{2:11}{3:11}   1 Mo     0.\".format(step.rjust(7),\n",
    "                                                                     str(\"%.3E\" % Decimal(sample_demag[n]['X'])).rjust(11),\n",
    "                                                                     str(\"%.3E\" % Decimal(sample_demag[n]['Y'])).rjust(11),\n",
    "                                                                     str(\"%.3E\" % Decimal(sample_demag[n]['Z'])).rjust(11))\n",
    "            print(step_line)\n",
    "            f.write(step_line)\n",
    "            f.write('\\n')"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "91360062",
   "metadata": {},
   "source": [
    "### Selecting a collection of files"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "08ed4b11",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(('C:/Users/mgarc/OneDrive - Universitat de Barcelona/GitHub/paleomagnetism/arqueomag/dat_arqueomag/fins420/1b-3b.dat',\n",
       "  'C:/Users/mgarc/OneDrive - Universitat de Barcelona/GitHub/paleomagnetism/arqueomag/dat_arqueomag/fins420/2a-2a.dat',\n",
       "  'C:/Users/mgarc/OneDrive - Universitat de Barcelona/GitHub/paleomagnetism/arqueomag/dat_arqueomag/fins420/3a-2a.dat',\n",
       "  'C:/Users/mgarc/OneDrive - Universitat de Barcelona/GitHub/paleomagnetism/arqueomag/dat_arqueomag/fins420/4b-3b.dat'),\n",
       " ('C:/Users/mgarc/OneDrive - Universitat de Barcelona/GitHub/paleomagnetism/arqueomag/dat_arqueomag/1b-3b.dat',\n",
       "  'C:/Users/mgarc/OneDrive - Universitat de Barcelona/GitHub/paleomagnetism/arqueomag/dat_arqueomag/2a-2a.dat',\n",
       "  'C:/Users/mgarc/OneDrive - Universitat de Barcelona/GitHub/paleomagnetism/arqueomag/dat_arqueomag/3a-2a.dat',\n",
       "  'C:/Users/mgarc/OneDrive - Universitat de Barcelona/GitHub/paleomagnetism/arqueomag/dat_arqueomag/4b-3b.dat'))"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "root = Tk()\n",
    "# Hide the main window\n",
    "root.withdraw()\n",
    "root.call('wm', 'attributes', '.', '-topmost', True)\n",
    "root.title(\"select first set of files\")\n",
    "filenames_1part = filedialog.askopenfilenames(filetypes = (('binary DAT','typeName {DAT}')))\n",
    "filenames_1part\n",
    "filenames_2part = filedialog.askopenfilenames(filetypes = (('binary DAT','typeName {DAT}')))\n",
    "filenames_1part, filenames_2part"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "149dddeb",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Site name: prueba_\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "'prueba_Tellier.TXT'"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "output = str(input(\"Site name: \"))\n",
    "outputfile = output+\"Tellier.TXT\"\n",
    "outputfile"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "4dad04ef",
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Id:1b-3b        in:004.0 az:90.0  dip:00.0  str:000.0 m:08.00 com:50\n",
      "\n",
      "    0N0  1.800E-06 -1.043E-06  2.521E-06   1 Mo     0.\n",
      "  100RA  1.638E-06 -1.060E-06  2.683E-06   1 Mo     0.\n",
      "  100VA  1.609E-06 -1.108E-06  2.490E-06   1 Mo     0.\n",
      "  150RB  1.529E-06 -9.378E-07  2.742E-06   1 Mo     0.\n",
      "  150VB  1.558E-06 -9.909E-07  2.348E-06   1 Mo     0.\n",
      "  200RC  1.506E-06 -9.290E-07  2.702E-06   1 Mo     0.\n",
      "  200VC  1.496E-06 -9.233E-07  2.090E-06   1 Mo     0.\n",
      "  100PC  1.460E-06 -1.004E-06  2.241E-06   1 Mo     0.\n",
      "  250RD  1.369E-06 -9.172E-07  2.682E-06   1 Mo     0.\n",
      "  250VD  1.253E-06 -8.165E-07  1.851E-06   1 Mo     0.\n",
      "  300RE  1.173E-06 -8.620E-07  2.726E-06   1 Mo     0.\n",
      "  300VE  1.247E-06 -7.416E-07  1.265E-06   1 Mo     0.\n",
      "  200PE  1.200E-06 -8.415E-07  1.909E-06   1 Mo     0.\n",
      "  330RF  1.070E-06 -7.575E-07  2.741E-06   1 Mo     0.\n",
      "  330VF  1.084E-06 -6.635E-07  1.011E-06   1 Mo     0.\n",
      "  360RG  7.793E-07 -5.843E-07  2.848E-06   1 Mo     0.\n",
      "  360VG  1.026E-06 -5.976E-07  5.574E-07   1 Mo     0.\n",
      "  300PG  9.851E-07 -6.893E-07  2.049E-06   1 Mo     0.\n",
      "  390RH  8.573E-07 -6.806E-07  2.766E-06   1 Mo     0.\n",
      "  390VH  9.339E-07 -5.217E-07  2.928E-08   1 Mo     0.\n",
      "  420RI  6.416E-07 -5.540E-07  2.827E-06   1 Mo     0.\n",
      "  420VI  8.367E-07 -4.225E-07 -5.485E-07   1 Mo     0.\n",
      "  360PI  6.690E-07 -4.767E-07  1.745E-06   1 Mo     0.\n",
      "  420X+  2.609E-06 -3.614E-07  1.185E-06   1 Mo     0.\n",
      "  420X- -1.233E-06 -3.829E-07  8.408E-07   1 Mo     0.\n",
      "  420Y+  7.482E-07  1.530E-06  1.033E-06   1 Mo     0.\n",
      "  420Y-  5.396E-07 -2.472E-06  9.270E-07   1 Mo     0.\n",
      "  420Z+  6.228E-07 -5.456E-07  2.745E-06   1 Mo     0.\n",
      "  450RJ  4.458E-07 -3.784E-07  2.846E-06   1 Mo     0.\n",
      "  450VJ  6.588E-07 -2.225E-07 -1.198E-06   1 Mo     0.\n",
      "  480RK  2.508E-07 -2.341E-07  2.852E-06   1 Mo     0.\n",
      "  480VK  5.775E-07 -1.254E-07 -1.841E-06   1 Mo     0.\n",
      "  420PK  3.021E-07 -1.496E-07  1.621E-06   1 Mo     0.\n",
      "  510RL  1.032E-07 -1.192E-07  2.879E-06   1 Mo     0.\n",
      "  510VL  2.897E-07  3.686E-08 -2.413E-06   1 Mo     0.\n",
      "  540RM  8.187E-08 -6.851E-08  2.894E-06   1 Mo     0.\n",
      "  540VM  2.498E-07  1.555E-07 -2.672E-06   1 Mo     0.\n",
      "  480RN  3.283E-08 -1.687E-07  2.183E-06   1 Mo     0.\n",
      "  540PN  1.207E-07 -5.140E-08  2.877E-06   1 Mo     0.\n",
      "  540X+  2.669E-07  7.979E-08 -2.735E-06   1 Mo     0.\n",
      "Id:2a-2a        in:004.0 az:90.0  dip:00.0  str:000.0 m:08.00 com:50\n",
      "\n",
      "    0N0  7.553E-07 -4.930E-07  1.073E-06   1 Mo     0.\n",
      "  100RJ  7.005E-07 -4.610E-07  1.117E-06   1 Mo     0.\n",
      "  100VJ  6.806E-07 -4.508E-07  1.118E-06   1 Mo     0.\n",
      "  150RK  6.896E-07 -4.334E-07  1.134E-06   1 Mo     0.\n",
      "  150VK  6.827E-07 -4.481E-07  1.101E-06   1 Mo     0.\n",
      "  200RL  6.869E-07 -4.636E-07  1.129E-06   1 Mo     0.\n",
      "  200VL  6.941E-07 -4.586E-07  1.064E-06   1 Mo     0.\n",
      "  100PL  6.661E-07 -4.880E-07  1.081E-06   1 Mo     0.\n",
      "  250RM  6.591E-07 -4.802E-07  1.137E-06   1 Mo     0.\n",
      "  250VM  7.270E-07 -4.354E-07  9.803E-07   1 Mo     0.\n",
      "  300RN  6.548E-07 -4.379E-07  1.143E-06   1 Mo     0.\n",
      "  300VN  6.333E-07 -4.238E-07  9.017E-07   1 Mo     0.\n",
      "  200PN  6.172E-07 -4.681E-07  9.575E-07   1 Mo     0.\n",
      "  330RO  6.039E-07 -3.962E-07  1.170E-06   1 Mo     0.\n",
      "  330VO  5.773E-07 -4.177E-07  8.146E-07   1 Mo     0.\n",
      "  360RP  5.522E-07 -3.726E-07  1.190E-06   1 Mo     0.\n",
      "  360VP  5.735E-07 -3.724E-07  6.397E-07   1 Mo     0.\n",
      "  300PP  5.697E-07 -3.701E-07  8.964E-07   1 Mo     0.\n",
      "  390RQ  4.855E-07 -3.524E-07  1.214E-06   1 Mo     0.\n",
      "  390VQ  5.299E-07 -3.283E-07  4.014E-07   1 Mo     0.\n",
      "  420RR  4.118E-07 -2.600E-07  1.230E-06   1 Mo     0.\n",
      "  420VR  4.326E-07 -2.602E-07 -2.486E-09   1 Mo     0.\n",
      "  360PR  4.033E-07 -2.697E-07  6.220E-07   1 Mo     0.\n",
      "  450RS  2.770E-07 -2.163E-07  1.214E-06   1 Mo     0.\n",
      "  450VS  2.920E-07 -1.377E-07 -4.202E-07   1 Mo     0.\n",
      "  450X+  1.199E-06 -2.105E-07  3.897E-07   1 Mo     0.\n",
      "  450X- -6.513E-07 -1.382E-07  3.136E-07   1 Mo     0.\n",
      "  450Y+  3.294E-07  7.989E-07  3.658E-07   1 Mo     0.\n",
      "  450Y-  3.268E-07 -1.098E-06  3.397E-07   1 Mo     0.\n",
      "  450Z+  2.136E-07 -2.287E-07  1.193E-06   1 Mo     0.\n",
      "  480RT -6.662E-08  6.944E-08 -6.094E-07   1 Mo     0.\n",
      "  480VT -5.156E-07  9.945E-08 -7.165E-07   1 Mo     0.\n",
      "  480VT  1.765E-07 -2.795E-08 -8.058E-07   1 Mo     0.\n",
      "  420PT  1.487E-07 -7.702E-08  4.923E-07   1 Mo     0.\n",
      "  510RU  4.062E-08 -5.355E-08  1.231E-06   1 Mo     0.\n",
      "  510VU  1.150E-07 -1.635E-09 -1.092E-06   1 Mo     0.\n",
      "  540RV -8.469E-08 -4.469E-09  1.232E-06   1 Mo     0.\n",
      "  540VV  1.030E-07 -7.958E-09 -1.196E-06   1 Mo     0.\n",
      "  480RW  3.157E-08 -4.349E-08  8.578E-07   1 Mo     0.\n",
      "  540PW  1.645E-08 -8.727E-08  1.226E-06   1 Mo     0.\n",
      "  540X+  8.815E-08  1.191E-08 -1.198E-06   1 Mo     0.\n",
      "Id:3a-2a        in:157.0 az:90.0  dip:00.0  str:000.0 m:08.00 com:50\n",
      "\n",
      "    0N0 -1.853E-06 -1.781E-07  2.061E-06   1 Mo     0.\n",
      "  100RS -1.866E-06 -1.638E-07  2.073E-06   1 Mo     0.\n",
      "  100VS -1.839E-06 -1.642E-07  1.891E-06   1 Mo     0.\n",
      "  150RT -1.811E-06 -2.229E-07  2.005E-06   1 Mo     0.\n",
      "  150VT -1.771E-06 -2.669E-08  1.807E-06   1 Mo     0.\n",
      "  200RU -1.672E-06 -2.580E-07  2.023E-06   1 Mo     0.\n",
      "  200VU -1.640E-06 -1.398E-07  1.599E-06   1 Mo     0.\n",
      "  100PU -1.670E-06 -1.634E-07  1.682E-06   1 Mo     0.\n",
      "  250RV -1.577E-06 -1.709E-08  2.032E-06   1 Mo     0.\n",
      "  250VV -1.513E-06 -1.500E-07  1.306E-06   1 Mo     0.\n",
      "  300RW -1.420E-06 -1.289E-07  2.044E-06   1 Mo     0.\n",
      "  300VW -1.374E-06 -4.386E-08  1.021E-06   1 Mo     0.\n",
      "  200PW -1.513E-06 -2.150E-08  1.401E-06   1 Mo     0.\n",
      "  330RX -1.351E-06 -1.755E-07  2.021E-06   1 Mo     0.\n",
      "  330VX -1.274E-06 -9.656E-08  8.251E-07   1 Mo     0.\n",
      "  360RY -1.232E-06 -2.404E-07  2.039E-06   1 Mo     0.\n",
      "  360VY -1.205E-06 -7.408E-09  5.234E-07   1 Mo     0.\n",
      "  300PY -1.231E-06 -9.681E-08  1.651E-06   1 Mo     0.\n",
      "  390RZ -1.093E-06 -2.499E-07  2.053E-06   1 Mo     0.\n",
      "  390VZ -9.899E-07 -1.572E-08  3.495E-07   1 Mo     0.\n",
      "  420R[ -1.002E-06 -2.023E-07  2.082E-06   1 Mo     0.\n",
      "  420V[ -8.421E-07 -4.470E-09  5.239E-08   1 Mo     0.\n",
      "  360P[ -1.063E-06 -1.365E-07  1.486E-06   1 Mo     0.\n",
      "  450R\\ -8.280E-07 -2.407E-07  2.091E-06   1 Mo     0.\n",
      "  450V\\ -7.446E-07  6.831E-08 -2.818E-07   1 Mo     0.\n",
      "  450X+  6.591E-07 -1.354E-07  8.523E-07   1 Mo     0.\n",
      "  450X- -2.245E-06 -6.240E-08  8.300E-07   1 Mo     0.\n",
      "  450Y+ -6.708E-07  1.476E-06  9.293E-07   1 Mo     0.\n",
      "  450Y- -9.226E-07 -1.449E-06  7.204E-07   1 Mo     0.\n",
      "  450Z+ -8.172E-07 -1.751E-08  2.059E-06   1 Mo     0.\n",
      "  480R] -6.744E-07 -1.110E-07  2.101E-06   1 Mo     0.\n",
      "  480V] -5.448E-07  7.333E-08 -6.973E-07   1 Mo     0.\n",
      "  420P] -7.253E-07 -1.231E-07  1.358E-06   1 Mo     0.\n",
      "  510R^ -4.430E-07 -2.478E-08  2.120E-06   1 Mo     0.\n",
      "  510V^ -2.769E-07 -1.098E-07 -1.184E-06   1 Mo     0.\n",
      "  540R_ -2.348E-07 -1.013E-07  2.112E-06   1 Mo     0.\n",
      "  540V_  3.556E-09  5.948E-08 -1.619E-06   1 Mo     0.\n",
      "  480R` -1.574E-07 -3.679E-08  1.309E-06   1 Mo     0.\n",
      "  540P` -1.897E-07 -1.951E-07  2.129E-06   1 Mo     0.\n",
      "  540X+ -4.972E-08 -1.014E-06 -1.422E-06   1 Mo     0.\n",
      "Id:4b-3b        in:158.0 az:90.0  dip:00.0  str:000.0 m:08.00 com:50\n",
      "\n",
      "    0N0 -4.375E-06  4.276E-07  6.752E-06   1 Mo     0.\n",
      "  100R\\ -4.302E-06  4.440E-07  6.774E-06   1 Mo     0.\n",
      "  100V\\ -4.285E-06  3.734E-07  6.638E-06   1 Mo     0.\n",
      "  150R] -4.823E-06  4.472E-07  6.389E-06   1 Mo     0.\n",
      "  150V] -4.097E-06  1.496E-07  6.541E-06   1 Mo     0.\n",
      "  200R^ -4.126E-06  4.688E-08  6.828E-06   1 Mo     0.\n",
      "  200V^ -4.235E-06  1.284E-07  6.099E-06   1 Mo     0.\n",
      "  100P^ -4.383E-06  3.166E-07  6.155E-06   1 Mo     0.\n",
      "  250R_ -4.047E-06  1.576E-07  6.798E-06   1 Mo     0.\n",
      "  250V_ -3.918E-06 -7.397E-09  5.685E-06   1 Mo     0.\n",
      "  300R` -4.286E-06  2.838E-07  6.614E-06   1 Mo     0.\n",
      "  300V` -3.956E-06  2.369E-07  5.003E-06   1 Mo     0.\n",
      "  200P` -3.900E-06  2.694E-07  5.768E-06   1 Mo     0.\n",
      "  330Ra -3.873E-06  1.492E-07  6.809E-06   1 Mo     0.\n",
      "  330Va -3.651E-06  2.008E-07  4.806E-06   1 Mo     0.\n",
      "  360Rb -3.621E-06 -4.956E-08  6.923E-06   1 Mo     0.\n",
      "  360Vb -3.528E-06  3.835E-07  4.309E-06   1 Mo     0.\n",
      "  300Pb -3.666E-06  9.230E-08  6.077E-06   1 Mo     0.\n",
      "  390Rc -3.341E-06  2.752E-08  7.026E-06   1 Mo     0.\n",
      "  390Vc -3.224E-06  3.010E-07  3.731E-06   1 Mo     0.\n",
      "  420Rd -3.669E-06  3.890E-07  6.815E-06   1 Mo     0.\n",
      "  420Vd -2.891E-06  2.080E-07  2.680E-06   1 Mo     0.\n",
      "  360Pd -3.011E-06  4.225E-07  5.444E-06   1 Mo     0.\n",
      "  450Re -2.702E-06  1.074E-07  7.119E-06   1 Mo     0.\n",
      "  450Ve -2.532E-06  3.400E-07  1.206E-06   1 Mo     0.\n",
      "  480Rf -2.104E-06  3.042E-07  7.145E-06   1 Mo     0.\n",
      "  480Vf -1.753E-06  3.727E-07 -1.086E-06   1 Mo     0.\n",
      "  420Pf -1.877E-06  1.565E-07  3.810E-06   1 Mo     0.\n",
      "  480X+  2.297E-06 -1.760E-07  3.090E-06   1 Mo     0.\n",
      "  480X- -6.167E-06  4.287E-07  2.583E-06   1 Mo     0.\n",
      "  480Y+ -1.659E-06  4.251E-06  3.224E-06   1 Mo     0.\n",
      "  480Y- -2.079E-06 -4.049E-06  2.671E-06   1 Mo     0.\n",
      "  480Z+ -1.944E-06  3.345E-08  7.070E-06   1 Mo     0.\n",
      "  510Rg -1.410E-06 -1.676E-07  7.251E-06   1 Mo     0.\n",
      "  510Vg -9.518E-07  3.289E-07 -2.931E-06   1 Mo     0.\n",
      "  540Rh -9.598E-07 -1.129E-08  7.313E-06   1 Mo     0.\n",
      "  540Vh -3.216E-07  3.255E-07 -4.722E-06   1 Mo     0.\n",
      "  480Ri -7.234E-07  2.821E-08  3.752E-06   1 Mo     0.\n",
      "  560Rj -5.318E-07 -9.432E-08  7.316E-06   1 Mo     0.\n",
      "  560Vj  7.863E-08  3.048E-07 -5.678E-06   1 Mo     0.\n",
      "  570Rk -4.702E-07 -3.888E-07  7.288E-06   1 Mo     0.\n",
      "  570Vk  1.845E-07  2.054E-07 -6.101E-06   1 Mo     0.\n",
      "  540Pk  5.638E-11  2.251E-11 -3.594E-12   1 Mo     0.\n",
      "  570Pk -4.424E-07 -2.437E-07  7.260E-06   1 Mo     0.\n",
      "  570X+  5.429E-08 -1.282E-07 -6.058E-06   1 Mo     0.\n"
     ]
    }
   ],
   "source": [
    "global step_index\n",
    "step = 0\n",
    "\n",
    "for file in range(len(filenames_1part)):\n",
    "    sample,sample_demag=readDATfile(filenames_1part[file])\n",
    "    WriteTXT_heading_Tellier(sample, outputfile)\n",
    "    WriteTXT_steps_Tellier(sample_demag, step, outputfile)\n",
    "    if filenames_2part:\n",
    "        step = step_index\n",
    "        sample,sample_demag=readDATfile(filenames_2part[file])\n",
    "        WriteTXT_steps_Tellier(sample_demag, step, outputfile)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "b91eeecd",
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'\\n# Create Tk root\\nroot = Tk()\\n# Hide the main window\\nroot.withdraw()\\nroot.call(\\'wm\\', \\'attributes\\', \\'.\\', \\'-topmost\\', True)\\nfilename = filedialog.askopenfilename(filetypes = ((\\'binary DAT\\',\\'typeName {DAT}\\')))\\n \\n\\ntry:\\n    with open(filename, \"rb\") as f:\\n        byte = f.read(1)\\n        sum_bytes = np.array(byte)\\n        n=1\\n        while byte:\\n            # Do stuff with byte.\\n            byte = f.read(1)\\n            sum_bytes = np.append(sum_bytes, byte)\\n            print(n, byte)\\n            n += 1\\nexcept IOError:\\n     print(\\'Error While Opening the file!\\')\\n'"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#UNMARK TO PRINT THE LIST OF BYTES\n",
    "'''\n",
    "# Create Tk root\n",
    "root = Tk()\n",
    "# Hide the main window\n",
    "root.withdraw()\n",
    "root.call('wm', 'attributes', '.', '-topmost', True)\n",
    "filename = filedialog.askopenfilename(filetypes = (('binary DAT','typeName {DAT}')))\n",
    " \n",
    "\n",
    "try:\n",
    "    with open(filename, \"rb\") as f:\n",
    "        byte = f.read(1)\n",
    "        sum_bytes = np.array(byte)\n",
    "        n=1\n",
    "        while byte:\n",
    "            # Do stuff with byte.\n",
    "            byte = f.read(1)\n",
    "            sum_bytes = np.append(sum_bytes, byte)\n",
    "            print(n, byte)\n",
    "            n += 1\n",
    "except IOError:\n",
    "     print('Error While Opening the file!')\n",
    "'''"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}