[!!!][WP] (issue7 refract): introduce new option -SN

Provide new option - still must be implemented.
This version does not compile.

src/ts/refract/include/refract_para.inc

@@ -30,6 +30,8 @@ c    14/11/2011 V1.4   remember whether minoff is forced on command line
 c    12/11/2012 V1.5   optionally do not time shift offset shifted traces
 c    12/11/2012 V1.6   added frame parameter
 c    24/10/2013 V1.7   added parameters for scaling mode -S3
+c    07/12/2015 V1.8   refering amplitude scaling to average in offest range
+c                      becomes the standard (without alternative)
 c
 c==============================================================================
 c
@@ -54,8 +56,8 @@ c minimum offset step to be taken as two different positions
 c true if plpar_minoff is forced by command line
       logical plpar_forceminoff
 c
-      logical plflag_m3avg
-      real plpar_m3avgxmin, plpar_m3avgxmax
+c offset range to calculate reference amplitude for amplitude scaling
+      real plpar_avgxmin, plpar_avgxmax
 c 
 c control modes
 c -------------
@@ -160,7 +162,7 @@ c common blocks
      &                       plflag_seistyle, plflag_ttstyle, plpar_radius,
      &                       plflag_hypoffset, plflag_tracenum,
      &                       plflag_tracename, plpar_forceminoff,
-     &                       plflag_m3avg, plpar_m3avgxmin, plpar_m3avgxmax
+     &                       plpar_avgrefxmin, plpar_avgrefxmax
       common /refract_elem/ elem_modbox, elem_filenames, elem_version,
      &                      elem_annot, elem_scales, elem_data,
      &                      elem_syntt, elem_picks, plflag_subscale,

src/ts/refract/sub/refract_cmdopt.f

@@ -39,6 +39,7 @@ c    24/10/2013   V1.12  - new option -So
 c                        - new option -S3
 c    18/11/2013   V1.13  new option -TF
 c    21/03/2104 thof:    new option -TR
+c    07/12/2015 thof:    new option -SN
 c
 c==============================================================================
 c
@@ -64,7 +65,7 @@ c declare local variables
       integer i
 c commandline
       integer maxopt
-      parameter(maxopt=68)
+      parameter(maxopt=69)
       character*3 optid(maxopt)
       character*120 optarg(maxopt)
       logical optset(maxopt), opthasarg(maxopt)
@@ -136,7 +137,10 @@ c seismograms scaling
       data optarg(65)              /'0'/
       data optid(66)               /'-S3'/
       data opthasarg(66)           /.TRUE./
-      data optarg(66)              /'0.,100.'/
+      data optarg(66)              /'0.,10.'/
+      data optid(69)               /'-SN'/
+      data opthasarg(69)           /.TRUE./
+      data optarg(69)              /'0.,10.'/
 c
 c additionals
       data (optid(i),     i=55,56) /'-Lt','-Ta'/
@@ -221,8 +225,10 @@ c seismogram scaling
       read(optarg(65), *) opt_Sordinate
       if ((opt_Sordinate.lt.0).or.(opt_Sordinate.gt.3))
      &  stop 'ERROR: argument to -So is out of range'
-      opt_Savgref=optset(66)
-      read(optarg(66), *), opt_Savgrefxmin, opt_Savgrefxmax
+      opt_Savgref=optset(66).or.optset(69)
+      read(optarg(69), *), opt_Savgrefxmin, opt_Savgrefxmax
+      if (optset(66)) 
+     &  read(optarg(66), *), opt_Savgrefxmin, opt_Savgrefxmax
       if ((opt_Savgrefxmin.ge.opt_Savgrefxmax)
      &    .or.(opt_Savgrefxmin.lt.0.)
      &    .or.(opt_Savgrefxmax.lt.0.))
@@ -265,7 +271,7 @@ c
       plpar_forceminoff=optset(53)
       if (plpar_forceminoff) plpar_minoff=opt_Sminoff
       if (optset(54)) plpar_radius=opt_Sradius
-      if (opt_Savgref) opt_Smode=3
+      if (optset(66)) opt_Smode=3
       plpar_mode=opt_Smode
       print *,plpar_mode
 c
@@ -286,9 +292,8 @@ c
       plflag_osnoreduce=opt_Sosnoreduce
       plpar_vred=opt_Svel
 c
-      plflag_m3avg=opt_Savgref
-      plpar_m3avgxmin=opt_Savgrefxmin
-      plpar_m3avgxmax=opt_Savgrefxmax
+      plpar_avgrefxmin=opt_Savgrefxmin
+      plpar_avgrefxmax=opt_Savgrefxmax
 c 
       plflag_grid=opt_Egrid
       plflag_vara=opt_Ewiggle

src/ts/refract/sub/refract_usage.f

@@ -38,6 +38,7 @@ c    24/10/2013   V1.9   - new option -So
 c                        - new option -S3
 c    18/11/2013   V1.10  new option -TF
 c    21/03/2014 thof:    new option -TR
+c    07/12/2015 htof:    new option -SN
 c
 c==============================================================================
 c
@@ -62,7 +63,7 @@ c go
       print *,'               [-Sx x1,x2] [-St t1,t2] [-Se exp] [-Sa lev]'
       print *,'               [-Sc lev] [-Sm mode] [-Sr vel] [-Si] [-SM]'
       print *,'               [-Sn] [-SR radius] [-SO minoff] [-Sh] [-So n]'
-      print *,'               [-S3 min,max]'
+      print *,'               [-S3 min,max] [-SN min,max]'
       print *,'               [-C] [-L] [-R radius] [-O minoff]'
       print *,'          file [t:n,n-n] [o:s] [s:i,s,w] [n:name]'
       print *,'               [h:h,l,s] [r:r,g,b] [f:format] [v:f]'
@@ -193,8 +194,10 @@ c
       print *,'                 2: scale all traces to first trace as reference'
       print *,'                 3: scale all traces to nearest offset'
       print *,'                    trace of each dataset as reference'
-      print *,'-S3 min,max    use scaling mode 3 but refer to average'
-      print *,'               amplitude in offset range from min to max'
+      print *,'-SN min,max    refer amplitude scaling to average'
+      print *,'               amplitude in range from min to max.'
+      print *,'-S3 min,max    equivalent to ''-Sm 3 -SN min,max'' '
+      print *,'               provided for compatibility only'
       print *,'-Sr vel        set traveltime reduction velocity'
       print *,'-Sn            do not align offset shifted traces along'
       print *,'               slope of reduced times'