fixscalevalue()

void datrw::su::fixscalevalue	(	short &	s,
		const bool &	strict = true
	)

fix a SeismicUn*x scale value

fix a SeismicUn*x scale value

some code out there produces SU data with scale values understood as exponents of ten; SEGY defines the scale value to be taken as numerator or denominator of the scaling factor with values -10000, -1000, -100, -10, 1, 10, 100, 1000, 10000 only; here we allow for the non-standard setting indicated by the absolute value of scale being larger than 1 and smaller than 10

Parameters

s	scale value to be fixed (input and output parameter)
strict	if true: abort if scale does not match original definition

Definition at line 498 of file sucomanager.cc.

References DATRW_assert, DATRW_report_assert, and datrw::su::subformat::key::strict.

Referenced by outputmodifiers(), datrw::su::SUheader::read(), datrw::su::ScalCoo::ScalCoo(), scalefactor(), scaletopower(), datrw::su::ScalCoo::set(), and datrw::su::SUheader::setdefaults().

    {
      /*
       * a scale value of 0 is not defined for the trace header in the SEGY
       * standard; since Geode data acquisitions produce SEGY data with scalel
       * set to 0, we have to deal with this value
       */
      DATRW_report_assert(s != 0,
                          "WARNING (ScalCoo::set): incorrect scale\n"
                          "The value for the scalco and scalel field\n"
                          "violates the SEG-Y trace header definition.\n"
                          "value: " << s << "; will be set to 1");
      DATRW_assert(!((strict) && (s==0)), 
                   "Violation of SeismicUn*x format definition!");
      // make a copy to allow modification
      if (s==0) { s=1; }
      int sabs=s > 0 ? s : -s;
      if ((sabs < 10) && (s != 1))
      {
        /*
         * some code out there produces SU data with scale values understood
         * as exponents of ten; SEGY defines the scale value to be taken as
         * numerator or denominator of the scaling factor with values -10000,
         * -1000, -100, -10, 1, 10, 100, 1000, 10000 only; here we allow for
         * the non-standard setting indicated by the absolute value of scale
         * being larger than 1 and smaller than 10
         */
        DATRW_report_assert(((sabs < 10) && (s != 1)),
                            "WARNING (ScalCoo::set): non-standard scale\n"
                            "The value for the scalco and scalel field\n"
                            "violates the SEG-Y trace header definition.\n"
                            "value: " << s);
        DATRW_assert(!strict,
                     "Violation of SeismicUn*x format definition!");
        /*
         * 10 to a power larger than 4 cannot be represented by a short int
         */
        DATRW_assert((sabs<5),
                      "ERROR (ScalCoo::set): "
                      "scale value does not fit in short int field");
        s= s > 0 ? std::pow(10,sabs) : -std::pow(10,sabs);
      }
    } // void fixscalevalue(short& s, const bool& strict=true)

Here is the caller graph for this function: