
 # subdomains  in x-direction nsx =     1
 # subdomains  in y-direction nsy =     1

 # control subdomain ncontrol =     0

 Global # elements  in x-direction ngx =   400
 Global # elements  in y-direction ngy =   400

 Local # elements  in x-direction nx =   400
 Local # elements  in y-direction ny =   400

  processor grid size in the x-dir  hpx =        200.00000    
  processor grid size in the y-dir  hpy  =        200.00000    

  domain size xsize =       0.20000000E+03 m

  domain size ysize =       0.20000000E+03 m

   maxlow (max.iter. in low freq range ) =   4000 perclow ( perc. of(low) freq. solved with maxlow =       0.50000000E-01

  rounded xsize =       0.20000000E+03 m ; 
 rounded ysize =       0.20000000E+03 m 


  mesh size h=       0.50000000E+00 m

  number of snapshots in freq. domain=     1

  freq. to compute  snapshots in freq.domain
    1         0.30000000E-01 kHz

  number of snapshots in time domain=     4

  times to compute  snapshots in time domain
    1         0.10000000E+02 msec

  times to compute  snapshots in time domain
    2         0.20000000E+02 msec

  times to compute  snapshots in time domain
    3         0.25000000E+02 msec

  times to compute  snapshots in time domain
    4         0.30000000E+02 msec

  quadrature choice=1,standard,=2,mid point
  keyqua =     2

  number of points to compute solution= 
 in the frequency domain    npf  =   200

   key==1 for type 1 source (tsang&rader paper) ;
       key==2 for type 2 source(gaussian derivative)             ;    2



   frequency  freq=       0.70000000E+00 khz



  w0 =       0.43982297E+01  alfa =       0.11060000E+01



  total simulation time  stime =        0.12000000E+03 msec


   parameters of  source band-pass filter

 set equal to 0. if not filter is desired

 f1  =       0.00000000E+00   khz 

 f2  =       0.00000000E+00   khz 

 f3  =       0.15000000E+01   khz 

 f4  =       0.16000000E+01   khz 



   maxiter =  max. # inner iterations  4000

   reduc =  relative error required        0.10000000E-03

  number of materials (nmat) =     5

  mat. number   density (ros)    bulk modulus (xkr)   Shear modulus (xmur)


     1          0.260000E+07     0.133000E+07          0.850000E+06

     2          0.260000E+07     0.133000E+07          0.850000E+06

     3          0.260000E+07     0.133000E+07          0.850000E+06

     4          0.260000E+07     0.133000E+07          0.850000E+06

     5          0.260000E+07     0.133000E+07          0.850000E+06


  mat. number  deltp (p-log.decr.)     delts (s-log.decr.)   


     1          0.000000E+00          0.000000E+00

     2          0.000000E+00          0.000000E+00

     3          0.000000E+00          0.000000E+00

     4          0.000000E+00          0.000000E+00

     5          0.000000E+00          0.000000E+00


  mat. number   density (ros)    bulk modulus (xkr)   Shear modulus (xmur)


     1          0.260000E+07     0.133000E+07          0.850000E+06

     2          0.260000E+07     0.133000E+07          0.850000E+06

     3          0.260000E+07     0.133000E+07          0.850000E+06

     4          0.260000E+07     0.133000E+07          0.850000E+06

     5          0.260000E+07     0.133000E+07          0.850000E+06


  bulmin-umin-bulmax-umax   

     0.13300E+07    0.85000E+06    0.13300E+07    0.85000E+06


 minimum  phase pp velocity   =  0.97336E+00 m/msec

 maximum  phase pp velocity  =  0.97336E+00 m/msec



 minimum  phase ss velocity   =  0.57177E+00 m/msec

 maximum  phase ss velocity  =  0.57177E+00 m/msec



 phase  velocity  pp1 =  0.97336E+00 m/msec

 phase velocity pp2 =  0.97336E+00 m/msec

 phase velocity ss1 =  0.57177E+00 m/msec

 phase velocity ss2 =  0.57177E+00 m/msec



   mat. number      fq1                     fq2    (Liu model)  


     1        0.10000000E-06 kHz       0.10000000E+08 kHz

     2        0.10000000E-06 kHz       0.10000000E+08 kHz

     3        0.10000000E-06 kHz       0.10000000E+08 kHz

     4        0.10000000E-06 kHz       0.10000000E+08 kHz

     5        0.10000000E-06 kHz       0.10000000E+08 kHz

  material number     qp(i)  (Liu model)


     1              0.00000000E+00

     2              0.00000000E+00

     3              0.00000000E+00

     4              0.00000000E+00

     5              0.00000000E+00

  material number     const(i)       tau1(i)(msec)          tau2 (i)(msec) 


     1             0.000000E+00      0.159155E+07           0.159155E-07

     2             0.000000E+00      0.159155E+07           0.159155E-07

     3             0.000000E+00      0.159155E+07           0.159155E-07

     4             0.000000E+00      0.159155E+07           0.159155E-07

     5             0.000000E+00      0.159155E+07           0.159155E-07

 Material  Number               Ref. Speed ve    relaxed  plane wave modulus xmr


            1                   0.97336E+00         0.24633E+07

            2                   0.97336E+00         0.24633E+07

            3                   0.97336E+00         0.24633E+07

            4                   0.97336E+00         0.24633E+07

            5                   0.97336E+00         0.24633E+07

  time step  dt =       0.51368287E-01 msec



  number of time steps  nt =   2336



  kdel=0 : delta, kdel = 1:  d/dx(delta). kdel =      1

   mute in msec  tmute =       0.00000000E+00
 mute in samples   ntmute =     0


  jsou =     5  ksou =     5

  xsou =       0.20000000E+01 m ysou =       0.20000000E+01 m

  number of receivers, nrec =    80
  xinic    2.0000000000000000     

  yinic    2.0000000000000000     

  xfin    190.00000000000000     

  yfin    2.0000000000000000     

  deltarecx    2.3797468354430378     

  deltarecy    0.0000000000000000     


 receiver i =    1   irx & iry =     5    5
               rx(i) & ry(i) =   2.0000000      2.0000000      m

 receiver i =    2   irx & iry =     9    5
               rx(i) & ry(i) =   4.3797468      2.0000000      m

 receiver i =    3   irx & iry =    14    5
               rx(i) & ry(i) =   6.7594937      2.0000000      m

 receiver i =    4   irx & iry =    19    5
               rx(i) & ry(i) =   9.1392405      2.0000000      m

 receiver i =    5   irx & iry =    24    5
               rx(i) & ry(i) =   11.518987      2.0000000      m

 receiver i =    6   irx & iry =    28    5
               rx(i) & ry(i) =   13.898734      2.0000000      m

 receiver i =    7   irx & iry =    33    5
               rx(i) & ry(i) =   16.278481      2.0000000      m

 receiver i =    8   irx & iry =    38    5
               rx(i) & ry(i) =   18.658228      2.0000000      m

 receiver i =    9   irx & iry =    43    5
               rx(i) & ry(i) =   21.037975      2.0000000      m

 receiver i =   10   irx & iry =    47    5
               rx(i) & ry(i) =   23.417722      2.0000000      m

 receiver i =   11   irx & iry =    52    5
               rx(i) & ry(i) =   25.797468      2.0000000      m

 receiver i =   12   irx & iry =    57    5
               rx(i) & ry(i) =   28.177215      2.0000000      m

 receiver i =   13   irx & iry =    62    5
               rx(i) & ry(i) =   30.556962      2.0000000      m

 receiver i =   14   irx & iry =    66    5
               rx(i) & ry(i) =   32.936709      2.0000000      m

 receiver i =   15   irx & iry =    71    5
               rx(i) & ry(i) =   35.316456      2.0000000      m

 receiver i =   16   irx & iry =    76    5
               rx(i) & ry(i) =   37.696203      2.0000000      m

 receiver i =   17   irx & iry =    81    5
               rx(i) & ry(i) =   40.075949      2.0000000      m

 receiver i =   18   irx & iry =    85    5
               rx(i) & ry(i) =   42.455696      2.0000000      m

 receiver i =   19   irx & iry =    90    5
               rx(i) & ry(i) =   44.835443      2.0000000      m

 receiver i =   20   irx & iry =    95    5
               rx(i) & ry(i) =   47.215190      2.0000000      m

 receiver i =   21   irx & iry =   100    5
               rx(i) & ry(i) =   49.594937      2.0000000      m

 receiver i =   22   irx & iry =   104    5
               rx(i) & ry(i) =   51.974684      2.0000000      m

 receiver i =   23   irx & iry =   109    5
               rx(i) & ry(i) =   54.354430      2.0000000      m

 receiver i =   24   irx & iry =   114    5
               rx(i) & ry(i) =   56.734177      2.0000000      m

 receiver i =   25   irx & iry =   119    5
               rx(i) & ry(i) =   59.113924      2.0000000      m

 receiver i =   26   irx & iry =   123    5
               rx(i) & ry(i) =   61.493671      2.0000000      m

 receiver i =   27   irx & iry =   128    5
               rx(i) & ry(i) =   63.873418      2.0000000      m

 receiver i =   28   irx & iry =   133    5
               rx(i) & ry(i) =   66.253165      2.0000000      m

 receiver i =   29   irx & iry =   138    5
               rx(i) & ry(i) =   68.632911      2.0000000      m

 receiver i =   30   irx & iry =   143    5
               rx(i) & ry(i) =   71.012658      2.0000000      m

 receiver i =   31   irx & iry =   147    5
               rx(i) & ry(i) =   73.392405      2.0000000      m

 receiver i =   32   irx & iry =   152    5
               rx(i) & ry(i) =   75.772152      2.0000000      m

 receiver i =   33   irx & iry =   157    5
               rx(i) & ry(i) =   78.151899      2.0000000      m

 receiver i =   34   irx & iry =   162    5
               rx(i) & ry(i) =   80.531646      2.0000000      m

 receiver i =   35   irx & iry =   166    5
               rx(i) & ry(i) =   82.911392      2.0000000      m

 receiver i =   36   irx & iry =   171    5
               rx(i) & ry(i) =   85.291139      2.0000000      m

 receiver i =   37   irx & iry =   176    5
               rx(i) & ry(i) =   87.670886      2.0000000      m

 receiver i =   38   irx & iry =   181    5
               rx(i) & ry(i) =   90.050633      2.0000000      m

 receiver i =   39   irx & iry =   185    5
               rx(i) & ry(i) =   92.430380      2.0000000      m

 receiver i =   40   irx & iry =   190    5
               rx(i) & ry(i) =   94.810127      2.0000000      m

 receiver i =   41   irx & iry =   195    5
               rx(i) & ry(i) =   97.189873      2.0000000      m

 receiver i =   42   irx & iry =   200    5
               rx(i) & ry(i) =   99.569620      2.0000000      m

 receiver i =   43   irx & iry =   204    5
               rx(i) & ry(i) =   101.94937      2.0000000      m

 receiver i =   44   irx & iry =   209    5
               rx(i) & ry(i) =   104.32911      2.0000000      m

 receiver i =   45   irx & iry =   214    5
               rx(i) & ry(i) =   106.70886      2.0000000      m

 receiver i =   46   irx & iry =   219    5
               rx(i) & ry(i) =   109.08861      2.0000000      m

 receiver i =   47   irx & iry =   223    5
               rx(i) & ry(i) =   111.46835      2.0000000      m

 receiver i =   48   irx & iry =   228    5
               rx(i) & ry(i) =   113.84810      2.0000000      m

 receiver i =   49   irx & iry =   233    5
               rx(i) & ry(i) =   116.22785      2.0000000      m

 receiver i =   50   irx & iry =   238    5
               rx(i) & ry(i) =   118.60759      2.0000000      m

 receiver i =   51   irx & iry =   242    5
               rx(i) & ry(i) =   120.98734      2.0000000      m

 receiver i =   52   irx & iry =   247    5
               rx(i) & ry(i) =   123.36709      2.0000000      m

 receiver i =   53   irx & iry =   252    5
               rx(i) & ry(i) =   125.74684      2.0000000      m

 receiver i =   54   irx & iry =   257    5
               rx(i) & ry(i) =   128.12658      2.0000000      m

 receiver i =   55   irx & iry =   262    5
               rx(i) & ry(i) =   130.50633      2.0000000      m

 receiver i =   56   irx & iry =   266    5
               rx(i) & ry(i) =   132.88608      2.0000000      m

 receiver i =   57   irx & iry =   271    5
               rx(i) & ry(i) =   135.26582      2.0000000      m

 receiver i =   58   irx & iry =   276    5
               rx(i) & ry(i) =   137.64557      2.0000000      m

 receiver i =   59   irx & iry =   281    5
               rx(i) & ry(i) =   140.02532      2.0000000      m

 receiver i =   60   irx & iry =   285    5
               rx(i) & ry(i) =   142.40506      2.0000000      m

 receiver i =   61   irx & iry =   290    5
               rx(i) & ry(i) =   144.78481      2.0000000      m

 receiver i =   62   irx & iry =   295    5
               rx(i) & ry(i) =   147.16456      2.0000000      m

 receiver i =   63   irx & iry =   300    5
               rx(i) & ry(i) =   149.54430      2.0000000      m

 receiver i =   64   irx & iry =   304    5
               rx(i) & ry(i) =   151.92405      2.0000000      m

 receiver i =   65   irx & iry =   309    5
               rx(i) & ry(i) =   154.30380      2.0000000      m

 receiver i =   66   irx & iry =   314    5
               rx(i) & ry(i) =   156.68354      2.0000000      m

 receiver i =   67   irx & iry =   319    5
               rx(i) & ry(i) =   159.06329      2.0000000      m

 receiver i =   68   irx & iry =   323    5
               rx(i) & ry(i) =   161.44304      2.0000000      m

 receiver i =   69   irx & iry =   328    5
               rx(i) & ry(i) =   163.82278      2.0000000      m

 receiver i =   70   irx & iry =   333    5
               rx(i) & ry(i) =   166.20253      2.0000000      m

 receiver i =   71   irx & iry =   338    5
               rx(i) & ry(i) =   168.58228      2.0000000      m

 receiver i =   72   irx & iry =   342    5
               rx(i) & ry(i) =   170.96203      2.0000000      m

 receiver i =   73   irx & iry =   347    5
               rx(i) & ry(i) =   173.34177      2.0000000      m

 receiver i =   74   irx & iry =   352    5
               rx(i) & ry(i) =   175.72152      2.0000000      m

 receiver i =   75   irx & iry =   357    5
               rx(i) & ry(i) =   178.10127      2.0000000      m

 receiver i =   76   irx & iry =   361    5
               rx(i) & ry(i) =   180.48101      2.0000000      m

 receiver i =   77   irx & iry =   366    5
               rx(i) & ry(i) =   182.86076      2.0000000      m

 receiver i =   78   irx & iry =   371    5
               rx(i) & ry(i) =   185.24051      2.0000000      m

 receiver i =   79   irx & iry =   376    5
               rx(i) & ry(i) =   187.62025      2.0000000      m

 receiver i =   80   irx & iry =   380    5
               rx(i) & ry(i) =   190.00000      2.0000000      m

  nyquist frequency =       0.97336319E+01 khz

  dt  =       0.51368287E-01  msec

 delf =       0.80000000E-02  khz

 delw =       0.50265482E-01  2 pi*khz



   parameters of  source band-pass filter

 set equal to 0. if not filter is desired

 f1  =       0.00000000E+00   khz 

 f2  =       0.00000000E+00   khz 

 f3  =       0.15000000E+01   khz 

 f4  =       0.16000000E+01   khz 



  freq. indices for snapshots 
    1      4

  time steps for snapshots 
    1    194

  time steps for snapshots 
    2    389

  time steps for snapshots 
    3    486

  time steps for snapshots 
    4    584

  if1 if2 if3 if4 =     1     1   188   201

 jw1 - jw2 - jwlow       1    200     10

  pulse length in samples     105

 psi=      0.39200000E+01width=      0.71428571E+00msec

 ts=      0.17857143E+01msec. totpul=      0.53571429E+01msec

  pulse length in samples    2336
  writes in 20 source in time
  writes in 21 amplitude source
  from  xi                     =       0.00000000E+00 m  
 up to xf                     =       0.20000000E+03 m 
 from  yi                     =       0.00000000E+00 m  
 up to yf                     =       0.35000000E+02 m 
 material type number  numat  =     1

  from  ixi (x-global grid #)    =     1
 up to ixf (x-global grid #)    =   400
 from  iyi (y-global grid #)    =     1
 up to yf  (y-global grid #)    =    70
 material type number  numat    =     1

  from  xi                     =       0.00000000E+00 m  
 up to xf                     =       0.20000000E+03 m 
 from  yi                     =       0.35000000E+02 m  
 up to yf                     =       0.65000000E+02 m 
 material type number  numat  =     3

  from  ixi (x-global grid #)    =     1
 up to ixf (x-global grid #)    =   400
 from  iyi (y-global grid #)    =    71
 up to yf  (y-global grid #)    =   130
 material type number  numat    =     3

  from  xi                     =       0.00000000E+00 m  
 up to xf                     =       0.20000000E+03 m 
 from  yi                     =       0.65000000E+02 m  
 up to yf                     =       0.71000000E+02 m 
 material type number  numat  =     4

  from  ixi (x-global grid #)    =     1
 up to ixf (x-global grid #)    =   400
 from  iyi (y-global grid #)    =   131
 up to yf  (y-global grid #)    =   142
 material type number  numat    =     4

  from  xi                     =       0.00000000E+00 m  
 up to xf                     =       0.20000000E+03 m 
 from  yi                     =       0.71000000E+02 m  
 up to yf                     =       0.12500000E+03 m 
 material type number  numat  =     3

  from  ixi (x-global grid #)    =     1
 up to ixf (x-global grid #)    =   400
 from  iyi (y-global grid #)    =   143
 up to yf  (y-global grid #)    =   250
 material type number  numat    =     3

  from  xi                     =       0.00000000E+00 m  
 up to xf                     =       0.20000000E+03 m 
 from  yi                     =       0.12500000E+03 m  
 up to yf                     =       0.13100000E+03 m 
 material type number  numat  =     5

  from  ixi (x-global grid #)    =     1
 up to ixf (x-global grid #)    =   400
 from  iyi (y-global grid #)    =   251
 up to yf  (y-global grid #)    =   262
 material type number  numat    =     5

  from  xi                     =       0.00000000E+00 m  
 up to xf                     =       0.20000000E+03 m 
 from  yi                     =       0.13100000E+03 m  
 up to yf                     =       0.16500000E+03 m 
 material type number  numat  =     3

  from  ixi (x-global grid #)    =     1
 up to ixf (x-global grid #)    =   400
 from  iyi (y-global grid #)    =   263
 up to yf  (y-global grid #)    =   330
 material type number  numat    =     3

  from  xi                     =       0.00000000E+00 m  
 up to xf                     =       0.20000000E+03 m 
 from  yi                     =       0.16500000E+03 m  
 up to yf                     =       0.20000000E+03 m 
 material type number  numat  =     2

  from  ixi (x-global grid #)    =     1
 up to ixf (x-global grid #)    =   400
 from  iyi (y-global grid #)    =   331
 up to yf  (y-global grid #)    =   400
 material type number  numat    =     2

  from  xi                     =      -0.10000000E+01 m  
 up to xf                     =      -0.10000000E+01 m 
 from  yi                     =      -0.10000000E+01 m  
 up to yf                     =      -0.10000000E+01 m 
 material type number  numat  =    -1


  jw1 , jw2  in main =     1      200

  frequency number   200 frequency :       0.15960000E+01
 jw  real sour, imag sour =    200     -0.53985E-03      0.39211E-03

  global error =        0.8448732E-02 iter. number =   100 denomin. =        0.3412115E-04

  global error =        0.6842986E-02 iter. number =   200 denomin. =        0.4716775E-04

  global error =        0.5427133E-02 iter. number =   300 denomin. =        0.6045024E-04

  global error =        0.4823721E-02 iter. number =   400 denomin. =        0.7273025E-04

  global error =        0.4169478E-02 iter. number =   500 denomin. =        0.8639193E-04

  global error =        0.2659167E-02 iter. number =   600 denomin. =        0.9588958E-04

  global error =        0.1745674E-02 iter. number =   700 denomin. =        0.9841321E-04

  global error =        0.6678606E-03 iter. number =   800 denomin. =        0.9981875E-04

  global error =        0.5885918E-03 iter. number =   900 denomin. =        0.1004619E-03

  global error =        0.5941031E-03 iter. number =  1000 denomin. =        0.1006395E-03

  global error =        0.5289645E-03 iter. number =  1100 denomin. =        0.1008089E-03

  global error =        0.3335453E-03 iter. number =  1200 denomin. =        0.1011929E-03

  global error =        0.1874430E-03 iter. number =  1300 denomin. =        0.1016489E-03

  global error =        0.1163631E-03 iter. number =  1400 denomin. =        0.1017399E-03

  global error =        0.1007112E-03 iter. number =  1500 denomin. =        0.1016467E-03

  convergence  for freq. #   200 achieved in  1505 inner iterations 
 relative error value =       0.99866329E-04


  frequency number   199 frequency :       0.15880000E+01
 jw  real sour, imag sour =    199     -0.12154E-02      0.72596E-03

  global error =        0.4656672E-02 iter. number =   100 denomin. =        0.2241465E-03

  global error =        0.3686921E-02 iter. number =   200 denomin. =        0.2685923E-03

  global error =        0.2996516E-02 iter. number =   300 denomin. =        0.3118308E-03

  global error =        0.2628593E-02 iter. number =   400 denomin. =        0.3545788E-03

  global error =        0.2269972E-02 iter. number =   500 denomin. =        0.4016694E-03

  global error =        0.1425812E-02 iter. number =   600 denomin. =        0.4437397E-03

  global error =        0.9211428E-03 iter. number =   700 denomin. =        0.4742278E-03

  global error =        0.3979722E-03 iter. number =   800 denomin. =        0.4969506E-03

  global error =        0.3433793E-03 iter. number =   900 denomin. =        0.5166181E-03

  global error =        0.3187141E-03 iter. number =  1000 denomin. =        0.5327511E-03

  global error =        0.2683547E-03 iter. number =  1100 denomin. =        0.5466909E-03

  global error =        0.1664571E-03 iter. number =  1200 denomin. =        0.5553350E-03

  global error =        0.1004046E-03 iter. number =  1300 denomin. =        0.5607215E-03

  convergence  for freq. #   199 achieved in  1301 inner iterations 
 relative error value =       0.99456960E-04


  frequency number   198 frequency :       0.15800000E+01
 jw  real sour, imag sour =    198     -0.20292E-02      0.97691E-03

  global error =        0.4488558E-02 iter. number =   100 denomin. =        0.8080931E-03

  global error =        0.3462452E-02 iter. number =   200 denomin. =        0.9021540E-03

  global error =        0.2917606E-02 iter. number =   300 denomin. =        0.9758201E-03

  global error =        0.2447628E-02 iter. number =   400 denomin. =        0.1036645E-02

  global error =        0.1975572E-02 iter. number =   500 denomin. =        0.1094620E-02

  global error =        0.1193529E-02 iter. number =   600 denomin. =        0.1122358E-02

  global error =        0.7542314E-03 iter. number =   700 denomin. =        0.1120771E-02

  global error =        0.4919524E-03 iter. number =   800 denomin. =        0.1082204E-02

  global error =        0.4298299E-03 iter. number =   900 denomin. =        0.1054192E-02

  global error =        0.3533353E-03 iter. number =  1000 denomin. =        0.1041070E-02

  global error =        0.2677501E-03 iter. number =  1100 denomin. =        0.1034360E-02

  global error =        0.1697622E-03 iter. number =  1200 denomin. =        0.1032235E-02

  global error =        0.1230427E-03 iter. number =  1300 denomin. =        0.1030826E-02

  global error =        0.1044217E-03 iter. number =  1400 denomin. =        0.1032055E-02

  convergence  for freq. #   198 achieved in  1431 inner iterations 
 relative error value =       0.99887297E-04


  frequency number   197 frequency :       0.15720000E+01
 jw  real sour, imag sour =    197     -0.29806E-02      0.11183E-02

  global error =        0.4121444E-02 iter. number =   100 denomin. =        0.1436400E-02

  global error =        0.3296099E-02 iter. number =   200 denomin. =        0.1628234E-02

  global error =        0.2687333E-02 iter. number =   300 denomin. =        0.1877268E-02

  global error =        0.2256409E-02 iter. number =   400 denomin. =        0.2150472E-02

  global error =        0.1848072E-02 iter. number =   500 denomin. =        0.2447420E-02

  global error =        0.1198802E-02 iter. number =   600 denomin. =        0.2679096E-02

  global error =        0.7721286E-03 iter. number =   700 denomin. =        0.2828130E-02

  global error =        0.3585814E-03 iter. number =   800 denomin. =        0.2877607E-02

  global error =        0.3174552E-03 iter. number =   900 denomin. =        0.2893055E-02

  global error =        0.2914166E-03 iter. number =  1000 denomin. =        0.2897977E-02

  global error =        0.2474553E-03 iter. number =  1100 denomin. =        0.2898970E-02

  global error =        0.1757503E-03 iter. number =  1200 denomin. =        0.2894510E-02

  global error =        0.1141540E-03 iter. number =  1300 denomin. =        0.2885276E-02

  convergence  for freq. #   197 achieved in  1320 inner iterations 
 relative error value =       0.99939035E-04


  frequency number   196 frequency :       0.15640000E+01
 jw  real sour, imag sour =    196     -0.40654E-02      0.11215E-02

  global error =        0.4598479E-02 iter. number =   100 denomin. =        0.3406689E-02

  global error =        0.3313648E-02 iter. number =   200 denomin. =        0.3420404E-02

  global error =        0.2848611E-02 iter. number =   300 denomin. =        0.3422237E-02

  global error =        0.2239468E-02 iter. number =   400 denomin. =        0.3451670E-02

  global error =        0.1468196E-02 iter. number =   500 denomin. =        0.3561453E-02

  global error =        0.8741928E-03 iter. number =   600 denomin. =        0.3743122E-02

  global error =        0.5528441E-03 iter. number =   700 denomin. =        0.3958081E-02

  global error =        0.4247519E-03 iter. number =   800 denomin. =        0.4155685E-02

  global error =        0.3615913E-03 iter. number =   900 denomin. =        0.4325382E-02

  global error =        0.2778408E-03 iter. number =  1000 denomin. =        0.4447687E-02

  global error =        0.1944945E-03 iter. number =  1100 denomin. =        0.4514598E-02

  global error =        0.1415747E-03 iter. number =  1200 denomin. =        0.4545496E-02

  global error =        0.1170378E-03 iter. number =  1300 denomin. =        0.4562816E-02

  global error =        0.1086614E-03 iter. number =  1400 denomin. =        0.4574376E-02

  global error =        0.1015848E-03 iter. number =  1500 denomin. =        0.4587761E-02

  convergence  for freq. #   196 achieved in  1520 inner iterations 
 relative error value =       0.99990359E-04


  frequency number   195 frequency :       0.15560000E+01
 jw  real sour, imag sour =    195     -0.52754E-02      0.95682E-03

  global error =        0.4503465E-02 iter. number =   100 denomin. =        0.5466485E-02

  global error =        0.3491889E-02 iter. number =   200 denomin. =        0.5899395E-02

  global error =        0.2939874E-02 iter. number =   300 denomin. =        0.6237854E-02

  global error =        0.2428040E-02 iter. number =   400 denomin. =        0.6443537E-02

  global error =        0.1887929E-02 iter. number =   500 denomin. =        0.6523023E-02

  global error =        0.1311040E-02 iter. number =   600 denomin. =        0.6410320E-02

  global error =        0.7960719E-03 iter. number =   700 denomin. =        0.6143483E-02

  global error =        0.5163444E-03 iter. number =   800 denomin. =        0.5816902E-02

  global error =        0.4530072E-03 iter. number =   900 denomin. =        0.5583405E-02

  global error =        0.3901890E-03 iter. number =  1000 denomin. =        0.5460043E-02

  global error =        0.3154997E-03 iter. number =  1100 denomin. =        0.5414634E-02

  global error =        0.2333257E-03 iter. number =  1200 denomin. =        0.5420432E-02

  global error =        0.1568317E-03 iter. number =  1300 denomin. =        0.5451339E-02

  global error =        0.1342374E-03 iter. number =  1400 denomin. =        0.5484188E-02

  global error =        0.1264557E-03 iter. number =  1500 denomin. =        0.5503790E-02

  global error =        0.1203642E-03 iter. number =  1600 denomin. =        0.5510854E-02

  global error =        0.1132463E-03 iter. number =  1700 denomin. =        0.5510608E-02

  global error =        0.1060341E-03 iter. number =  1800 denomin. =        0.5506602E-02

  global error =        0.1001622E-03 iter. number =  1900 denomin. =        0.5501683E-02

  convergence  for freq. #   195 achieved in  1905 inner iterations 
 relative error value =       0.99965630E-04


  frequency number   194 frequency :       0.15480000E+01
 jw  real sour, imag sour =    194     -0.65978E-02      0.59325E-03

  global error =        0.3923576E-02 iter. number =   100 denomin. =        0.6609949E-02

  global error =        0.3066225E-02 iter. number =   200 denomin. =        0.7200339E-02

  global error =        0.2494778E-02 iter. number =   300 denomin. =        0.7967388E-02

  global error =        0.1990372E-02 iter. number =   400 denomin. =        0.8779890E-02

  global error =        0.1489438E-02 iter. number =   500 denomin. =        0.9578761E-02

  global error =        0.9405575E-03 iter. number =   600 denomin. =        0.1008795E-01

  global error =        0.6355850E-03 iter. number =   700 denomin. =        0.1025793E-01

  global error =        0.3865899E-03 iter. number =   800 denomin. =        0.1015204E-01

  global error =        0.3446108E-03 iter. number =   900 denomin. =        0.9985321E-02

  global error =        0.2902944E-03 iter. number =  1000 denomin. =        0.9841026E-02

  global error =        0.2250109E-03 iter. number =  1100 denomin. =        0.9727143E-02

  global error =        0.1650180E-03 iter. number =  1200 denomin. =        0.9636919E-02

  global error =        0.1280363E-03 iter. number =  1300 denomin. =        0.9560170E-02

  global error =        0.1135948E-03 iter. number =  1400 denomin. =        0.9477024E-02

  global error =        0.1070439E-03 iter. number =  1500 denomin. =        0.9406478E-02

  global error =        0.1006452E-03 iter. number =  1600 denomin. =        0.9355805E-02

  convergence  for freq. #   194 achieved in  1612 inner iterations 
 relative error value =       0.99984313E-04


  frequency number   193 frequency :       0.15400000E+01
 jw  real sour, imag sour =    193     -0.80150E-02     -0.59604E-06

  global error =        0.4098109E-02 iter. number =   100 denomin. =        0.1073522E-01

  global error =        0.3011217E-02 iter. number =   200 denomin. =        0.1090137E-01

  global error =        0.2425078E-02 iter. number =   300 denomin. =        0.1115756E-01

  global error =        0.1802022E-02 iter. number =   400 denomin. =        0.1158001E-01

  global error =        0.1242134E-02 iter. number =   500 denomin. =        0.1227921E-01

  global error =        0.7861434E-03 iter. number =   600 denomin. =        0.1314621E-01

  global error =        0.5918135E-03 iter. number =   700 denomin. =        0.1405680E-01

  global error =        0.4297396E-03 iter. number =   800 denomin. =        0.1484586E-01

  global error =        0.3183964E-03 iter. number =   900 denomin. =        0.1542009E-01

  global error =        0.2379754E-03 iter. number =  1000 denomin. =        0.1577511E-01

  global error =        0.1785089E-03 iter. number =  1100 denomin. =        0.1598738E-01

  global error =        0.1447860E-03 iter. number =  1200 denomin. =        0.1611394E-01

  global error =        0.1301213E-03 iter. number =  1300 denomin. =        0.1622185E-01

  global error =        0.1171978E-03 iter. number =  1400 denomin. =        0.1631118E-01

  global error =        0.1065464E-03 iter. number =  1500 denomin. =        0.1637782E-01

  convergence  for freq. #   193 achieved in  1588 inner iterations 
 relative error value =       0.99995613E-04


  frequency number   192 frequency :       0.15320000E+01
 jw  real sour, imag sour =    192     -0.95041E-02     -0.85600E-03

  global error =        0.4388729E-02 iter. number =   100 denomin. =        0.1862690E-01

  global error =        0.3531732E-02 iter. number =   200 denomin. =        0.1955770E-01

  global error =        0.3036931E-02 iter. number =   300 denomin. =        0.2030817E-01

  global error =        0.2544737E-02 iter. number =   400 denomin. =        0.2080773E-01

  global error =        0.1967329E-02 iter. number =   500 denomin. =        0.2095370E-01

  global error =        0.1336519E-02 iter. number =   600 denomin. =        0.2053332E-01

  global error =        0.8280547E-03 iter. number =   700 denomin. =        0.1977076E-01

  global error =        0.5647256E-03 iter. number =   800 denomin. =        0.1883316E-01

  global error =        0.5154968E-03 iter. number =   900 denomin. =        0.1811022E-01

  global error =        0.4482321E-03 iter. number =  1000 denomin. =        0.1771654E-01

  global error =        0.3590149E-03 iter. number =  1100 denomin. =        0.1761414E-01

  global error =        0.2761618E-03 iter. number =  1200 denomin. =        0.1768308E-01

  global error =        0.2134640E-03 iter. number =  1300 denomin. =        0.1779307E-01

  global error =        0.1916603E-03 iter. number =  1400 denomin. =        0.1789169E-01

  global error =        0.1819816E-03 iter. number =  1500 denomin. =        0.1794994E-01

  global error =        0.1743979E-03 iter. number =  1600 denomin. =        0.1796639E-01

  global error =        0.1689392E-03 iter. number =  1700 denomin. =        0.1795289E-01

  global error =        0.1650167E-03 iter. number =  1800 denomin. =        0.1792904E-01

  global error =        0.1603611E-03 iter. number =  1900 denomin. =        0.1790980E-01

  global error =        0.1553842E-03 iter. number =  2000 denomin. =        0.1789771E-01

  global error =        0.1510240E-03 iter. number =  2100 denomin. =        0.1789307E-01

  global error =        0.1486174E-03 iter. number =  2200 denomin. =        0.1789460E-01

  global error =        0.1481575E-03 iter. number =  2300 denomin. =        0.1789950E-01

  global error =        0.1479438E-03 iter. number =  2400 denomin. =        0.1790435E-01

  global error =        0.1461964E-03 iter. number =  2500 denomin. =        0.1790734E-01

  global error =        0.1430039E-03 iter. number =  2600 denomin. =        0.1790854E-01

  global error =        0.1396205E-03 iter. number =  2700 denomin. =        0.1790831E-01

  global error =        0.1374142E-03 iter. number =  2800 denomin. =        0.1790715E-01

  global error =        0.1369419E-03 iter. number =  2900 denomin. =        0.1790574E-01

  global error =        0.1369230E-03 iter. number =  3000 denomin. =        0.1790467E-01

  global error =        0.1356769E-03 iter. number =  3100 denomin. =        0.1790413E-01

  global error =        0.1333349E-03 iter. number =  3200 denomin. =        0.1790396E-01

  global error =        0.1305576E-03 iter. number =  3300 denomin. =        0.1790403E-01

  global error =        0.1285343E-03 iter. number =  3400 denomin. =        0.1790425E-01

  global error =        0.1278594E-03 iter. number =  3500 denomin. =        0.1790454E-01

  global error =        0.1278377E-03 iter. number =  3600 denomin. =        0.1790478E-01

  global error =        0.1268370E-03 iter. number =  3700 denomin. =        0.1790486E-01

  global error =        0.1249645E-03 iter. number =  3800 denomin. =        0.1790471E-01

  global error =        0.1226395E-03 iter. number =  3900 denomin. =        0.1790437E-01

  global error =        0.1207797E-03 iter. number =  4000 denomin. =        0.1790400E-01

  NO convergence  for freq. #   192 after   4001 inner iterations 
 relative error value reduction =       0.12077971E-03


  frequency number   191 frequency :       0.15240000E+01
 jw  real sour, imag sour =    191     -0.11037E-01     -0.20035E-02

  global error =        0.4025350E-02 iter. number =   100 denomin. =        0.2084405E-01

  global error =        0.3063234E-02 iter. number =   200 denomin. =        0.2270187E-01

  global error =        0.2495376E-02 iter. number =   300 denomin. =        0.2487284E-01

  global error =        0.1963109E-02 iter. number =   400 denomin. =        0.2724881E-01

  global error =        0.1520226E-02 iter. number =   500 denomin. =        0.2967940E-01

  global error =        0.1020893E-02 iter. number =   600 denomin. =        0.3155968E-01

  global error =        0.7037189E-03 iter. number =   700 denomin. =        0.3253726E-01

  global error =        0.3983223E-03 iter. number =   800 denomin. =        0.3260182E-01

  global error =        0.3347435E-03 iter. number =   900 denomin. =        0.3229053E-01

  global error =        0.2873058E-03 iter. number =  1000 denomin. =        0.3187780E-01

  global error =        0.2479730E-03 iter. number =  1100 denomin. =        0.3153291E-01

  global error =        0.2022965E-03 iter. number =  1200 denomin. =        0.3128568E-01

  global error =        0.1715764E-03 iter. number =  1300 denomin. =        0.3107813E-01

  global error =        0.1560288E-03 iter. number =  1400 denomin. =        0.3087072E-01

  global error =        0.1509076E-03 iter. number =  1500 denomin. =        0.3068402E-01

  global error =        0.1472200E-03 iter. number =  1600 denomin. =        0.3053350E-01

  global error =        0.1440469E-03 iter. number =  1700 denomin. =        0.3042540E-01

  global error =        0.1408741E-03 iter. number =  1800 denomin. =        0.3036081E-01

  global error =        0.1381295E-03 iter. number =  1900 denomin. =        0.3032772E-01

  global error =        0.1355814E-03 iter. number =  2000 denomin. =        0.3031353E-01

  global error =        0.1331774E-03 iter. number =  2100 denomin. =        0.3030735E-01

  global error =        0.1312605E-03 iter. number =  2200 denomin. =        0.3030516E-01

  global error =        0.1295375E-03 iter. number =  2300 denomin. =        0.3030438E-01

  global error =        0.1277665E-03 iter. number =  2400 denomin. =        0.3030491E-01

  global error =        0.1261134E-03 iter. number =  2500 denomin. =        0.3030722E-01

  global error =        0.1241231E-03 iter. number =  2600 denomin. =        0.3031057E-01

  global error =        0.1219803E-03 iter. number =  2700 denomin. =        0.3031386E-01

  global error =        0.1201751E-03 iter. number =  2800 denomin. =        0.3031661E-01

  global error =        0.1185634E-03 iter. number =  2900 denomin. =        0.3031846E-01

  global error =        0.1170801E-03 iter. number =  3000 denomin. =        0.3031952E-01

  global error =        0.1156646E-03 iter. number =  3100 denomin. =        0.3032019E-01

  global error =        0.1140038E-03 iter. number =  3200 denomin. =        0.3032069E-01

  global error =        0.1120930E-03 iter. number =  3300 denomin. =        0.3032103E-01

  global error =        0.1103433E-03 iter. number =  3400 denomin. =        0.3032121E-01

  global error =        0.1087535E-03 iter. number =  3500 denomin. =        0.3032129E-01

  global error =        0.1075088E-03 iter. number =  3600 denomin. =        0.3032128E-01

  global error =        0.1063038E-03 iter. number =  3700 denomin. =        0.3032125E-01

  global error =        0.1048501E-03 iter. number =  3800 denomin. =        0.3032120E-01

  global error =        0.1031663E-03 iter. number =  3900 denomin. =        0.3032114E-01

  global error =        0.1014677E-03 iter. number =  4000 denomin. =        0.3032110E-01

  NO convergence  for freq. #   191 after   4001 inner iterations 
 relative error value reduction =       0.10146766E-03


  frequency number   190 frequency :       0.15160000E+01
 jw  real sour, imag sour =    190     -0.12579E-01     -0.34721E-02

  global error =        0.4486560E-02 iter. number =   100 denomin. =        0.3347932E-01

  global error =        0.3251751E-02 iter. number =   200 denomin. =        0.3410941E-01

  global error =        0.2736396E-02 iter. number =   300 denomin. =        0.3446189E-01

  global error =        0.2098432E-02 iter. number =   400 denomin. =        0.3517365E-01

  global error =        0.1446402E-02 iter. number =   500 denomin. =        0.3662750E-01

  global error =        0.7993193E-03 iter. number =   600 denomin. =        0.3879472E-01

  global error =        0.5159083E-03 iter. number =   700 denomin. =        0.4124076E-01

  global error =        0.4202245E-03 iter. number =   800 denomin. =        0.4347332E-01

  global error =        0.3547396E-03 iter. number =   900 denomin. =        0.4530778E-01

  global error =        0.2784315E-03 iter. number =  1000 denomin. =        0.4659771E-01

  global error =        0.2071065E-03 iter. number =  1100 denomin. =        0.4739399E-01

  global error =        0.1469847E-03 iter. number =  1200 denomin. =        0.4784065E-01

  global error =        0.1275718E-03 iter. number =  1300 denomin. =        0.4815189E-01
