33#include <wx/datetime.h>
35#include "Station_Data.h"
37#include "TC_Error_Code.h"
38#include "TCDataSource.h"
53#define M_PI 3.141592654
56#define TIDE_MAX_DERIV (2)
60#define TIDE_TIME_PREC (15)
66#define TIDE_BLEND_TIME (3600)
72#define TIDE_TIME_STEP (TIDE_TIME_PREC)
73#define TIDE_BAD_TIME ((time_t) - 1)
96 TC_Error_Code LoadDataSources(std::vector<std::string> &sources);
97 std::vector<std::string> GetDataSet(
void) {
return m_sourcefile_array; }
99 bool IsReady(
void) {
return bTCMReady; }
101 bool GetTideOrCurrent(time_t t,
int idx,
float &value,
float &dir);
102 bool GetTideOrCurrent15(time_t t,
int idx,
float &tcvalue,
float &dir,
104 bool GetTideFlowSens(time_t t,
int sch_step,
int idx,
float &tcvalue_now,
105 float &tcvalue_prev,
bool &w_t);
106 void GetHightOrLowTide(time_t t,
int sch_step_1,
int sch_step_2,
107 float tide_val,
bool w_t,
int idx,
float &tcvalue,
110 int GetStationTimeOffset(
IDX_entry *pIDX);
111 int GetNextBigEvent(time_t *tm,
int idx);
112 std::wstring GetTidalEventStr(
int station_id, wxDateTime ref_dt,
double lat,
113 double lon,
int dt_type);
117 const IDX_entry *GetIDX_entry(
int index)
const;
119 int Get_max_IDX()
const {
return m_Combined_IDX_array.size() - 1; }
121 std::map<double, const IDX_entry *> GetStationsForLL(
double xlat,
124 int GetStationIDXbyName(
const wxString &prefix,
double xlat,
126 int GetStationIDXbyNameType(
const wxString &prefix,
double xlat,
double xlon,
128 void ScrubCurrentDepths();
139 wxString pmru_file_name;
141 ArrayOfTCDSources m_source_array;
142 std::vector<std::string> m_sourcefile_array;
144 std::vector<IDX_entry *> m_Combined_IDX_array;
150#ifndef __OAML_TIDES_H__
151#define __OAML_TIDES_H__
175#include <sys/types.h>
180#define NV_BYTE int8_t
181#define NV_INT16 int16_t
182#define NV_INT32 int32_t
183#define NV_INT64 int64_t
184#define NV_U_BYTE uint8_t
185#define NV_U_INT16 uint16_t
186#define NV_U_INT32 uint32_t
187#define NV_U_INT64 uint64_t
189#define NV_BYTE unsigned char
190#define NV_INT16 short
193#define NV_U_BYTE unsigned char
194#define NV_U_INT16 unsigned short
195#define NV_U_INT32 unsigned int
196#define NV_U_INT64 unsigned long
201#define NV_BOOL unsigned char
203#define NV_U_CHAR unsigned char
204#define NV_FLOAT32 float
205#define NV_FLOAT64 double
210#define NV_U_INT32_MAX 4294967295
211#define NV_INT32_MAX 2147483647
212#define NV_U_INT16_MAX 65535
213#define NV_INT16_MAX 32767
215#define LIBTCD_VERSION "PFM Software - libtcd v2.2.5 - 2010-08-17"
216#define LIBTCD_MAJOR_REV 2
217#define LIBTCD_MINOR_REV 2
230#undef USE_PRAGMA_MESSAGE
233#ifdef USE_PRAGMA_MESSAGE
234#pragma message("WARNING: COMPAT114 is enabled! See libtcd.html.")
236#warning COMPAT114 is enabled! See libtcd.html.
243#define ONELINER_LENGTH 90
245#define MONOLOGUE_LENGTH 10000
246#define MAX_CONSTITUENTS 255
256 NV_CHAR version[ONELINER_LENGTH];
257 NV_U_INT32 major_rev;
258 NV_U_INT32 minor_rev;
259 NV_CHAR last_modified[ONELINER_LENGTH];
260 NV_U_INT32 number_of_records;
262 NV_U_INT32 number_of_years;
263 NV_U_INT32 constituents;
264 NV_U_INT32 level_unit_types;
265 NV_U_INT32 dir_unit_types;
266 NV_U_INT32 restriction_types;
267 NV_U_INT32 datum_types;
268 NV_U_INT32 countries;
270 NV_U_INT32 legaleses;
273 NV_U_INT32 pedigree_types;
279enum TIDE_RECORD_TYPE { REFERENCE_STATION = 1, SUBORDINATE_STATION = 2 };
282 NV_INT32 record_number;
283 NV_U_INT32 record_size;
284 NV_U_BYTE record_type;
286 NV_FLOAT64 longitude;
287 NV_INT32 reference_station;
289 NV_CHAR name[ONELINER_LENGTH];
299 NV_CHAR source[ONELINER_LENGTH];
300 NV_U_BYTE restriction;
301 NV_CHAR comments[MONOLOGUE_LENGTH];
302 NV_CHAR notes[MONOLOGUE_LENGTH];
304 NV_CHAR station_id_context[ONELINER_LENGTH];
305 NV_CHAR station_id[ONELINER_LENGTH];
306 NV_U_INT32 date_imported;
307 NV_CHAR xfields[MONOLOGUE_LENGTH];
308 NV_U_BYTE direction_units;
309 NV_INT32 min_direction;
310 NV_INT32 max_direction;
311 NV_U_BYTE level_units;
314 NV_FLOAT32 datum_offset;
316 NV_INT32 zone_offset;
317 NV_U_INT32 expiration_date;
318 NV_U_INT16 months_on_station;
319 NV_U_INT32 last_date_on_station;
320 NV_U_BYTE confidence;
321 NV_FLOAT32 amplitude[MAX_CONSTITUENTS];
322 NV_FLOAT32 epoch[MAX_CONSTITUENTS];
325 NV_INT32 min_time_add;
326 NV_FLOAT32 min_level_add;
327 NV_FLOAT32 min_level_multiply;
328 NV_INT32 max_time_add;
329 NV_FLOAT32 max_level_add;
330 NV_FLOAT32 max_level_multiply;
331 NV_INT32 flood_begins;
338 NV_U_BYTE avg_level_units;
339 NV_FLOAT32 min_avg_level;
340 NV_FLOAT32 max_avg_level;
351#define NULLSLACKOFFSET 0xA00
355#define AMPLITUDE_EPSILON 0.00005
358#define NINT(a) ((a) < 0.0 ? (NV_INT32)((a) - 0.5) : (NV_INT32)((a) + 0.5))
370const NV_CHAR *get_country(NV_INT32 num);
371const NV_CHAR *get_tzfile(NV_INT32 num);
372const NV_CHAR *get_level_units(NV_INT32 num);
373const NV_CHAR *get_dir_units(NV_INT32 num);
374const NV_CHAR *get_restriction(NV_INT32 num);
375const NV_CHAR *get_datum(NV_INT32 num);
376const NV_CHAR *get_legalese(NV_INT32 num);
381const NV_CHAR *get_constituent(NV_INT32 num);
386const NV_CHAR *get_station(NV_INT32 num);
390NV_FLOAT64 get_speed(NV_INT32 num);
395NV_FLOAT32 get_equilibrium(NV_INT32 num, NV_INT32 year);
396NV_FLOAT32 get_node_factor(NV_INT32 num, NV_INT32 year);
403NV_FLOAT32 *get_equilibriums(NV_INT32 num);
404NV_FLOAT32 *get_node_factors(NV_INT32 num);
411NV_INT32 get_time(
const NV_CHAR *
string);
412NV_CHAR *ret_time(NV_INT32 time);
413NV_CHAR *ret_time_neat(NV_INT32 time);
420NV_CHAR *ret_date(NV_U_INT32 date);
426NV_INT32 search_station(
const NV_CHAR *
string);
430NV_INT32 find_station(
const NV_CHAR *name);
431NV_INT32 find_tzfile(
const NV_CHAR *name);
432NV_INT32 find_country(
const NV_CHAR *name);
433NV_INT32 find_level_units(
const NV_CHAR *name);
434NV_INT32 find_dir_units(
const NV_CHAR *name);
435NV_INT32 find_restriction(
const NV_CHAR *name);
436NV_INT32 find_datum(
const NV_CHAR *name);
437NV_INT32 find_constituent(
const NV_CHAR *name);
438NV_INT32 find_legalese(
const NV_CHAR *name);
461void set_speed(NV_INT32 num, NV_FLOAT64 value);
466void set_equilibrium(NV_INT32 num, NV_INT32 year, NV_FLOAT32 value);
467void set_node_factor(NV_INT32 num, NV_INT32 year, NV_FLOAT32 value);
473NV_BOOL open_tide_db(
const NV_CHAR *file);
481NV_BOOL create_tide_db(
const NV_CHAR *file, NV_U_INT32 constituents,
482 NV_CHAR
const *
const constituent[],
483 const NV_FLOAT64 *speed, NV_INT32 start_year,
484 NV_U_INT32 num_years,
485 NV_FLOAT32
const *
const equilibrium[],
486 NV_FLOAT32
const *
const node_factor[]);
504NV_INT32 get_nearest_partial_tide_record(NV_FLOAT64 lat, NV_FLOAT64 lon,
511NV_INT32 read_tide_record(NV_INT32 num,
TIDE_RECORD *rec);
527NV_BOOL update_tide_record(NV_INT32 num,
TIDE_RECORD *rec);
529NV_BOOL update_tide_record(NV_INT32 num,
TIDE_RECORD *rec,
544#define NAME_LENGTH ONELINER_LENGTH
545#define SOURCE_LENGTH ONELINER_LENGTH
546#define COMMENTS_LENGTH MONOLOGUE_LENGTH
547NV_CHAR *get_pedigree(NV_INT32 num);
548NV_INT32 find_pedigree(
const NV_CHAR *name);
605#ifndef __TIDE_HEADER__
606#define __TIDE_HEADER__
611 const NV_CHAR *keyphrase;
612 const NV_CHAR *datatype;
622 NV_CHAR **constituent;
624 NV_FLOAT32 **equilibrium;
625 NV_FLOAT32 **node_factor;
626 NV_CHAR **level_unit;
628 NV_CHAR **restriction;
633 NV_U_INT32 header_size;
634 NV_U_INT32 speed_bits;
635 NV_U_INT32 speed_scale;
636 NV_INT32 speed_offset;
637 NV_U_INT32 equilibrium_bits;
638 NV_U_INT32 equilibrium_scale;
639 NV_INT32 equilibrium_offset;
640 NV_U_INT32 node_bits;
641 NV_U_INT32 node_scale;
642 NV_INT32 node_offset;
643 NV_U_INT32 amplitude_bits;
644 NV_U_INT32 amplitude_scale;
645 NV_U_INT32 epoch_bits;
646 NV_U_INT32 epoch_scale;
647 NV_U_INT32 constituent_bits;
648 NV_U_INT32 record_type_bits;
649 NV_U_INT32 latitude_scale;
650 NV_U_INT32 latitude_bits;
651 NV_U_INT32 longitude_scale;
652 NV_U_INT32 longitude_bits;
653 NV_U_INT32 record_size_bits;
654 NV_U_INT32 station_bits;
655 NV_U_INT32 level_unit_bits;
656 NV_U_INT32 dir_unit_bits;
657 NV_U_INT32 restriction_bits;
658 NV_U_INT32 max_restriction_types;
659 NV_U_INT32 tzfile_bits;
660 NV_U_INT32 max_tzfiles;
661 NV_U_INT32 country_bits;
662 NV_U_INT32 max_countries;
663 NV_U_INT32 datum_bits;
664 NV_U_INT32 max_datum_types;
665 NV_U_INT32 legalese_bits;
666 NV_U_INT32 max_legaleses;
667 NV_U_INT32 datum_offset_bits;
668 NV_U_INT32 datum_offset_scale;
669 NV_U_INT32 date_bits;
670 NV_U_INT32 months_on_station_bits;
671 NV_U_INT32 confidence_value_bits;
672 NV_U_INT32 time_bits;
673 NV_U_INT32 level_add_bits;
674 NV_U_INT32 level_add_scale;
675 NV_U_INT32 level_multiply_bits;
676 NV_U_INT32 level_multiply_scale;
677 NV_U_INT32 direction_bits;
678 NV_U_INT32 constituent_size;
679 NV_U_INT32 level_unit_size;
680 NV_U_INT32 dir_unit_size;
681 NV_U_INT32 restriction_size;
682 NV_U_INT32 tzfile_size;
683 NV_U_INT32 country_size;
684 NV_U_INT32 datum_size;
685 NV_U_INT32 legalese_size;
686 NV_U_INT32 end_of_file;
689 NV_U_INT32 pedigree_bits;
690 NV_U_INT32 pedigree_size;
700static KEY keys[] = {{
703 {(NV_CHAR *)hd.pub.version},
708 {(NV_CHAR *)&hd.pub.major_rev},
713 {(NV_CHAR *)&hd.pub.minor_rev},
718 {(NV_CHAR *)hd.pub.last_modified},
721 "[NUMBER OF RECORDS]",
723 {(NV_CHAR *)&hd.pub.number_of_records},
728 {(NV_CHAR *)&hd.pub.start_year},
733 {(NV_CHAR *)&hd.pub.number_of_years},
738 {(NV_CHAR *)&hd.pub.constituents},
741 "[LEVEL UNIT TYPES]",
743 {(NV_CHAR *)&hd.pub.level_unit_types},
746 "[DIRECTION UNIT TYPES]",
748 {(NV_CHAR *)&hd.pub.dir_unit_types},
751 "[RESTRICTION TYPES]",
753 {(NV_CHAR *)&hd.pub.restriction_types},
758 {(NV_CHAR *)&hd.pub.pedigree_types},
763 {(NV_CHAR *)&hd.pub.tzfiles},
768 {(NV_CHAR *)&hd.pub.countries},
773 {(NV_CHAR *)&hd.pub.datum_types},
778 {(NV_CHAR *)&hd.pub.legaleses},
783 {(NV_CHAR *)&hd.header_size},
788 {(NV_CHAR *)&hd.speed_bits},
793 {(NV_CHAR *)&hd.speed_scale},
798 {(NV_CHAR *)&hd.speed_offset},
801 "[EQUILIBRIUM BITS]",
803 {(NV_CHAR *)&hd.equilibrium_bits},
806 "[EQUILIBRIUM SCALE]",
808 {(NV_CHAR *)&hd.equilibrium_scale},
811 "[EQUILIBRIUM OFFSET]",
813 {(NV_CHAR *)&hd.equilibrium_offset},
818 {(NV_CHAR *)&hd.node_bits},
823 {(NV_CHAR *)&hd.node_scale},
828 {(NV_CHAR *)&hd.node_offset},
833 {(NV_CHAR *)&hd.amplitude_bits},
838 {(NV_CHAR *)&hd.amplitude_scale},
843 {(NV_CHAR *)&hd.epoch_bits},
848 {(NV_CHAR *)&hd.epoch_scale},
851 "[CONSTITUENT BITS]",
853 {(NV_CHAR *)&hd.constituent_bits},
858 {(NV_CHAR *)&hd.level_unit_bits},
861 "[DIRECTION UNIT BITS]",
863 {(NV_CHAR *)&hd.dir_unit_bits},
866 "[RESTRICTION BITS]",
868 {(NV_CHAR *)&hd.restriction_bits},
873 {(NV_CHAR *)&hd.pedigree_bits},
878 {(NV_CHAR *)&hd.tzfile_bits},
883 {(NV_CHAR *)&hd.country_bits},
888 {(NV_CHAR *)&hd.datum_bits},
893 {(NV_CHAR *)&hd.legalese_bits},
896 "[RECORD TYPE BITS]",
898 {(NV_CHAR *)&hd.record_type_bits},
903 {(NV_CHAR *)&hd.latitude_scale},
908 {(NV_CHAR *)&hd.latitude_bits},
913 {(NV_CHAR *)&hd.longitude_scale},
918 {(NV_CHAR *)&hd.longitude_bits},
921 "[RECORD SIZE BITS]",
923 {(NV_CHAR *)&hd.record_size_bits},
928 {(NV_CHAR *)&hd.station_bits},
931 "[DATUM OFFSET BITS]",
933 {(NV_CHAR *)&hd.datum_offset_bits},
936 "[DATUM OFFSET SCALE]",
938 {(NV_CHAR *)&hd.datum_offset_scale},
943 {(NV_CHAR *)&hd.date_bits},
946 "[MONTHS ON STATION BITS]",
948 {(NV_CHAR *)&hd.months_on_station_bits},
951 "[CONFIDENCE VALUE BITS]",
953 {(NV_CHAR *)&hd.confidence_value_bits},
958 {(NV_CHAR *)&hd.time_bits},
963 {(NV_CHAR *)&hd.level_add_bits},
968 {(NV_CHAR *)&hd.level_add_scale},
971 "[LEVEL MULTIPLY BITS]",
973 {(NV_CHAR *)&hd.level_multiply_bits},
976 "[LEVEL MULTIPLY SCALE]",
978 {(NV_CHAR *)&hd.level_multiply_scale},
983 {(NV_CHAR *)&hd.direction_bits},
986 "[CONSTITUENT SIZE]",
988 {(NV_CHAR *)&hd.constituent_size},
993 {(NV_CHAR *)&hd.level_unit_size},
996 "[DIRECTION UNIT SIZE]",
998 {(NV_CHAR *)&hd.dir_unit_size},
1001 "[RESTRICTION SIZE]",
1003 {(NV_CHAR *)&hd.restriction_size},
1008 {(NV_CHAR *)&hd.pedigree_size},
1013 {(NV_CHAR *)&hd.tzfile_size},
1018 {(NV_CHAR *)&hd.country_size},
1023 {(NV_CHAR *)&hd.datum_size},
1028 {(NV_CHAR *)&hd.legalese_size},
1033 {(NV_CHAR *)&hd.end_of_file},
Represents an index entry for tidal and current data.