[Shake-dev] ShakeMap v3.2 grind question

[Vince Quitoriano]

Hi,

It looks like there's something wrong with the output of the 'surface'
command. I don't know what is different between Linux and SunOS grdfiles
that would cause this problem. It may be a roundoff error.

Unfortunately the output of the blockmean|surface step is not saved, so
you will have to do some editing to capture the result.

Try this:

Change the line
/usr/local/GMT4.2.1/bin//surface -G=1
-R123.916667/131.100000/32.900000/39.100000
-I0.0166666666666667/0.0166666666666667 -T0.9 -Ll0

to
/usr/local/GMT4.2.1/bin//surface -Gtest.grd=1
-R123.916667/131.100000/32.900000/39.100000
-I0.0166666666666667/0.0166666666666667 -T0.9 -Ll0

and run grind again. This should stop with an error, but create
'test.grd'. Run grdinfo on it to see if the lat-lon coordinates make sense.

If that doesn't work, you could check if your OS has a different size
for floating point numbers. Look at the filesize of test.grd. It should
be nx * ny * 4 + 892 bytes for the header.

Let me know if that doesn't help.

- Vince

Dong-Hoon Sheen wrote:
> Dear Shake-dev,
>  
> I just started working with a ShakeMap v3.2 installation on SunOS 
> 5.9 with GMT4.2.1.
> When running grind, I get an error (see below).
> I've tested this on Linux with GMT4.2.1, but it was fine.
>  
> Thank you for your help.
>  
> Dong-Hoon Sheen.
>  
> ---------------------------------------------------
> bin$./grind -event 20070120 -xml -qtm -verbose
> Using Regression::Small for regression
> grind: Using fixed boundary 124/33/131/39
> grind: 39.1:32.9:131.1:123.9
> grind:
>  
> Big bound: 123.916667/131.100000/32.900000/39.100000
> Map bound: 124.016667/131.000000/33.000000/39.000000
> Bias: original stations: 76, close stations: 16
> grind: Type: pga; bias 2.44375; absolute deviation was 1.472168 is 0.737168
> grind: Type: pgv; bias 1.9921875; absolute deviation was 0.000000 is 
> 0.000000
> grind: Bias 2.444 1.992
> Outlier standard deviations: 3.
> grind: 0 stations flagged out this iteration
> Done initializing grid, now iterating (4599).
> Done with point 01000.
> Done with point 02000.
> Done with point 03000.
> Done with point 04000.
> Done computing Richter grid estimates.
> Done with richter.
> Culling estimates.
> grind: Killing centroid (37.688900 128.584100) for type acc because of 
> station DGY (37.690300 128.674200)
>  
> Making pga rock grid.
> Running: /usr/local/GMT4.2.1/bin//blockmean 
> -R123.916667/131.100000/32.900000/39.100000 
> -I0.0166666666666667/0.0166666666666667 | 
> /usr/local/GMT4.2.1/bin//surface -G=1 
> -R123.916667/131.100000/32.900000/39.100000 
> -I0.0166666666666667/0.0166666666666667 -T0.9 -Ll0 
> grind: Couldn't read the grd data.
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> Shake-dev mailing list
> Shake-dev at geohazards.usgs.gov
> https://geohazards.usgs.gov/mailman/listinfo/shake-dev