[Shake-dev] small bug

[Bruce Worden]

Hi Renate,

It appears that you are using GMT 4.5.2, which is known to cause some  
problems. See the attached email exchange for an earlier discussion of  
a problem. The easiest solution is probably to downgrade your GMT to  
4.4.X (or even 4.5.0). I don't have 4.5.2 locally, but I'll try to get  
it and see if I can make it work with SM v3.5. In the meantime,  
though, you may want to avoid 4.5.2.

Bruce

-------------- next part --------------
An embedded message was scrubbed...
From: Kuo-Wan Lin <klin at usgs.gov>
Subject: Re: [Shake-dev] problems with grind
Date: Tue, 02 Feb 2010 11:21:56 -0700
Size: 7297
URL: <http://geohazards.usgs.gov/pipermail/shake-dev/attachments/20100514/4cc8723c/attachment.eml>
-------------- next part --------------


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