mherman09 / hdef Goto Github PK

Currently, distaz2lola and lola2distaz require a command line argument to read standard input. It would be better if they could automatically check for standard input.

Polyfit_special is a subprogram that windows data before fitting it. It should be an option in polyfit, not a stand alone program. Merge it into polyfit.

Please note that your utmgeo() routine provided to convert long/lat to UTM coordinates contains two 2 mistakes:

source file: Hdef/src/map_module.f90

line 147:

wrong: f3 = 15.d0*e4/256.d0*45.d0*e6/1024.d0
correct would be: f3 = 15.d0*e4/256.d0 + 45.d0*e6/1024.d0

line 209:

wrong: f3 = 5.d0*3.d0*t1 + 10.d0*c1 - 4.d0*c1**2 - 9.d0*ep2
correct would be: f3 = 5.d0 + 3.d0*t1 + 10.d0*c1 - 4.d0*c1**2 - 9.d0*ep2

the above corrections could be verified by looking at the original publication:
Snyder, J.P., Map projections - a working manual, USGS professional paper 1395, 1987
https://doi.org/10.3133/pp1395

line 147 mistake/correction: see page 61, eq. 3-21 for M
line 209 mistake/correction: see page 63, eq. 8-17 for phi

the above error will lead to a positional error of ~7 m when converting lon/lat to UTM.

best regards,
daniel

Can I use this code to modeeling afterslip in time series?

Implement multiple empirical relations into a new program. Already partially in o92util, but not easy to use independently of computing deformation.

Add options in fltinv for:

1. Grid search
2. Evolutionary algorithm - fix simulated annealing...
3. SVD-based threshold inversion

Add volume change source to o92util