[Shake-dev] shakemap topography

[Hagos, Lijam]

Jordi,

 

In line to what Bruce suggested, you may also want to include the "strictbound" parameter in grind.conf (found towards the bottom of the file)  to force your map to have specific lat/lon within the boundary of your topo grid:

 

strictbound:  West/South/East/North     (in degrees)

 

At least this worked for me in the past.

 

Lijam

 

 

 

 

 

 

 

 

From: shake-dev-bounces at geohazards.usgs.gov [mailto:shake-dev-bounces at geohazards.usgs.gov] On Behalf Of Bruce Worden
Sent: Monday, March 12, 2012 11:21 AM
To: "Roviró Pou, Jordi"
Cc: shake-dev at geohazards.usgs.gov
Subject: Re: [Shake-dev] shakemap topography

 

The problem seems to be that your map's dimensions exceed the range of the file. The map goes from 1.55 to 4.05 degrees longitude, but the file's maximum longitude is 3.85 degrees.

 

You'll need to 1) reduce the size of the map, or 2) shift the map so that the right edge is 3.85 degrees or less, or 3) increase the size of your topo file.

 

Hope this helps,

Bruce

 

On Mar 12, 2012, at 1:59 AM, Roviró Pou, Jordi wrote:





Hello ,

We made two files of topography, one of these has more resolution (sispyr_topo.grd) than the other one (sispyr_m.grd).

The idea to create both files was to try to identify the problem, and we changed the topography to see if the problem was our GRD file but we had the same problem as well as above.

 

I attach both grdinfo of the two topography files.

 

Grdinfo sispyr_m.grd

 

xsismica at nobel grdinfo sispyr_m.grd

sispyr_m.grd: Title: sispyr_m.grd

sispyr_m.grd: Command: xyz2grd -V -R-2.51/3.85/40.83/44.33 -I1m Sispyr_m-8780.xyz -Gsispyr_m.grd

sispyr_m.grd: Remark:

sispyr_m.grd: Gridline node registration used

sispyr_m.grd: Grid file format: nf (# 18) GMT netCDF format (float)  (COARDS-compliant) [DEFAULT]

sispyr_m.grd: x_min: -2.51 x_max: 3.85 x_inc: 0.0166492146597 name: x nx: 383

sispyr_m.grd: y_min: 40.83 y_max: 44.33 y_inc: 0.0166666666667 name: y ny: 211

sispyr_m.grd: z_min: -2569 z_max: 3044 name: z

sispyr_m.grd: scale_factor: 1 add_offset: 0

 

grdinfo sispyr_topo.grd

 

xsismica at nobel grdinfo sispyr_topo.grd

sispyr_topo.grd: Title: sispyr_topo.grd

sispyr_topo.grd: Command: xyz2grd -V -R-2.51/3.85/40.83/44.33 -I23c/16c topo_sispyr.txt -Gsispyr_topo.grd

sispyr_topo.grd: Remark:

sispyr_topo.grd: Gridline node registration used

sispyr_topo.grd: Grid file format: nf (# 18) GMT netCDF format (float)  (COARDS-compliant) [DEFAULT]

sispyr_topo.grd: x_min: -2.51 x_max: 3.85 x_inc: 0.00639195979899 name: x nx: 996

sispyr_topo.grd: y_min: 40.83 y_max: 44.33 y_inc: 0.00444162436548 name: y ny: 789

sispyr_topo.grd: z_min: -3.83333325386 z_max: 3223.77783203 name: z

sispyr_topo.grd: scale_factor: 1 add_offset: 0

 

Thank you.

 

Best Regards

_______________________________________________
Shake-dev mailing list
Shake-dev at geohazards.usgs.gov
https://geohazards.usgs.gov/mailman/listinfo/shake-dev

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geohazards.usgs.gov/pipermail/shake-dev/attachments/20120312/34c540c1/attachment.html>