[CIG-SEISMO] specfem2d: snapshots of domain subset

[Dimitri Komatitsch]

Hi Phil,

Perfect!

If you use MPI i.e. if you use more than one processor core, the right 
change should probably be:

    ! horizontal min and max coordinates of the image
    xmin_color_image_loc = max(minval(coord(1,:)), 
yyyyyy_put_here_the_min_X_that_you_want_yyyyyyyyyy)
    xmax_color_image_loc = min(maxval(coord(1,:), 
yyyyyy_put_here_the_max_X_that_you_want_yyyyyyyyyy)

    ! vertical min and max coordinates of the image, slightly increase it
to go beyond maximum topography
    zmin_color_image_loc = max(minval(coord(2,:), 
yyyyyy_put_here_the_min_Z_that_you_want_yyyyyyyyyy)
    zmax_color_image_loc = min(maxval(coord(2,:), 
yyyyyy_put_here_the_max_Z_that_you_want_yyyyyyyyyy)


If you use a single slice (no MPI), this should do it:

    ! horizontal min and max coordinates of the image
    xmin_color_image_loc = yyyyy_put_here_the_min_X_that_you_want_yyyyyyyyyy
    xmax_color_image_loc = 
yyyyyy_put_here_the_max_X_that_you_want_yyyyyyyyyy

    ! vertical min and max coordinates of the image, slightly increase it
to go beyond maximum topography
    zmin_color_image_loc = 
yyyyyy_put_here_the_min_Z_that_you_want_yyyyyyyyyy
    zmax_color_image_loc = 
yyyyyy_put_here_the_max_Z_that_you_want_yyyyyyyyyy

Best regards,
Dimitri.

On 10/11/2017 04:47 PM, Phil Cummins wrote:
> Hi Dimitri,
> 
> Well try it and let yo know. Thanks!
> 
> - Phil
> 
> ------------------------------------------------------------------------
> *From:* Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
> *Sent:* Wednesday, October 11, 2017 10:12:15 PM
> *To:* cig-seismo at geodynamics.org; Phil Cummins
> *Subject:* Re: [CIG-SEISMO] specfem2d: snapshots of domain subset
> 
> Hi Phil,
> 
> It is not implemented, but in src/specfem2D/prepare_color_image.F90 if
> you replace
> 
>     ! horizontal min and max coordinates of the image
>     xmin_color_image_loc = minval(coord(1,:))
>     xmax_color_image_loc = maxval(coord(1,:))
> 
>     ! vertical min and max coordinates of the image, slightly increase it
> to go beyond maximum topography
>     zmin_color_image_loc = minval(coord(2,:))
>     zmax_color_image_loc = maxval(coord(2,:))
> 
> with the min and max that you want it will probably work (please let us
> know).
> 
> Thanks,
> Dimitri.
> 
> On 10/11/2017 10:27 AM, Phil Cummins wrote:
>> Hi,
>> 
>> I'm using spefem2d to compute seismic waves n a sedimentary basin.
>> Does anyone know how I can output snapshots of a subset of the 
>> computational domain only, not the whole thing? I want to view the 
>> seismic waves in a shallow basin that is barely visible in a snapshot of 
>> the whole computational domain. I.e., I want to zoom in on the basin. 
>> Many thanks!
>> 
> 
> -- 
> Dimitri Komatitsch, CNRS Research Director (DR CNRS)
> Laboratory of Mechanics and Acoustics, Marseille, France
> http://komatitsch.free.fr

-- 
Dimitri Komatitsch, CNRS Research Director (DR CNRS)
Laboratory of Mechanics and Acoustics, Marseille, France
http://komatitsch.free.fr