Interface: map_amp

  The map_amp function computes the transformation from n_tr_inp input traces into n_tr_out output traces. The transformation is defined by the five arrays map, amp, first_samp, last_samp, and max_freq.

The map array is indexed in the ouput time/depth sampling, and contains the corresponding time coordinates of the input trace. These input-time coordinates are normalized so to be in samples, and not physical units.

The amp array is indexed in the ouput time/depth sampling, and contains the amplitude scaling factor of the transformation.

The first_samp and last_samp arrays define the valid range of the transformation, as the first and last output samples.

The max_freq array is indexed in the ouput time/depth sampling, and contains the high cut-off frequency of the bandpass filter applied to the input trace. The maximum frequency must be expressed as fraction of the Nyquist frequency.

The shapes of the output arrays, and whether each input array is allocated and initialized, depend on the spatial mapping method that is used to apply the operator.

• Trace-to-trace
  • n_tr_inp=1
  • n_tr_out=1
  • max_freq: not allocated
  • inp_cent_coord: not allocated
  • inp_geom%reg: not allocated
  • inp_geom%irr: initialized by Genkir3D
  • out_geom%reg: not allocated
  • out_geom%irr: initialized by Genkir3D
• Gathering
  • n_tr_inp=n_tr_inp
  • n_tr_out=1
  • max_freq: allocated
  • inp_cent_coord: not allocated
  • inp_geom%reg: not allocated
  • inp_geom%irr: initialized by Genkir3D
  • out_geom%reg: not allocated
  • out_geom%irr: initialized by Genkir3D
• Spraying
  • n_tr_inp=1
  • n_tr_out=(2*nh_z1+1)*(2*nh_z2+1)
  • max_freq: allocated
  • inp_cent_coord=initialized by oper_domain
  • inp_geom%reg: not allocated
  • inp_geom%irr: initialized by Genkir3D
  • out_geom%reg: initialized by Genkir3D
  • out_geom%irr: initialized by oper_domain

function map_amp (want_info,inv_map,inv_amp,inp_geom,out_geom, &
      first_samp, last_samp, map, amp, inp_cent_coord) result (info)

type(genkir3d_info_type) :: info
logical, intent(in)      :: want_info,inv_map,inv_amp
type (seis3d_geom_type), intent(in) :: inp_geom
type (seis3d_geom_type), intent(in) :: out_geom
integer, dimension(:,:), intent(out) :: first_samp,last_samp
real,  dimension (:,:,:), intent(out)  :: map,amp,max_freq
real, pointer,  optional,  dimension (:,:,:)  :: inp_cent_coord

Name Description

want_info If .true. return memory requirements in info

inv_map Compute mapping of inverse operator

inv_amp Compute amplitude of inverse operator

Name Description

info Information on operator

inp_geom Geometry of input traces

out_geom Geometry of output traces

Name Description Shape

first_samp First valid output sample (1:n_tr_inp,1:n_tr_out)

last_samp Last valid output sample (1:n_tr_inp,1:n_tr_out)

map Mapping function (1:n_tz_out,1:n_tr_inp,1:n_tr_out)

amp Amplitude function (1:n_tz_out,1:n_tr_inp,1:n_tr_out)

max_freq Maximum frequency (1:n_tz_out,1:n_tr_inp,1:n_tr_out)

inp_cent_coord Input centered coordinates (-nh_z1:nh_z1,-nh_z2:nh_z2,2)