   /******************************************************************
    *                        File earlybirdlib.h                     *
    *                                                                *
    *  Include file for WC/ATWC EarlyBird library functions.         *
    *                                                                *
    *  August, 2010: Added Dave Nyland's display sorts               *
    *  2001: Paul Whitmore, NOAA-WCATWC - paul.whitmore@noaa.gov     *
    *                                                                *
    ******************************************************************/
                                                 
#if defined(_WINNT)
#include <windows.h>
#else
#define WORD short
typedef struct _SYSTEMTIME {
  WORD wYear;
  WORD wMonth;
  WORD wDayOfWeek;
  WORD wDay;
  WORD wHour;
  WORD wMinute;
  WORD wSecond;
  WORD wMilliseconds;
} SYSTEMTIME;
#define min(a,b) (((a)<(b))?(a):(b))
#define max(a,b) (((a)>(b))?(a):(b))
#endif

#ifndef EARLYBIRDLIB_H                               
#define EARLYBIRDLIB_H                              
                                  
/* Necessary include files */                      
#include <math.h>
#include <time.h>
#include <trace_buf.h>
#include <transport.h>     
                                      
/* Definitions */
#define PI      3.14159265358979323846
#define TWO_PI  6.28318530717958647692
#define TWOPI   6.28318530717958647692
#define RAD     0.017453292519943   /* Convert degrees to radians */
#define DEG     57.295779513083     /* Convert radians to degrees */
#define DEGREE  57.295779513083     
#define DEGTOKM 111.194927          /* km per degree */

#define MAX_TEMP       30000  /* Diskrw - Max size of temp data array */

#define MAX_SCREENS       10  /* Maximum # screen subsets */
#define DISPLAY_BUFFER_SIZE 65536 /* Trace display buffer size in samps. */
#define BROADBAND_SCALE  4.0  /* Broadband trace scaling decrement */
#define MAX_NUM_WWA_SEGS   2  /* Max # of w/w/a segments in Threat db */
#define MAX_NUM_INFO_SEGS  6  /* Max # of info segments in Threat db */
#define MAX_NUM_MODELS  1000  /* Max. # of pre-computed models on disk */

#define PHASE_LENGTH       8  /* Filters - Maximum phase length for PICKTWC */
#define MAX_FILT_SECTIONS  5  /* Maximum # of 2nd order IIR filters to apply */
#define FILTER_ORDER       3  /* Number of poles (filter order) */
// Added more stations.
#define MAX_STATIONS    1000  /* Max # stations in pick station file */
#define MAX_PHASES        20  /* Max # phases to overlay on trace */
#define MAX_ALARM_STN    100  /* Max number of stations per region */
#define NUM_AT_ZERO        6  /* # samps at 0 before assuming dead trace */
#define SPECTRO_SIZE   1024  /* Max Array size for spectrograms */
 
#define MAXREGS          100  /* Maximum depth regions */
#define CIRC_BUFFER_SIZE 24000/* # of samples in main buffer, 10 min @ 40 Hz */
#define LGSECONDS        150  /* # secs after P to look for max LG for Ml */
#define PRE_P_TIME        30  /* Secs of pre-event data to include in P_ARRAY */
#define PRE_LP_TIME      120  /* Amt of Pre-P data to determine bckgrnd noise */
#define P_ARRAY ((PRE_P_TIME+LGSECONDS)*100)  /* Size of P-wave buffer */
#define MAX_NN_SAMPS    4096  /* # Samples to save for Neural Net checker */
#define AIC_NUM_SEC       20  /* Number of seconds to re-analyze pick */
#define NN_NUM_SEC        12   /* # secs for neural net check - use whole # */

#define MAX_HYPO_SIZE    512  /* Maximum size of TYPE_HYPOTWC */
#define MAX_QUAKES       100  /* Max # quakes to show on display */

#define MAX_PICKTWC_SIZE 512  /* Max size (bytes) of PickTWC messages */ 
#define FIRST_MOTION_SAMPS 3  /* # samples after trigger which must be in same
                                 direction for a 1st motion to be declared */
#define RECENT_ALARM_TIME 900 /* # seconds late data can be and still run thru
                                 alarm */							   
#define FREQ_MIN         1.5  /* Minimum frequency to declare local event */
#define DELTA_TELE       45.  /* Minimum distance (degrees) for teleseismic
                                 declaration */

#define NUM_MAJOR_CITIES    10  /* At beginning of cities.dat */
#define NUM_CITIES          85  /* Total reference cities (including above) */
#define NUM_MAJOR_CITIES_EC 17  /* At beginning of cities-ec.dat */
#define NUM_CITIES_EC      131  /* Total reference cities (including above) */

// Added more stations.
#define MAX_STATION_DATA  2500  /* Max number of stations with metadata */
#define DEPTH_LEVELS_IASP   76  /* # depth levels in iaspei91 table */
#define IASP_DEPTH_INC     10.  /* Depth increment (km) in table */
#define IASP_DIST_INC      0.5  /* Distance increment (deg.) in table */
#define IASP_NUM_PER_DEP   361  /* # p-times per depth increment */
#define QUAKE_ITER          16	/* Iterations in Solution loop */
#define	COEFFSIZE    MAX_STATION_DATA	/* Coefficient array in QuakeSolve */
#define DEPTH         0.006278  /* 40km depth normalized with radius 6371km */
#define DEPTHKM             40  /* Depth in km specified by DEPTH */
#define EARTHRAD         6371.0 /* Earth radius in km */
#define MAX_TO_KO            3  /* Max # stations to remove from hypo by res. -
                                   Presently (2/2001) max is 3. */
#define DEFAULT_DEPTH       20  /* Default depth for location in km */
#define DEFAULT_MAXDEPTH    50  /* Default max depth for location in km */
#define EQSIZE           10000  /* Number of quake depth data */

#define MB_TIME             30  /* # seconds after P  to allow MBs */
#define MM_TIME            800  /* # seconds after R  to allow Mms (make same as MM_BACK...) */
#define MS_TIME           1800  /* # seconds after R  to allow Mss */
#define MS_BACKGROUND_TIME 240  /* Time to compute background MDF (seconds) */
#define MM_BACKGROUND_TIME 800  /* Time to compute background Mm noise (sec) */
#define MAXMWPARRAY      20000  /* Max size of Mwp buffers */
#define MAX_MWP_TIME       300  /* Max. time (secs) needed for Mwp comps */
#define NORM_MWP_TIME      120  /* Normal window (secs) used in Mwp comps */
#define MWP_BACKGROUND_TIME 30  /* Time (seconds) to evaluate background Noise */
#define MAX_ZP              50  /* Maximum number of poles/zeros */
#define GF_TRACELENGTH     200  /* # samples to save for each GF */
#define GF_PREEVENT         10  /* # GF samples to save prior to P */
#define DATA_PREEVENT      100  /* # samples to save prior to P in data */
#define DATA_POSTEVENT     100  /* # samples to save beyond EndTime in data */
#define MTBUFFER_SIZE (GF_TRACELENGTH+DATA_PREEVENT+DATA_POSTEVENT+1)
                                /* Max # data samples */

/* Mm Stuff */
#define MAX_SPECTRA        512  /* Max # spectral amplitudes */
#define NUM_PATH_PER        28  /* # different periods specified in regions */
#define NUM_ATTEN            7  /* # crustal types */
#define NUM_LAT_BOUNS       18  /* # lat grids for crust definition */
#define NUM_LON_BOUNS       36  /* # lon grids for crust definition */
#define MMBUFFER_SIZE    16384  /* Max # samples to send to Mm functions */
#define MAX_FFT_PER        100
#define MAX_FFT          16384  /* Max samples in FFT */
#define MM_MIN_TIME         60. /* Minimum # seconds necessary for Mm */

/* Structure defnitions */
typedef struct FCOMPLEX {float r, i;} fcomplex;  /* Complex definition */

typedef struct {
   char     szRegionName[32];  /* Name this alarm region */
   int      iAlarmThresh;      /* Number of picks needed to call alarm */
   int      iNumPicksCnt;      /* Running total of picks */
   double   dMaxTime;          /* Max time in s allowed between all picks */
   double   dLastTime;         /* Time (1/1/70 s) of pick #1 */
   double   dThresh;           /* PStrength Threshold to exceed */
   int      iNumStnInReg;      /* Number of stations to use in region */
   char     szStation[MAX_ALARM_STN][TRACE_STA_LEN];  /* Station Name */
   char     szStnAlarm[MAX_ALARM_STN][TRACE_STA_LEN]; /* Stations alarmed */
} ALARMSTRUCT;

typedef struct {               /* AZIDELT - Distance / azimuth structure */
   double      dDelta;         /* Distance in degrees between 2 pts. on earth */
   double      dAzimuth;       /* Azimuth in degrees between two points */
} AZIDELT;

typedef struct {           /* CHNLHEADER - If this must change, add to bottom */
   char        szStation[6]; 	/* Station name */
   char        szChannel[6];    /* Channel type */
   char        szNetID[4];      /* Seismic Network name */
   SYSTEMTIME  stStartTime;     /* Start Time (exact) for each trace */
   double      dSampRate;       /* Samples per second */
   long        lNumSamps;       /* # of samples in this file for this chn */
   int         iBytePerSamp;    /* Number of bytes / sample (2 or 4) */
   int         iTrigger;        /* 1=triggered data, 0=continuous */
   int         iSignalToNoise;  /* S:N ratio used in P-picker */
   int         iPickStatus;     /* 1=initialize, 0=don't pick, 2=pick 
                                   it, 3=already picked */
   int         iStationType;    /* Model of seismometer (see STATION) */
   double      dLat;            /* Latitude in geographic coords(after 2/2001)*/
   double      dLon;            /* Longitude in geographic coords */
   double      dElevation;      /* Station Elevation (in meters-prior to 2/2001
                                                      normalized) */
   double      dGain;           /* Gain (bb in counts/m/s) */
   double      dGainCalibration;/* WC&ATWC analog gain factors */
   double      dClipLevel;      /* Max. number of counts for signal */
   double      dTimeCorrection; /* Analog signal time delays */
   double      dScaleFactor;    /* Empirical factor for screen trace scaling */
   long        lUnused1;        /* For future expansion */
   int         iUnused1;        /* For future expansion */
   int         iUnused2;        /* For future expansion */
   double      dUnused1;        /* For future expansion */
   double      dUnused2;        /* For future expansion */
   double      dUnused3;        /* For future expansion */
   char        szUnunsed1[32];  /* For future expansion */
   char        szLocation[4];  /* Location */
} CHNLHEADER;                       

typedef struct {           /* CITY - Structure with city locations */
   double  dLat;           /* City lat. (geocentric) */
   double  dLon;           /* City long. (geocentric) */
   char    szLoc[32];      /* City name */
} CITY;

typedef struct {           /* CITYDIS - Struct. with loc... from nearest city */  
   int     iDis[2];        /* distance in miles */
   char    *pszDir[2];     /* direction (SW, NE, etc.) */
   char    *pszLoc[2];     /* city name */
} CITYDIS;

typedef struct {              /* DISKHEADER - Disk Header Structure */
   SYSTEMTIME  stStartTime;   /* Start Time of file (nominal) */
   int         iNumChans;     /* Number of channels of data / file */
   int         iChnHdrSize;   /* Number of bytes / channel header */
} DISKHEADER;

typedef struct {        /* Structure of earthquake avg and max depths */
   int iLat;            /* Geographic (+/-) latitude of point with depth data */
   int iLon;            /* Geographic (+/-) longitude of point with depth data*/
   int iAveDepth;       /* Average shallow (<50km) eq depth for this area */
   int iMaxDepth;       /* Maximum likely depth for this area (km) */
} EQDEPTHDATA;

typedef struct {              /* LATLON - Location structure */
   double  dLat;              /* Lat (geocentric) */
   double  dLon;              /* Lon (geocentric) */
   double  dCoslat;           /* Cosine of Latitude */
   double  dSinlat;           /* Sine of latitude */
   double  dCoslon;           /* Cosine of Longitude */
   double  dSinlon;           /* Sine of longitude */
} LATLON;

typedef struct {            /* HYPO - hypocenter location structure */
   double  dLat;            /* Geocentric epicentral latitude */
   double  dLon;            /* Geocentric epicentral longitude */ 
   double  dCoslat;         /* Cosine of geocentric epicentral latitude */
   double  dSinlat;         /* Sine of geocentric epicentral latitude */
   double  dCoslon;         /* Cosine of geocentric epicentral longitude */
   double  dSinlon;         /* Sine of geocentric epicentral longitude */
   double  dOriginTime;     /* Origin time in 1/1/70 seconds */
   SYSTEMTIME stOTime;      /* Origin Time in SYSTEMTIME structure */
   SYSTEMTIME stOTimeRnd;   /* Origin Time in SYSTEMTIME structure rounded to
                               nearest minute */
   char    szOTimeRnd[5];   /* Origin time (string) rounded to min. */
   double  dDepth;          /* Hypocenter depth (in km) */
   double  dPreferredMag;   /* Magnitude to use for this event */
   char    szPMagType[3];   /* Magnitude type of dPreferredMag (l, b, etc.) */
   LATLON  llEpiGG;         /* Location in geographic coords. (+/-) */
   char    szLat[12];       /* Latitude in string form (with N,S...) */
   char    szLon[12];       /* Longitude in string form (with E,W..) */
   int     iNumPMags;       /* Number of stations used in dPreferredMag */
   int     iAlarmIssued;    /* 0-no alarm, 1-RESPOND sent, 2-good loc made */
   char    szQuakeID[32];   /* Event ID */
   int     iVersion;        /* Version of location (1, 2, ...) */
   int     iBullNo;         /* Bulletin number */
   int     iDepthControl;   /* 1->Fix at 40; 2->Float; 3->User fix; */
   int     iNumPs;          /* Number of P's used in quake solution */
   int     iNumBadPs;       /* Number of P's excluded from solution */
   double  dAvgRes;         /* Residual average */
   int     iAzm;            /* Azimuthal coverage in degrees */
   double  dNearestDist;    /* Closest station epicentral distance in deg. */
   double  dFirstPTime;     /* 1/1/70 time (seconds) of earliest P */
   int     iGoodSoln;       /* 0-bad, 1-soso, 2-good, 3-very good */
   double  dMbAvg;          /* Modified average Mb */
   int     iNumMb;          /* Number of Mb's in modified average */
   int     iNumMbClip;      /* Number of "clipped" Mb's */
   double  dMlAvg;          /* Modified average Ml */
   int     iNumMl;          /* Number of Ml's in modified average */
   int     iNumMlClip;      /* Number of "clipped" Ml's */
   double  dMSAvg;          /* Modified average MS */
   int     iNumMS;          /* Number of MS's in modified average */
   int     iNumMSClip;      /* Number of "clipped" MS's */
   double  dMwpAvg;         /* Modified average Mwp */
   int     iNumMwp;         /* Number of Mwp's in modified average */
   int     iNumMwpClip;     /* Number of "clipped" Mwp's */
   double  dMwAvg;          /* Modified average Mw */
   int     iNumMw;          /* Number of Mw's in modified average */
   int     iNumMwClip;      /* Number of "clipped" Mw's */
   double  dTheta;          /* Theta (energy/moment) ratio */
   double  dThetaSD;        /* Standard deviation in Theta computation */
   int     iNumTheta;       /* Number of stations used in theta computation */
   int     iMagOnly;        /* 0->New location, 1->updated magnitudes only */
   int     iUpdateMap;      /* 1->Update map in Earthvu */
} HYPO;

typedef struct {              /* MMSTUFF - Mm information */
   char    szStation[TRACE_STA_LEN];  /* Station name */
   char    szChannel[TRACE_CHAN_LEN]; /* Channel identifier (SEED notation) */
   char    szNetID[TRACE_NET_LEN];    /* Network ID */
   char    szLocation[TRACE_STA_LEN];  /* Location */
   double  dMmTravTime;   /* Rayleigh wave travel time for Mm (in seconds) */
   double  dMmStartTime;  /* Start time of Mm Rayleigh wave window (1/1/70)*/
   double  dMmEndTime;    /* End time of Mm Rayleigh wave window (1/1/70) */
   double  dMmMax;        /* Maximum Mm computed for this station */
   double  dPerMax;       /* Period of the maximum Mm */
   AZIDELT azdelt;        /* Epicentral distance (deg) and azimuth epi-sta */
} MMSTUFF;

typedef struct {                      /* MTSTUFF - mtinver information */
   char    szStation[TRACE_STA_LEN];  /* Station name */
   char    szChannel[TRACE_CHAN_LEN]; /* Channel identifier (SEED notation) */
   char    szNetID[TRACE_NET_LEN];    /* Network ID */
   char    szLocation[TRACE_STA_LEN];  /* Location */
   int     iNPts;                     /* # samples to use; Minimum GF and data*/
   double  dWt;                       /* Station weighting */
   int     iWvType;                   /* Wave Type (1=p vert, 2=p rad, 3=s Hor*/
   int     iUse;                      /* 0=don't use station, 1=use */
   int     iUseOrig;                  /* 0=don't use station, 1=use */
   double  dDt;                       /* Sample interval (s) */
   AZIDELT azdelt;        /* Epicentral distance (deg) and azimuth epi-sta */
   double  dMtStartTime;              /* Trace start time */
   double  dMtEndTime;                /* Trace end time */
   double  dSumSq;                    /* Signal Sum of squares */
   double  dData[MTBUFFER_SIZE];      /* Original velocity data */
   double  dDataFilt[GF_TRACELENGTH]; /* Filtered/deconvolved velocity data */
   double  dSynth[GF_TRACELENGTH];    /* Synthetic seismogram */
   double  dGFFiltVss[GF_TRACELENGTH];/* Filtered green's function */
   double  dGFFiltVds[GF_TRACELENGTH];/* Filtered green's function */
   double  dGFFiltFfds[GF_TRACELENGTH];/* Filtered green's function */
} MTSTUFF;

typedef struct {                      /* MTRESULTS - mtinver solution */
   double  dM0;                       /* in dyne-cm */
   double  dScalarMom;                /* Scalar moment in dyne-cm */
   double  dMag;                      /* Mw derived from dScalarMom */
   double  dMisFit[MAX_STATIONS];
   double  dStrike1;                  /* Fault plane solution params */
   double  dDip1;
   double  dRake1;
   double  dStrike2;
   double  dDip2;
   double  dRake2;
   double  dError;
   double  dMisFitMd;
   double  dMisFitAvg;
   double  dTWindow;                  /* Trace length in seconds */
   double  dMaxTLag;                  /* Max shift to traces (s) */
   int     iAlign;
   int     iPass;                     /* Pass 1, 2, or 3 */
   int     iAzGap;
} MTRESULTS;

typedef struct {           /* List of previously located quakes by loc_wcatwc */
   double  dOTime;         /* 1/1/70 origin time of quake */
   double  dLat;           /* Geographic (+/-) latitude */
   double  dLon;           /* Geographic (+/-) longitude */
   double  d1stPTime;      /* First good P-time in auto-loc (1/1/70 sec.) */
   int     iDepth;         /* Hypocenter depth (km) */
   double  dPreferredMag;  /* Magnitude to use for this event */
   int     iNumPMags;      /* Number of stations used in dPreferredMag */
   int     iGoodSoln;      /* 0-bad, 1-soso, 2-good, 3-very good */
   char    szQuakeID[32];  /* Event ID */
   int     iVersion;       /* Version of location (1, 2, ...) */
   int     iNumPs;         /* Number of P's used in quake solution */
   double  dAvgRes;        /* Sum of absolute values of residuals */
   double  dAzm;           /* Azimuthal coverage in degrees */
   double  dMbAvg;         /* Modified average Mb */
   int     iNumMb;         /* Number of Mb's in modified average */
   double  dMlAvg;         /* Modified average Ml */
   int     iNumMl;         /* Number of Ml's in modified average */
   double  dMSAvg;         /* Modified average MS */
   int     iNumMS;         /* Number of MS's in modified average */
   double  dMwpAvg;        /* Modified average Mwp */
   int     iNumMwp;        /* Number of Mwp's in modified average */
   double  dMwAvg;         /* Modified average Mw */
   int     iNumMw;         /* Number of Mw's in modified average */
   char    szPMagType[4];  /* Magnitude type of dPreferredMag (l, b, etc.) */
   double  dTheta;          /* Theta (energy/moment) ratio */
   double  dThetaSD;        /* Standard deviation in Theta computation */
   int     iNumTheta;       /* Number of stations used in theta computation */
} OLDQUAKE;

typedef struct {             /* Phase structure (interacts with NEICs taulib) */
   int     iNumPhases;                /* # phases at this station */
   char    szPhase[MAX_PHASES][PHASE_LENGTH]; /* Phase name */
   double  dPhaseTT[MAX_PHASES];      /* Phase travel time (since o-time) */
   double  dPhaseTime[MAX_PHASES];    /* 1/1/70 seconds time of phase */
} PHASE;

typedef struct                        /* Phase start and end times */
{
   double  dPStart;                   /* P start time in 1/1/70 s */
   double  dPEnd;                     /* P end time in 1/1/70 s */
   double  dPTravTime;                /* P travel time in seconds */
   double  dRStart;                   /* R start time in 1/1/70 s */
   double  dREnd;                     /* R end time in 1/1/70 s */
   double  dRTravTime;                /* R travel time in seconds */
} PHASE_WINDOWS;

typedef struct                      /* Structure for P-pick trace data */
{
   SYSTEMTIME  stStartTime;         /* Time at P-pick Index */
   double      dStartTime;          /* Time (in 1/1/70 secs) at P-pick index */
   long        lPIndex;             /* Index of P-pick */
   long        lNumSamps;           /* Number of samples to this point */
   long        lData[P_ARRAY];      /* Broadband data array */
   long        lFiltData[P_ARRAY];  /* Filtered data array */
   long        lDCOffset;           /* DC Offset for this trace */
   long        lFiltDCOffset;       /* Filtered DC Offset for this trace */
   long        lNoise;              /* Peak noise in BB signal */
} PSTRUCT;

typedef struct {
   int     iNumPers;                /* Number of frequencies with amps */
   double  dAmp[MAX_SPECTRA];       /* Spectral amplitude */
   double  dPer[MAX_SPECTRA];       /* Spectral period */
} SPECTRA;

typedef struct {              /* STATION - station parameters and variables */
   double  dLat;              /* Station geographic latitude (+=N, -=S) */
   double  dLon;              /* Station geographic longitude (+=E, -=W) */
   double  dCoslat;           /* These cos/sines are only used where the above*/
   double  dSinlat;           /*  are in geocentric, which is only in the */
   double  dCoslon;           /*  near station lookup table determination */
   double  dSinlon;  
   char    cFirstMotion;      /* ?=unknown, U=up, D=down */
   double  dAlarmAmp;         /* Signal amplitude to exceedin m/s */
   double  dAlarmDur;         /* Duration (sec) signal must exceed dAlarmDur */
   double  dAlarmLastSamp;    /* Last samp to exceed threshold in Alarm */
   double  dAlarmLastTriggerTime;/* Time when alarm was last triggered for
                                    this station */
   double  dAlarmMinFreq;     /* In hertz; low frequency limit condition */
   double  dAlarmSamp;        /* Alarm sample value in m/s */
   double  dAmp0;             /* Response scaling factor from calibs */
   double  dAvAmp;            /* Average signal amp per 12 cycle */
   double  dAveFiltNoise;     /* Moving average of RMS (filtered) */
   double  dAveLDC;           /* Moving average of filtered DC offset */
   double  dAveLDCRaw;        /* Moving average of unfilteredDC offset */
   double  dAveLDCRawOrig;    /* Moving average of unfilteredDC offset when
                                 Phase 1 passed */
   double  dAveLTA;           /* Moving average of average signal amplitude */
   double  dAveMDF;           /* Moving average of MDF */
   double  dAveMDFOrig;       /* Moving avg of MDF when Phase 1 was passed */
   double  dAveRawNoise;      /* Moving average of RMS */
   double  dAveRawNoiseOrig;  /* Moving avg of RMS when Phase 1 was passed */
   double  dAzimuth;          /* Epicenter/station azimuth */
   double  dClipLevel;        /* Max Counts which signal can attain */
   double  dCooze;            /* Variable used in locations */
   double  dDataEndTime;      /* 1/1/70s of last non-zero data */
   double  dDataStartTime;    /* 1/1/70s of 1st data != 0 */
   double  dDelta;            /* Epicentral distance in degrees */
   double  dElevation;        /* Station elevation in meters */
   double  dEndTime;          /* Time at end of last packet (1/1/70 seconds) */
   double  dExpectedPTime;    /* Expected P-time (1/1/70) for hypo in dummy */
   double  dFiltX1[MAX_FILT_SECTIONS];   /* Saved data for filter */ 
   double  dFiltX2[MAX_FILT_SECTIONS];   /* Saved data for filter */ 
   double  dFiltY1[MAX_FILT_SECTIONS];   /* Saved data for filter */ 
   double  dFiltY2[MAX_FILT_SECTIONS];   /* Saved data for filter */ 
   double  dFiltX1LP[MAX_FILT_SECTIONS]; /* Saved data for filter */ 
   double  dFiltX2LP[MAX_FILT_SECTIONS]; /* Saved data for filter */ 
   double  dFiltY1LP[MAX_FILT_SECTIONS]; /* Saved data for filter */ 
   double  dFiltY2LP[MAX_FILT_SECTIONS]; /* Saved data for filter */ 
   double  dFracDelta;        /* Fractional part of epi. distance */
   double  dFreq;             /* Dominant frequency of this P-pick */
   double  dGainCalibration;  /* Factor to converts SW cal amplitude to gain */
   double  dLTAThresh;        /* Phase 2 ampltude threshold */
   double  dLTAThreshOrig;    /* dLTAThresh at time Phase 1 first passed */
   double  dMaxD;             /* Maximum displacement value in signal */
   double  dMaxID;            /* Max integrated displacement value in signal */
   double  dMaxV;             /* Maximum velocity value in signal */
   double  dMbAmpGM;          /* Mb amplitude (ground motion in nm) */
   double  dMbMag;            /* Mb magnitude */
   double  dMbPer;            /* Mb Per data, per of lMbAmp (sec) */
   double  dMbTime;           /* 1/1/70 time (sec) at end of Mb T/A */
   double  dMDFThresh;        /* MDF to exceed to pass Phase 1 */
   double  dMDFThreshOrig;    /* dMDFThresh at time Phase 1 first passed */
   double  dMlAmpGM;          /* Ml amplitude (ground motion in nm) */
   double  dMlMag;            /* Ml magnitude */
   double  dMlPer;            /* Ml Per data, per of lMlAmp (sec) */
   double  dMlTime;           /* 1/1/70 time (sec) at end of Ml T/A */
   double  dMSAmpGM;          /* MS amplitude (ground motion in um) */
   double  dMSMag;            /* MS magnitude */
   double  dMSPer;            /* MS Per data, per of lMSAmp (sec) */
   double  dMSTime;           /* 1/1/70 time (sec) at end of MS T/A */
   double  dMwpIntDisp;       /* Maximum integrated disp. peak-to-peak amp */
   double  dMwpMag;           /* Mwp magnitude */
   double  dMwpTime;          /* Mwp window time in seconds */
   double  dMwAmpGM;          /* Mw amplitude (ground motion in um) */
   double  dMwMag;            /* Mw magnitude */
   double  dMwMagBG;          /* Mw magnitude based on background data */
   double  dMwPer;            /* Mw Per data */
   double  dMwTime;           /* 1/1/70 time (sec) at end of Mw T/A */
   double  dMwAmpSp[MAX_SPECTRA];   /* Spectral amplitude */
   double  dMwAmpSpBG[MAX_SPECTRA]; /* Background Spectral amplitude */
   double  dMwMagSp[MAX_SPECTRA];   /* Spectral magnitude */
   double  dMwMagSpBG[MAX_SPECTRA]; /* Background Spectral magnitude */
   double  dMwPerSp[MAX_SPECTRA];   /* Spectral period */
   double  dPEndTime;         /* End time of P wave window (1/1/70) */
   double  dPerMax;           /* Period of the maximum Mm */
   double  dPhaseTime;        /* Time (1/1/70 seconds) of picked phase */
   double  dPhaseTT[MAX_PHASES]; /* IASPEI91 Phase travel time (since o-time) */
   double  dPhaseTimeIasp[MAX_PHASES];/* 1/1/70 seconds time of Iaspei phase */
   double  dPStartTime;       /* Start time of P wave window (1/1/70)*/
   double  dPTravTime;        /* P wave travel time (in seconds) */
   double  dPStrength;        /* Ratio of P motion to background */
   double  dPTime;            /* P-time in seconds from 1/1/70 */
   double  dRes;              /* Location residual for this station */
   double  dResidualWeights;  /* Weighting factors for residual */
   double  dREndTime;         /* End time of Rayleigh wave window (1/1/70) */
   double  dRStartTime;       /* Start time of Rayleigh wave window (1/1/70) */
   double  dRTravTime;        /* Rayleigh wave travel time (in seconds) */
   double  dSampRate;         /* Sample rate in samps/sec */
   double  dScaleFactor;      /* Empirical factor for screen trace scaling */
   double  dScreenLat[MAX_SCREENS][2];/* Lat bounds of screen's stns. */
   double  dScreenLon[MAX_SCREENS][2];/* Lon bounds of screen's stns. */
   double  dScreenStart;      /* Time (1/1/70s) at trace start on screen */
   double  dSens;             /* Station sensitivity at 1 Hz in cts/m/s */
   double  dSnooze;           /* Variable used in locations */
   double  dStartTime;        /* 1/1/70 second at index = 0 (in ANALYZE: oldest
                                 time in buffer) */
   double  dSumLDC;           /* Accumulator for average DC amplitude */
   double  dSumLDCRaw;        /* Accumulator for average, unfiltered DC amp */
   double  dSumLTA;           /* Accumulator for average signal amplitude */
   double  dSumMDF;           /* Accumulator for MDF summation */
   double  dSumRawNoise;      /* Accumulator for RMS summation */
   double  dThetaEnergy;      /* Energy release calculated from this station */   
   double  dThetaMoment;      /* Moment calculated from this stn for Theta */   
   double  dTheta;            /* Theta for this station */   
   double  dTimeCorrection;   /* Transmission time delay (sec) to subtract
                                 from data */
   double  dTrigTime;         /* 1/1/70 time (sec) that Phase1 was passed */
   double  dVScale;           /* Vertical scaling factor */
   int     iAgency;           /* Agency indicator (not used after 5/08) */
   int     iAhead;            /* 1->data ahead of present; 0->data ok */
   int     iAlarmPage;        /* 0->no pager message, 1->pager alarm*/								 
   int     iAlarmSpeak;       /* 0->no spoken message, 1->speak alarm*/								 
   int     iAlarmStatus;      /* 1=Initialize digital alarm variables, 2=Process
                                 data in alarm, 3=Alarm declared, 0=No alarms*/
   int     iBin;              /* Associator bin for this pick */
   int     iBytePerSamp;      /* # bytes/sample - normally 4 */
   int     iCal;              /* 1 -> pick looks like calibration */
   int     iClipIt;           /* 1 -> Clip display amplitude; 0 -> don't */
   int     iComputeMwp;       /* 1=use this stn for Mwp, 0=don't */
   int     iDisplayStatus;    /* 0=toggled off, 1=on, 2=always off */
   int     iFirst;            /* 1=First packet for this station */
   int     iFiltStatus;       /* 1=Run data through SP filter, 0=don't */
   int     iHasWrapped;       /* 1 -> enough data recorded to fill buffer */
   int     iLPAlarmSent;      /* 1 -> LP Alarm for this station recently sent */ 
   int     iMbClip;           /* 1 if Clipped; 0 if not */       
   int     iMlClip;           /* 1 if Clipped; 0 if not */      
   int     iMSClip;           /* 1 if Clipped; 0 if not */
   int     iMwClip;           /* 1 if Clipped; 0 if not */
   int     iMwNumPers;        /* Number of frequencies with amps (Mm) */
   int     iNearbyStnArray[5];/* Array of 5 closest stations */
   int     iNPole;            /* # poles in response function */
   int     iNumPhases;        /* # phases at this station */
   int     iNZero;            /* # zeros in response function */
   int     iPickStatus;       /* 0=don't pick, 1=initialize, 2=pick it
                                 3=it's been picked, get mags */
   int     iSignalToNoise;    /* S:N ratio to exceed for P-pick */
   int     iStationDisp[MAX_SCREENS]; /* Flags to show if this stn on screen */
   int     iStationSortIndex; /* Used by the qsort to obtain the original
                                 station lineup dln 2/13/10 */
   int     iStationType;      /* Model of seismometer:
                                 1 = STS1     360s
                                 2 = STS2     130s
                                 3 = CMG-3NSN 30s
                                 4 = CMG-3T   100s
                                 5 = KS360i   360s
                                 6 = KS5400   350s
                                 7 = CMG-3    30s
                                 8 = CMG-40T  30s
                                 9 = CMG3TNSN 30s 
                                 10 = KS-10   20s
                                 11 = CMG3ESP_30  30s
                                 12 = CMG3ESP_60  60s
				 13 = Trillium 40 40s 
				 14 = CMG3ESP_120 (100) 120s 
				 15 = CMG40T_20   20s 
				 16 = CMG3T_360  360s 
				 17 = KS2000_120 120s 
				 18 = CMG 6TD    30s 
				 19 = Trillium 120 120s 
				 20 = Trillium 240 240s 
				 29 = unknown broadband (no cal) 
				 30 = EpiSensor FBA ES-T 
				 31 = Guralp-5T* 
				 40 = GT_S13 
				 41 = MP_L4 
				 42 = Ranger SS1 
                                 50 = generic LP (no cal)
                                 51 = ATWC LP Response (Hi gain)
                                 52 = ATWC LP Response (Low gain)
                                 100 = Generic SP (no cal)
                                 101 = ATWC SP Response (Hi gain)
                                 102 = ATWC SP Response (Medium gain)
                                 103 = ATWC SP Response (Low gain)
                                 999 = Not determined */
   int     iTrigger;          /* 1=triggered data, 0=continuous */
   int     iUseMe;            /* 2->This pick will not be removed by FindBadPs,
                                 1->Use this P in location; 0->Don't (auto KO)
                                 -1->Don't (manually knocked out) */
   long    lAlarmCycs;        /* # samples per half cycle */
   long    lAlarmP1;          /* Alarm Phase 1 pass flag */
   long    lAlarmSamps;       /* # samples ctr while alarm threshold exceeded */
   long    lCurSign;          /* Sign of current MDF for Phase 3 */
   long    lCycCnt;           /* Cycle ctr (if T/A in first MbCycles, 
                                 this is associated with Mb magnitude) */
   long    lCycCntLTA;        /* Cycle counter for LTAs */
   long    lEndData;          /* Data at last point of previous buffer */								 
   long    lFiltSamps;        /* Number of samples processed by filter per 
                                 sequence */
   int     lFirstMotionCtr;   /* Number of samples checked so far */
   long    lHit;              /* Number of hits counter for Phases 2 & 3 */
   long    lIndex;            /* Next index to write sample (on read, index of
				 dStartTime) - ANALYZE */
   long    lIndexToStartWrite;/* Used in ReadDiskNew - Index to write to */
   long    lLastSign;         /* Sign of last MDF for Phase 3 */
   long    lLTACtr;           /* Long term averages counter */
   long    lMagAmp;           /* Summation of 1/2 cycle amplitudes */
   double  dMaxPk;            /* Amp (p-p nm) for use in magnitude comp. */
   long    lMbPer;            /* Mb Per data, per of dMbAmpGM doubled in sec*10 */
   long    lMDFCnt;           /* Counter of cycles used in MDFTotal */
   long    lMDFNew;           /* Present MDF value */
   long    lMDFOld;           /* Last MDF value */
   long    lMDFRunning;       /* Running total of sample differences */
   long    lMDFRunningLTA;    /* Running total of sample differences for LTAs */
   long    lMDFTotal;         /* Total of MDFs over several cycles */
   long    lMis;              /* Num. of misses count for Phases 2 & 3 */
   long    lMlPer;            /* Ml Per data, per of dMlAmpGM doubled in sec*10 */
   long    lMwpCtr;           /* Index which counts samples from P for Mwp */
   long    lNumOsc;           /* # of osc. counter for Phase 3 */
   long    lPer;              /* Temporary period array - in seconds*10 */
   long    lPhase1;           /* Phase 1 passed flag */
   long    lPhase2;           /* Phase 2 passed flag */
   long    lPhase3;           /* Phase 3 passed flag */
   long    lPhase2Cnt;        /* Sample counter for timing Phase 2 */
   long    lPhase3Cnt;        /* Sample counter for timing Phase 3 */
   long    lPickBufIndex;     /* Present index of P adjustment buffer */
   long    lPickBufRawIndex;  /* Present index of Neural net buffer */
   long    lPickIndex;        /* Pick index; 0-10000->pick_wcatwc, 10000-20000
                                 ->develo; 20000-30000->hypo_display source*/
   long    lRawCircCtr;       /* Index counter for raw-data, circ buff */
   long    lRawCircSize;      /* Size (# samples) in plRawCircBuff */
   long    lRawNoise;         /* Max peak/trough signal difference */
   long    lRawNoiseOrig;     /* Max peak/trough signal difference when Phase1
                                 passed */
   long    lSampIndexF;       /* Next index to write in filt data circ. buff */
   long    lSampIndexR;       /* Next index to write in raw data circular buff */
   long    lSampNew;          /* Present sample */
   long    lSampOld;          /* Last sample */  
   long    lSampRaw;          /* Un-filter present sample */
   long    lSampsInLastPacket;/* # samps in last packet read in - ANALYZE */
   long    lSampsPerCyc;      /* Number of samples per half cycle */
   long    lSampsPerCycLTA;   /* Number of samples per half cycle in LTA*/
   long    lSWSim;            /* Count of similarities to sin wave cal */
   long    lTest1;            /* Phase 2 Test 1 passed */
   long    lTest2;            /* Phase 2 Test 2 passed */
   long    lTrigFlag;         /* Flag -> samp has passed MDF threshold */
   long    l13sCnt;           /* 13s counter; this delays p-picks 
                                 by 13 seconds after pick for Neural net */
   double  *pdRawDispData;    /* Pointer to Mwp buffer of displacement signal */
   double  *pdRawIDispData;   /* Pointer to integrated displacement signal */
   long    *plFiltCircBuff;   /* Pointer to buffer of filtered signal */
   long    *plPickBuf;        /* Pick buffer for P adjustment */
   long    *plPickBufRaw;     /* Pick buffer for Neural net checker */
   long    *plRawData;        /* Pointer to Mwp buffer of unfiltered signal */
   long    *plRawCircBuff;    /* Pointer to buffer of unfiltered signal */
   double  *pdTimeCircBuff;   /* Pointer to the buffer time -- JFP */
   SYSTEMTIME stPTime;        /* Time at P-pick Index in structure */
   SYSTEMTIME stStartTime;    /* Time of data at index 0 */
   char    szChannel[TRACE_CHAN_LEN]; /* Channel identifier (SEED notation) */
   char    szHypoID[32];      /* If pick made in hypo_display, this is the
                                 associated hypocenter ID */
   char    szNetID[TRACE_NET_LEN];    /* Network ID */
   char    szPhase[PHASE_LENGTH];     /* Pick Phase */
   char    szPhaseIasp[MAX_PHASES][PHASE_LENGTH]; /* Iaspei91 Phase name */
   char    szStation[TRACE_STA_LEN];  /* Station name */
   char    szLocation[TRACE_STA_LEN]; /* Location */
   char    szStationName[64]; /* Seismometer name */
   fcomplex zPoles[MAX_ZP];   /* Poles of response function */
   fcomplex zZeros[MAX_ZP];   /* Zeros of response function */
} STATION;               

typedef struct _THREAT      /* Threat structure read in from threat data file */
{                        
   char   szWarnUGC[MAX_NUM_WWA_SEGS][512];       /* Warning zone UGCs */
   char   szWarnHeadLineStd[MAX_NUM_WWA_SEGS][256];/* Warning zone Headlines (cap) */
   char   szWarnHeadLine[MAX_NUM_WWA_SEGS][512];  /* Warning zone Headlines */
   char   szWarnSegArea[MAX_NUM_WWA_SEGS][256];   /* Warning zone segment areas */
   char   szWarnPager[MAX_NUM_WWA_SEGS][128];     /* Warning zone segment areas */
   char   szAdvUGC[MAX_NUM_WWA_SEGS][512];        /* Advisory zone UGCs */
   char   szAdvHeadLineStd[MAX_NUM_WWA_SEGS][256];/* Adv zone Headlines (cap) */
   char   szAdvHeadLine[MAX_NUM_WWA_SEGS][512];   /* Advisory zone Headlines */
   char   szAdvSegArea[MAX_NUM_WWA_SEGS][256];    /* Advisory zone segment areas */
   char   szAdvPager[MAX_NUM_WWA_SEGS][128];      /* Warning zone segment areas */
   char   szWatchUGC[MAX_NUM_WWA_SEGS][512];      /* Watch zone UGCs */
   char   szWatchHeadLineStd[MAX_NUM_WWA_SEGS][256];/* Watch zone Headlines (cap) */
   char   szWatchHeadLine[MAX_NUM_WWA_SEGS][512]; /* Watch zone Headlines */
   char   szWatchSegArea[MAX_NUM_WWA_SEGS][256];  /* Watch zone segment areas */
   char   szWatchPager[MAX_NUM_WWA_SEGS][128];    /* Warning zone segment areas */
   char   szInfoAreaStd[MAX_NUM_INFO_SEGS][256];  /* Info only headlines (cap) */
   char   szInfoArea[MAX_NUM_INFO_SEGS][256];     /* Info only headlines */
   char   szCancelUGC[512];                       /* Cancel zone UGCs */
   char   szCancelHeadLineStd[256];               /* Cancel Headlines (cap) */
   char   szCancelHeadLine[256];                  /* Cancel Headlines */
   char   szCancelSegArea[256];                   /* Cancel segment areas */
   char   szSpecArea[128];                        /* String to place in Summary */
} THREAT;
typedef THREAT *PTHREAT;

typedef struct {       /* TSUTRAVTIME - Everything you need to know about TTT */
   short      sTravTimeTotalMin;     /* Tsunami travel times in minutes */
   short      sTravTimeHrMin;        /* Tsunami travel times in hhmm */
   time_t     lArrivalTimeEpochSec;  /* Arrival time in 1/1/70 seconds */
   char       szSiteNameLowLong[64]; /* Tsunami travel location (long/lower) */
   char       szSiteNameCapLong[64]; /* Tsunami travel location (long/caps) */
   char       szGeoCodes[128];       /* UGCs (onshore) tied to this site */
   char       szMarineZones[64];     /* UGCs (offshore) tied to this site */
   float      fPlaceRelBkPts;        /* Indicator showing relative loc to bkpt*/
   short      sSiteType;             /* Site Type - 0=other, 1=ETA pt,
                                        2=ETA and short list */
   short      sBreakPoint;           /* 1->site used as break point; 0->no */
   char       szSiteTGAbbr[8];       /* If tide gage site, abbreviation */
   float      fSiteLatLon[2];        /* Site lan/lon (geographic +/-) */
   char       szSiteCode[8];         /* Site name abbreviation */
   char       szSiteCountry[32];     /* Site country */
   char       szSiteCountryCode[8];  /* iso Site country code */
   char       szSiteName[32];        /* Site name city */
   char       szSiteState[32];       /* Site state or province */
   char       szSiteStateCode[8];    /* Site state abbr */
   char       cOcean;                /* Ocean indicator */
   char       cSiteNSEW[2];          /* Site lat/lon North, South,  indicator */
} TSUTRAVTIME;

/* Function prototypes */
/* complex function prototypes */
fcomplex Cadd( fcomplex a, fcomplex b );
fcomplex Csub( fcomplex a, fcomplex b );
fcomplex Cmul( fcomplex a, fcomplex b );
fcomplex Complex( double re, double im );
fcomplex Conjg( fcomplex z );
fcomplex Cdiv( fcomplex a, fcomplex b );
float    Cabs( fcomplex z );
fcomplex RCmul( double x, fcomplex a );
fcomplex Csqrt( fcomplex z );

/* diskrw function prototypes */
char    *CreateFileName( double, int, char *, char * );
int      GetNumStnsInFile( char * );
double   GetTimeFromFileName( char * );
int      ReadDiskData( char *, STATION [], int *, int, int, int, STATION [] );
int      ReadDiskDataForHypo( int, int, char *, char *, double, int, 
                              STATION [] );
int      ReadDiskDataForMTSolo( int, int, char *, char *, double, int,
                                STATION [], int * );
int      ReadDiskDataNew( char *, STATION [], int, int, int );
int      ReadDiskHeader( char *, STATION [], int );
int      ReadLineupFile( char *, STATION **, int );

#ifdef _WINNT
/* display function prototypes */
void     DisplayAlignPTime( HDC, int, int, double, int, double );
void     DisplayBodyPhase( HDC, STATION [], int, int, double, long, long, int, 
                           int, int, int, int, int, double, int, double );
void     DisplayChannel( HDC, STATION [], int, int, long, long, int,
                         int, int, int, int );
void     DisplayChannelID( HDC , STATION [], int, int, long, long,
                           int, int, int, int, int, int );
void     DisplayDCOffset( HDC , STATION [], int, int, long, long,
                          int, int, int, int, int );
void     DisplayEpiDist( HDC , STATION [], int, int, long, long,
                         int, int, int, int, int );
void     DisplayExpectedP( HDC, STATION [], int, int, int, int, int, int, int,
                           int, double, double, long );
void     DisplayKOPs( HDC, STATION [], int, int, int, int, int, 
                      long, long, int );
void     DisplayMagBox( HDC, int, int, int, int, int, int, int, double, double,
                        long, STATION [] );
void     DisplayMagCycle( HDC, int, int, int, int, int, int, int, double, 
                          double, long, STATION [], int, int, int, double );
void     DisplayMs( HDC, int, STATION [], long, long, int, int, int, int, int );
void     DisplayMwp( HDC, STATION [], int, int, int, int, int, int, int,
                     int, double, double, int, double );
void     DisplayNetwork( HDC, STATION [], int, int, long, long,
                         int, int, int, int, int );
void     DisplayOriginTime( HDC, double, double, int, int, HYPO *, int );
void     DisplayPPicks( HDC, STATION [], int, int, int, int, int, int, int,
                        int, double, double, long, int, double );
void     DisplayRTimes( HDC, int, int, int, int, int, int, int,
                        double, double, long, STATION [], int, int, double );
void     DisplayTimeLines( HDC, long, long, int, int, int, double, double,
                           double, int ); 
void     DisplayTitle( HDC, long, long, HYPO *, int, int, int ); 
void     DisplayTraces( HDC, STATION [], int, int, int, int, int, int, int, int,
                        int, int, double, double, int *, int, int, int, int, 
                        int, double );
int      ReadScreenIni( int *, int, int, char [][32], STATION [], char * );
#endif
int      struct_cmp_by_EpicentralDistance( const void *, const void * );
int      struct_cmp_by_ID( const void *, const void * );  
int      struct_cmp_by_name( const void *, const void * );
int      struct_cmp_by_StationSortIndex( const void *, const void * ); 

/* dummy function prototypes */   
int      CopyDummyData( char *, char * );
int      ReadDummyData( HYPO *, char * );
int      ReadPTimeFile( STATION *, char *, int );
int      ReadPTimeFileForMS( STATION *, char *, int );
int      WriteDummyData( HYPO *, char *, int, int );
void     WriteLPFile( int, STATION [], char *, HYPO * );
void     WriteMwFile( int, STATION [], char *, HYPO * );
void     WritePTimeFile( int, STATION [], char * );
int      WriteMwSignalFile( int, STATION [], char *, HYPO * );
void     WriteThetaFile( int, STATION [], char *, HYPO * );
			                         
/* fereg function prototypes */
char    *getnam( int, char * );                              
char    *getnamLC (int, char * );
int      getnum( double, double, char *, char * );
char    *namnum( double, double, int *, char *, char *, char * );
char    *namnumLC (double, double, int *, char *, char *, char *); 

/* filter function prototypes */
void     apiir( long *, long, int, STATION *, double [], double [], int );
void 	 apiirG( long[], long, int, long, double[], double[], int, int, int,
                 long, int, STATION * );
void     bilin2( double [], double [], int );
void     bupoles( fcomplex [], char [], int *, int );
void     FilterPacket( long *, STATION *, long, double, double, double,
                       double );
void     FilterPacketAIC( long *, long, double, double, double, double );
int      FilterPacketG( STATION *, int, double, double, double, int );
void     GetFilterCoefs( double, double, double, double [], double [], int *);
void     InitVar( STATION * );
void     lptbpa( fcomplex [], char [], int, double, 
                 double, double [], double [], int * );
void     ResetFilter( STATION * ); 
void     secord( double *, double *, double *, double, double, double, double *,
                 double *, double, double, double );
void     Taper( long *, STATION *, long, int, long );
void     TaperAIC( long *, long, int, long );
void     taperG( long[], long, long, int, long, long, int );
double   warp( double, double );
void     zero( double[], int );

/* geotools function prototypes */
void  	 ConvertLoc( LATLON *, LATLON *, char *, char * );
double   DepthCheck( double, double, char * );
void     GeoCent( LATLON * );                           
void     GeoGraphic( LATLON *, LATLON * );
void     GetDistanceAz( LATLON *, LATLON *, AZIDELT * );
void     GetLatLonTrig( LATLON * );
char    *GetOpAgency( char * );
int      GetPhaseTime( STATION *, HYPO *, PHASE *, int, char *, char * );
void     GetQuakeID( double, char * );
int      GetRegion( double, double );                     
char    *GetSeisInfo( int, double * );
int      IsItDangerous( double, double );
int      IsItInHawaii( double, double );
char    *itoaX( int, char * );
double   LLConv( char * );
void     PadBlanksL( int, char * );
void     PadZeroes( int, char * );
void     PadZeroesR( int, char * );
void     PointToEpi( LATLON *, AZIDELT *, LATLON * );
int      ReadThreatFile( char *, HYPO *, THREAT *, int *, int *, int *, int[],
                         int[], int[], int[], int[], int[] );
int      WarnOrAdvisory( double, double );

/* get_hypo function prototypes */
int      CopyHypo( HYPO *, HYPO [] );
int      HypoStruct( char *, HYPO * );
void     InitHypo( HYPO * );
void     LoadHypo( char *, HYPO [], int );	  

/* get_pick function prototypes */
long     AICPick( double[], long, double );
void     CheckForAlarm( STATION *, ALARMSTRUCT [], int, unsigned char, SHM_INFO,
                        unsigned char, unsigned char );
void     CheckForLPAlarm( STATION *, int, int NBStations[MAX_STATIONS][15],
                          unsigned char, SHM_INFO, unsigned char,
                          unsigned char );
void     CheckForLPAlarmSWD( STATION *, int, unsigned char, SHM_INFO, 
                             unsigned char, unsigned char );
void     CopyPBuf( STATION *, STATION * );
void     FindDataEndSWD( STATION * );
void     InitP( STATION * );
void     MovingAvg( long, STATION *, long, long, long );  
int      NNClassifier( double[], long );
void     PickV( STATION *, double, int, int, double, double, double, int,
                double, int, int, double, double, unsigned char, SHM_INFO,
                SHM_INFO, unsigned char, unsigned char, unsigned char, long [],
                long *, ALARMSTRUCT [], int, int, int *, int );
int      PPickStruct( char *, STATION *, unsigned char );
void     Reset( STATION * );
void     SeismicAlarm( STATION *, long, long *, int, unsigned char,
                       SHM_INFO, unsigned char, unsigned char, double, double,
                       int );
void     ShortCopyPBuf( STATION *, STATION * );
void     VeryShortCopyPBuf( STATION *, STATION * );
                                                      
/* littoral function prototypes */
void     CityAziDelt( AZIDELT *, LATLON *, CITY *, int );
int      LoadCities( CITY *, int, char *, char * );
int      LoadCitiesEC( CITY *, int, char *, char * );
void     NearestCities( LATLON *, CITY *, CITYDIS * );
void     NearestCitiesEC( LATLON *, CITY *, CITYDIS * );

/* locquake function prototypes */
void     ComputePRTimeWindows( HYPO *, STATION [], int, int );
void     ComputeTimeWindows( double, int, STATION[], int, LATLON *,
                             PHASE_WINDOWS[] );
void     FindBadPs( int, int, STATION [], HYPO *, double, double, int,
                    EQDEPTHDATA[] );
int      FindDepth( double, double, EQDEPTHDATA [], int *, int * );
void     ForceLocIasp( int, STATION [], HYPO * );
void     GetEpiAzDelta( int, STATION *, HYPO * );
void     GetLatLonTrig( LATLON * );
void     GetPTimes( int, STATION *, HYPO * );
void     InitialLocator( int, int, int, STATION *, HYPO *, double, double );
void     IsItGoodSoln( int, STATION *, HYPO *, int );
int      IsItSameQuake( HYPO *, HYPO * );
int      IsItSameQuakeTight( HYPO *, HYPO * );
int      LoadEQData( char *, int, EQDEPTHDATA * );
int      QuakeAzimuthSort( int, STATION * );
double   QuakeDeta( int, double [][4] );
void     QuakeDets( int, double [][4], double [], double [] );
void     QuakeSolveIasp( int, STATION *, HYPO *, EQDEPTHDATA[], int );
int      Round( double );

/* logquake function prototypes */
void     LogToMySQL( HYPO *, char *, char *, char *, char *, char *,
                     char *, char *, int, long, int, char *, char * );
void     MagnitudeLog( int, STATION [], HYPO * );       
void     QuakeLog( int, STATION [], HYPO *, CITY *, CITY *, char *, char *,
                   char *, char *, int, char * );
void     QuakeLog2( char *, HYPO * );

/* mags.h function prototypes */
void     AutoMwp( STATION *, double, int, int, int );
void     CalcAverageTheta( HYPO *, STATION [], int );
double   ComputeAverageMm( int, STATION [], int * );
double   ComputeAvgMm( int, MMSTUFF [], int * );
double   ComputeAvgMS( int, STATION [], int * );
void     ComputeMagnitudes( int, STATION [], HYPO * );
double   ComputeMbMag( char *, double, double, double, double, double );
void     ComputeMbMl( STATION *, int, int );
double   ComputeMlMag( char *, double, double, double, double );
double   ComputeMSMag( char *, double, double, double, double );
double   ComputeMwMag( double, double, double, double );
double   ComputeMwpMag( double, double );
int      detrend( double, double, double [], long, double, double [],
                  double [], int, double );
int      GetMbMl( STATION *, int, unsigned char, SHM_INFO,
                  unsigned char, unsigned char, int, int, double );
void     GetMDFFilt( long, long, STATION * );
void     GetNoise( long, long, STATION * );
void     GetPreferredMag( HYPO * );
int      integrate (double *, double *, int *, int *, long,
                    double [], double, double);
int      LoadBVals( char * );
long     MbMlCtsFromGM( char *, double, long, double );
double   MbMlGroundMotion( char *, double, long, long );
long     MsCtsFromGM( char *, double, long, double );
double   MsGroundMotion( char *, double, long, long );
double   Mwp_adjustment( double );
double   wavelet_decomp( double, double [], long, double );
void     ZeroMagnitudes( STATION *, int );
                             
/* mjstime function prototypes */
void     ConvertTM2ST( struct tm *, SYSTEMTIME * );
void     CopyDate (SYSTEMTIME *, SYSTEMTIME *);
int      DateToDay (SYSTEMTIME *);
void     DateToModJulianSec (SYSTEMTIME *, double *);
void     DayToDate (long, SYSTEMTIME *);
int      IsDayLightSavings (SYSTEMTIME *, int);
int      IsDayLightSavings2( double, int );
struct tm *TWCgmtime( time_t );
void     NewDateFromModSec (SYSTEMTIME *, double);
void     NewDateFromModSecRounded (SYSTEMTIME *, double);
void     MSDayToDate( long, struct tm *);

/* Mm function prototypes */
#ifdef _WINNT
void     DisplayLittoral( HDC, long, long, int, int, HYPO * );
void     DisplaySpectra( HDC, long, long, int, int, MMSTUFF [], SPECTRA [],
                         int [] );
void     DisplayTitles( HDC, long, long, int, int );
#endif
double   Mm( long, long *, MMSTUFF *, STATION *, double, double,
             SPECTRA *, double );           
int      RdRegion( char * );
                           
/* ReadStationData function prototypes */
int      IsComment( char [] );
int      LoadResponseData( STATION *, char * );
int      LoadResponseDataAll( STATION [], char *, int );
int      LoadStationData( STATION *, char * );
void     LogStaList( STATION *, int );
int      ReadStationData( char *, char *, STATION [], int );
int      ReadStationDataAlloc( char *, char *, STATION **, int, int * );
int      ReadStationList( STATION **, int *, char *, char *, char *, int, int );

/* report function prototypes */
void     ReportAlarm( STATION *, unsigned char, SHM_INFO, unsigned char,
                      unsigned char, int, char *, int );
void     ReportHypo( char *, char *, unsigned char,
                     SHM_INFO, unsigned char, unsigned char, unsigned char,
                     int, int );
void     ReportPick( STATION *, unsigned char, SHM_INFO, unsigned char,
                     unsigned char, int );
void     SendAlarm2Earthvu( char *, char *, char * );
void     SendPPick2Earthvu( char * );

/* Added new function to get meta data from MySQL database */
/* ******************************************************* */
int ReadStationList_DB( STATION **, int *,  int, char *, char *, char *, char * );
int ReadStationData_DB(  STATION [], int , char *, char *, char *, char * );
int LoadResponseData_DB( STATION *, char *, char *, char *, char * );
int LoadStationData_DB( STATION *, char *, char *, char *, char * );

#endif
