forked from geodynamics/hc
-
Notifications
You must be signed in to change notification settings - Fork 0
/
ggrd_grdtrack_util.h
156 lines (122 loc) · 4.6 KB
/
ggrd_grdtrack_util.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
/*
header files for modified GMT grd interpolation routines dealing
with grd interpolation
*/
#ifndef __GGRD_READ_GMT__
#include "gmt.h"
void GMT_grdio_init (void);
#define __GGRD_READ_GMT__
#endif
#include "ggrd_base.h"
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#ifndef USE_GMT3 /* new GMT */
#define GGRD_GMT_GLOBAL_STRING "-fg"
#define GGRD_GMT_GEOGRAPHIC_STRING "-fg"
#define GGRD_GMT_XPERIODIC_STRING "-Lx"
#else /* old GMT */
#define GGRD_GMT_GLOBAL_STRING "-Lg"
#define GGRD_GMT_GEOGRAPHIC_STRING ""
#define GGRD_GMT_XPERIODIC_STRING ""
#endif
/*
wrappers
*/
int ggrd_grdtrack_init_general(ggrd_boolean ,char *, char *,char *,
struct ggrd_gt *,ggrd_boolean ,
ggrd_boolean, ggrd_boolean);
ggrd_boolean ggrd_grdtrack_interpolate_lonlatz(double ,double ,double ,
struct ggrd_gt *,double *,
ggrd_boolean );
ggrd_boolean ggrd_grdtrack_interpolate_rtp(double ,double ,double ,
struct ggrd_gt *,double *,
ggrd_boolean,ggrd_boolean,double);
ggrd_boolean ggrd_grdtrack_interpolate_xyz(double ,double ,double ,
struct ggrd_gt *,double *,
ggrd_boolean);
ggrd_boolean ggrd_grdtrack_interpolate_xy(double ,double ,
struct ggrd_gt *,
double *,
ggrd_boolean );
ggrd_boolean ggrd_grdtrack_interpolate_tp(double ,double ,
struct ggrd_gt *,
double *,
ggrd_boolean ,ggrd_boolean);
void ggrd_grdtrack_free_gstruc(struct ggrd_gt *);
int ggrd_grdtrack_rescale(struct ggrd_gt *,ggrd_boolean , ggrd_boolean ,
ggrd_boolean ,double);
void ggrd_init_vstruc(struct ggrd_master *);
/*
moderately external
*/
void ggrd_init_master(struct ggrd_master *);
int ggrd_init_thist_from_file(struct ggrd_t *,char *,ggrd_boolean ,ggrd_boolean);
int ggrd_read_vel_grids(struct ggrd_master *, double, unsigned short, unsigned short, char *,ggrd_boolean);
#ifndef USE_GMT3
/* GMT >4.1.2 */
ggrd_boolean ggrd_grdtrack_interpolate(double *, ggrd_boolean , struct GRD_HEADER *, float *,
struct GMT_EDGEINFO *, int, float *, int , double *,ggrd_boolean,
struct GMT_BCR *);
// GMT < 4.5.1
//int ggrd_grdtrack_init(double *, double *, double *, double *, float **, int *, char *, struct GRD_HEADER **, struct GMT_EDGEINFO **, char *, ggrd_boolean *, int *, ggrd_boolean, char *, float **, int *, ggrd_boolean, ggrd_boolean, ggrd_boolean, struct GMT_BCR *);
// GMT >= 4.5.1
int ggrd_grdtrack_init(double *, double *, double *, double *, float **, int *, char *,
struct GRD_HEADER **, struct GMT_EDGEINFO **, char *,
ggrd_boolean *, GMT_LONG *, ggrd_boolean, char *,
float **, int *,GMT_LONG, ggrd_boolean,
ggrd_boolean, struct GMT_BCR *);
/* < 4.5.1 */
//void ggrd_print_layer_avg(float *,float *,int , int ,int, FILE *,int *);
/* >= 4.5.1 */
void ggrd_print_layer_avg(float *,float *,int , int ,int, FILE *,GMT_LONG *);
#else
void ggrd_print_layer_avg(float *,float *,int , int ,int, FILE *,int *);
ggrd_boolean ggrd_grdtrack_interpolate(double *, ggrd_boolean , struct GRD_HEADER *, float *,
struct GMT_EDGEINFO *, int,
float *, int ,
double *,ggrd_boolean,
struct BCR *);
int ggrd_grdtrack_init(double *, double *,double *, double *, /* geographic bounds,
set all to zero to
get the whole range from the
input grid files
*/
float **, /* data */
int *, /* size of data */
char *, /* name, or prefix, of grd file with scalars */
struct GRD_HEADER **,
struct GMT_EDGEINFO **,
char *,ggrd_boolean *,
int *, /* [4] array with padding (output) */
ggrd_boolean _d, char *, /* depth file name */
float **, /* layers, pass as NULL */
int *, /* number of layers */
ggrd_boolean , /* linear/cubic? */
ggrd_boolean ,ggrd_boolean,
struct BCR *);
void ggrd_global_bcr_assign(struct BCR *);
void my_GMT_bcr_init (struct GRD_HEADER *, int *,
int ,struct BCR *);
#endif
/*
local
*/
void ggrd_gt_interpolate_z(double,float *,int ,
int *, int *, double *, double *,ggrd_boolean,
ggrd_boolean *); /* */
float ggrd_gt_rms(float *,int );
float ggrd_gt_mean(float *,int );
FILE *ggrd_open(char *, char *, char *);
void ggrd_find_spherical_vel_from_rigid_cart_rot(double *,
double *,
double *,
double *,
double *);
void ggrd_vecalloc(double **,int,char *);
void ggrd_vecrealloc(double **,int,char *);
void ggrd_indexx(int ,GGRD_CPREC *, int *);
void ggrd_calc_mean_and_stddev(GGRD_CPREC *, GGRD_CPREC *,int ,GGRD_CPREC *,
GGRD_CPREC *,GGRD_CPREC *,
ggrd_boolean , ggrd_boolean,GGRD_CPREC *);