@@ -556,8 +556,6 @@ def read_segy_trace_header(trace_index: int) -> bytes:
556556        flat_mask  =  trace_mask .ravel ()
557557        flat_raw_headers  =  raw_headers_data .ravel ()  # Flatten to 1D array of 240-byte header records 
558558
559-         # operation = "w" 
560- 
561559        for  grid_idx  in  range (flat_mask .size ):
562560            if  not  flat_mask [grid_idx ]:
563561                print (f"Skipping trace { grid_idx }   because it is masked" )
@@ -573,58 +571,4 @@ def read_segy_trace_header(trace_index: int) -> bytes:
573571
574572            assert_array_equal (mdio_header_bytes , segy_header_bytes )
575573
576-             # Compare byte-by-byte 
577-             # Write hexdumps to separate files for analysis 
578-             # def hexdump_to_string(data: bytes, title: str) -> str: 
579-             #     """Create hexdump string.""" 
580-             #     lines = [f"{title}", "=" * len(title), ""] 
581- 
582-             #     for i in range(0, len(data), 16): 
583-             #         # Address 
584-             #         addr = i 
585-             #         hex_part = "" 
586-             #         ascii_part = "" 
587- 
588-             #         # Process 16 bytes at a time 
589-             #         for j in range(16): 
590-             #             if i + j < len(data): 
591-             #                 byte_val = data[i + j] 
592-             #                 hex_part += f"{byte_val:02x} " 
593-             #                 ascii_part += chr(byte_val) if 32 <= byte_val <= 126 else "." 
594-             #             else: 
595-             #                 hex_part += "   " 
596-             #                 ascii_part += " " 
597- 
598-             #         lines.append(f"{addr:08x}: {hex_part} |{ascii_part}|") 
599- 
600-             #     return "\n".join(lines) 
601- 
602-             # # Generate filenames for this test case 
603-             # segy_filename = f"segy_headers_{grid_conf.name}.txt" 
604-             # mdio_filename = f"mdio_headers_{grid_conf.name}.txt" 
605- 
606-             # # Append SEG-Y hexdump to file 
607-             # with open(segy_filename, operation) as f: 
608-             #     if segy_trace_idx == 0: 
609-             #         f.write("")  # Start fresh for first trace 
610-             #     else: 
611-             #         f.write("\n\n")  # Add spacing between traces 
612-             #     f.write(hexdump_to_string(segy_header_bytes, 
613-             #                             f"SEG-Y Header - {grid_conf.name} Trace {segy_trace_idx} (240 bytes)")) 
614- 
615-             # # Append MDIO hexdump to file 
616-             # with open(mdio_filename, operation) as f: 
617-             #     if segy_trace_idx == 0: 
618-             #         f.write("")  # Start fresh for first trace 
619-             #     else: 
620-             #         f.write("\n\n")  # Add spacing between traces 
621-             #     f.write(hexdump_to_string(mdio_header_bytes, 
622-             #                             f"MDIO Raw Header - {grid_conf.name} Trace {segy_trace_idx} (240 bytes)")) 
623-             # operation = 'a' 
624- 
625-             # if segy_trace_idx == 0: 
626-             #     print(f"\nHeader hexdumps being written for {grid_conf.name}:") 
627-             #     print(f"  SEG-Y: {segy_filename}") 
628-             #     print(f"  MDIO:  {mdio_filename}") 
629- 
630574            segy_trace_idx  +=  1 
0 commit comments