view release on metacpan or  search on metacpan

lib/App/SeismicUnixGui/big_streams/iSuvelan.pm  view on Meta::CPAN

		$suvelan->smute( quotemeta( ( $iSuvelan->{_smute} ) ) );
		$suvelan->dtratio( quotemeta( ( $iSuvelan->{_dtratio} ) ) );
		$suvelan->nsmooth( quotemeta( ( $iSuvelan->{_nsmooth} ) ) );

		#				$suvelan->verbose(
		#			quotemeta(
		#				( $iSuvelan->{_verbose} )
		#			)
		#		);
		$suvelan->pwr( quotemeta( ( $iSuvelan->{_pwr} ) ) );

		$suvelan[1] = suvelan->Step();

=head2 DISPLAY DATA
 
 new dx_major_divisions 1-14-20

=cut

		my $time_inc_major = $iSuvelan->{_tmax_s} / 10;

		my $number_minor_time_divisions = 5;

		# to make Jorge Reyes' data work
		# my $new_dt_s            = $iSuvelan->{_tmax_s} * $iSuvelan->{_dt_s};

		# to make Daniel Lopez's data work
		#	my $s2ms     = 1000;
		#	my $new_dt_s = $s2ms * $iSuvelan->{_dt_s};

		my $dx_major_divisions = (
			$iSuvelan->{_first_velocity} + (
				$iSuvelan->{_number_of_velocities} *
				  $iSuvelan->{_velocity_increment}
			)
		) / 10;

#	print("iSuvelan,dx_major_divisions:$dx_major_divisions\n");
#	print("iSuvelan,time_inc_major: $time_inc_major\n");
#	print("iSuvelan,number_minor_time_divisions:$number_minor_time_divisions \n");
#	print("iSuvelan,tmax_s $iSuvelan->{_tmax_s}\n");
#	print("iSuvelan,new_dt_s $new_dt_s (ms)\n");

		$N = 2;
		$windowtitle[1] = '\('
		  . $N . '\)\ '
		  . $iSuvelan->{_sufile_in}
		  . '\ CDP=\ '
		  . $iSuvelan->{_cdp_num};

		#units=Semblance					\\
		$suximage->clear();
		$suximage->box_width( quotemeta(300) );
		$suximage->box_height( quotemeta(450) );
		$suximage->box_X0( quotemeta(600) );
		$suximage->box_Y0( quotemeta(0) );
		$suximage->title( quotemeta( $iSuvelan->{_cdp_num_suffix} ) );
		$suximage->windowtitle( quotemeta('Semblance') );
		$suximage->ylabel( quotemeta('TWTT (s)') );
		$suximage->xlabel( quotemeta('Velocity (m/s)') );
		$suximage->legend( quotemeta($on) );
		$suximage->cmap( quotemeta('hsv1') );
		$suximage->first_x( quotemeta( $iSuvelan->{_first_velocity} ) );
		$suximage->dx( quotemeta( $iSuvelan->{_velocity_increment} ) );
		$suximage->dt_s(quotemeta($new_dt_s));
		$suximage->loclip( quotemeta( $iSuvelan->{_min_semblance} ) );
		$suximage->hiclip( quotemeta( $iSuvelan->{_max_semblance} ) );
		$suximage->verbose( quotemeta($off) );

		#print ("d2num is ($iSuvelan->{_velocity_increment})\n\n");
		$suximage->dx_major_divisions( quotemeta($dx_major_divisions) );
		$suximage->dy_minor_divisions(
			quotemeta($number_minor_time_divisions) );
		$suximage->dy_major_divisions( quotemeta($time_inc_major) );

		$suximage->tend_s( quotemeta( $iSuvelan->{_tmax_s} ) );

		#$suximage -> percent4clip(quotemeta(95.0);
		$suximage->first_tick_number_x(
			quotemeta( $iSuvelan->{_first_velocity} ) );
		$suximage->picks( $iSuvelan->{_Tvel_outbound} );

print("iSuvelan, calcNdisplay: Writing picks to $iSuvelan->{_Tvel_outbound}\n\n");

=head2 conditions
 
 when number_of_tries is >=2 
 there should be a pre-exisiting digitized
 overlay curve to plot as well

=cut

		if ( $iSuvelan->{_number_of_tries} >= 2 ) {

			print("using a curve file:\n");
			print("\t$iSuvelan->{_Tvel_inbound}\n\n");
			$suximage->curvefile( quotemeta( $iSuvelan->{_Tvel_inbound} ) )
			  ;
			my ( $ref_T_nmo, $ref_Vnmo, $num_tvel_pairs ) =
			$manage_files_by->read_2cols( \$iSuvelan->{_Tvel_inbound} );
			$suximage->npair( quotemeta($num_tvel_pairs) );
			$suximage->curvecolor( quotemeta(2) );
		}

		$suximage[1] = $suximage->Step();

		# print("$suximage[1]\n\n");

		$N = 1;

=head2  set suxwigb parameters

  In the perl module for suxwigb we should
  have (but we do not yet) an explanation of each of
  these parameters

=cut

		$suxwigb->clear();
		$suxwigb->title( quotemeta( $iSuvelan->{_sufile_in} ) );
		$suxwigb->ylabel( quotemeta('Time s') );