/*-------------------------------------------------------------------- * The MB-system: mb_coor_scale.c 1/21/93 * $Id: mb_coor_scale.c 1917 2012-01-10 19:25:33Z caress $ * * Copyright (c) 1993-2012 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. *--------------------------------------------------------------------*/ /* * The function mb_coor_scale.c returns scaling factors * to turn longitude and latitude differences into distances in meters. * This function is based on code by James Charters (Scripps Institution * of Oceanography). * * Author: D. W. Caress * Date: January 21, 1993 * * $Log: mb_coor_scale.c,v $ * Revision 5.4 2008/09/27 03:27:10 caress * Working towards release 5.1.1beta24 * * Revision 5.3 2006/09/11 18:55:52 caress * Changes during Western Flyer and Thomas Thompson cruises, August-September * 2006. * * Revision 5.2 2003/04/17 21:05:23 caress * Release 5.0.beta30 * * Revision 5.1 2002/09/18 23:32:59 caress * Release 5.0.beta23 * * Revision 5.0 2000/12/01 22:48:41 caress * First cut at Version 5.0. * * Revision 4.9 2000/10/11 01:02:30 caress * Convert to ANSI C * * Revision 4.8 2000/09/30 06:26:58 caress * Snapshot for Dale. * * Revision 4.7 1998/10/05 18:32:27 caress * MB-System version 4.6beta * * Revision 4.6 1997/04/21 17:02:07 caress * MB-System 4.5 Beta Release. * * Revision 4.6 1997/04/21 17:02:07 caress * MB-System 4.5 Beta Release. * * Revision 4.5 1996/04/22 13:21:19 caress * Now have DTR and MIN/MAX defines in mb_define.h * * Revision 4.5 1996/04/22 13:21:19 caress * Now have DTR and MIN/MAX defines in mb_define.h * * Revision 4.4 1996/04/22 10:57:09 caress * DTR define now in mb_io.h * * Revision 4.3 1995/03/06 19:38:54 caress * Changed include strings.h to string.h for POSIX compliance. * * Revision 4.2 1994/10/21 12:11:53 caress * Release V4.0 * * Revision 4.1 1994/07/29 18:46:51 caress * Changes associated with supporting Lynx OS (byte swapped) and * using unix second time base (for time_d values). * * Revision 4.0 1994/03/05 23:55:38 caress * First cut at version 4.0 * * Revision 4.1 1994/03/03 03:39:43 caress * Fixed copyright message. * * Revision 4.0 1994/02/21 04:03:10 caress * First cut at new version. No changes. * * Revision 3.1 1993/05/14 22:32:39 sohara * fixed rcs_id message * * Revision 3.1 1993/05/14 22:32:39 sohara * fixed rcs_id message * * Revision 3.0 1993/04/23 15:45:07 dale * Initial version * */ /* standard include files */ #include #include #include /* mbio include files */ #include "../../include/mb_status.h" #include "../../include/mb_io.h" #include "../../include/mb_define.h" /* ellipsoid coefficients from World Geodetic System Ellipsoid of 1972 * - see Bowditch (H.O. 9 -- American Practical Navigator). */ #define C1 111412.84 #define C2 -93.5 #define C3 0.118 #define C4 111132.92 #define C5 -559.82 #define C6 1.175 #define C7 0.0023 static char rcs_id[]="$Id: mb_coor_scale.c 1917 2012-01-10 19:25:33Z caress $"; /*--------------------------------------------------------------------*/ int mb_coor_scale(int verbose, double latitude, double *mtodeglon, double *mtodeglat) { char *function_name = "mb_coor_scale"; int status; double radlat; /* 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 latitude: %f\n",latitude); } /* check that the latitude value is sensible */ if (fabs(latitude) <= 90.0) { radlat = DTR*latitude; *mtodeglon = 1./fabs(C1*cos(radlat) + C2*cos(3*radlat) + C3*cos(5*radlat)); *mtodeglat = 1./fabs(C4 + C5*cos(2*radlat) + C6*cos(4*radlat) + C7*cos(6*radlat)); status = MB_SUCCESS; } /* set error flag if needed */ else { status = MB_FAILURE; } /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n",function_name); fprintf(stderr,"dbg2 Revision id: %s\n",rcs_id); fprintf(stderr,"dbg2 Return arguments:\n"); fprintf(stderr,"dbg2 mtodeglon: %g\n",*mtodeglon); fprintf(stderr,"dbg2 mtodeglat: %g\n",*mtodeglat); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return status */ return(status); } /*--------------------------------------------------------------------*/ int mb_apply_lonflip(int verbose, int lonflip, double *longitude) { char *function_name = "mb_apply_lonflip"; int status; /* 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 lonflip: %d\n",lonflip); fprintf(stderr,"dbg2 longitude: %f\n",*longitude); } /* apply lonflip */ if (lonflip < 0) { if (*longitude > 0.) *longitude = *longitude - 360.; else if (*longitude < -360.) *longitude = *longitude + 360.; } else if (lonflip == 0) { if (*longitude > 180.) *longitude = *longitude - 360.; else if (*longitude < -180.) *longitude = *longitude + 360.; } else { if (*longitude > 360.) *longitude = *longitude - 360.; else if (*longitude < 0.) *longitude = *longitude + 360.; } /* print output debug statements */ if (verbose >= 2) { fprintf(stderr,"\ndbg2 MBIO function <%s> completed\n",function_name); fprintf(stderr,"dbg2 Revision id: %s\n",rcs_id); fprintf(stderr,"dbg2 Return arguments:\n"); fprintf(stderr,"dbg2 longitude: %f\n",*longitude); fprintf(stderr,"dbg2 Return status:\n"); fprintf(stderr,"dbg2 status: %d\n",status); } /* return status */ return(status); } /*--------------------------------------------------------------------*/