Skip to content

LOCTOPO_SURFACE 0/0 Z-level ASCII data points read. #76

@ProfeSergioM

Description

@ProfeSergioM

I sucessfull make to read .grd.asc file, but i think is doing nothing. At some part of the log, i get the line "0/0 Z-level ASCII data points read"

This is part of the file:

Maule.grd: Title: Grid imported via GDAL
Maule.grd: Command:
Maule.grd: Remark:
Maule.grd: Gridline node registration used [Cartesian grid]
Maule.grd: Grid file format: gd = Import/export through GDAL
Maule.grd: x_min: -70.87375 x_max: -70.1370833334 x_inc: 0.000277777777778 name: x n_columns: 2653
Maule.grd: y_min: -36.2890277778 y_max: -35.7693055555 y_inc: 0.000277777777778 name: y n_rows: 1872
Maule.grd: v_min: 784 v_max: 3966 name: z
Maule.grd: scale_factor: 1 add_offset: 0
Maule.grd: Default CPT:
1383
1372
1391
1404
1412
1425
1438
1451

And the log:

NLLoc (NonLinLoc v7.1.04 05Sep2024)
GRID: {x, y, z}
Num: {101, 101, 101}
Orig: {-50, -50, -3.092}
LenSide: {1, 1, 1}
Type: PROB_DENSITY
NLLoc LOCQUAL2ERR:
0 -> 0.100000
1 -> 0.500000
2 -> 1.000000
3 -> 2.000000
4 -> 3.000000
5 -> 4.000000
6 -> 5.000000
7 -> 6.000000
8 -> 7.000000
9 -> 99999.900000

GMT grd file header: ./static/Maule.grd.asc
phline: "Maule.grd: Title: Grid imported via GDAL
"
" Grid imported via GDAL" /* Descriptive title of the data set /
phline: "Maule.grd: Command:
"
" " /
Command line that produced the grdfile /
phline: "Maule.grd: Remark:
"
" " /
Any additional comments /
phline: "Maule.grd: Gridline node registration used [Cartesian grid]
"
sscanf istat=2
1 /
0 for grid line reg, 1 for pixel reg /
phline: "Maule.grd: Grid file format: gd = Import/export through GDAL
"
phline: "Maule.grd: x_min: -70.87375 x_max: -70.1370833334 x_inc: 0.000277777777778 name: x n_columns: 2653
"
sscanf istat=5
-70.873750 /
Minimum x-value of region /
-70.137083 /
Maximum x-value of region /
0.000278 /
Node spacing in x-dimension /
x /
Units of the x-dimension /
0 /
Number of nodes in the x-dimension /
phline: "Maule.grd: y_min: -36.2890277778 y_max: -35.7693055555 y_inc: 0.000277777777778 name: y n_rows: 1872
"
sscanf istat=5
-36.289028 /
Minimum y-value of region /
-35.769306 /
Maximum y-value of region /
0.000278 /
Node spacing in y-dimension /
y /
Units of the y-dimension /
0 /
Number of nodes in the y-dimension /
/
X/Y grid is lat/lon. /
phline: "Maule.grd: v_min: 784 v_max: 3966 name: z
"
sscanf istat=1
0.000000 /
Minimum z-value in data set /
0.000000 /
Maximum z-value in data set /
/
Units of the z-dimension /
phline: "Maule.grd: scale_factor: 1 add_offset: 0
"
sscanf istat=3
1.000000 /
Factor to multiply z-values after read /
0.000000 /
Offset to add to scaled z-values */
0/0 Z-level ASCII data points read.
INFO: no Search Prior (LOCPRIOR) values read.
INFO: no Search Prior (LOCPOSTERIOR) values read.
INFO: no Component Descirption (LOCCMP) params read.
INFO: no Alias (LOCALIAS) params read.
INFO: no Exclude (LOCEXCLUDE) params read.
INFO: no Include (LOCINCLUDE) params read.
INFO: no Time Delay (LOCDELAY) params read.
INFO: no Time Delay Surface (LOCDELAY_SURFACE) params read.
INFO: no Elevation Correction (LOCELEVCORR) params read.
INFO: no Fixed Origin Time (LOCFIXOTIME) params read.
INFO: no Station (LOCSRCE or GTSRCE) params read.

... Reading observation file ./obs/LagMaule_2025-09-02T03-53-33.630.obs

Reading next set of observations (Files open: Tot:4 Buf:0 Hdr:0 Alloc: 0) ...

... 6 observations read, 6 will be used for location (./loc/LagMaule_2025-09-02T03-53-33.630.20250902.035333).
Station Dist Weight: Ave Station Distance: 8.198116
Station Dist Weight: SPUE 1.147739 (4.791013,1.303355,-2.380000)
Station Dist Weight: SMAU 0.846050 (-2.214560,-0.804563,-2.215000)
Station Dist Weight: SNI2 1.006211 (-5.808418,-4.532965,-2.326000)
Station Dist Weight: SPUE 1.147739 (4.791013,1.303355,-2.380000)
Station Dist Weight: SMAU 0.846050 (-2.214560,-0.804563,-2.215000)
Station Dist Weight: SNI2 1.006211 (-5.808418,-4.532965,-2.326000)
LOCGAU param CorrLen is zero, will not be used: 0.000000
Locating... (Files open: Tot:4 Buf:0 Hdr:0 Alloc: 12 3DMem: used:0/avail:0/load:0) ...

Applying Octtree search within Grid 0:
INFO: EDT_otime_weight activated, OT_WT exceeds EDT_OT_WT_FLOOR.
OctTree num samples = 50000 / 50000
Octree oct_node_value_max= 3.128423e+01 oct_tree_integral= 4.011183e+00
OCTREE nInitial 1000 nEvaluated 50024 smallestNodeSide 0.078125/0.078125/0.078125 oct_tree_integral 4.011183e+00
INFO: EDT_otime_weight: ot_ml_std 0.008329

ellipsoid_volume = 5.675377e+03

Finished location: ./loc/LagMaule_2025-09-02T03-53-33.630.20250902.035333.grid0.loc.hyp

...end of observation file detected.

No more observation files. 1 events read, 1 events located, 1 locations completed.

Thanks in advance!

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions