[Shake-dev] small bug

[Bruce Worden]

Hi Renate,

Thanks for the bug report. grdcontour definitely has some  
idiosyncrasies. I'll see what I can do with it.

Bruce


On May 11, 2010, at 2:32 PM, Renate Hartog wrote:

>
> 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});
> <pga_std.grd>_______________________________________________
> Shake-dev mailing list
> Shake-dev at geohazards.usgs.gov
> https://geohazards.usgs.gov/mailman/listinfo/shake-dev