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)
93 TC_Error_Code LoadDataSources(std::vector<std::string> &sources);
94 std::vector<std::string> GetDataSet(
void) {
return m_sourcefile_array; }
96 bool IsReady(
void) {
return bTCMReady; }
98 bool GetTideOrCurrent(time_t t,
int idx,
float &value,
float &dir);
99 bool GetTideOrCurrent15(time_t t,
int idx,
float &tcvalue,
float &dir,
101 bool GetTideFlowSens(time_t t,
int sch_step,
int idx,
float &tcvalue_now,
102 float &tcvalue_prev,
bool &w_t);
103 void GetHightOrLowTide(time_t t,
int sch_step_1,
int sch_step_2,
104 float tide_val,
bool w_t,
int idx,
float &tcvalue,
107 int GetStationTimeOffset(
IDX_entry *pIDX);
108 int GetNextBigEvent(time_t *tm,
int idx);
109 std::wstring GetTidalEventStr(
int station_id, wxDateTime ref_dt,
double lat,
110 double lon,
int dt_type);
114 const IDX_entry *GetIDX_entry(
int index)
const;
116 int Get_max_IDX()
const {
return m_Combined_IDX_array.size() - 1; }
118 std::map<double, const IDX_entry *> GetStationsForLL(
double xlat,
121 int GetStationIDXbyName(
const wxString &prefix,
double xlat,
123 int GetStationIDXbyNameType(
const wxString &prefix,
double xlat,
double xlon,
125 void ScrubCurrentDepths();
136 wxString pmru_file_name;
138 ArrayOfTCDSources m_source_array;
139 std::vector<std::string> m_sourcefile_array;
141 std::vector<IDX_entry *> m_Combined_IDX_array;
147#ifndef __OAML_TIDES_H__
148#define __OAML_TIDES_H__
172#include <sys/types.h>
177#define NV_BYTE int8_t
178#define NV_INT16 int16_t
179#define NV_INT32 int32_t
180#define NV_INT64 int64_t
181#define NV_U_BYTE uint8_t
182#define NV_U_INT16 uint16_t
183#define NV_U_INT32 uint32_t
184#define NV_U_INT64 uint64_t
186#define NV_BYTE unsigned char
187#define NV_INT16 short
190#define NV_U_BYTE unsigned char
191#define NV_U_INT16 unsigned short
192#define NV_U_INT32 unsigned int
193#define NV_U_INT64 unsigned long
198#define NV_BOOL unsigned char
200#define NV_U_CHAR unsigned char
201#define NV_FLOAT32 float
202#define NV_FLOAT64 double
207#define NV_U_INT32_MAX 4294967295
208#define NV_INT32_MAX 2147483647
209#define NV_U_INT16_MAX 65535
210#define NV_INT16_MAX 32767
212#define LIBTCD_VERSION "PFM Software - libtcd v2.2.5 - 2010-08-17"
213#define LIBTCD_MAJOR_REV 2
214#define LIBTCD_MINOR_REV 2
227#undef USE_PRAGMA_MESSAGE
230#ifdef USE_PRAGMA_MESSAGE
231#pragma message("WARNING: COMPAT114 is enabled! See libtcd.html.")
233#warning COMPAT114 is enabled! See libtcd.html.
240#define ONELINER_LENGTH 90
242#define MONOLOGUE_LENGTH 10000
243#define MAX_CONSTITUENTS 255
253 NV_CHAR version[ONELINER_LENGTH];
254 NV_U_INT32 major_rev;
255 NV_U_INT32 minor_rev;
256 NV_CHAR last_modified[ONELINER_LENGTH];
257 NV_U_INT32 number_of_records;
259 NV_U_INT32 number_of_years;
260 NV_U_INT32 constituents;
261 NV_U_INT32 level_unit_types;
262 NV_U_INT32 dir_unit_types;
263 NV_U_INT32 restriction_types;
264 NV_U_INT32 datum_types;
265 NV_U_INT32 countries;
267 NV_U_INT32 legaleses;
270 NV_U_INT32 pedigree_types;
276enum TIDE_RECORD_TYPE { REFERENCE_STATION = 1, SUBORDINATE_STATION = 2 };
279 NV_INT32 record_number;
280 NV_U_INT32 record_size;
281 NV_U_BYTE record_type;
283 NV_FLOAT64 longitude;
284 NV_INT32 reference_station;
286 NV_CHAR name[ONELINER_LENGTH];
296 NV_CHAR source[ONELINER_LENGTH];
297 NV_U_BYTE restriction;
298 NV_CHAR comments[MONOLOGUE_LENGTH];
299 NV_CHAR notes[MONOLOGUE_LENGTH];
301 NV_CHAR station_id_context[ONELINER_LENGTH];
302 NV_CHAR station_id[ONELINER_LENGTH];
303 NV_U_INT32 date_imported;
304 NV_CHAR xfields[MONOLOGUE_LENGTH];
305 NV_U_BYTE direction_units;
306 NV_INT32 min_direction;
307 NV_INT32 max_direction;
308 NV_U_BYTE level_units;
311 NV_FLOAT32 datum_offset;
313 NV_INT32 zone_offset;
314 NV_U_INT32 expiration_date;
315 NV_U_INT16 months_on_station;
316 NV_U_INT32 last_date_on_station;
317 NV_U_BYTE confidence;
318 NV_FLOAT32 amplitude[MAX_CONSTITUENTS];
319 NV_FLOAT32 epoch[MAX_CONSTITUENTS];
322 NV_INT32 min_time_add;
323 NV_FLOAT32 min_level_add;
324 NV_FLOAT32 min_level_multiply;
325 NV_INT32 max_time_add;
326 NV_FLOAT32 max_level_add;
327 NV_FLOAT32 max_level_multiply;
328 NV_INT32 flood_begins;
335 NV_U_BYTE avg_level_units;
336 NV_FLOAT32 min_avg_level;
337 NV_FLOAT32 max_avg_level;
348#define NULLSLACKOFFSET 0xA00
352#define AMPLITUDE_EPSILON 0.00005
355#define NINT(a) ((a) < 0.0 ? (NV_INT32)((a) - 0.5) : (NV_INT32)((a) + 0.5))
367const NV_CHAR *get_country(NV_INT32 num);
368const NV_CHAR *get_tzfile(NV_INT32 num);
369const NV_CHAR *get_level_units(NV_INT32 num);
370const NV_CHAR *get_dir_units(NV_INT32 num);
371const NV_CHAR *get_restriction(NV_INT32 num);
372const NV_CHAR *get_datum(NV_INT32 num);
373const NV_CHAR *get_legalese(NV_INT32 num);
378const NV_CHAR *get_constituent(NV_INT32 num);
383const NV_CHAR *get_station(NV_INT32 num);
387NV_FLOAT64 get_speed(NV_INT32 num);
392NV_FLOAT32 get_equilibrium(NV_INT32 num, NV_INT32 year);
393NV_FLOAT32 get_node_factor(NV_INT32 num, NV_INT32 year);
400NV_FLOAT32 *get_equilibriums(NV_INT32 num);
401NV_FLOAT32 *get_node_factors(NV_INT32 num);
408NV_INT32 get_time(
const NV_CHAR *
string);
409NV_CHAR *ret_time(NV_INT32 time);
410NV_CHAR *ret_time_neat(NV_INT32 time);
417NV_CHAR *ret_date(NV_U_INT32 date);
423NV_INT32 search_station(
const NV_CHAR *
string);
427NV_INT32 find_station(
const NV_CHAR *name);
428NV_INT32 find_tzfile(
const NV_CHAR *name);
429NV_INT32 find_country(
const NV_CHAR *name);
430NV_INT32 find_level_units(
const NV_CHAR *name);
431NV_INT32 find_dir_units(
const NV_CHAR *name);
432NV_INT32 find_restriction(
const NV_CHAR *name);
433NV_INT32 find_datum(
const NV_CHAR *name);
434NV_INT32 find_constituent(
const NV_CHAR *name);
435NV_INT32 find_legalese(
const NV_CHAR *name);
458void set_speed(NV_INT32 num, NV_FLOAT64 value);
463void set_equilibrium(NV_INT32 num, NV_INT32 year, NV_FLOAT32 value);
464void set_node_factor(NV_INT32 num, NV_INT32 year, NV_FLOAT32 value);
470NV_BOOL open_tide_db(
const NV_CHAR *file);
478NV_BOOL create_tide_db(
const NV_CHAR *file, NV_U_INT32 constituents,
479 NV_CHAR
const *
const constituent[],
480 const NV_FLOAT64 *speed, NV_INT32 start_year,
481 NV_U_INT32 num_years,
482 NV_FLOAT32
const *
const equilibrium[],
483 NV_FLOAT32
const *
const node_factor[]);
501NV_INT32 get_nearest_partial_tide_record(NV_FLOAT64 lat, NV_FLOAT64 lon,
508NV_INT32 read_tide_record(NV_INT32 num,
TIDE_RECORD *rec);
524NV_BOOL update_tide_record(NV_INT32 num,
TIDE_RECORD *rec);
526NV_BOOL update_tide_record(NV_INT32 num,
TIDE_RECORD *rec,
541#define NAME_LENGTH ONELINER_LENGTH
542#define SOURCE_LENGTH ONELINER_LENGTH
543#define COMMENTS_LENGTH MONOLOGUE_LENGTH
544NV_CHAR *get_pedigree(NV_INT32 num);
545NV_INT32 find_pedigree(
const NV_CHAR *name);
602#ifndef __TIDE_HEADER__
603#define __TIDE_HEADER__
608 const NV_CHAR *keyphrase;
609 const NV_CHAR *datatype;
619 NV_CHAR **constituent;
621 NV_FLOAT32 **equilibrium;
622 NV_FLOAT32 **node_factor;
623 NV_CHAR **level_unit;
625 NV_CHAR **restriction;
630 NV_U_INT32 header_size;
631 NV_U_INT32 speed_bits;
632 NV_U_INT32 speed_scale;
633 NV_INT32 speed_offset;
634 NV_U_INT32 equilibrium_bits;
635 NV_U_INT32 equilibrium_scale;
636 NV_INT32 equilibrium_offset;
637 NV_U_INT32 node_bits;
638 NV_U_INT32 node_scale;
639 NV_INT32 node_offset;
640 NV_U_INT32 amplitude_bits;
641 NV_U_INT32 amplitude_scale;
642 NV_U_INT32 epoch_bits;
643 NV_U_INT32 epoch_scale;
644 NV_U_INT32 constituent_bits;
645 NV_U_INT32 record_type_bits;
646 NV_U_INT32 latitude_scale;
647 NV_U_INT32 latitude_bits;
648 NV_U_INT32 longitude_scale;
649 NV_U_INT32 longitude_bits;
650 NV_U_INT32 record_size_bits;
651 NV_U_INT32 station_bits;
652 NV_U_INT32 level_unit_bits;
653 NV_U_INT32 dir_unit_bits;
654 NV_U_INT32 restriction_bits;
655 NV_U_INT32 max_restriction_types;
656 NV_U_INT32 tzfile_bits;
657 NV_U_INT32 max_tzfiles;
658 NV_U_INT32 country_bits;
659 NV_U_INT32 max_countries;
660 NV_U_INT32 datum_bits;
661 NV_U_INT32 max_datum_types;
662 NV_U_INT32 legalese_bits;
663 NV_U_INT32 max_legaleses;
664 NV_U_INT32 datum_offset_bits;
665 NV_U_INT32 datum_offset_scale;
666 NV_U_INT32 date_bits;
667 NV_U_INT32 months_on_station_bits;
668 NV_U_INT32 confidence_value_bits;
669 NV_U_INT32 time_bits;
670 NV_U_INT32 level_add_bits;
671 NV_U_INT32 level_add_scale;
672 NV_U_INT32 level_multiply_bits;
673 NV_U_INT32 level_multiply_scale;
674 NV_U_INT32 direction_bits;
675 NV_U_INT32 constituent_size;
676 NV_U_INT32 level_unit_size;
677 NV_U_INT32 dir_unit_size;
678 NV_U_INT32 restriction_size;
679 NV_U_INT32 tzfile_size;
680 NV_U_INT32 country_size;
681 NV_U_INT32 datum_size;
682 NV_U_INT32 legalese_size;
683 NV_U_INT32 end_of_file;
686 NV_U_INT32 pedigree_bits;
687 NV_U_INT32 pedigree_size;
697static KEY keys[] = {{
700 {(NV_CHAR *)hd.pub.version},
705 {(NV_CHAR *)&hd.pub.major_rev},
710 {(NV_CHAR *)&hd.pub.minor_rev},
715 {(NV_CHAR *)hd.pub.last_modified},
718 "[NUMBER OF RECORDS]",
720 {(NV_CHAR *)&hd.pub.number_of_records},
725 {(NV_CHAR *)&hd.pub.start_year},
730 {(NV_CHAR *)&hd.pub.number_of_years},
735 {(NV_CHAR *)&hd.pub.constituents},
738 "[LEVEL UNIT TYPES]",
740 {(NV_CHAR *)&hd.pub.level_unit_types},
743 "[DIRECTION UNIT TYPES]",
745 {(NV_CHAR *)&hd.pub.dir_unit_types},
748 "[RESTRICTION TYPES]",
750 {(NV_CHAR *)&hd.pub.restriction_types},
755 {(NV_CHAR *)&hd.pub.pedigree_types},
760 {(NV_CHAR *)&hd.pub.tzfiles},
765 {(NV_CHAR *)&hd.pub.countries},
770 {(NV_CHAR *)&hd.pub.datum_types},
775 {(NV_CHAR *)&hd.pub.legaleses},
780 {(NV_CHAR *)&hd.header_size},
785 {(NV_CHAR *)&hd.speed_bits},
790 {(NV_CHAR *)&hd.speed_scale},
795 {(NV_CHAR *)&hd.speed_offset},
798 "[EQUILIBRIUM BITS]",
800 {(NV_CHAR *)&hd.equilibrium_bits},
803 "[EQUILIBRIUM SCALE]",
805 {(NV_CHAR *)&hd.equilibrium_scale},
808 "[EQUILIBRIUM OFFSET]",
810 {(NV_CHAR *)&hd.equilibrium_offset},
815 {(NV_CHAR *)&hd.node_bits},
820 {(NV_CHAR *)&hd.node_scale},
825 {(NV_CHAR *)&hd.node_offset},
830 {(NV_CHAR *)&hd.amplitude_bits},
835 {(NV_CHAR *)&hd.amplitude_scale},
840 {(NV_CHAR *)&hd.epoch_bits},
845 {(NV_CHAR *)&hd.epoch_scale},
848 "[CONSTITUENT BITS]",
850 {(NV_CHAR *)&hd.constituent_bits},
855 {(NV_CHAR *)&hd.level_unit_bits},
858 "[DIRECTION UNIT BITS]",
860 {(NV_CHAR *)&hd.dir_unit_bits},
863 "[RESTRICTION BITS]",
865 {(NV_CHAR *)&hd.restriction_bits},
870 {(NV_CHAR *)&hd.pedigree_bits},
875 {(NV_CHAR *)&hd.tzfile_bits},
880 {(NV_CHAR *)&hd.country_bits},
885 {(NV_CHAR *)&hd.datum_bits},
890 {(NV_CHAR *)&hd.legalese_bits},
893 "[RECORD TYPE BITS]",
895 {(NV_CHAR *)&hd.record_type_bits},
900 {(NV_CHAR *)&hd.latitude_scale},
905 {(NV_CHAR *)&hd.latitude_bits},
910 {(NV_CHAR *)&hd.longitude_scale},
915 {(NV_CHAR *)&hd.longitude_bits},
918 "[RECORD SIZE BITS]",
920 {(NV_CHAR *)&hd.record_size_bits},
925 {(NV_CHAR *)&hd.station_bits},
928 "[DATUM OFFSET BITS]",
930 {(NV_CHAR *)&hd.datum_offset_bits},
933 "[DATUM OFFSET SCALE]",
935 {(NV_CHAR *)&hd.datum_offset_scale},
940 {(NV_CHAR *)&hd.date_bits},
943 "[MONTHS ON STATION BITS]",
945 {(NV_CHAR *)&hd.months_on_station_bits},
948 "[CONFIDENCE VALUE BITS]",
950 {(NV_CHAR *)&hd.confidence_value_bits},
955 {(NV_CHAR *)&hd.time_bits},
960 {(NV_CHAR *)&hd.level_add_bits},
965 {(NV_CHAR *)&hd.level_add_scale},
968 "[LEVEL MULTIPLY BITS]",
970 {(NV_CHAR *)&hd.level_multiply_bits},
973 "[LEVEL MULTIPLY SCALE]",
975 {(NV_CHAR *)&hd.level_multiply_scale},
980 {(NV_CHAR *)&hd.direction_bits},
983 "[CONSTITUENT SIZE]",
985 {(NV_CHAR *)&hd.constituent_size},
990 {(NV_CHAR *)&hd.level_unit_size},
993 "[DIRECTION UNIT SIZE]",
995 {(NV_CHAR *)&hd.dir_unit_size},
998 "[RESTRICTION SIZE]",
1000 {(NV_CHAR *)&hd.restriction_size},
1005 {(NV_CHAR *)&hd.pedigree_size},
1010 {(NV_CHAR *)&hd.tzfile_size},
1015 {(NV_CHAR *)&hd.country_size},
1020 {(NV_CHAR *)&hd.datum_size},
1025 {(NV_CHAR *)&hd.legalese_size},
1030 {(NV_CHAR *)&hd.end_of_file},
Represents an index entry for tidal and current data.