[Shake-dev] small bug

[Renate Hartog]

Hi Bruce et al.,

I found a small bug with the grdcontour part in mapping. It seg-faults 
when trying to make a contour map of the attached pga_std.grd file, which 
has z-values between 0.71 and 0.79. It tries to use a contour interval of 
0.1, which bombs. I ran the identical command with C0.01 and it's fine.


Seems like changing the cont_scale function a bit should do it.

cheers, Renate


  mapping: Error running command:
/usr/local/GMT4.5.2/bin//grdcontour 
/home/SHAKE/bin/../data/02281854/output/pga_std.grd -P 
-R-124.7266000000/-120.7266000000/45.7890000000/48.5090000000 -JM16.51 
-C0.1 -W5/255/255/255 -Af12/255/255/255 
-D/home/SHAKE/bin/../data/02281854/mapping/cont_pga_std.dat -Q8 -m 
-L0.05/9999 -O -K


#-----------------------------------------------------------------------
       # Make contours for everything except Intensity (i.e. PGA, PGV, 
SP_??)

#-----------------------------------------------------------------------
       my $cont_int = cont_scale($ghd->zmax());
       my $surf_ll  = $cont_int / 2;
       $command = "$gmt_bin/grdcontour ${file} $gmtflags -C$cont_int "
                . "-W$cont_width/$cont_color -Af12/$cont_color "
                . "-D$outputdir/cont_${display}.dat -Q8 $mflag "
                . "-L$surf_ll/9999 -O -K";
       $params{$display} = [];
       slurp($command, $params{$display});
-------------- next part --------------
A non-text attachment was scrubbed...
Name: pga_std.grd
Type: application/octet-stream
Size: 73972 bytes
Desc: 
URL: <http://geohazards.usgs.gov/pipermail/shake-dev/attachments/20100511/49da974e/attachment-0001.obj>