/*-------------------------------------------------------------------- * The MB-system: mbr_mgd77txt.c 5/18/99 * $Id: mbr_mgd77txt.c 2227 2015-02-05 21:12:31Z caress $ * * Copyright (c) 2014-2015 by * David W. Caress (caress@mbari.org) * Monterey Bay Aquarium Research Institute * Moss Landing, CA 95039 * and Dale N. Chayes (dale@ldeo.columbia.edu) * Lamont-Doherty Earth Observatory * Palisades, NY 10964 * * See README file for copying and redistribution conditions. *--------------------------------------------------------------------*/ /* * mbr_mgd77txt.c contains the functions for reading and writing * multibeam data in the MGD77TXT format. * These functions include: * mbr_alm_mgd77txt - allocate read/write memory * mbr_dem_mgd77txt - deallocate read/write memory * mbr_rt_mgd77txt - read and translate data * mbr_wt_mgd77txt - translate and write data * * Author: D. W. Caress * Date: February 25, 2014 * * $Log: mbr_mgd77txt.c,v $ * */ /* standard include files */ #include #include #include /* mbio include files */ #include "mb_status.h" #include "mb_format.h" #include "mb_io.h" #include "mb_define.h" #include "mbsys_singlebeam.h" /* * Notes on the MBF_MGD77TXT data format: * 1. The MGD77 format is is an exchange format for marine * geophysical data (bathymetry, magnetics, and gravity). * The format standard is maintained by the National * Geophysical Data Center of NOAA. * 2. The standard MGD77 format includes a 1920 byte header * followed by 120 byte data records. The header consists * of 24 80 byte records. The first character of the first * header record is either 1 (pre-Y2K fix) or 4 (post-Y2K fix). * MB-System treats the header as 16 120 byte records and * provides no means of modifying the header. * 3. The data records are each 120 bytes long. The first * character of each data record is either * 3 (pre-Y2K fix) or 5 (post-Y2K fix). * 4. The MB-System implementation includes the support of * an arbitrary number of comment records at the beginning * of each file. The comment records are 120 bytes each * and begin with the character '#'. * */ /* header and data record in bytes */ #define MBF_MGD77TXT_HEADER_NUM 16 #define MBF_MGD77TXT_DATA_LEN 128 struct mbf_mgd77txt_struct { /* type of data record */ int kind; /* survey id */ char survey_id[8]; /* Identifier supplied by the contributing organization, else given by NGDC in a manner which represents the data. */ /* time stamp */ double time_d; int time_i[7]; int timezone; /* Corrects time (in characters 13-27) to GMT when added: equals zero when time is GMT. Timezone normally falls between -13 and +12 inclusively. */ /* navigation */ double longitude; double latitude; double heading; /* degrees */ double speed; /* km/hr */ int nav_type; /* Indicates how lat/lon was obtained: 1 = Observed fix 3 = Interpolated 9 = Unspecified */ int nav_quality; /* QUALITY CODE FOR NAVIGATION - 5 - Suspected, by the originating institution 6 - Suspected, by the data center 9 - No identifiable problem found */ /* motion sensor data */ double roll; double pitch; double heave; /* bathymetry */ int flag; /* MB-System style beamflag */ double tt; /* two way travel time in sec */ double bath; /* corrected depth in m */ int bath_corr; /* BATHYMETRIC CORRECTION CODE This code details the procedure used for determining the sound velocity correction to depth: 01-55 Matthews' Zones with zone 59 Matthews' Zones, no zone 60 S. Kuwahara Formula 61 Wilson Formula 62 Del Grosso Formula 63 Carter's Tables 88 Other (see Add. Doc.) 99 Unspecified */ int bath_type; /* BATHYMETRIC TYPE CODE Indicates how the data record's bathymetric value was obtained: 1 = Observed 3 = Interpolated 9 = Unspecified */ /* magnetics */ double mag_tot_1; /* MAGNETICS TOTAL FIELD, 1ST SENSOR In tenths of nanoteslas (gammas). For leading sensor. Use this field for single sensor. */ double mag_tot_2; /* MAGNETICS TOTAL FIELD, 2ND SENSOR In tenths of nanoteslas (gammas). For trailing sensor. */ double mag_res; /* MAGNETICS RESIDUAL FIELD In tenths of nanoteslas (gammas). */ int mag_res_sensor; /* SENSOR FOR RESIDUAL FIELD 1 = 1st or leading sensor 2 = 2nd or trailing sensor 9 = Unspecified */ double mag_diurnal; /* MAGNETICS DIURNAL CORRECTION - In tenths of nanoteslas (gammas). (In nanoteslas) if 9-filled (i.e., set to "+9999"), total and residual fields are assumed to be uncorrected; if used, total and residuals are assumed to have been already corrected. */ double mag_altitude; /* DEPTH OR ALTITUDE OF MAGNETICS SENSOR In meters. + = Below sealevel - = Above sealevel */ /* gravity */ double gravity; /* OBSERVED GRAVITY In milligals. Corrected for Eotvos, drift, and tares */ double eotvos; /* EOTVOS CORRECTION In milligals. E = 7.5 V cos phi sin alpha + 0.0042 V*V */ double free_air; /* FREE-AIR ANOMALY In milligals. Free-air Anomaly = G(observed) - G(theoretical) */ /* seismic */ int seismic_line; /* SEISMIC LINE NUMBER Used for cross referencing with seismic data. */ int seismic_shot; /* SEISMIC SHOT-POINT NUMBER */ /* comment */ char comment[MB_COMMENT_MAXLINE]; }; /* essential function prototypes */ int mbr_register_mgd77txt(int verbose, void *mbio_ptr, int *error); int mbr_info_mgd77txt(int verbose, int *system, int *beams_bath_max, int *beams_amp_max, int *pixels_ss_max, char *format_name, char *system_name, char *format_description, int *numfile, int *filetype, int *variable_beams, int *traveltime, int *beam_flagging, int *nav_source, int *heading_source, int *vru_source, int *svp_source, double *beamwidth_xtrack, double *beamwidth_ltrack, int *error); int mbr_alm_mgd77txt(int verbose, void *mbio_ptr, int *error); int mbr_dem_mgd77txt(int verbose, void *mbio_ptr, int *error); int mbr_zero_mgd77txt(int verbose, char *data_ptr, int *error); int mbr_rt_mgd77txt(int verbose, void *mbio_ptr, void *store_ptr, int *error); int mbr_wt_mgd77txt(int verbose, void *mbio_ptr, void *store_ptr, int *error); int mbr_mgd77txt_rd_data(int verbose, void *mbio_ptr, int *error); int mbr_mgd77txt_wr_data(int verbose, void *mbio_ptr, void *data_ptr, int *error); static char rcs_id[]="$Id: mbr_mgd77txt.c 2227 2015-02-05 21:12:31Z caress $"; /*--------------------------------------------------------------------*/ int mbr_register_mgd77txt(int verbose, void *mbio_ptr, int *error) { char *function_name = "mbr_register_mgd77txt"; int status = MB_SUCCESS; struct mb_io_struct *mb_io_ptr; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n",function_name); fprintf(stderr,"dbg2 Revision id: %s\n",rcs_id); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 verbose: %d\n",verbose); } /* get mb_io_ptr */ mb_io_ptr = (struct mb_io_struct *) mbio_ptr; /* set format info parameters */ status = mbr_info_mgd77txt(verbose, &mb_io_ptr->system, &mb_io_ptr->beams_bath_max, &mb_io_ptr->beams_amp_max, &mb_io_ptr->pixels_ss_max, mb_io_ptr->format_name, mb_io_ptr->system_name, mb_io_ptr->format_description, &mb_io_ptr->numfile, &mb_io_ptr->filetype, &mb_io_ptr->variable_beams, &mb_io_ptr->traveltime, &mb_io_ptr->beam_flagging, &mb_io_ptr->nav_source, &mb_io_ptr->heading_source, &mb_io_ptr->vru_source, &mb_io_ptr->svp_source, &mb_io_ptr->beamwidth_xtrack, &mb_io_ptr->beamwidth_ltrack, error); /* set format and system specific function pointers */ mb_io_ptr->mb_io_format_alloc = &mbr_alm_mgd77txt; mb_io_ptr->mb_io_format_free = &mbr_dem_mgd77txt; mb_io_ptr->mb_io_store_alloc = &mbsys_singlebeam_alloc; mb_io_ptr->mb_io_store_free = &mbsys_singlebeam_deall; mb_io_ptr->mb_io_read_ping = &mbr_rt_mgd77txt; mb_io_ptr->mb_io_write_ping = &mbr_wt_mgd77txt; mb_io_ptr->mb_io_dimensions = &mbsys_singlebeam_dimensions; mb_io_ptr->mb_io_extract = &mbsys_singlebeam_extract; mb_io_ptr->mb_io_insert = &mbsys_singlebeam_insert; mb_io_ptr->mb_io_extract_nav = &mbsys_singlebeam_extract_nav; mb_io_ptr->mb_io_insert_nav = &mbsys_singlebeam_insert_nav; mb_io_ptr->mb_io_extract_altitude = &mbsys_singlebeam_extract_altitude; mb_io_ptr->mb_io_insert_altitude = NULL; mb_io_ptr->mb_io_extract_svp = NULL; mb_io_ptr->mb_io_insert_svp = NULL; mb_io_ptr->mb_io_ttimes = &mbsys_singlebeam_ttimes; mb_io_ptr->mb_io_detects = &mbsys_singlebeam_detects; mb_io_ptr->mb_io_copyrecord = &mbsys_singlebeam_copy; mb_io_ptr->mb_io_extract_rawss = NULL; mb_io_ptr->mb_io_insert_rawss = NULL; /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n",function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 system: %d\n",mb_io_ptr->system); fprintf(stderr,"dbg2 beams_bath_max: %d\n",mb_io_ptr->beams_bath_max); fprintf(stderr,"dbg2 beams_amp_max: %d\n",mb_io_ptr->beams_amp_max); fprintf(stderr,"dbg2 pixels_ss_max: %d\n",mb_io_ptr->pixels_ss_max); fprintf(stderr,"dbg2 format_name: %s\n",mb_io_ptr->format_name); fprintf(stderr,"dbg2 system_name: %s\n",mb_io_ptr->system_name); fprintf(stderr,"dbg2 format_description: %s\n",mb_io_ptr->format_description); fprintf(stderr,"dbg2 numfile: %d\n",mb_io_ptr->numfile); fprintf(stderr,"dbg2 filetype: %d\n",mb_io_ptr->filetype); fprintf(stderr,"dbg2 variable_beams: %d\n",mb_io_ptr->variable_beams); fprintf(stderr,"dbg2 traveltime: %d\n",mb_io_ptr->traveltime); fprintf(stderr,"dbg2 beam_flagging: %d\n",mb_io_ptr->beam_flagging); fprintf(stderr,"dbg2 nav_source: %d\n",mb_io_ptr->nav_source); fprintf(stderr,"dbg2 heading_source: %d\n",mb_io_ptr->heading_source); fprintf(stderr,"dbg2 vru_source: %d\n",mb_io_ptr->vru_source); fprintf(stderr,"dbg2 svp_source: %d\n",mb_io_ptr->svp_source); fprintf(stderr,"dbg2 beamwidth_xtrack: %f\n",mb_io_ptr->beamwidth_xtrack); fprintf(stderr,"dbg2 beamwidth_ltrack: %f\n",mb_io_ptr->beamwidth_ltrack); fprintf(stderr,"dbg2 format_alloc: %p\n",(void *)mb_io_ptr->mb_io_format_alloc); fprintf(stderr,"dbg2 format_free: %p\n",(void *)mb_io_ptr->mb_io_format_free); fprintf(stderr,"dbg2 store_alloc: %p\n",(void *)mb_io_ptr->mb_io_store_alloc); fprintf(stderr,"dbg2 store_free: %p\n",(void *)mb_io_ptr->mb_io_store_free); fprintf(stderr,"dbg2 read_ping: %p\n",(void *)mb_io_ptr->mb_io_read_ping); fprintf(stderr,"dbg2 write_ping: %p\n",(void *)mb_io_ptr->mb_io_write_ping); fprintf(stderr,"dbg2 extract: %p\n",(void *)mb_io_ptr->mb_io_extract); fprintf(stderr,"dbg2 insert: %p\n",(void *)mb_io_ptr->mb_io_insert); fprintf(stderr,"dbg2 extract_nav: %p\n",(void *)mb_io_ptr->mb_io_extract_nav); fprintf(stderr,"dbg2 insert_nav: %p\n",(void *)mb_io_ptr->mb_io_insert_nav); fprintf(stderr,"dbg2 extract_altitude: %p\n",(void *)mb_io_ptr->mb_io_extract_altitude); fprintf(stderr,"dbg2 insert_altitude: %p\n",(void *)mb_io_ptr->mb_io_insert_altitude); fprintf(stderr,"dbg2 extract_svp: %p\n",(void *)mb_io_ptr->mb_io_extract_svp); fprintf(stderr,"dbg2 insert_svp: %p\n",(void *)mb_io_ptr->mb_io_insert_svp); fprintf(stderr,"dbg2 ttimes: %p\n",(void *)mb_io_ptr->mb_io_ttimes); fprintf(stderr,"dbg2 detects: %p\n",(void *)mb_io_ptr->mb_io_detects); fprintf(stderr,"dbg2 extract_rawss: %p\n",(void *)mb_io_ptr->mb_io_extract_rawss); fprintf(stderr,"dbg2 insert_rawss: %p\n",(void *)mb_io_ptr->mb_io_insert_rawss); fprintf(stderr,"dbg2 copyrecord: %p\n",(void *)mb_io_ptr->mb_io_copyrecord); fprintf(stderr,"dbg2 error: %d\n",*error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return status */ return(status); } /*--------------------------------------------------------------------*/ int mbr_info_mgd77txt(int verbose, int *system, int *beams_bath_max, int *beams_amp_max, int *pixels_ss_max, char *format_name, char *system_name, char *format_description, int *numfile, int *filetype, int *variable_beams, int *traveltime, int *beam_flagging, int *nav_source, int *heading_source, int *vru_source, int *svp_source, double *beamwidth_xtrack, double *beamwidth_ltrack, int *error) { char *function_name = "mbr_info_mgd77txt"; int status = MB_SUCCESS; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n",function_name); fprintf(stderr,"dbg2 Revision id: %s\n",rcs_id); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 verbose: %d\n",verbose); } /* set format info parameters */ status = MB_SUCCESS; *error = MB_ERROR_NO_ERROR; *system = MB_SYS_SINGLEBEAM; *beams_bath_max = 1; *beams_amp_max = 0; *pixels_ss_max = 0; strncpy(format_name, "MGD77TXT", MB_NAME_LENGTH); strncpy(system_name, "SINGLEBEAM", MB_NAME_LENGTH); strncpy(format_description, "Format name: MBF_MGD77TXT\nInformal Description: NGDC MGD77 underway geophysics format\nAttributes: single beam bathymetry, nav, magnetics, gravity,\n 122 byte ascii records with CRLF line breaks, NOAA NGDC\n", MB_DESCRIPTION_LENGTH); *numfile = 1; *filetype = MB_FILETYPE_NORMAL; *variable_beams = MB_NO; *traveltime = MB_YES; *beam_flagging = MB_NO; *nav_source = MB_DATA_DATA; *heading_source = MB_DATA_DATA; *vru_source = MB_DATA_DATA; *svp_source = MB_DATA_NONE; *beamwidth_xtrack = 0.0; *beamwidth_ltrack = 0.0; /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n",function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 system: %d\n",*system); fprintf(stderr,"dbg2 beams_bath_max: %d\n",*beams_bath_max); fprintf(stderr,"dbg2 beams_amp_max: %d\n",*beams_amp_max); fprintf(stderr,"dbg2 pixels_ss_max: %d\n",*pixels_ss_max); fprintf(stderr,"dbg2 format_name: %s\n",format_name); fprintf(stderr,"dbg2 system_name: %s\n",system_name); fprintf(stderr,"dbg2 format_description: %s\n",format_description); fprintf(stderr,"dbg2 numfile: %d\n",*numfile); fprintf(stderr,"dbg2 filetype: %d\n",*filetype); fprintf(stderr,"dbg2 variable_beams: %d\n",*variable_beams); fprintf(stderr,"dbg2 traveltime: %d\n",*traveltime); fprintf(stderr,"dbg2 beam_flagging: %d\n",*beam_flagging); fprintf(stderr,"dbg2 nav_source: %d\n",*nav_source); fprintf(stderr,"dbg2 heading_source: %d\n",*heading_source); fprintf(stderr,"dbg2 vru_source: %d\n",*vru_source); fprintf(stderr,"dbg2 svp_source: %d\n",*svp_source); fprintf(stderr,"dbg2 beamwidth_xtrack: %f\n",*beamwidth_xtrack); fprintf(stderr,"dbg2 beamwidth_ltrack: %f\n",*beamwidth_ltrack); fprintf(stderr,"dbg2 error: %d\n",*error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return status */ return(status); } /*--------------------------------------------------------------------*/ int mbr_alm_mgd77txt(int verbose, void *mbio_ptr, int *error) { char *function_name = "mbr_alm_mgd77txt"; int status = MB_SUCCESS; struct mb_io_struct *mb_io_ptr; struct mbf_mgd77txt_struct *data; char *data_ptr; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n",function_name); fprintf(stderr,"dbg2 Revision id: %s\n",rcs_id); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 verbose: %d\n",verbose); fprintf(stderr,"dbg2 mbio_ptr: %p\n",(void *)mbio_ptr); } /* get pointer to mbio descriptor */ mb_io_ptr = (struct mb_io_struct *) mbio_ptr; /* set initial status */ status = MB_SUCCESS; /* allocate memory for data structure */ mb_io_ptr->structure_size = sizeof(struct mbf_mgd77txt_struct); mb_io_ptr->data_structure_size = 0; status = mb_malloc(verbose,mb_io_ptr->structure_size, &mb_io_ptr->raw_data,error); status = mb_malloc(verbose,sizeof(struct mbsys_singlebeam_struct), &mb_io_ptr->store_data,error); /* get pointer to mbio descriptor */ mb_io_ptr = (struct mb_io_struct *) mbio_ptr; data = (struct mbf_mgd77txt_struct *) mb_io_ptr->raw_data; data_ptr = (char *) data; /* set number of header records read to zero */ mb_io_ptr->save1 = 0; /* initialize everything to zeros */ mbr_zero_mgd77txt(verbose,data_ptr,error); /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n",function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",*error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return status */ return(status); } /*--------------------------------------------------------------------*/ int mbr_dem_mgd77txt(int verbose, void *mbio_ptr, int *error) { char *function_name = "mbr_dem_mgd77txt"; int status = MB_SUCCESS; struct mb_io_struct *mb_io_ptr; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n",function_name); fprintf(stderr,"dbg2 Revision id: %s\n",rcs_id); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 verbose: %d\n",verbose); fprintf(stderr,"dbg2 mbio_ptr: %p\n",(void *)mbio_ptr); } /* get pointer to mbio descriptor */ mb_io_ptr = (struct mb_io_struct *) mbio_ptr; /* deallocate memory for data descriptor */ status = mb_free(verbose,&mb_io_ptr->raw_data,error); status = mb_free(verbose,&mb_io_ptr->store_data,error); /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n",function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",*error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return status */ return(status); } /*--------------------------------------------------------------------*/ int mbr_zero_mgd77txt(int verbose, char *data_ptr, int *error) { char *function_name = "mbr_zero_mgd77txt"; int status = MB_SUCCESS; struct mbf_mgd77txt_struct *data; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n",function_name); fprintf(stderr,"dbg2 Revision id: %s\n",rcs_id); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 verbose: %d\n",verbose); fprintf(stderr,"dbg2 data_ptr: %p\n",(void *)data_ptr); } /* get pointer to data descriptor */ data = (struct mbf_mgd77txt_struct *) data_ptr; /* initialize everything to zeros */ if (data != NULL) { data->kind = MB_DATA_NONE; for (i=0;i<8;i++) data->survey_id[i] = 0; data->time_d = 0.0; for (i=0;i<7;i++) data->time_i[i] = 0; data->timezone = 0; data->longitude = 0.0; data->latitude = 0.0; data->heading = 0.0; data->speed = 0.0; data->nav_type = 9; data->nav_quality = 9; data->roll = 0.0; data->pitch = 0.0; data->heave = 0.0; data->tt = 0.0; data->flag = MB_FLAG_NULL; data->bath = 0.0; data->bath_corr = 99; data->bath_type = 9; data->mag_tot_1 = 0.0; data->mag_tot_2 = 0.0; data->mag_res = 0.0; data->mag_res_sensor = 9; data->mag_diurnal = 0.0; data->mag_altitude = 0.0; data->gravity = 0.0; data->eotvos = 0.0; data->free_air = 0.0; data->seismic_line = 0; data->seismic_shot = 0; for (i=0;icomment[i] = 0; } /* assume success */ status = MB_SUCCESS; *error = MB_ERROR_NO_ERROR; /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n",function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",*error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return status */ return(status); } /*--------------------------------------------------------------------*/ int mbr_rt_mgd77txt(int verbose, void *mbio_ptr, void *store_ptr, int *error) { char *function_name = "mbr_rt_mgd77txt"; int status = MB_SUCCESS; struct mb_io_struct *mb_io_ptr; struct mbf_mgd77txt_struct *data; struct mbsys_singlebeam_struct *store; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n",function_name); fprintf(stderr,"dbg2 Revision id: %s\n",rcs_id); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 verbose: %d\n",verbose); fprintf(stderr,"dbg2 mbio_ptr: %p\n",(void *)mbio_ptr); fprintf(stderr,"dbg2 store_ptr: %p\n",(void *)store_ptr); } /* get pointers to mbio descriptor and data structures */ mb_io_ptr = (struct mb_io_struct *) mbio_ptr; data = (struct mbf_mgd77txt_struct *) mb_io_ptr->raw_data; store = (struct mbsys_singlebeam_struct *) store_ptr; /* read next data from file */ status = mbr_mgd77txt_rd_data(verbose,mbio_ptr,error); /* set error and kind in mb_io_ptr */ mb_io_ptr->new_error = *error; mb_io_ptr->new_kind = data->kind; /* translate values to data storage structure */ if (status == MB_SUCCESS && store != NULL) { store->kind = data->kind; for (i=0;i<8;i++) store->survey_id[i] = data->survey_id[i]; store->time_d = data->time_d; for (i=0;i<7;i++) store->time_i[i] = data->time_i[i]; store->timezone = data->timezone; store->longitude = data->longitude; store->latitude = data->latitude; store->heading = data->heading; store->speed = data->speed; store->nav_type = data->nav_type; store->nav_quality = data->nav_quality; store->roll = data->roll; store->pitch = data->pitch; store->heave = data->heave; store->flag = data->flag; store->tt = data->tt; store->bath = data->bath; store->bath_corr = data->bath_corr; store->bath_type = data->bath_type; store->mag_tot_1 = data->mag_tot_1; store->mag_tot_2 = data->mag_tot_2; store->mag_res = data->mag_res; store->mag_res_sensor = data->mag_res_sensor; store->mag_diurnal = data->mag_diurnal; store->mag_altitude = data->mag_altitude; store->gravity = data->gravity; store->eotvos = data->eotvos; store->free_air = data->free_air; store->seismic_line = data->seismic_line; store->seismic_shot = data->seismic_shot; for (i=0;icomment[i] = data->comment[i]; } /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n",function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",*error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return status */ return(status); } /*--------------------------------------------------------------------*/ int mbr_wt_mgd77txt(int verbose, void *mbio_ptr, void *store_ptr, int *error) { char *function_name = "mbr_wt_mgd77txt"; int status = MB_SUCCESS; struct mb_io_struct *mb_io_ptr; struct mbf_mgd77txt_struct *data; struct mbsys_singlebeam_struct *store; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n",function_name); fprintf(stderr,"dbg2 Revision id: %s\n",rcs_id); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 verbose: %d\n",verbose); fprintf(stderr,"dbg2 mbio_ptr: %p\n",(void *)mbio_ptr); fprintf(stderr,"dbg2 store_ptr: %p\n",(void *)store_ptr); } /* get pointer to mbio descriptor */ mb_io_ptr = (struct mb_io_struct *) mbio_ptr; /* get pointer to raw data structure */ data = (struct mbf_mgd77txt_struct *) mb_io_ptr->raw_data; store = (struct mbsys_singlebeam_struct *) store_ptr; /* first translate values from data storage structure */ if (store != NULL) { data->kind = store->kind; for (i=0;i<8;i++) data->survey_id[i] = store->survey_id[i]; data->time_d = store->time_d; for (i=0;i<7;i++) data->time_i[i] = store->time_i[i]; data->timezone = store->timezone; data->longitude = store->longitude; data->latitude = store->latitude; data->heading = store->heading; data->speed = store->speed; data->nav_type = store->nav_type; data->nav_quality = store->nav_quality; data->roll = store->roll; data->pitch = store->pitch; data->heave = store->heave; data->flag = store->flag; data->tt = store->tt; data->bath = store->bath; data->bath_corr = store->bath_corr; data->bath_type = store->bath_type; data->mag_tot_1 = store->mag_tot_1; data->mag_tot_2 = store->mag_tot_2; data->mag_res = store->mag_res; data->mag_res_sensor = store->mag_res_sensor; data->mag_diurnal = store->mag_diurnal; data->mag_altitude = store->mag_altitude; data->gravity = store->gravity; data->eotvos = store->eotvos; data->free_air = store->free_air; data->seismic_line = store->seismic_line; data->seismic_shot = store->seismic_shot; for (i=0;icomment[i] = store->comment[i]; } /* write next data to file */ status = mbr_mgd77txt_wr_data(verbose,mbio_ptr,(void *)data,error); /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n",function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",*error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return status */ return(status); } /*--------------------------------------------------------------------*/ int mbr_mgd77txt_rd_data(int verbose, void *mbio_ptr, int *error) { char *function_name = "mbr_mgd77txt_rd_data"; int status = MB_SUCCESS; struct mb_io_struct *mb_io_ptr; struct mbf_mgd77txt_struct *data; int *header_read; char line[MB_COMMENT_MAXLINE]; char *read_ptr; int shift; int neg_unit; int itmp; double dtmp; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n",function_name); fprintf(stderr,"dbg2 Revision id: %s\n",rcs_id); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 verbose: %d\n",verbose); fprintf(stderr,"dbg2 mbio_ptr: %p\n",(void *)mbio_ptr); } /* get pointer to mbio descriptor */ mb_io_ptr = (struct mb_io_struct *) mbio_ptr; /* get pointer to raw data structure */ data = (struct mbf_mgd77txt_struct *) mb_io_ptr->raw_data; header_read = (int *) &mb_io_ptr->save1; /* initialize everything to zeros */ mbr_zero_mgd77txt(verbose,mb_io_ptr->raw_data,error); /* set file position */ mb_io_ptr->file_bytes = ftell(mb_io_ptr->mbfp); mb_io_ptr->file_pos = mb_io_ptr->file_bytes; /* read next record */ if ((read_ptr = fgets(line, MB_PATH_MAXLINE, mb_io_ptr->mbfp)) != NULL) { mb_io_ptr->file_bytes += strlen(line); status = MB_SUCCESS; *error = MB_ERROR_NO_ERROR; } else { status = MB_FAILURE; *error = MB_ERROR_EOF; } mb_io_ptr->file_bytes = ftell(mb_io_ptr->mbfp); /* handle the data */ if (status == MB_SUCCESS && *header_read > 0 && *header_read < MBF_MGD77TXT_HEADER_NUM) { data->kind = MB_DATA_HEADER; (*header_read)++; strncpy(data->comment, line, strlen(line)-2); } else if (status == MB_SUCCESS && (line[0] == '1' || line[0] == '4')) { data->kind = MB_DATA_HEADER; (*header_read) = 1; strncpy(data->comment, line, strlen(line)-2); } else if (status == MB_SUCCESS && line[0] == '#') { data->kind = MB_DATA_COMMENT; strncpy(data->comment,&line[1],strlen(line)-3); } else if (status == MB_SUCCESS && line[0] == '3') { data->kind = MB_DATA_DATA; /* get survey id */ shift = 1; for (i=0;i<8;i++) data->survey_id[i] = line[i+shift]; /* get time */ shift += 8; mb_get_int(&data->timezone, &line[shift], 5); shift += 5; data->timezone = data->timezone / 100; mb_get_int(&itmp, &line[shift], 2); shift += 2; mb_fix_y2k(verbose, itmp, &data->time_i[0]); mb_get_int(&data->time_i[1], &line[shift], 2); shift += 2; mb_get_int(&data->time_i[2], &line[shift], 2); shift += 2; mb_get_int(&data->time_i[3], &line[shift], 2); shift += 2; mb_get_int(&itmp, &line[shift], 5); shift += 5; data->time_i[4] = 0.001 * itmp; dtmp = (itmp - 1000 * data->time_i[4]) * 0.06; data->time_i[5] = (int) dtmp; data->time_i[6] = 1000000 * (dtmp - data->time_i[5]); mb_get_time(verbose,data->time_i,&data->time_d); /* get nav */ neg_unit = 8; if (line[shift] == '-') { shift += 1; neg_unit = 7; } mb_get_int(&itmp, &line[shift], neg_unit); shift += neg_unit; data->latitude = 0.00001 * itmp; if (neg_unit == 7) data->latitude = -data->latitude; neg_unit = 9; if (line[shift] == '-') { shift += 1; neg_unit = 8; } mb_get_int(&itmp, &line[shift], neg_unit); shift += neg_unit; data->longitude = 0.00001 * itmp; if (neg_unit == 8) data->longitude = -data->longitude; mb_get_int(&data->nav_type, &line[shift], 1); shift += 1; /* get bath */ mb_get_int(&itmp, &line[shift], 6); shift += 6; data->tt = 0.0001 * itmp; mb_get_int(&itmp, &line[shift], 6); shift += 6; data->bath = 0.1 * itmp; mb_get_int(&data->bath_corr, &line[shift], 2); shift += 2; mb_get_int(&data->bath_type, &line[shift], 1); shift += 1; if (data->bath > 0.0 && data->bath < 99999.9) { data->flag = MB_FLAG_NONE; } else { data->flag = MB_FLAG_NULL; } /* get magnetics */ mb_get_int(&itmp, &line[shift], 6); shift += 6; data->mag_tot_1 = 0.1 * itmp; mb_get_int(&itmp, &line[shift], 6); shift += 6; data->mag_tot_2 = 0.1 * itmp; mb_get_int(&itmp, &line[shift], 6); shift += 6; data->mag_res = 0.1 * itmp; mb_get_int(&data->mag_res_sensor, &line[shift], 1); shift += 1; mb_get_int(&itmp, &line[shift], 5); shift += 5; data->mag_diurnal = 0.1 * itmp; mb_get_int(&itmp, &line[shift], 6); shift += 6; data->mag_altitude = itmp; /* get gravity */ mb_get_int(&itmp, &line[shift], 7); shift += 7; data->gravity = 0.1 * itmp; mb_get_int(&itmp, &line[shift], 6); shift += 6; data->eotvos = 0.1 * itmp; mb_get_int(&itmp, &line[shift], 5); shift += 5; data->free_air = 0.1 * itmp; mb_get_int(&data->seismic_line, &line[shift], 5); shift += 5; mb_get_int(&data->seismic_shot, &line[shift], 6); shift += 6; /* get nav quality */ mb_get_int(&data->nav_quality, &line[shift], 1); shift += 1; } else if (status == MB_SUCCESS && line[0] == '5') { data->kind = MB_DATA_DATA; /* get survey id */ shift = 1; for (i=0;i<8;i++) data->survey_id[i] = line[i+shift]; /* get time */ shift += 8; mb_get_int(&data->timezone, &line[shift], 3); shift += 3; mb_get_int(&data->time_i[0], &line[shift], 4); shift += 4; mb_get_int(&data->time_i[1], &line[shift], 2); shift += 2; mb_get_int(&data->time_i[2], &line[shift], 2); shift += 2; mb_get_int(&data->time_i[3], &line[shift], 2); shift += 2; mb_get_int(&itmp, &line[shift], 5); shift += 5; data->time_i[4] = 0.001 * itmp; dtmp = (itmp - 1000 * data->time_i[4]) * 0.06; data->time_i[5] = (int) dtmp; data->time_i[6] = 1000000 * (dtmp - data->time_i[5]); mb_get_time(verbose,data->time_i,&data->time_d); /* get nav */ neg_unit = 8; if (line[shift] == '-') { shift += 1; neg_unit = 7; } mb_get_int(&itmp, &line[shift], neg_unit); shift += neg_unit; data->latitude = 0.00001 * itmp; if (neg_unit == 7) data->latitude = -data->latitude; neg_unit = 9; if (line[shift] == '-') { shift += 1; neg_unit = 8; } mb_get_int(&itmp, &line[shift], neg_unit); shift += neg_unit; data->longitude = 0.00001 * itmp; if (neg_unit == 8) data->longitude = -data->longitude; mb_get_int(&data->nav_type, &line[shift], 1); shift += 1; /* get bath */ mb_get_int(&itmp, &line[shift], 6); shift += 6; data->tt = 0.0001 * itmp; mb_get_int(&itmp, &line[shift], 6); shift += 6; data->bath = 0.1 * itmp; mb_get_int(&data->bath_corr, &line[shift], 2); shift += 2; mb_get_int(&data->bath_type, &line[shift], 1); shift += 1; if (data->bath > 0.0 && data->bath < 99999.9) { data->flag = MB_FLAG_NONE; } else { data->flag = MB_FLAG_NULL; } /* get magnetics */ mb_get_int(&itmp, &line[shift], 6); shift += 6; data->mag_tot_1 = 0.1 * itmp; mb_get_int(&itmp, &line[shift], 6); shift += 6; data->mag_tot_2 = 0.1 * itmp; mb_get_int(&itmp, &line[shift], 6); shift += 6; data->mag_res = 0.1 * itmp; mb_get_int(&data->mag_res_sensor, &line[shift], 1); shift += 1; mb_get_int(&itmp, &line[shift], 5); shift += 5; data->mag_diurnal = 0.1 * itmp; mb_get_int(&itmp, &line[shift], 6); shift += 6; data->mag_altitude = itmp; /* get gravity */ mb_get_int(&itmp, &line[shift], 7); shift += 7; data->gravity = 0.1 * itmp; mb_get_int(&itmp, &line[shift], 6); shift += 6; data->eotvos = 0.1 * itmp; mb_get_int(&itmp, &line[shift], 5); shift += 5; data->free_air = 0.1 * itmp; mb_get_int(&data->seismic_line, &line[shift], 5); shift += 5; mb_get_int(&data->seismic_shot, &line[shift], 6); shift += 6; /* get nav quality */ mb_get_int(&data->nav_quality, &line[shift], 1); shift += 1; } /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n",function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",*error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return status */ return(status); } /*--------------------------------------------------------------------*/ int mbr_mgd77txt_wr_data(int verbose, void *mbio_ptr, void *data_ptr, int *error) { char *function_name = "mbr_mgd77txt_wr_data"; int status = MB_SUCCESS; struct mb_io_struct *mb_io_ptr; struct mbf_mgd77txt_struct *data; char line[MB_COMMENT_MAXLINE]; int itmp; int write_status; int shift; int i; /* print input debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> called\n",function_name); fprintf(stderr,"dbg2 Revision id: %s\n",rcs_id); fprintf(stderr,"dbg2 Input arguments:\n"); fprintf(stderr,"dbg2 verbose: %d\n",verbose); fprintf(stderr,"dbg2 mbio_ptr: %p\n",(void *)mbio_ptr); fprintf(stderr,"dbg2 data_ptr: %p\n",(void *)data_ptr); } /* get pointer to mbio descriptor */ mb_io_ptr = (struct mb_io_struct *) mbio_ptr; /* get pointer to raw data structure */ data = (struct mbf_mgd77txt_struct *) data_ptr; /* handle the data */ if (data->kind == MB_DATA_HEADER) { sprintf(line,"%s\r\n",data->comment); } else if (data->kind == MB_DATA_COMMENT) { sprintf(line, "#%s\r\n", data->comment); } else if (data->kind == MB_DATA_DATA) { /* set data record id */ shift = 0; line[0] = '5'; shift += 1; /* get survey id */ for (i=0;i<8;i++) line[i+shift] = data->survey_id[i]; shift += 8; /* get time */ sprintf(&line[shift], "%3.3d", data->timezone); shift += 3; sprintf(&line[shift], "%4.4d", data->time_i[0]); shift += 4; sprintf(&line[shift], "%2.2d", data->time_i[1]); shift += 2; sprintf(&line[shift], "%2.2d", data->time_i[2]); shift += 2; sprintf(&line[shift], "%2.2d", data->time_i[3]); shift += 2; itmp = (1000.0 * data->time_i[4]) + (1000.0 * (data->time_i[5]/60.0)) + (1000.0 * ((data->time_i[6]/1000000.0)/60.0)); sprintf(&line[shift], "%5.5d", itmp); shift += 5; /* get nav */ itmp = 100000 * data->latitude; if (itmp < 0) { sprintf(&line[shift], "-"); shift += 1; sprintf(&line[shift], "%7.7d", itmp*-1); shift += 7; } else { sprintf(&line[shift], "%8.8d", itmp); shift += 8; } itmp = 100000 * data->longitude; if (itmp < 0) { sprintf(&line[shift], "-"); shift += 1; sprintf(&line[shift], "%8.8d", itmp*-1); shift += 8; } else { sprintf(&line[shift], "%9.9d", itmp); shift += 9; } sprintf(&line[shift], "%1.1d", data->nav_type); shift += 1; /* get bath */ if (data->flag == MB_FLAG_NONE) { itmp = 10000 * data->tt; sprintf(&line[shift], "%6.6d", itmp); shift += 6; itmp = 10 * data->bath; sprintf(&line[shift], "%6.6d", itmp); shift += 6; } else { itmp = 999999; sprintf(&line[shift], "%6.6d", itmp); shift += 6; sprintf(&line[shift], "%6.6d", itmp); shift += 6; } sprintf(&line[shift], "%2.2d", data->bath_corr); shift += 2; sprintf(&line[shift], "%1.1d", data->bath_type); shift += 1; /* get magnetics */ itmp = 10 * data->mag_tot_1; sprintf(&line[shift], "%6.6d", itmp); shift += 6; itmp = 10 * data->mag_tot_2; sprintf(&line[shift], "%6.6d", itmp); shift += 6; itmp = 10 * data->mag_res; sprintf(&line[shift], "%6.6d", itmp); shift += 6; sprintf(&line[shift], "%1.1d", data->mag_res_sensor); shift += 1; itmp = 10 * data->mag_diurnal; sprintf(&line[shift], "%5.5d", itmp); shift += 5; itmp = data->mag_altitude; sprintf(&line[shift], "%6.6d", itmp); shift += 6; /* get gravity */ itmp = 10 * data->gravity; sprintf(&line[shift], "%7.7d", itmp); shift += 7; itmp = 10 * data->eotvos; sprintf(&line[shift], "%6.6d", itmp); shift += 6; itmp = 10 * data->free_air; sprintf(&line[shift], "%5.5d", itmp); shift += 5; sprintf(&line[shift], "%5.5d", data->seismic_line); shift += 5; sprintf(&line[shift], "%6.6d", data->seismic_shot); shift += 6; /* get nav quality */ sprintf(&line[shift], "%1.1d", data->nav_quality); shift += 1; /* set end of line */ sprintf(&line[shift], "\r\n"); shift += 2; line[shift] = '\0'; } if ((write_status = fputs(line, mb_io_ptr->mbfp)) > 0) { *error = MB_ERROR_NO_ERROR; status = MB_SUCCESS; } else { *error = MB_ERROR_WRITE_FAIL; status = MB_FAILURE; } /* print output debug statements */ if (verbose >= 5) { fprintf(stderr,"\ndbg5 Data record kind in MBIO function <%s>\n",function_name); fprintf(stderr,"dbg5 kind: %d\n",data->kind); } /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n",function_name); fprintf(stderr,"dbg2 Return values:\n"); fprintf(stderr,"dbg2 error: %d\n",*error); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return status */ return(status); } /*--------------------------------------------------------------------*/