Home

Context Navigation

← Previous Change
Next Change →

Changeset 1257 for palm/trunk

Timestamp:

Nov 8, 2013 3:18:40 PM (11 years ago)

Author:

raasch

Message:

New:
---

openACC porting of timestep calculation
(modules, timestep, time_integration)

Changed:

openACC loop directives and vector clauses removed (because they do not give any performance improvement with PGI
compiler versions > 13.6)
(advec_ws, buoyancy, coriolis, diffusion_e, diffusion_s, diffusion_u, diffusion_v, diffusion_w, diffusivities, exchange_horiz, fft_xy, pres, production_e, transpose, tridia_solver, wall_fluxes)

openACC loop independent clauses added
(boundary_conds, prandtl_fluxes, pres)

openACC declare create statements moved after FORTRAN declaration statement
(diffusion_u, diffusion_v, diffusion_w, fft_xy, poisfft, production_e, tridia_solver)

openACC end parallel replaced by end parallel loop
(flow_statistics, pres)

openACC "kernels do" replaced by "kernels loop"
(prandtl_fluxes)

output format for theta* changed to avoid output of *
(run_control)

Errors:

bugfix for calculation of advective timestep (old version may cause wrong timesteps in case of
vertixcally stretched grids)
Attention: standard run-control output has changed!
(timestep)

Location:

palm/trunk

Files:

: 24 edited

INSTALL/example_cbl_rc (modified) (2 diffs)
SOURCE/advec_ws.f90 (modified) (5 diffs)
SOURCE/boundary_conds.f90 (modified) (4 diffs)
SOURCE/buoyancy.f90 (modified) (2 diffs)
SOURCE/coriolis.f90 (modified) (4 diffs)
SOURCE/diffusion_e.f90 (modified) (4 diffs)
SOURCE/diffusion_s.f90 (modified) (5 diffs)
SOURCE/diffusion_u.f90 (modified) (7 diffs)
SOURCE/diffusion_v.f90 (modified) (7 diffs)
SOURCE/diffusion_w.f90 (modified) (4 diffs)
SOURCE/exchange_horiz.f90 (modified) (3 diffs)
SOURCE/fft_xy.f90 (modified) (9 diffs)
SOURCE/flow_statistics.f90 (modified) (30 diffs)
SOURCE/modules.f90 (modified) (2 diffs)
SOURCE/prandtl_fluxes.f90 (modified) (13 diffs)
SOURCE/pres.f90 (modified) (10 diffs)
SOURCE/production_e.f90 (modified) (13 diffs)
SOURCE/prognostic_equations.f90 (modified) (6 diffs)
SOURCE/run_control.f90 (modified) (2 diffs)
SOURCE/time_integration.f90 (modified) (2 diffs)
SOURCE/timestep.f90 (modified) (7 diffs)
SOURCE/transpose.f90 (modified) (19 diffs)
SOURCE/tridia_solver.f90 (modified) (10 diffs)
SOURCE/wall_fluxes.f90 (modified) (3 diffs)

Legend:

: Unmodified
: Added
: Removed

palm/trunk/INSTALL/example_cbl_rc

-                      r1182
+                      r1257
  ******************************      ------------------------------------------
  * PALM 3.9  Rev: 1179        *      atmosphere - 3D - run without 1D - prerun
+ * PALM 3.9  Rev: 1255M       *      atmosphere - 3D - run without 1D - prerun
  ******************************      ------------------------------------------
  Date:                 14-06-13      Run:       example_cbl
  Time:                 09:15:07      Run-No.:   00
  Run on host:            lcsgih
+ Date:                 08-11-13      Run:       example_cbl
+ Time:                 16:08:04      Run-No.:   00
+ Run on host:           lccrayb
  Number of PEs:               8      Processor grid (x,y): (  2,  4) calculated
  ------------------------------------------------------------------------------
 …
 ------------------
+RUN  ITER. HH:MM:SS.SS    DT(E)     UMAX     VMAX     WMAX     U*    W*   THETA*    Z_I     ENERG.   DISTENERG    DIVOLD     DIVNEW     UMAX(KJI)    VMAX(KJI)    WMAX(KJI)   ADVECX   ADVECY   MGCYC
+-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+      0 00:00:00.00  20.0000A  -0.2192D -0.2169D -0.1214  0.000  1.38   0.000    800.  0.106E-02  0.000E+00  0.000E+00  0.000E+00   11  32   5    5  37  29    4  26   0    0.000    0.000       0
+      1 00:00:20.00  20.0000A  -0.2192  -0.2169  -0.1214  0.001  0.69  ******    100.  0.105E-02  0.948E-03  0.491E-05  0.355E-15   11  32   5    5  37  29    4  26   0    0.000    0.000       0
+      2 00:00:40.00  20.0000D  -0.2151  -0.2176  -0.1209  0.001  0.69  ******    100.  0.104E-02  0.940E-03  0.802E-05  0.359E-15   11  32   5    5  37  29    4  26   0    0.000    0.000       0
+      3 00:01:00.00  20.0000D  -0.2107  -0.2178  -0.1204  0.001  0.69  ******    100.  0.103E-02  0.932E-03  0.111E-04  0.363E-15   11  32   5    5  37  29    4  26   0    0.000    0.000       0
+      4 00:01:20.00  20.0000D  -0.2060  -0.2177  -0.1199  0.001  0.69  ******    100.  0.102E-02  0.924E-03  0.142E-04  0.366E-15   11  32   5    4   5   2    4  26   0    0.000    0.000       0
+      5 00:01:40.00  20.0000D  -0.2011  -0.2184  -0.1193  0.001  0.93  ******    250.  0.101E-02  0.916E-03  0.173E-04  0.368E-15   11  32   5    4   5   2    4  26   0    0.000    0.000       0
+      6 00:02:00.00  20.0000D  -0.1959  -0.2186  -0.1188  0.001  0.93  ******    250.  0.100E-02  0.908E-03  0.204E-04  0.370E-15   11  32   5    4   5   2    4  26   0    0.000    0.000       0
+      7 00:02:20.00  20.0000D  -0.1942  -0.2184  -0.1182  0.001  1.22  ******    550.  0.994E-03  0.900E-03  0.235E-04  0.371E-15   11  32   4    4   5   2    4  26   0    0.000    0.000       0
+      8 00:02:40.00  20.0000D  -0.1948D -0.2177D -0.1175  0.002  1.22  ******    550.  0.204E-02  0.892E-03  0.516E-03  0.735E-13   11  32   4    4   5   2    4  26   0    0.000    0.000       0
+      9 00:03:00.00  20.0000D   0.3635  -0.2778  -0.1533  0.002  1.22  ******    550.  0.201E-02  0.182E-02  0.321E-04  0.742E-15    4  35   4    7  14  12    5  10  39    0.000    0.000       0
+     10 00:03:20.00  20.0000D   0.3563  -0.2763  -0.1517  0.002  1.22  ******    550.  0.198E-02  0.180E-02  0.351E-04  0.739E-15    4  35   4    7  14  12    5  10  39    0.000    0.000       0
+     11 00:03:40.00  20.0000D   0.3475  -0.2747  -0.1498  0.002  1.22  ******    550.  0.195E-02  0.177E-02  0.380E-04  0.739E-15    4  35   4    6  29  30    5  10  39    0.000    0.000       0
+     12 00:04:00.00  20.0000D   0.3374  -0.2747  -0.1476  0.002  1.22  ******    550.  0.193E-02  0.175E-02  0.410E-04  0.741E-15    4  35   4    6  29  30    5  10  39    0.000    0.000       0
+     13 00:04:20.00  20.0000D   0.3262  -0.2735  -0.1451  0.002  1.22  ******    550.  0.190E-02  0.173E-02  0.439E-04  0.744E-15    4  35   4    6  29  30    5  10  39    0.000    0.000       0
+     14 00:04:40.00  20.0000D   0.3249  -0.2712  -0.1455  0.002  1.22  ******    550.  0.188E-02  0.171E-02  0.468E-04  0.749E-15    4  35   5    6  29  30    6  23   4    0.000    0.000       0
+     15 00:05:00.00  20.0000D   0.3275D -0.2694D -0.1470  0.002  1.22  ******    550.  0.293E-02  0.169E-02  0.529E-03  0.720E-13    4  35   5    7  13  38    6  23   4    0.000    0.000       0
+     16 00:05:20.00  20.0000D   0.3782  -0.3667  -0.2151  0.003  1.22  ******    550.  0.288E-02  0.262E-02  0.546E-04  0.110E-14   11  16  39    7  13  12    8  21  26    0.000    0.000       0
+     17 00:05:40.00  20.0000D   0.3741  -0.3587  -0.2156  0.003  1.22  ******    550.  0.284E-02  0.258E-02  0.573E-04  0.110E-14   11  16  39    7  15  26    8  21  26    0.000    0.000       0
+     18 00:06:00.00  20.0000D   0.3717  -0.3525  -0.2156  0.003  1.22  ******    550.  0.279E-02  0.255E-02  0.600E-04  0.111E-14    6   3  26    7  15  26    8  21  26    0.000    0.000       0
+     19 00:06:20.00  20.0000D  -0.3700  -0.3448  -0.2153  0.003  1.22  ******    550.  0.276E-02  0.251E-02  0.626E-04  0.112E-14    6  12  31    7  15  26    8  21  26    0.000    0.000       0
+     20 00:06:40.00  20.0000D  -0.3674  -0.3439  -0.2146  0.004  1.22  ******    550.  0.273E-02  0.249E-02  0.652E-04  0.114E-14    6  12  31    8  34  19    8  21  26    0.000    0.000       0
+     21 00:07:00.00  20.0000D  -0.3632  -0.3493  -0.2134  0.004  1.22  ******    550.  0.270E-02  0.246E-02  0.677E-04  0.116E-14    6  12  31    7  14  26    8  21  26    0.000    0.000       0
+     22 00:07:20.00  20.0000D  -0.3577  -0.3559   0.2292  0.005  1.22  ******    550.  0.269E-02  0.245E-02  0.701E-04  0.120E-14    6  12  31    7  14  26    2  15  28    0.000    0.000       0
+     23 00:07:40.00  20.0000D  -0.3509D -0.3603D  0.2560  0.005  1.22  ******    550.  0.377E-02  0.244E-02  0.552E-03  0.750E-13    6  12  31    7  14  26    2  15  28    0.000    0.000       0
+     24 00:08:00.00  20.0000D   0.4111   0.4185   0.3167  0.006  1.18  ******    500.  0.374E-02  0.341E-02  0.770E-04  0.163E-14    8  36  11   10  34  34    2  15  28    0.000    0.000       0
+     25 00:08:20.00  20.0000D  -0.4110   0.4103   0.3500  0.007  1.18  ******    500.  0.374E-02  0.340E-02  0.792E-04  0.168E-14    5  15  27   10  34  34    2  15  28    0.000    0.000       0
+     26 00:08:40.00  20.0000D  -0.4162   0.3992   0.3893  0.008  1.18  ******    500.  0.375E-02  0.342E-02  0.812E-04  0.175E-14    5  15  27   10  34  34    2  15  28    0.000    0.000       0
+     27 00:09:00.00  20.0000D  -0.4207   0.3858   0.4347  0.009  1.18  ******    500.  0.381E-02  0.347E-02  0.832E-04  0.184E-14    5  15  27   10  34  34    2  15  28    0.000    0.000       0
+     28 00:09:20.00  20.0000D  -0.4246  -0.3825   0.4864  0.010  1.18  ******    500.  0.390E-02  0.355E-02  0.850E-04  0.195E-14    5  15  27    9  21  22    2  15  28    0.000    0.000       0
+     29 00:09:40.00  20.0000D  -0.4283  -0.3813   0.5443  0.011  1.18  ******    500.  0.405E-02  0.368E-02  0.866E-04  0.209E-14    5  15  27    9  21  22    2  15  28    0.000    0.000       0
+     30 00:10:00.00  20.0000D  -0.4322D -0.3772D  0.6082  0.013  1.18  ******    500.  0.530E-02  0.387E-02  0.555E-03  0.730E-13    5  15  27    9  21  22    2  15  28    0.000    0.000       0
+     31 00:10:20.00  20.0000D  -0.4451   0.4333   0.6835  0.014  1.18  ******    500.  0.557E-02  0.505E-02  0.911E-04  0.276E-14    6  11  30    9  20  25    3  14  27    0.000    0.000       0
+     32 00:10:40.00  20.0000D  -0.4498   0.4329   0.7743  0.016  1.18  -9.947    500.  0.598E-02  0.541E-02  0.920E-04  0.300E-14    6  11  30    9  20  25    3  14  27    0.000    0.000       0
+     33 00:11:00.00  20.0000D  -0.4499   0.4284   0.8794  0.018  1.18  -8.653    500.  0.655E-02  0.591E-02  0.925E-04  0.329E-14    6  11  30    9  20  25    3  15  28    0.000    0.000       0
+     34 00:11:20.00  20.0000D  -0.4686   0.4203   1.0077  0.020  1.18  -7.700    500.  0.733E-02  0.660E-02  0.926E-04  0.365E-14    1  15  30    9  20  25    3  15  28    0.000    0.000       0
+     35 00:11:40.00  20.0000D  -0.5131   0.4390   1.1411  0.023  1.18  -6.766    500.  0.839E-02  0.754E-02  0.924E-04  0.407E-14    5  15  26    1  13  27    3  15  28    0.000    0.000       0
+     36 00:12:00.00  20.0000A  -0.5579   0.4874   1.2860  0.026  1.18  -6.023    500.  0.980E-02  0.879E-02  0.919E-04  0.459E-14    1  14  30    1  17  36    3  20  36    0.000    0.000       0
+     37 00:12:20.00  20.0000A  -0.6069   0.5356   1.4354  0.028  1.18  -5.484    500.  0.117E-01  0.104E-01  0.913E-04  0.519E-14    1  14  30    1  17  36    3  19  36    0.000    0.000       0
+     38 00:12:40.00  20.0000A  -0.6537   0.5816   1.6319  0.031  1.05  -5.117    350.  0.141E-01  0.126E-01  0.907E-04  0.591E-14    1  14  30    1  17  36    4  20  36    0.000    0.000       0
+     39 00:13:00.00  20.0000A  -0.6962   0.6231   1.8563  0.034  1.18  -4.484    500.  0.171E-01  0.153E-01  0.907E-04  0.675E-14    1  14  30    1  17  36    4  20  36    0.000    0.000       0
+     40 00:13:20.00  20.0000A  -0.7336   0.6662   2.0403  0.038  1.09  -4.142    400.  0.210E-01  0.188E-01  0.916E-04  0.772E-14    1  16  29    1  29  28    4  19  36    0.000    0.000       0
+     41 00:13:40.00  20.0000A  -0.7826   0.7271   2.2419  0.041  1.09  -3.852    400.  0.259E-01  0.232E-01  0.939E-04  0.884E-14    1  16  29    1  20  22    5  19  36    0.000    0.000       0
+     42 00:13:58.00  18.0000A  -0.8274   0.7915   2.4938  0.044  1.14  -3.465    450.  0.312E-01  0.283E-01  0.874E-04  0.902E-14    1  15  29    1  20  22    5  19  36    0.000    0.000       0
+     43 00:14:14.70  16.7000A  -0.8825   0.8465   2.6922  0.046  1.14  -3.273    450.  0.370E-01  0.339E-01  0.846E-04  0.939E-14    1   2  23    1  20  22    5  19  36    0.000    0.000       0
+     44 00:14:30.20  15.5000A  -0.9373   0.8934   2.8940  0.048  1.18  -3.114    500.  0.432E-01  0.399E-01  0.822E-04  0.968E-14    1   2  23    1  20  22    6  19  36    0.000    0.000       0
+     45 00:14:44.70  14.5000A  -0.9840   0.9319   3.1005  0.050  1.18  -3.023    500.  0.497E-01  0.462E-01  0.810E-04  0.997E-14    1   2  23    1  20  22    6  19  36    0.000    0.000       0
+     46 00:14:58.50  13.8000A  -1.0209   0.9629   3.2577  0.052  1.22  -3.169    550.  0.566E-01  0.529E-01  0.813E-04  0.104E-13    1   2  23    1  20  22    6  19  36    0.000    0.000       0
+     47 00:15:11.50  13.0000A   1.0590   0.9874   3.4608  0.053  1.22  -2.944    550.  0.637E-01  0.599E-01  0.809E-04  0.106E-13    1  31  18    1  20  22    7  19  36    0.000    0.000       0
+     48 00:15:23.90  12.4000A   1.1029   1.0076   3.6147  0.055  1.28  -2.874    650.  0.711E-01  0.671E-01  0.815E-04  0.109E-13    1  31  18    1  29  28    7  19  36    0.000    0.000       0
+     49 00:15:35.90  12.0000A   1.1393   1.0350   3.7578  0.056  1.28  -2.769    650.  0.787E-01  0.746E-01  0.832E-04  0.113E-13    1  31  18    1  20  23    8  19  36    0.000    0.000       0
+     50 00:15:47.40  11.5000A   1.1689   1.0620   3.9096  0.057  1.29  -2.648    650.  0.865E-01  0.822E-01  0.842E-04  0.116E-13    1  31  18    1  20  23    8  19  36    0.000    0.000       0
+     51 00:15:58.70  11.3000A   1.1914   1.0868   3.9957  0.058  1.29  -2.546    650.  0.945E-01  0.900E-01  0.870E-04  0.121E-13    1  31  18    1  20  23    9  19  36    0.000    0.000       0
+     52 00:16:09.60  10.9000A   1.2072   1.1100   4.1423  0.059  1.29  -2.467    650.  0.103E+00  0.980E-01  0.882E-04  0.124E-13    1  31  18    1  20  23    9  19  36    0.000    0.000       0
+     53 00:16:20.30  10.7000A   1.2163   1.1311   4.2098  0.059  1.38  -2.410    800.  0.111E+00  0.106E+00  0.907E-04  0.129E-13    1  31  18    1  20  23    9  19  36    0.000    0.000       0
+     54 00:16:30.70  10.4000A   1.2230   1.1501   4.3170  0.060  1.38  -2.387    800.  0.119E+00  0.114E+00  0.923E-04  0.132E-13    1  30  18    1  20  23   10  19  36    0.000    0.000       0
+     55 00:16:40***  10.3000A   1.2285   1.1668   4.3818  0.061  1.38  -2.372    800.  0.127E+00  0.122E+00  0.954E-04  0.137E-13    1  30  18    1  20  23   10  19  36    0.000    0.000       0
+     56 00:16:51.10  10.1000A   1.2378   1.2243   4.4452  0.061  1.38  -2.372    800.  0.136E+00  0.131E+00  0.977E-04  0.140E-13    1  37  15    7  12   8   11  19  36    0.000    0.000       0
+     57 00:17:01.09   9.9900A   1.2531   1.2263   4.5048  0.062  1.38  -2.413    800.  0.144E+00  0.139E+00  0.101E-03  0.145E-13    1  37  15    6  20  34   11  19  36    0.000    0.000       0
+     58 00:17:11.04   9.9500A  -1.3217   1.2613   4.5240  0.062  1.38  -2.415    800.  0.153E+00  0.147E+00  0.104E-03  0.150E-13    4  35   3   10  17   2   12  19  36    0.000    0.000       0
+     59 00:17:20.87   9.8300A  -1.3776   1.3448   4.5770  0.063  1.38  -2.386    800.  0.161E+00  0.156E+00  0.107E-03  0.155E-13    4  35   3   10  17   2   12  19  36    0.000    0.000       0
+     60 00:17:30.73   9.8600A  -1.4186   1.2882   4.5631  0.063  1.38  -2.396    800.  0.170E+00  0.164E+00  0.111E-03  0.161E-13    4  35   3    9  16   3   13  19  36    0.000    0.000       0
+     61 00:17:40.50   9.7700A  -1.4466   1.4232   4.6056  0.064  1.41  -2.426    850.  0.178E+00  0.173E+00  0.114E-03  0.166E-13    4  35   3   11  17   2   13  19  36    0.000    0.000       0
+     62 00:17:50.32   9.8200A  -1.4630   1.4757   4.5817  0.064  1.41  -2.506    850.  0.186E+00  0.181E+00  0.119E-03  0.173E-13    4  35   3   11  17   2   13  19  36    0.000    0.000       0
+     63 00:18:00.12   9.8000A  -1.4689  -1.4758   4.5924  0.064  1.41  -2.441    850.  0.195E+00  0.189E+00  0.122E-03  0.178E-13    4  35   3    9   4  39   14  19  36    0.000    0.000       0
+     64 00:18:09.98   9.8600A  -1.4900   1.5484   4.5623  0.065  1.41  -2.537    850.  0.203E+00  0.198E+00  0.127E-03  0.185E-13    4  35   2   12  17   2   14  19  36    0.000    0.000       0
+     65 00:18:19.93   9.9500A  -1.5149   1.5649   4.5234  0.065  1.41  -2.404    850.  0.211E+00  0.206E+00  0.132E-03  0.193E-13    5  35   2   12  17   2   15  19  36    0.000    0.000       0
+     66 00:18:29.90   9.9700A  -1.5562  -1.6149   4.5130  0.066  1.41  -2.297    850.  0.219E+00  0.214E+00  0.136E-03  0.199E-13    5  35   2   10   4  39   15  15  27    0.000    0.000       0
+     67 00:18:39.84   9.9400A  -1.5803  -1.6213   4.5279  0.066  1.41  -2.265    850.  0.227E+00  0.221E+00  0.140E-03  0.205E-13    5  35   2   10   4  39   12  15   1    0.000    0.000       0
+     68 00:18:49.80   9.9600A  -1.5944  -1.5975   4.5161  0.067  1.41  -2.220    850.  0.234E+00  0.229E+00  0.144E-03  0.211E-13    5  35   2    9   5  38   12  15   1    0.000    0.000       0
+     69 00:18:59.73   9.9300A  -1.6046  -1.6886   4.5325  0.067  1.41  -2.209    850.  0.241E+00  0.236E+00  0.148E-03  0.216E-13    5  35   2    9   5  38   13  15   1    0.000    0.000       0
+     70 00:19:09.63   9.9000A  -1.6147  -1.8090   4.5451  0.067  1.41  -2.253    850.  0.247E+00  0.242E+00  0.151E-03  0.221E-13    5  35   2   10   4  38   13  15   1    0.000    0.000       0
+     71 00:19:19.63  10.0000A  -1.6261  -1.8785   4.4972  0.068  1.41  -2.321    850.  0.253E+00  0.249E+00  0.156E-03  0.229E-13    5  35   2   10   4  38   13  15   1    0.000    0.000       0
+     72 00:19:29.63  10.0000A  -1.6387  -1.8839   4.4921  0.068  1.41  -2.413    850.  0.259E+00  0.254E+00  0.160E-03  0.234E-13    5  35   2   10   4  38   14  15   1    0.000    0.000       0
+     73 00:19:39.73  10.1000A  -1.6503  -1.9429   4.4673  0.069  1.41  -2.182    850.  0.264E+00  0.260E+00  0.166E-03  0.242E-13    5  35   2   11   4  38   14  15   1    0.000    0.000       0
+     74 00:19:49.93  10.2000A  -1.6583  -2.0326   4.3937  0.069  1.41  -2.154    850.  0.269E+00  0.265E+00  0.171E-03  0.249E-13    5  35   2   11   4  38   14  15   1    0.000    0.000       0
+     75 00:20:00.23  10.3000A  -1.7171  -2.0769   4.3835  0.070  1.41  -2.083    850.  0.273E+00  0.269E+00  0.177E-03  0.257E-13   12  36  39   12   3  38   12  35  25    0.000    0.000       0
+     76 00:20:10.43  10.2000A  -1.7567  -2.1408   4.4178  0.070  1.41  -2.097    850.  0.277E+00  0.273E+00  0.179E-03  0.259E-13   12  36  39   12   3  38   13  35  25    0.000    0.000       0
+     77 00:20:20.53  10.1000A  -1.7527  -2.1271   4.4670  0.071  1.41  -2.080    850.  0.280E+00  0.277E+00  0.180E-03  0.261E-13   12  36  39   12   3  38   13  35  25    0.000    0.000       0
+     78 00:20:30.73  10.2000A  -1.8161  -2.1682   4.4292  0.071  1.41  -2.213    850.  0.283E+00  0.280E+00  0.185E-03  0.267E-13   13  36  38   13   3  38   14  35  25    0.000    0.000       0
+     79 00:20:40.73  10.0000A  -1.8889  -2.1900   4.4840  0.072  1.41  -1.983    850.  0.285E+00  0.282E+00  0.185E-03  0.266E-13   13  36  38   13   3  38   14  35  25    0.000    0.000       0
+     80 00:20:50.83  10.1000A  -1.9297  -2.1973   4.4481  0.073  1.41  -1.942    850.  0.286E+00  0.284E+00  0.190E-03  0.272E-13   13  36  38   14   2  38   14  35  25    0.000    0.000       0
+     81 00:21:01.13  10.3000A  -1.9368  -2.2598   4.3901  0.073  1.41  -1.922    850.  0.288E+00  0.285E+00  0.196E-03  0.280E-13   13  36  38   14   2  38   15  35  25    0.000    0.000       0
+     82 00:21:11.43  10.3000A  -1.9076  -2.1711   4.3540  0.074  1.41  -1.909    850.  0.288E+00  0.286E+00  0.198E-03  0.284E-13   13  36  38   14   2  38   15  35  25    0.000    0.000       0
+     83 00:21:22.03  10.6000A  -1.9439  -2.1446   4.2332  0.075  1.41  -1.899    850.  0.288E+00  0.286E+00  0.206E-03  0.295E-13   14  36  38   15   2  38   15  35  25    0.000    0.000       0
+     84 00:21:33.13  11.1000A  -1.9810  -2.1166   4.0656  0.076  1.41  -1.899    850.  0.288E+00  0.286E+00  0.218E-03  0.311E-13   14  36  38   15   2  38   16  35  25    0.000    0.000       0
+     85 00:21:44.53  11.4000A  -1.9978  -2.0805   3.9358  0.076  1.41  -1.879    850.  0.288E+00  0.286E+00  0.226E-03  0.322E-13   14  36  38   15   1  38   16  35  25    0.000    0.000       0
+     86 00:21:56.53  12.0000A  -1.9911  -1.9954   3.7411  0.077  1.38  -1.876    800.  0.287E+00  0.285E+00  0.240E-03  0.342E-13   14  36  38   15   2  38   16  35  25    0.000    0.000       0
+     87 00:22:09.33  12.8000A  -1.9587  -2.0083   3.5118  0.078  1.38  -1.843    800.  0.286E+00  0.285E+00  0.258E-03  0.367E-13   14  36  38   17  38   9   12   6  38    0.000    0.000       0
+     88 00:22:22.23  12.9000A  -1.8968  -2.0019   3.4815  0.079  1.38  -1.817    800.  0.285E+00  0.283E+00  0.263E-03  0.372E-13   14  36  38   17  38   9   13  20   6    0.000    0.000       0
+     89 00:22:35.43  13.2000A  -1.8835  -1.9266   3.4165  0.081  1.38  -1.817    800.  0.283E+00  0.282E+00  0.271E-03  0.382E-13   15  35  38   17  38   9   14  36  16    0.000    0.000       0
+     90 00:22:48.73  13.3000A  -1.8964  -1.9624   3.3714  0.082  1.38  -1.773    800.  0.281E+00  0.281E+00  0.275E-03  0.386E-13   16  35  37   10  14  36   14  36  16    0.000    0.000       0
+     91 00:23:02.53  13.8000A  -1.9119  -1.7922   3.2564  0.083  1.38  -1.762    800.  0.280E+00  0.279E+00  0.286E-03  0.401E-13   16  35  36   10  14  36   14  36  14    0.000    0.000       0
+     92 00:23:16.73  14.2000A  -1.9743  -1.7683   3.1765  0.084  1.35  -1.706    750.  0.277E+00  0.277E+00  0.293E-03  0.412E-13   16  35  36   11  14  36   14  19   6    0.000    0.000       0
+     93 00:23:31.43  14.7000A  -1.9862  -1.8340   3.0522  0.086  1.35  -1.664    750.  0.275E+00  0.275E+00  0.303E-03  0.426E-13   16  35  36   11  14  36   14  19   6    0.000    0.000       0
+     94 00:23:46.43  15.0000A  -1.9412  -1.6691   2.9938  0.087  1.35  -1.628    750.  0.272E+00  0.273E+00  0.308E-03  0.433E-13   16  35  36   11  14  36   16  18   6    0.000    0.000       0
+     95 00:24:01.33  14.9000A  -1.8521  -1.7158   3.0103  0.088  1.35  -1.592    750.  0.269E+00  0.270E+00  0.303E-03  0.428E-13   16  35  35   11  14  32   11  22  10    0.000    0.000       0
+     96 00:24:16.13  14.8000A  -1.8021  -1.7470   3.0382  0.090  1.35  -1.565    750.  0.266E+00  0.267E+00  0.299E-03  0.421E-13    7  20  10   11  14  32   11  22  10    0.000    0.000       0
+     97 00:24:31.03  14.9000A  -1.8736  -1.7532   3.0227  0.091  1.32  -1.528    700.  0.263E+00  0.264E+00  0.297E-03  0.419E-13    8  20  10    1  22  29   12  22  10    0.000    0.000       0
+     98 00:24:45.93  14.9000A  -1.8942  -1.7566   3.0293  0.092  1.32  -1.486    700.  0.259E+00  0.261E+00  0.295E-03  0.414E-13    8  20  10    1  22  29   12  22  10    0.000    0.000       0
+     99 00:25:00.93  15.0000A  -1.9088   1.7638   3.0033  0.094  1.32  -1.459    700.  0.256E+00  0.258E+00  0.294E-03  0.411E-13    8  20   9    1   7  14   12  22  10    0.000    0.000       0
+    100 00:25:15.93  15.0000A  -1.9789   1.8127   3.0007  0.095  1.32  -1.474    700.  0.252E+00  0.255E+00  0.292E-03  0.406E-13    8  20   9    1   7  14   13  22  10    0.000    0.000       0
+    101 00:25:31.03  15.1000A  -1.9401   1.8502   2.9834  0.096  1.32  -1.500    700.  0.249E+00  0.252E+00  0.292E-03  0.404E-13    8  20   9    1   7  14   13  22  10    0.000    0.000       0
+    102 00:25:46.43  15.4000A  -1.8588   1.8909   2.9152  0.097  1.32  -1.452    700.  0.245E+00  0.249E+00  0.295E-03  0.406E-13    9  20   9    1   8  13   13  22  10    0.000    0.000       0
+    103 00:26:02.03  15.6000A   1.8657   1.9566   2.8848  0.098  1.32  -1.442    700.  0.242E+00  0.245E+00  0.297E-03  0.405E-13    8   0  24    1   8  13   14  22  10    0.000    0.000       0
+    104 00:26:17.93  15.9000A   1.8880   2.0113   2.8223  0.099  1.32  -1.405    700.  0.238E+00  0.242E+00  0.300E-03  0.408E-13    9   0  24    1   8  14   14  22  10    0.000    0.000       0
+    105 00:26:34.23  16.3000A   1.8023   2.0555   2.7619  0.101  1.32  -1.369    700.  0.234E+00  0.239E+00  0.305E-03  0.412E-13    9   0  24    1   8  14    8  10  12    0.000    0.000       0
+    106 00:26:50.63  16.4000A   1.7683   2.0830   2.7515  0.101  1.29  -1.369    650.  0.230E+00  0.235E+00  0.304E-03  0.409E-13   10   0  24    1   8  14   12   8  17    0.000    0.000       0
+    107 00:27:07.23  16.6000A   1.6965   2.0957   2.7134  0.102  1.29  -1.395    650.  0.227E+00  0.232E+00  0.304E-03  0.408E-13   10   0  24    1   8  14   12   8  17    0.000    0.000       0
+    108 00:27:24.23  17.0000A  -1.6847   2.1450   2.6428  0.103  1.29  -1.555    650.  0.223E+00  0.229E+00  0.308E-03  0.412E-13    8  37  33    1   9  14   12   8  17    0.000    0.000       0
+    109 00:27:41.23  17.0000A   1.7150   2.1839   2.6486  0.103  1.32  -1.381    700.  0.219E+00  0.225E+00  0.304E-03  0.405E-13   10  10   4    1   9  14    9  10  11    0.000    0.000       0
+    110 00:27:58.13  16.9000A  -1.7298   2.2080   2.6616  0.103  1.32  -1.379    700.  0.216E+00  0.222E+00  0.299E-03  0.397E-13    5  36  26    1   9  14    9  10  11    0.000    0.000       0
+    111 00:28:15.33  17.2000A  -1.7799   2.2247   2.6237  0.103  1.35  -1.410    750.  0.213E+00  0.219E+00  0.300E-03  0.397E-13    5  36  26    1   9  14    9   9  18    0.000    0.000       0
+    112 00:28:32.13  16.8000A  -1.7880   2.2310   2.6717  0.103  1.35  -1.375    750.  0.210E+00  0.216E+00  0.290E-03  0.382E-13    7  35  25    1   9  14    9   9  18    0.000    0.000       0
+    113 00:28:49.33  17.2000A  -1.8509   2.2187   2.6171  0.103  1.35  -1.382    750.  0.207E+00  0.213E+00  0.293E-03  0.385E-13    7  35  25    1  10  14   10  11  11    0.000    0.000       0
+    114 00:29:06.73  17.4000A  -1.8885   2.2092   2.5798  0.102  1.35  -1.375    750.  0.204E+00  0.210E+00  0.292E-03  0.383E-13    6  35  25    1  10  14   10  11  11    0.000    0.000       0
+    115 00:29:24.83  18.1000A  -1.9577   2.1622   2.4896  0.102  1.35  -1.370    750.  0.201E+00  0.207E+00  0.300E-03  0.391E-13    6  35  25    1  10  14   10  11  11    0.000    0.000       0
+    116 00:29:42.83  18.0000A  -1.9567   2.1160   2.4985  0.101  1.35  -1.352    750.  0.199E+00  0.205E+00  0.295E-03  0.382E-13    6  35  25    1  10  15   13  35   7    0.000    0.000       0
+    117 00:30:01.13  18.3000A  -1.8608   2.0695   2.4612  0.100  1.35  -1.382    750.  0.197E+00  0.202E+00  0.296E-03  0.381E-13    7  35  24    1  10  16    9  11   0    0.000    0.000       0
+    118 00:30:19.23  18.1000A  -1.8672   2.0382   2.4803  0.099  1.35  -1.383    750.  0.195E+00  0.200E+00  0.289E-03  0.370E-13    5  35  25    1  10  16    9   0  35    0.000    0.000       0
+    119 00:30:36.73  17.5000A  -1.8222   1.9954   2.5678  0.099  1.38  -1.372    800.  0.193E+00  0.198E+00  0.276E-03  0.351E-13    5  35  25    1  10  16    9   0  35    0.000    0.000       0
+    120 00:30:53.93  17.2000A   1.8055   1.9528   2.6115  0.098  1.38  -1.369    800.  0.192E+00  0.197E+00  0.268E-03  0.339E-13    1  22   8    1  10  17    9   0  35    0.000    0.000       0
+    121 00:31:10.63  16.7000A   1.8094   1.9314   2.6918  0.097  1.38  -1.355    800.  0.191E+00  0.195E+00  0.258E-03  0.324E-13    1  22   8    1  10  17   10   0  35    0.000    0.000       0
+    122 00:31:27.03  16.4000A   1.7928   1.8929   2.7468  0.096  1.38  -1.374    800.  0.190E+00  0.194E+00  0.251E-03  0.313E-13    1  22   8    1  10  17   10   0  35    0.000    0.000       0
+    123 00:31:43.33  16.3000A   1.7763   1.9108   2.7602  0.095  1.38  -1.400    800.  0.190E+00  0.194E+00  0.248E-03  0.308E-13    1  22   9    1  10  11   10   0  35    0.000    0.000       0
+    124 00:31:59.63  16.3000A   1.7786   1.9160   2.7662  0.095  1.38  -1.421    800.  0.189E+00  0.193E+00  0.247E-03  0.305E-13    1  22   9    1  10  11   11   0  35    0.000    0.000       0
+    125 00:32:15.93  16.3000A   1.7607   1.9153   2.7645  0.094  1.38  -1.456    800.  0.189E+00  0.193E+00  0.247E-03  0.303E-13    1  22   9    1  11  11   11   0  35    0.000    0.000       0
+    126 00:32:32.23  16.3000A   1.7319   1.9441   2.7583  0.094  1.38  -1.469    800.  0.189E+00  0.193E+00  0.247E-03  0.302E-13    1  22   9    1  11  11   11  39  35    0.000    0.000       0
+    127 00:32:48.23  16.0000A   1.7104   1.9380   2.8098  0.093  1.38  -1.478    800.  0.189E+00  0.193E+00  0.243E-03  0.296E-13    1  21   9    1  11  11   12  39  35    0.000    0.000       0
+    128 00:33:04.03  15.8000A   1.6848   1.9166   2.8413  0.093  1.38  -1.473    800.  0.190E+00  0.193E+00  0.240E-03  0.291E-13    1  21   9    1  12  11   12  39  35    0.000    0.000       0
+    129 00:33:19.93  15.9000A   1.6595   1.9415   2.8385  0.092  1.38  -1.481    800.  0.190E+00  0.194E+00  0.242E-03  0.293E-13    1  22   9    1  12  11   12  39  35    0.000    0.000       0
+    130 00:33:35.93  16.0000A   1.6488   1.9224   2.8124  0.092  1.38  -1.487    800.  0.191E+00  0.194E+00  0.244E-03  0.295E-13    1  22   9    1  12  11   13  39  35    0.000    0.000       0
+    131 00:33:52.03  16.1000A   1.6413   1.8629   2.8018  0.091  1.35  -1.507    750.  0.192E+00  0.195E+00  0.246E-03  0.297E-13    1  22   9    1  12  11   12  39  38    0.000    0.000       0
+    132 00:34:07.93  15.9000A   1.6318   1.8526   2.8335  0.091  1.35  -1.539    750.  0.192E+00  0.196E+00  0.243E-03  0.294E-13    1  22   9    1  13  11   13  39  38    0.000    0.000       0
+    133 00:34:23.23  15.3000A   1.6230   1.8376   2.9328  0.091  1.35  -1.505    750.  0.193E+00  0.196E+00  0.235E-03  0.283E-13    7  28  15    1  13  11   10  25  21    0.000    0.000       0
+    134 00:34:38.33  15.1000A   1.6344   1.8336   2.9819  0.090  1.35  -1.509    750.  0.194E+00  0.197E+00  0.233E-03  0.280E-13    7  28  15    1  13  10   11  25  21    0.000    0.000       0
+    135 00:34:53.53  15.2000A   1.5790   1.8780   2.9665  0.090  1.35  -1.519    750.  0.194E+00  0.197E+00  0.235E-03  0.283E-13    8  28  15    5  28  13   11  25  21    0.000    0.000       0
+    136 00:35:08.93  15.4000A   1.6237   1.8794   2.9163  0.089  1.41  -1.516    850.  0.195E+00  0.198E+00  0.239E-03  0.287E-13    7  29  15    5  28  13   12  25  22    0.000    0.000       0
+    137 00:35:23.93  15.0000A   1.6603   1.8668   3.0042  0.089  1.41  -1.506    850.  0.196E+00  0.199E+00  0.234E-03  0.281E-13    7  29  15    6  29  14   13  24  21    0.000    0.000       0
+    138 00:35:38.63  14.7000A   1.6174   1.8698   3.0587  0.089  1.41  -1.485    850.  0.196E+00  0.199E+00  0.230E-03  0.276E-13    7  29  15    6  29  14   14  24  21    0.000    0.000       0
+    139 00:35:53.13  14.5000A   1.5861   1.8785   3.0967  0.089  1.41  -1.482    850.  0.197E+00  0.200E+00  0.228E-03  0.274E-13    7  30  16    6  29  13   14  24  21    0.000    0.000       0
+    140 00:36:07.83  14.7000A   1.6144   1.9207   3.0703  0.088  1.41  -1.543    850.  0.197E+00  0.200E+00  0.232E-03  0.278E-13    7  30  16    5  29  12   14  24  22    0.000    0.000       0
+    141 00:36:22.83  15.0000A  -1.5908   1.9489   2.9967  0.088  1.41  -1.507    850.  0.198E+00  0.201E+00  0.238E-03  0.285E-13    5  10  18    5  29  12   15  24  22    0.000    0.000       0
+    142 00:36:38.13  15.3000A  -1.6124   1.9034   2.9318  0.088  1.41  -1.554    850.  0.199E+00  0.201E+00  0.244E-03  0.292E-13    5  10  18    5  29  12   15  24  22    0.000    0.000       0
+    143 00:36:53.23  15.1000A   1.6289   1.8495   2.9798  0.088  1.41  -1.558    850.  0.199E+00  0.202E+00  0.242E-03  0.289E-13   10   2   1    6  29  13   14  24  23    0.000    0.000       0
+    144 00:37:08.03  14.8000A   1.6782   1.8285   3.0490  0.088  1.41  -1.577    850.  0.200E+00  0.202E+00  0.238E-03  0.283E-13   10   2   1    6  29  12   15  24  23    0.000    0.000       0
+    145 00:37:23.23  15.2000A   1.6806   1.8179   2.9608  0.088  1.41  -1.571    850.  0.201E+00  0.203E+00  0.245E-03  0.291E-13   10   2   1    6  29  12   15  24  23    0.000    0.000       0
+    146 00:37:38.53  15.3000A   1.7344   1.8072   2.9365  0.088  1.41  -1.517    850.  0.201E+00  0.204E+00  0.248E-03  0.293E-13   10   2   2    6  30  12   12  26  16    0.000    0.000       0
+    147 00:37:53.83  15.3000A   1.7740   1.7768   2.9393  0.088  1.41  -1.500    850.  0.202E+00  0.204E+00  0.249E-03  0.294E-13   10   2   2    6  30  12   12  26  16    0.000    0.000       0
+    148 00:38:09.23  15.4000A   1.7807   1.6732   2.9181  0.088  1.41  -1.498    850.  0.203E+00  0.205E+00  0.251E-03  0.296E-13   10   2   2   15  31  16   12  26  16    0.000    0.000       0
+    149 00:38:24.73  15.5000A  -1.7681   1.7137   2.8972  0.088  1.41  -1.514    850.  0.203E+00  0.205E+00  0.254E-03  0.298E-13    6  10  17   15  31  16   13  26  16    0.000    0.000       0
+    150 00:38:40.43  15.7000A  -1.7616   1.7355   2.8623  0.088  1.38  -1.614    800.  0.204E+00  0.206E+00  0.259E-03  0.303E-13    6  11  17   15  31  16   13  26  16    0.000    0.000       0
+    151 00:38:56.53  16.1000A  -1.7349   1.7326   2.7955  0.089  1.38  -1.560    800.  0.205E+00  0.207E+00  0.267E-03  0.311E-13    6  11  17   15  31  16   13  26  16    0.000    0.000       0
+    152 00:39:12.73  16.2000A  -1.7278   1.6948   2.7766  0.089  1.38  -1.548    800.  0.206E+00  0.208E+00  0.269E-03  0.313E-13    6  10  16   16  32  16   12  19  11    0.000    0.000       0
+    153 00:39:28.73  16.0000A  -1.6749   1.7246   2.8056  0.089  1.38  -1.511    800.  0.207E+00  0.209E+00  0.267E-03  0.309E-13    6  10  16   16  32  16    6  38  37    0.000    0.000       0
+    154 00:39:44.23  15.5000A   1.6902   1.7261   2.8993  0.090  1.38  -1.509    800.  0.208E+00  0.209E+00  0.259E-03  0.301E-13    7  13  11   16  32  16    7  38  37    0.000    0.000       0
+    155 00:39:59.63  15.4000A   1.6935   1.7031   2.9163  0.090  1.38  -1.515    800.  0.209E+00  0.210E+00  0.259E-03  0.299E-13    7  13  11   15  32  17    7  38  37    0.000    0.000       0
+    156 00:40:14.73  15.1000A  -1.7797   1.6675   2.9802  0.090  1.38  -1.514    800.  0.210E+00  0.211E+00  0.255E-03  0.295E-13    7  10  18   15  33  17    8  38  36    0.000    0.000       0
+    157 00:40:29.33  14.6000A  -1.8277   1.7187   3.0862  0.091  1.38  -1.504    800.  0.211E+00  0.212E+00  0.248E-03  0.287E-13    7  10  18   15  33  17    9  38  36    0.000    0.000       0
+    158 00:40:43.63  14.3000A  -1.8209   1.7025   3.1454  0.091  1.38  -1.520    800.  0.212E+00  0.214E+00  0.244E-03  0.283E-13    7  10  18   15  33  17    9  38  36    0.000    0.000       0
+    159 00:40:57.33  13.7000A  -1.8584   1.6482   3.2818  0.091  1.38  -1.479    800.  0.213E+00  0.215E+00  0.235E-03  0.273E-13    7  10  17   16  34  17   10  38  36    0.000    0.000       0
+    160 00:41:10.93  13.6000A  -1.8838   1.6713   3.3117  0.092  1.38  -1.454    800.  0.215E+00  0.216E+00  0.234E-03  0.272E-13    7  10  17   16  34  17   10  38  36    0.000    0.000       0
+    161 00:41:24.33  13.4000A  -1.8844   1.6480   3.3512  0.092  1.38  -1.443    800.  0.216E+00  0.217E+00  0.231E-03  0.269E-13    9   9  17   16  34  17   11  38  36    0.000    0.000       0
+    162 00:41:37.63  13.3000A  -1.9152   1.6080   3.3815  0.092  1.38  -1.452    800.  0.217E+00  0.218E+00  0.230E-03  0.267E-13    9   9  17   11  17  11   11  38  36    0.000    0.000       0
+    163 00:41:51.03  13.4000A  -1.9299   1.6385   3.3695  0.093  1.38  -1.480    800.  0.219E+00  0.220E+00  0.232E-03  0.269E-13    9   9  17   11  17  11   12  37  36    0.000    0.000       0
+    164 00:42:04.13  13.1000A  -1.9304   1.6459   3.4379  0.093  1.38  -1.566    800.  0.221E+00  0.222E+00  0.227E-03  0.263E-13    9   9  17   11  17  11   12  37  36    0.000    0.000       0
+    165 00:42:17.23  13.1000A  -1.9441   1.6894   3.4372  0.093  1.38  -1.531    800.  0.222E+00  0.223E+00  0.227E-03  0.262E-13    9   9  16    6  13  23   12  37  36    0.000    0.000       0
+    166 00:42:30.33  13.1000A  -1.9605   1.7224   3.4349  0.094  1.38  -1.517    800.  0.224E+00  0.225E+00  0.227E-03  0.261E-13    9   9  16    6  13  23   13  38  35    0.000    0.000       0
+    167 00:42:43.13  12.8000A  -1.9517   1.6409   3.5062  0.094  1.38  -1.480    800.  0.226E+00  0.227E+00  0.222E-03  0.255E-13    9   9  16    6  13  23   13  37  35    0.000    0.000       0
+    168 00:42:55.63  12.5000A  -1.9214   1.6858   3.6044  0.094  1.41  -1.498    850.  0.228E+00  0.229E+00  0.217E-03  0.249E-13    9   9  16    7  13  23   14  37  35    0.000    0.000       0
+    169 00:43:07.73  12.1000A  -1.9821   1.8072   3.7222  0.094  1.41  -1.483    850.  0.230E+00  0.230E+00  0.210E-03  0.241E-13   16  38  33    7  13  23   14  37  35    0.000    0.000       0
+    170 00:43:19.73  12.0000A  -1.9788   1.7642   3.7600  0.094  1.41  -1.437    850.  0.231E+00  0.232E+00  0.208E-03  0.238E-13   16  38  32    7  13  23   14  37  35    0.000    0.000       0
+    171 00:43:31.73  12.0000A  -1.9337   1.6513   3.7552  0.095  1.41  -1.432    850.  0.233E+00  0.234E+00  0.207E-03  0.238E-13   16  38  32    6  13  24   13  35  37    0.000    0.000       0
+    172 00:43:43.63  11.9000A  -1.9126   1.6512   3.7694  0.095  1.43  -1.442    900.  0.235E+00  0.236E+00  0.206E-03  0.236E-13   17  37  33    7  13  24   13  35  37    0.000    0.000       0
+    173 00:43:55.63  12.0000A  -1.9645   1.7137   3.7580  0.095  1.43  -1.579    900.  0.237E+00  0.238E+00  0.207E-03  0.238E-13   17  37  33    7  13  24   14  35  37    0.000    0.000       0
+    174 00:44:07.53  11.9000A  -1.9699   1.7579   3.7681  0.095  1.43  -1.479    900.  0.239E+00  0.240E+00  0.206E-03  0.236E-13   17  36  34   13  30  13   14  35  37    0.000    0.000       0
+    175 00:44:19.53  12.0000A  -1.9899  -1.7872   3.7442  0.095  1.43  -1.497    900.  0.241E+00  0.242E+00  0.207E-03  0.237E-13   17  36  34    7   3   0   14  35  37    0.000    0.000       0
+    176 00:44:31.73  12.2000A  -1.9782  -1.7649   3.6914  0.095  1.43  -1.508    900.  0.243E+00  0.244E+00  0.211E-03  0.241E-13   17  35  34    7   3   0   14  35  37    0.000    0.000       0
+    177 00:44:44.13  12.4000A  -2.1187  -1.7856   3.6148  0.095  1.43  -1.501    900.  0.246E+00  0.246E+00  0.215E-03  0.245E-13   18  36  33    7   2   0   14  35  37    0.000    0.000       0
+    178 00:44:56.93  12.8000A  -2.1643   1.7763   3.5217  0.095  1.43  -1.505    900.  0.248E+00  0.249E+00  0.222E-03  0.253E-13   18  36  33   14  30  13   14  35  37    0.000    0.000       0
+    179 00:45:10.13  13.2000A  -2.1899   1.8103   3.4188  0.095  1.43  -1.498    900.  0.250E+00  0.251E+00  0.230E-03  0.260E-13   18  36  32   14  30  13   14  35  37    0.000    0.000       0
+    180 00:45:23.23  13.1000A  -2.1199   1.8669   3.4290  0.095  1.43  -1.481    900.  0.252E+00  0.253E+00  0.229E-03  0.259E-13   18  35  33   14  31  13   12  27  19    0.000    0.000       0
+    181 00:45:36.13  12.9000A  -2.1411   1.8771   3.4781  0.094  1.43  -1.481    900.  0.255E+00  0.256E+00  0.226E-03  0.256E-13   18  35  33   14  31  13   12  27  19    0.000    0.000       0
+    182 00:45:48.93  12.8000A  -2.0012   1.8359   3.5123  0.094  1.41  -1.473    850.  0.257E+00  0.258E+00  0.226E-03  0.255E-13   18  35  33   14  31  13   13  27  19    0.000    0.000       0
+    183 00:46:01.63  12.7000A  -1.9681   1.7665   3.5453  0.094  1.41  -1.468    850.  0.259E+00  0.260E+00  0.225E-03  0.253E-13   19  36  32   15  31  12   13  27  19    0.000    0.000       0
+    184 00:46:14.33  12.7000A  -2.0050   1.8226   3.5420  0.095  1.41  -1.478    850.  0.262E+00  0.263E+00  0.226E-03  0.254E-13   19  36  31   15  31  13   13  27  19    0.000    0.000       0
+    185 00:46:27.03  12.7000A  -2.0360   1.8702   3.5337  0.095  1.41  -1.467    850.  0.264E+00  0.265E+00  0.227E-03  0.255E-13   19  36  31   15  32  12   11  27  17    0.000    0.000       0
+    186 00:46:39.43  12.4000A  -1.9794   1.8889   3.6395  0.095  1.41  -1.461    850.  0.266E+00  0.268E+00  0.223E-03  0.250E-13   19  34  33   15  32  12   12  27  17    0.000    0.000       0
+    187 00:46:51.53  12.1000A  -1.9829   1.8800   3.7238  0.095  1.41  -1.474    850.  0.268E+00  0.270E+00  0.218E-03  0.245E-13   19  34  35   15  32  13   12  27  17    0.000    0.000       0
+    188 00:47:03.33  11.8000A  -1.9563   1.8534   3.8072  0.095  1.43  -1.571    900.  0.270E+00  0.272E+00  0.214E-03  0.241E-13   19  34  35   15  32  13   13  27  17    0.000    0.000       0
+    189 00:47:14.93  11.6000A  -1.9058   1.8136   3.8844  0.095  1.43  -1.515    900.  0.272E+00  0.274E+00  0.210E-03  0.238E-13   19  34  35   15  33  12   13  27  17    0.000    0.000       0
+    190 00:47:26.43  11.5000A  -1.9185   1.7941   3.9300  0.095  1.43  -1.498    900.  0.274E+00  0.276E+00  0.209E-03  0.237E-13   19  33  33   15  33  12   14  27  17    0.000    0.000       0
+    191 00:47:37.73  11.3000A  -1.9295   1.7738   3.9966  0.095  1.43  -1.514    900.  0.275E+00  0.278E+00  0.205E-03  0.234E-13   19  33  34   17  31  16   14  27  17    0.000    0.000       0
+    192 00:47:48.93  11.2000A  -1.9254   1.7376   4.0156  0.095  1.43  -1.662    900.  0.277E+00  0.279E+00  0.203E-03  0.233E-13   19  33  33   17  32  15   14  27  17    0.000    0.000       0
+    193 00:48:00.03  11.1000A  -1.9236   1.8334   4.0558  0.096  1.43  -1.545    900.  0.278E+00  0.281E+00  0.202E-03  0.232E-13   19  33  33   18  31  16   15  27  17    0.000    0.000       0
+    194 00:48:11.13  11.1000A   1.9224   1.8792   4.0518  0.096  1.43  -1.519    900.  0.279E+00  0.282E+00  0.202E-03  0.233E-13   19  26  23   18  31  16   15  27  17    0.000    0.000       0
+    195 00:48:22.33  11.2000A   1.9515   1.8232   4.0092  0.096  1.43  -1.599    900.  0.279E+00  0.283E+00  0.204E-03  0.237E-13   19  26  23   10  18   8   16  27  17    0.000    0.000       0
+    196 00:48:33.63  11.3000A   1.9459   1.8319   3.9875  0.096  1.43  -1.508    900.  0.280E+00  0.283E+00  0.206E-03  0.240E-13   19  26  23   18  32  16   16  27  17    0.000    0.000       0
+    197 00:48:45.13  11.5000A   2.0032   1.8133   3.9024  0.097  1.43  -1.464    900.  0.280E+00  0.284E+00  0.210E-03  0.245E-13   19  25  23   18  32  16   16  27  17    0.000    0.000       0
+    198 00:48:57.03  11.9000A   2.0484  -1.8480   3.7806  0.097  1.43  -1.477    900.  0.281E+00  0.284E+00  0.218E-03  0.255E-13   19  25  23   11  39   1   16  26  17    0.000    0.000       0
+    199 00:49:09.03  12.0000A   2.0522  -1.8870   3.7610  0.097  1.43  -1.559    900.  0.281E+00  0.284E+00  0.221E-03  0.258E-13   19  25  23   11  39   1   16  26  17    0.000    0.000       0
+    200 00:49:21.13  12.1000A   2.0229  -1.8906   3.7099  0.098  1.43  -1.446    900.  0.281E+00  0.284E+00  0.224E-03  0.261E-13   19  25  23   11  39   1    9  14  22    0.000    0.000       0
+    201 00:49:33.13  12.0000A   2.0743  -1.8688   3.7569  0.098  1.43  -1.426    900.  0.281E+00  0.284E+00  0.223E-03  0.260E-13   19  24  23   11  39   1   10  14  22    0.000    0.000       0
+    202 00:49:45.13  12.0000A   2.0950  -1.8804   3.7464  0.098  1.43  -1.431    900.  0.281E+00  0.284E+00  0.225E-03  0.260E-13   19  24  23   12  38   1   10  14  22    0.000    0.000       0
+    203 00:49:57.23  12.1000A  -2.0884  -1.9138   3.7273  0.099  1.41  -1.409    850.  0.281E+00  0.284E+00  0.228E-03  0.263E-13   17  37  32   12  38   1   11  14  22    0.000    0.000       0
+    204 00:50:09.43  12.2000A  -2.0887  -1.9119   3.6887  0.099  1.41  -1.363    850.  0.281E+00  0.284E+00  0.231E-03  0.266E-13   17  37  32   12  38   1   11  14  22    0.000    0.000       0
+    205 00:50:21.93  12.5000A  -2.1896  -2.0015   3.6051  0.100  1.41  -1.350    850.  0.281E+00  0.284E+00  0.238E-03  0.274E-13   17  36  31   16  11  21   12  14  22    0.000    0.000       0
+    206 00:50:34.63  12.7000A  -2.2651  -2.0636   3.5441  0.100  1.41  -1.353    850.  0.281E+00  0.285E+00  0.244E-03  0.279E-13   17  36  31   16  11  21   12  14  22    0.000    0.000       0
+    207 00:50:47.63  13.0000A  -2.2971  -2.0814   3.4575  0.100  1.41  -1.345    850.  0.281E+00  0.285E+00  0.250E-03  0.287E-13   17  36  31   16  11  21   11  14  21    0.000    0.000       0
+    208 00:51:00.53  12.9000A  -2.2936  -2.0890   3.4958  0.101  1.41  -1.344    850.  0.281E+00  0.285E+00  0.249E-03  0.285E-13   17  36  31   17  10  21   13  13  21    0.000    0.000       0
+    209 00:51:13.53  13.0000A  -2.2664  -2.2144   3.4629  0.101  1.41  -1.330    850.  0.282E+00  0.285E+00  0.252E-03  0.287E-13   17  36  30   17  10  21   13  13  21    0.000    0.000       0
+    210 00:51:27.03  13.5000A  -2.2679  -2.2833   3.3395  0.101  1.38  -1.332    800.  0.282E+00  0.286E+00  0.263E-03  0.298E-13   16  36  29   17  10  21   14  13  21    0.000    0.000       0
+    211 00:51:40.13  13.1000A  -2.3404  -2.2995   3.4477  0.101  1.38  -1.317    800.  0.282E+00  0.286E+00  0.255E-03  0.289E-13   16  36  29   17  10  21   10  29  22    0.000    0.000       0
+    212 00:51:53.33  13.2000A  -2.3039  -2.2830   3.4199  0.101  1.38  -1.305    800.  0.283E+00  0.286E+00  0.257E-03  0.291E-13   16  36  29   17  10  20   10  29  22    0.000    0.000       0
+    213 00:52:06.53  13.2000A  -2.3426  -2.2478   3.4072  0.102  1.38  -1.295    800.  0.283E+00  0.287E+00  0.257E-03  0.291E-13   16  36  28   17  10  20   11  15  20    0.000    0.000       0
+    214 00:52:19.93  13.4000A  -2.3017  -2.2220   3.3705  0.102  1.38  -1.290    800.  0.283E+00  0.287E+00  0.260E-03  0.295E-13   16  36  28   18   9  21   12  15  20    0.000    0.000       0
+    215 00:52:33.13  13.2000A  -2.1993  -2.2564   3.4128  0.102  1.38  -1.286    800.  0.284E+00  0.287E+00  0.256E-03  0.291E-13   16  35  28   18   9  21   12  15  20    0.000    0.000       0
+    216 00:52:46.53  13.4000A  -2.1743  -2.2587   3.3666  0.101  1.38  -1.289    800.  0.284E+00  0.288E+00  0.260E-03  0.295E-13   16  35  28   18   9  21   13  15  20    0.000    0.000       0
+    217 00:52:59.83  13.3000A  -2.0894  -2.2654   3.3951  0.101  1.38  -1.299    800.  0.285E+00  0.288E+00  0.258E-03  0.292E-13   16  35  27   18   8  21   13  15  20    0.000    0.000       0
+    218 00:53:13.33  13.5000A  -2.0527  -2.2657   3.3361  0.101  1.38  -1.319    800.  0.285E+00  0.289E+00  0.261E-03  0.296E-13   15  37  27   18   8  21   14  15  20    0.000    0.000       0
+    219 00:53:26.73  13.4000A  -2.0443  -2.2331   3.3648  0.101  1.38  -1.363    800.  0.286E+00  0.289E+00  0.260E-03  0.293E-13   15  36  26   18   8  21   14  15  20    0.000    0.000       0
+    220 00:53:40.43  13.7000A  -2.0110   2.2497   3.2756  0.101  1.38  -1.343    800.  0.287E+00  0.290E+00  0.266E-03  0.300E-13   15  36  26   11   1  33   14  15  20    0.000    0.000       0
+    221 00:53:54.23  13.8000A  -1.9943   2.2539   3.2716  0.100  1.38  -1.343    800.  0.288E+00  0.291E+00  0.269E-03  0.302E-13   15  35  26   11   1  33   15  15  20    0.000    0.000       0
+    222 00:54:08.43  14.2000A  -1.9126   2.2149   3.1769  0.100  1.43  -1.352    900.  0.289E+00  0.292E+00  0.278E-03  0.310E-13   15  35  26   11   1  33   15  15  20    0.000    0.000       0
+    223 00:54:22.53  14.1000A   1.8463   2.1740   3.1980  0.099  1.43  -1.374    900.  0.290E+00  0.293E+00  0.277E-03  0.308E-13   18  37   3   12   2  33    8  24  19    0.000    0.000       0
+    224 00:54:36.33  13.8000A   1.8217   2.2035   3.2562  0.099  1.43  -1.375    900.  0.291E+00  0.293E+00  0.272E-03  0.302E-13   18  37   4   13   2  32    9  24  19    0.000    0.000       0
+    225 00:54:50.03  13.7000A  -1.8307   2.2091   3.2891  0.099  1.43  -1.374    900.  0.291E+00  0.294E+00  0.271E-03  0.300E-13   14  34  34   12   2  32    7  24  19    0.000    0.000       0
+    226 00:55:03.53  13.5000A  -1.9921   2.2362   3.3245  0.098  1.43  -1.449    900.  0.292E+00  0.294E+00  0.268E-03  0.296E-13   14  35  33   12   2  32    8  24  19    0.000    0.000       0
+    227 00:55:16.63  13.1000A  -2.0976   2.2104   3.4318  0.098  1.43  -1.353    900.  0.292E+00  0.295E+00  0.260E-03  0.289E-13   14  35  33   12   2  32    8  24  19    0.000    0.000       0
+    228 00:55:29.63  13.0000A  -2.1127   2.1442   3.4703  0.098  1.43  -1.357    900.  0.293E+00  0.295E+00  0.259E-03  0.287E-13   14  35  33   12   2  32    8  24  19    0.000    0.000       0
+    229 00:55:42.63  13.0000A  -2.1757   2.1861   3.4724  0.097  1.43  -1.379    900.  0.293E+00  0.296E+00  0.260E-03  0.288E-13   14  36  32   12   3  32    8  24  19    0.000    0.000       0
+    230 00:55:55.33  12.7000A  -2.2178   2.1991   3.5402  0.097  1.43  -1.353    900.  0.294E+00  0.296E+00  0.254E-03  0.282E-13   14  36  32   12   3  32    9  24  19    0.000    0.000       0
+    231 00:56:07.83  12.5000A  -2.2191   2.1754   3.6009  0.097  1.43  -1.340    900.  0.294E+00  0.297E+00  0.251E-03  0.279E-13   14  36  32   12   3  32    9  24  19    0.000    0.000       0
+    232 00:56:20.23  12.4000A  -2.1862   2.1233   3.6316  0.096  1.43  -1.348    900.  0.295E+00  0.297E+00  0.250E-03  0.277E-13   14  36  32   12   3  32    9  24  19    0.000    0.000       0
+    233 00:56:32.33  12.1000A  -2.1534   2.0497   3.7145  0.096  1.43  -1.375    900.  0.295E+00  0.297E+00  0.244E-03  0.271E-13   14  37  31   12   3  32   10  24  20    0.000    0.000       0
+    234 00:56:44.23  11.9000A   2.3527  -1.9985   3.7876  0.096  1.40  -1.434    850.  0.296E+00  0.298E+00  0.241E-03  0.266E-13   13  23  23   15  14  18   10  24  20    0.000    0.000       0
+    235 00:56:56.13  11.9000A   2.5119  -2.0811   3.7821  0.095  1.40  -2.509    850.  0.296E+00  0.298E+00  0.241E-03  0.266E-13   13  23  23   15  14  18   11  24  20    0.000    0.000       0
+    236 00:57:07.73  11.6000A   2.5215  -2.0743   3.8687  0.095  1.40  -1.380    850.  0.296E+00  0.298E+00  0.235E-03  0.260E-13   13  23  23   15  14  18   11  24  20    0.000    0.000       0
+    237 00:57:19.33  11.6000A   2.4780  -1.9862   3.8704  0.095  1.40  -1.361    850.  0.297E+00  0.299E+00  0.235E-03  0.259E-13   14  23  23   15  14  18   11  24  20    0.000    0.000       0
+    238 00:57:31.13  11.8000A   2.6413  -1.9429   3.8264  0.095  1.43  -1.348    900.  0.297E+00  0.299E+00  0.239E-03  0.263E-13   14  23  24   15  14  17   11  24  20    0.000    0.000       0
+    239 00:57:43.03  11.9000A   2.8071  -1.9566   3.7711  0.094  1.43  -1.350    900.  0.297E+00  0.299E+00  0.241E-03  0.265E-13   14  23  24   15  13  17   12  24  20    0.000    0.000       0
+    240 00:57:55.03  12.0000A   2.7559  -1.9826   3.7504  0.094  1.43  -1.358    900.  0.297E+00  0.299E+00  0.243E-03  0.267E-13   14  23  24   15  13  17   12  24  20    0.000    0.000       0
+    241 00:58:07.13  12.1000A   2.6541  -1.9789   3.7270  0.094  1.43  -1.369    900.  0.296E+00  0.298E+00  0.245E-03  0.269E-13   14  23  25   17  38   6   12  24  20    0.000    0.000       0
+    242 00:58:19.33  12.2000A   2.5772   1.9401   3.6792  0.093  1.43  -1.379    900.  0.296E+00  0.298E+00  0.247E-03  0.271E-13   15  23  24    9  16  24   12  24  20    0.000    0.000       0
+    243 00:58:31.53  12.2000A   2.6764   2.0023   3.6993  0.093  1.43  -1.388    900.  0.295E+00  0.298E+00  0.247E-03  0.272E-13   15  23  25    9  17  24   11  23  18    0.000    0.000       0
+    244 00:58:43.53  12.0000A   2.5236   2.0522   3.7468  0.092  1.43  -1.402    900.  0.295E+00  0.297E+00  0.243E-03  0.267E-13   15  23  25    9  17  24   11  23  18    0.000    0.000       0
+    245 00:58:55.53  12.0000A   2.4766   2.0650   3.7573  0.092  1.43  -1.424    900.  0.294E+00  0.296E+00  0.243E-03  0.268E-13   15  23  26    9  17  24   11  23  18    0.000    0.000       0
+    246 00:59:07.43  11.9000A  -2.2941   2.0526   3.7810  0.092  1.43  -1.427    900.  0.293E+00  0.296E+00  0.241E-03  0.266E-13   15  23  12    9  17  24   12  23  18    0.000    0.000       0
+    247 00:59:19.33  11.9000A  -2.2936   2.0268   3.7910  0.091  1.43  -1.450    900.  0.293E+00  0.295E+00  0.242E-03  0.266E-13   17  18  13    9  17  24   12  23  18    0.000    0.000       0
+    248 00:59:30.93  11.6000A  -2.3306  -2.0074   3.8738  0.091  1.43  -1.429    900.  0.292E+00  0.294E+00  0.236E-03  0.260E-13   17  18  13   16  36   7   11  22  18    0.000    0.000       0
+    249 00:59:42.23  11.3000A  -2.2921   2.0080   3.9762  0.091  1.43  -1.438    900.  0.291E+00  0.293E+00  0.230E-03  0.254E-13   17  18  13    9  18  24   12  22  18    0.000    0.000       0
+    250 00:59:53.33  11.1000A  -2.2240   2.0011   4.0447  0.091  1.43  -1.539    900.  0.290E+00  0.292E+00  0.227E-03  0.250E-13   16  23  11    9  18  24   13  22  18    0.000    0.000       0
+    251 01:00:04.33  11.0000A  -2.1937   1.9824   4.0855  0.090  1.43  -1.466    900.  0.289E+00  0.291E+00  0.225E-03  0.248E-13   15  21  11    9  18  24   13  22  18    0.000    0.000       0
+RUN  ITER. HH:MM:SS.SS    DT(E)     UMAX     VMAX     WMAX     U*    W*      THETA*     Z_I     ENERG.   DISTENERG    DIVOLD     DIVNEW     UMAX(KJI)    VMAX(KJI)    WMAX(KJI)   ADVECX   ADVECY   MGCYC
+---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+      0 00:00:00.00  20.0000A  -0.2192D -0.2169D -0.1214  0.000  1.38   0.000E+00    800.  0.106E-02  0.000E+00  0.000E+00  0.000E+00   11  32   5    5  37  29    4  26   0    0.000    0.000       0
+      1 00:00:20.00  20.0000A  -0.2192  -0.2169  -0.1214  0.001  0.69  -0.124E+03    100.  0.105E-02  0.948E-03  0.491E-05  0.355E-15   11  32   5    5  37  29    4  26   0    0.000    0.000       0
+      2 00:00:40.00  20.0000D  -0.2151  -0.2176  -0.1209  0.001  0.69  -0.125E+03    100.  0.104E-02  0.940E-03  0.802E-05  0.359E-15   11  32   5    5  37  29    4  26   0    0.000    0.000       0
+      3 00:01:00.00  20.0000D  -0.2107  -0.2178  -0.1204  0.001  0.69  -0.125E+03    100.  0.103E-02  0.932E-03  0.111E-04  0.363E-15   11  32   5    5  37  29    4  26   0    0.000    0.000       0
+      4 00:01:20.00  20.0000D  -0.2060  -0.2177  -0.1199  0.001  0.69  -0.126E+03    100.  0.102E-02  0.924E-03  0.142E-04  0.366E-15   11  32   5    4   5   2    4  26   0    0.000    0.000       0
+      5 00:01:40.00  20.0000D  -0.2011  -0.2184  -0.1193  0.001  0.93  -0.125E+03    250.  0.101E-02  0.916E-03  0.173E-04  0.368E-15   11  32   5    4   5   2    4  26   0    0.000    0.000       0
+      6 00:02:00.00  20.0000D  -0.1959  -0.2186  -0.1188  0.001  0.93  -0.124E+03    250.  0.100E-02  0.908E-03  0.204E-04  0.370E-15   11  32   5    4   5   2    4  26   0    0.000    0.000       0
+      7 00:02:20.00  20.0000D  -0.1942  -0.2184  -0.1182  0.001  1.22  -0.123E+03    550.  0.994E-03  0.900E-03  0.235E-04  0.371E-15   11  32   4    4   5   2    4  26   0    0.000    0.000       0
+      8 00:02:40.00  20.0000D  -0.1948D -0.2177D -0.1175  0.002  1.22  -0.126E+03    550.  0.204E-02  0.892E-03  0.516E-03  0.735E-13   11  32   4    4   5   2    4  26   0    0.000    0.000       0
+      9 00:03:00.00  20.0000D   0.3635  -0.2778  -0.1533  0.002  1.22  -0.969E+02    550.  0.201E-02  0.182E-02  0.321E-04  0.742E-15    4  35   4    7  14  12    5  10  39    0.000    0.000       0
+     10 00:03:20.00  20.0000D   0.3563  -0.2763  -0.1517  0.002  1.22  -0.925E+02    550.  0.198E-02  0.180E-02  0.351E-04  0.739E-15    4  35   4    7  14  12    5  10  39    0.000    0.000       0
+     11 00:03:40.00  20.0000D   0.3475  -0.2747  -0.1498  0.002  1.22  -0.887E+02    550.  0.195E-02  0.177E-02  0.380E-04  0.739E-15    4  35   4    6  29  30    5  10  39    0.000    0.000       0
+     12 00:04:00.00  20.0000D   0.3374  -0.2747  -0.1476  0.002  1.22  -0.856E+02    550.  0.193E-02  0.175E-02  0.410E-04  0.741E-15    4  35   4    6  29  30    5  10  39    0.000    0.000       0
+     13 00:04:20.00  20.0000D   0.3262  -0.2735  -0.1451  0.002  1.22  -0.926E+02    550.  0.190E-02  0.173E-02  0.439E-04  0.744E-15    4  35   4    6  29  30    5  10  39    0.000    0.000       0
+     14 00:04:40.00  20.0000D   0.3249  -0.2712  -0.1455  0.002  1.22  -0.766E+02    550.  0.188E-02  0.171E-02  0.468E-04  0.749E-15    4  35   5    6  29  30    6  23   4    0.000    0.000       0
+     15 00:05:00.00  20.0000D   0.3275D -0.2694D -0.1470  0.002  1.22  -0.703E+02    550.  0.293E-02  0.169E-02  0.529E-03  0.720E-13    4  35   5    7  13  38    6  23   4    0.000    0.000       0
+     16 00:05:20.00  20.0000D   0.3782  -0.3667  -0.2151  0.003  1.22  -0.600E+02    550.  0.288E-02  0.262E-02  0.546E-04  0.110E-14   11  16  39    7  13  12    8  21  26    0.000    0.000       0
+     17 00:05:40.00  20.0000D   0.3741  -0.3587  -0.2156  0.003  1.22  -0.562E+02    550.  0.284E-02  0.258E-02  0.573E-04  0.110E-14   11  16  39    7  15  26    8  21  26    0.000    0.000       0
+     18 00:06:00.00  20.0000D   0.3717  -0.3525  -0.2156  0.003  1.22  -0.520E+02    550.  0.279E-02  0.255E-02  0.600E-04  0.111E-14    6   3  26    7  15  26    8  21  26    0.000    0.000       0
+     19 00:06:20.00  20.0000D  -0.3700  -0.3448  -0.2153  0.003  1.22  -0.482E+02    550.  0.276E-02  0.251E-02  0.626E-04  0.112E-14    6  12  31    7  15  26    8  21  26    0.000    0.000       0
+     20 00:06:40.00  20.0000D  -0.3674  -0.3439  -0.2146  0.004  1.22  -0.433E+02    550.  0.273E-02  0.249E-02  0.652E-04  0.114E-14    6  12  31    8  34  19    8  21  26    0.000    0.000       0
+     21 00:07:00.00  20.0000D  -0.3632  -0.3493  -0.2134  0.004  1.22  -0.399E+02    550.  0.270E-02  0.246E-02  0.677E-04  0.116E-14    6  12  31    7  14  26    8  21  26    0.000    0.000       0
+     22 00:07:20.00  20.0000D  -0.3577  -0.3559   0.2292  0.005  1.22  -0.374E+02    550.  0.269E-02  0.245E-02  0.701E-04  0.120E-14    6  12  31    7  14  26    2  15  28    0.000    0.000       0
+     23 00:07:40.00  20.0000D  -0.3509D -0.3603D  0.2560  0.005  1.22  -0.331E+02    550.  0.377E-02  0.244E-02  0.552E-03  0.750E-13    6  12  31    7  14  26    2  15  28    0.000    0.000       0
+     24 00:08:00.00  20.0000D   0.4111   0.4185   0.3167  0.006  1.18  -0.251E+02    500.  0.374E-02  0.341E-02  0.770E-04  0.163E-14    8  36  11   10  34  34    2  15  28    0.000    0.000       0
+     25 00:08:20.00  20.0000D  -0.4110   0.4103   0.3500  0.007  1.18  -0.224E+02    500.  0.374E-02  0.340E-02  0.792E-04  0.168E-14    5  15  27   10  34  34    2  15  28    0.000    0.000       0
+     26 00:08:40.00  20.0000D  -0.4162   0.3992   0.3893  0.008  1.18  -0.200E+02    500.  0.375E-02  0.342E-02  0.812E-04  0.175E-14    5  15  27   10  34  34    2  15  28    0.000    0.000       0
+     27 00:09:00.00  20.0000D  -0.4207   0.3858   0.4347  0.009  1.18  -0.177E+02    500.  0.381E-02  0.347E-02  0.832E-04  0.184E-14    5  15  27   10  34  34    2  15  28    0.000    0.000       0
+     28 00:09:20.00  20.0000D  -0.4246  -0.3825   0.4864  0.010  1.18  -0.157E+02    500.  0.390E-02  0.355E-02  0.850E-04  0.195E-14    5  15  27    9  21  22    2  15  28    0.000    0.000       0
+     29 00:09:40.00  20.0000D  -0.4283  -0.3813   0.5443  0.011  1.18  -0.143E+02    500.  0.405E-02  0.368E-02  0.866E-04  0.209E-14    5  15  27    9  21  22    2  15  28    0.000    0.000       0
+     30 00:10:00.00  20.0000D  -0.4322D -0.3772D  0.6082  0.013  1.18  -0.125E+02    500.  0.530E-02  0.387E-02  0.555E-03  0.730E-13    5  15  27    9  21  22    2  15  28    0.000    0.000       0
+     31 00:10:20.00  20.0000D  -0.4451   0.4333   0.6835  0.014  1.18  -0.125E+02    500.  0.557E-02  0.505E-02  0.911E-04  0.276E-14    6  11  30    9  20  25    3  14  27    0.000    0.000       0
+     32 00:10:40.00  20.0000D  -0.4498   0.4329   0.7743  0.016  1.18  -0.995E+01    500.  0.598E-02  0.541E-02  0.920E-04  0.300E-14    6  11  30    9  20  25    3  14  27    0.000    0.000       0
+     33 00:11:00.00  20.0000D  -0.4499   0.4284   0.8794  0.018  1.18  -0.865E+01    500.  0.655E-02  0.591E-02  0.925E-04  0.329E-14    6  11  30    9  20  25    3  15  28    0.000    0.000       0
+     34 00:11:20.00  20.0000D  -0.4686   0.4203   1.0077  0.020  1.18  -0.770E+01    500.  0.733E-02  0.660E-02  0.926E-04  0.365E-14    1  15  30    9  20  25    3  15  28    0.000    0.000       0
+     35 00:11:40.00  20.0000A  -0.5131   0.4390   1.1411  0.023  1.18  -0.677E+01    500.  0.839E-02  0.754E-02  0.924E-04  0.407E-14    5  15  26    1  13  27    3  15  28    0.000    0.000       0
+     36 00:12:00.00  20.0000A  -0.5579   0.4874   1.2860  0.026  1.18  -0.602E+01    500.  0.980E-02  0.879E-02  0.919E-04  0.459E-14    1  14  30    1  17  36    3  20  36    0.000    0.000       0
+     37 00:12:20.00  20.0000A  -0.6069   0.5356   1.4354  0.028  1.18  -0.548E+01    500.  0.117E-01  0.104E-01  0.913E-04  0.519E-14    1  14  30    1  17  36    3  19  36    0.000    0.000       0
+     38 00:12:40.00  20.0000A  -0.6537   0.5816   1.6319  0.031  1.05  -0.512E+01    350.  0.141E-01  0.126E-01  0.907E-04  0.591E-14    1  14  30    1  17  36    4  20  36    0.000    0.000       0
+     39 00:13:00.00  20.0000A  -0.6962   0.6231   1.8563  0.034  1.18  -0.448E+01    500.  0.171E-01  0.153E-01  0.907E-04  0.675E-14    1  14  30    1  17  36    4  20  36    0.000    0.000       0
+     40 00:13:20.00  20.0000A  -0.7336   0.6662   2.0403  0.038  1.09  -0.414E+01    400.  0.210E-01  0.188E-01  0.916E-04  0.772E-14    1  16  29    1  29  28    4  19  36    0.000    0.000       0
+     41 00:13:40.00  20.0000A  -0.7826   0.7271   2.2419  0.041  1.09  -0.385E+01    400.  0.259E-01  0.232E-01  0.939E-04  0.884E-14    1  16  29    1  20  22    5  19  36    0.000    0.000       0
+     42 00:13:58.00  18.0000A  -0.8274   0.7915   2.4938  0.044  1.14  -0.346E+01    450.  0.312E-01  0.283E-01  0.874E-04  0.902E-14    1  15  29    1  20  22    5  19  36    0.000    0.000       0
+     43 00:14:14.70  16.7000A  -0.8825   0.8465   2.6922  0.046  1.14  -0.327E+01    450.  0.370E-01  0.339E-01  0.846E-04  0.939E-14    1   2  23    1  20  22    5  19  36    0.000    0.000       0
+     44 00:14:30.20  15.5000A  -0.9373   0.8934   2.8940  0.048  1.18  -0.311E+01    500.  0.432E-01  0.399E-01  0.822E-04  0.968E-14    1   2  23    1  20  22    6  19  36    0.000    0.000       0
+     45 00:14:44.70  14.5000A  -0.9840   0.9319   3.1005  0.050  1.18  -0.302E+01    500.  0.497E-01  0.462E-01  0.810E-04  0.997E-14    1   2  23    1  20  22    6  19  36    0.000    0.000       0
+     46 00:14:58.50  13.8000A  -1.0209   0.9629   3.2577  0.052  1.22  -0.317E+01    550.  0.566E-01  0.529E-01  0.813E-04  0.104E-13    1   2  23    1  20  22    6  19  36    0.000    0.000       0
+     47 00:15:11.50  13.0000A   1.0590   0.9874   3.4608  0.053  1.22  -0.294E+01    550.  0.637E-01  0.599E-01  0.809E-04  0.106E-13    1  31  18    1  20  22    7  19  36    0.000    0.000       0
+     48 00:15:23.90  12.4000A   1.1029   1.0076   3.6147  0.055  1.28  -0.287E+01    650.  0.711E-01  0.671E-01  0.815E-04  0.109E-13    1  31  18    1  29  28    7  19  36    0.000    0.000       0
+     49 00:15:35.90  12.0000A   1.1393   1.0350   3.7578  0.056  1.28  -0.277E+01    650.  0.787E-01  0.746E-01  0.832E-04  0.113E-13    1  31  18    1  20  23    8  19  36    0.000    0.000       0
+     50 00:15:47.40  11.5000A   1.1689   1.0620   3.9096  0.057  1.29  -0.265E+01    650.  0.865E-01  0.822E-01  0.842E-04  0.116E-13    1  31  18    1  20  23    8  19  36    0.000    0.000       0
+     51 00:15:58.70  11.3000A   1.1914   1.0868   3.9957  0.058  1.29  -0.255E+01    650.  0.945E-01  0.900E-01  0.870E-04  0.121E-13    1  31  18    1  20  23    9  19  36    0.000    0.000       0
+     52 00:16:09.60  10.9000A   1.2072   1.1100   4.1423  0.059  1.29  -0.247E+01    650.  0.103E+00  0.980E-01  0.882E-04  0.124E-13    1  31  18    1  20  23    9  19  36    0.000    0.000       0
+     53 00:16:20.30  10.7000A   1.2163   1.1311   4.2098  0.059  1.38  -0.241E+01    800.  0.111E+00  0.106E+00  0.907E-04  0.129E-13    1  31  18    1  20  23    9  19  36    0.000    0.000       0
+     54 00:16:30.70  10.4000A   1.2230   1.1501   4.3170  0.060  1.38  -0.239E+01    800.  0.119E+00  0.114E+00  0.923E-04  0.132E-13    1  30  18    1  20  23   10  19  36    0.000    0.000       0
+     55 00:16:40***  10.3000A   1.2285   1.1668   4.3818  0.061  1.38  -0.237E+01    800.  0.127E+00  0.122E+00  0.954E-04  0.137E-13    1  30  18    1  20  23   10  19  36    0.000    0.000       0
+     56 00:16:51.10  10.1000A   1.2378   1.2243   4.4452  0.061  1.38  -0.237E+01    800.  0.136E+00  0.131E+00  0.977E-04  0.140E-13    1  37  15    7  12   8   11  19  36    0.000    0.000       0
+     57 00:17:01.09   9.9900A   1.2531   1.2263   4.5048  0.062  1.38  -0.241E+01    800.  0.144E+00  0.139E+00  0.101E-03  0.145E-13    1  37  15    6  20  34   11  19  36    0.000    0.000       0
+     58 00:17:11.04   9.9500A  -1.3217   1.2613   4.5240  0.062  1.38  -0.242E+01    800.  0.153E+00  0.147E+00  0.104E-03  0.150E-13    4  35   3   10  17   2   12  19  36    0.000    0.000       0
+     59 00:17:20.87   9.8300A  -1.3776   1.3448   4.5770  0.063  1.38  -0.239E+01    800.  0.161E+00  0.156E+00  0.107E-03  0.155E-13    4  35   3   10  17   2   12  19  36    0.000    0.000       0
+     60 00:17:30.73   9.8600A  -1.4186   1.2882   4.5631  0.063  1.38  -0.240E+01    800.  0.170E+00  0.164E+00  0.111E-03  0.161E-13    4  35   3    9  16   3   13  19  36    0.000    0.000       0
+     61 00:17:40.50   9.7700A  -1.4466   1.4232   4.6056  0.064  1.41  -0.243E+01    850.  0.178E+00  0.173E+00  0.114E-03  0.166E-13    4  35   3   11  17   2   13  19  36    0.000    0.000       0
+     62 00:17:50.32   9.8200A  -1.4630   1.4757   4.5817  0.064  1.41  -0.251E+01    850.  0.186E+00  0.181E+00  0.119E-03  0.173E-13    4  35   3   11  17   2   13  19  36    0.000    0.000       0
+     63 00:18:00.12   9.8000A  -1.4689  -1.4758   4.5924  0.064  1.41  -0.244E+01    850.  0.195E+00  0.189E+00  0.122E-03  0.178E-13    4  35   3    9   4  39   14  19  36    0.000    0.000       0
+     64 00:18:09.98   9.8600A  -1.4900   1.5484   4.5623  0.065  1.41  -0.254E+01    850.  0.203E+00  0.198E+00  0.127E-03  0.185E-13    4  35   2   12  17   2   14  19  36    0.000    0.000       0
+     65 00:18:19.93   9.9500A  -1.5149   1.5649   4.5234  0.065  1.41  -0.240E+01    850.  0.211E+00  0.206E+00  0.132E-03  0.193E-13    5  35   2   12  17   2   15  19  36    0.000    0.000       0
+     66 00:18:29.90   9.9700A  -1.5562  -1.6149   4.5130  0.066  1.41  -0.230E+01    850.  0.219E+00  0.214E+00  0.136E-03  0.199E-13    5  35   2   10   4  39   15  15  27    0.000    0.000       0
+     67 00:18:39.84   9.9400A  -1.5803  -1.6213   4.5279  0.066  1.41  -0.227E+01    850.  0.227E+00  0.221E+00  0.140E-03  0.205E-13    5  35   2   10   4  39   12  15   1    0.000    0.000       0
+     68 00:18:49.80   9.9600A  -1.5944  -1.5975   4.5161  0.067  1.41  -0.222E+01    850.  0.234E+00  0.229E+00  0.144E-03  0.211E-13    5  35   2    9   5  38   12  15   1    0.000    0.000       0
+     69 00:18:59.73   9.9300A  -1.6046  -1.6886   4.5325  0.067  1.41  -0.221E+01    850.  0.241E+00  0.236E+00  0.148E-03  0.216E-13    5  35   2    9   5  38   13  15   1    0.000    0.000       0
+     70 00:19:09.63   9.9000A  -1.6147  -1.8090   4.5451  0.067  1.41  -0.225E+01    850.  0.247E+00  0.242E+00  0.151E-03  0.221E-13    5  35   2   10   4  38   13  15   1    0.000    0.000       0
+     71 00:19:19.63  10.0000A  -1.6261  -1.8785   4.4972  0.068  1.41  -0.232E+01    850.  0.253E+00  0.249E+00  0.156E-03  0.229E-13    5  35   2   10   4  38   13  15   1    0.000    0.000       0
+     72 00:19:29.63  10.0000A  -1.6387  -1.8839   4.4921  0.068  1.41  -0.241E+01    850.  0.259E+00  0.254E+00  0.160E-03  0.234E-13    5  35   2   10   4  38   14  15   1    0.000    0.000       0
+     73 00:19:39.73  10.1000A  -1.6503  -1.9429   4.4673  0.069  1.41  -0.218E+01    850.  0.264E+00  0.260E+00  0.166E-03  0.242E-13    5  35   2   11   4  38   14  15   1    0.000    0.000       0
+     74 00:19:49.93  10.2000A  -1.6583  -2.0326   4.3937  0.069  1.41  -0.215E+01    850.  0.269E+00  0.265E+00  0.171E-03  0.249E-13    5  35   2   11   4  38   14  15   1    0.000    0.000       0
+     75 00:20:00.23  10.3000A  -1.7171  -2.0769   4.3835  0.070  1.41  -0.208E+01    850.  0.273E+00  0.269E+00  0.177E-03  0.257E-13   12  36  39   12   3  38   12  35  25    0.000    0.000       0
+     76 00:20:10.43  10.2000A  -1.7567  -2.1408   4.4178  0.070  1.41  -0.210E+01    850.  0.277E+00  0.273E+00  0.179E-03  0.259E-13   12  36  39   12   3  38   13  35  25    0.000    0.000       0
+     77 00:20:20.53  10.1000A  -1.7527  -2.1271   4.4670  0.071  1.41  -0.208E+01    850.  0.280E+00  0.277E+00  0.180E-03  0.261E-13   12  36  39   12   3  38   13  35  25    0.000    0.000       0
+     78 00:20:30.73  10.2000A  -1.8161  -2.1682   4.4292  0.071  1.41  -0.221E+01    850.  0.283E+00  0.280E+00  0.185E-03  0.267E-13   13  36  38   13   3  38   14  35  25    0.000    0.000       0
+     79 00:20:40.73  10.0000A  -1.8889  -2.1900   4.4840  0.072  1.41  -0.198E+01    850.  0.285E+00  0.282E+00  0.185E-03  0.266E-13   13  36  38   13   3  38   14  35  25    0.000    0.000       0
+     80 00:20:50.83  10.1000A  -1.9297  -2.1973   4.4481  0.073  1.41  -0.194E+01    850.  0.286E+00  0.284E+00  0.190E-03  0.272E-13   13  36  38   14   2  38   14  35  25    0.000    0.000       0
+     81 00:21:01.13  10.3000A  -1.9368  -2.2598   4.3901  0.073  1.41  -0.192E+01    850.  0.288E+00  0.285E+00  0.196E-03  0.280E-13   13  36  38   14   2  38   15  35  25    0.000    0.000       0
+     82 00:21:11.43  10.3000A  -1.9076  -2.1711   4.3540  0.074  1.41  -0.191E+01    850.  0.288E+00  0.286E+00  0.198E-03  0.284E-13   13  36  38   14   2  38   15  35  25    0.000    0.000       0
+     83 00:21:22.03  10.6000A  -1.9439  -2.1446   4.2332  0.075  1.41  -0.190E+01    850.  0.288E+00  0.286E+00  0.206E-03  0.295E-13   14  36  38   15   2  38   15  35  25    0.000    0.000       0
+     84 00:21:33.13  11.1000A  -1.9810  -2.1166   4.0656  0.076  1.41  -0.190E+01    850.  0.288E+00  0.286E+00  0.218E-03  0.311E-13   14  36  38   15   2  38   16  35  25    0.000    0.000       0
+     85 00:21:44.53  11.4000A  -1.9978  -2.0805   3.9358  0.076  1.41  -0.188E+01    850.  0.288E+00  0.286E+00  0.226E-03  0.322E-13   14  36  38   15   1  38   16  35  25    0.000    0.000       0
+     86 00:21:56.53  12.0000A  -1.9911  -1.9954   3.7411  0.077  1.38  -0.188E+01    800.  0.287E+00  0.285E+00  0.240E-03  0.342E-13   14  36  38   15   2  38   16  35  25    0.000    0.000       0
+     87 00:22:09.33  12.8000A  -1.9587  -2.0083   3.5118  0.078  1.38  -0.184E+01    800.  0.286E+00  0.285E+00  0.258E-03  0.367E-13   14  36  38   17  38   9   12   6  38    0.000    0.000       0
+     88 00:22:22.23  12.9000A  -1.8968  -2.0019   3.4815  0.079  1.38  -0.182E+01    800.  0.285E+00  0.283E+00  0.263E-03  0.372E-13   14  36  38   17  38   9   13  20   6    0.000    0.000       0
+     89 00:22:35.43  13.2000A  -1.8835  -1.9266   3.4165  0.081  1.38  -0.182E+01    800.  0.283E+00  0.282E+00  0.271E-03  0.382E-13   15  35  38   17  38   9   14  36  16    0.000    0.000       0
+     90 00:22:48.73  13.3000A  -1.8964  -1.9624   3.3714  0.082  1.38  -0.177E+01    800.  0.281E+00  0.281E+00  0.275E-03  0.386E-13   16  35  37   10  14  36   14  36  16    0.000    0.000       0
+     91 00:23:02.53  13.8000A  -1.9119  -1.7922   3.2564  0.083  1.38  -0.176E+01    800.  0.280E+00  0.279E+00  0.286E-03  0.401E-13   16  35  36   10  14  36   14  36  14    0.000    0.000       0
+     92 00:23:16.73  14.2000A  -1.9743  -1.7683   3.1765  0.084  1.35  -0.171E+01    750.  0.277E+00  0.277E+00  0.293E-03  0.412E-13   16  35  36   11  14  36   14  19   6    0.000    0.000       0
+     93 00:23:31.43  14.7000A  -1.9862  -1.8340   3.0522  0.086  1.35  -0.166E+01    750.  0.275E+00  0.275E+00  0.303E-03  0.426E-13   16  35  36   11  14  36   14  19   6    0.000    0.000       0
+     94 00:23:46.43  15.0000A  -1.9412  -1.6691   2.9938  0.087  1.35  -0.163E+01    750.  0.272E+00  0.273E+00  0.308E-03  0.433E-13   16  35  36   11  14  36   16  18   6    0.000    0.000       0
+     95 00:24:01.33  14.9000A  -1.8521  -1.7158   3.0103  0.088  1.35  -0.159E+01    750.  0.269E+00  0.270E+00  0.303E-03  0.428E-13   16  35  35   11  14  32   11  22  10    0.000    0.000       0
+     96 00:24:16.13  14.8000A  -1.8021  -1.7470   3.0382  0.090  1.35  -0.156E+01    750.  0.266E+00  0.267E+00  0.299E-03  0.421E-13    7  20  10   11  14  32   11  22  10    0.000    0.000       0
+     97 00:24:31.03  14.9000A  -1.8736  -1.7532   3.0227  0.091  1.32  -0.153E+01    700.  0.263E+00  0.264E+00  0.297E-03  0.419E-13    8  20  10    1  22  29   12  22  10    0.000    0.000       0
+     98 00:24:45.93  14.9000A  -1.8942  -1.7566   3.0293  0.092  1.32  -0.149E+01    700.  0.259E+00  0.261E+00  0.295E-03  0.414E-13    8  20  10    1  22  29   12  22  10    0.000    0.000       0
+     99 00:25:00.93  15.0000A  -1.9088   1.7638   3.0033  0.094  1.32  -0.146E+01    700.  0.256E+00  0.258E+00  0.294E-03  0.411E-13    8  20   9    1   7  14   12  22  10    0.000    0.000       0
+    100 00:25:15.93  15.0000A  -1.9789   1.8127   3.0007  0.095  1.32  -0.147E+01    700.  0.252E+00  0.255E+00  0.292E-03  0.406E-13    8  20   9    1   7  14   13  22  10    0.000    0.000       0
+    101 00:25:31.03  15.1000A  -1.9401   1.8502   2.9834  0.096  1.32  -0.150E+01    700.  0.249E+00  0.252E+00  0.292E-03  0.404E-13    8  20   9    1   7  14   13  22  10    0.000    0.000       0
+    102 00:25:46.43  15.4000A  -1.8588   1.8909   2.9152  0.097  1.32  -0.145E+01    700.  0.245E+00  0.249E+00  0.295E-03  0.406E-13    9  20   9    1   8  13   13  22  10    0.000    0.000       0
+    103 00:26:01.93  15.5000A   1.8657   1.9566   2.8848  0.098  1.32  -0.144E+01    700.  0.242E+00  0.245E+00  0.295E-03  0.403E-13    8   0  24    1   8  13   14  22  10    0.000    0.000       0
+    104 00:26:17.53  15.6000A   1.8883   2.0109   2.8229  0.099  1.32  -0.141E+01    700.  0.238E+00  0.242E+00  0.294E-03  0.400E-13    9   0  24    1   8  14   14  22  10    0.000    0.000       0
+    105 00:26:33.33  15.8000A   1.8054   2.0547   2.7633  0.100  1.32  -0.137E+01    700.  0.234E+00  0.239E+00  0.296E-03  0.400E-13    9   0  24    1   8  14    8  10  12    0.000    0.000       0
+    106 00:26:49.33  16.0000A   1.7700   2.0820   2.7515  0.101  1.29  -0.137E+01    650.  0.231E+00  0.236E+00  0.297E-03  0.400E-13   10   0  24    1   8  14   12   8  17    0.000    0.000       0
+    107 00:27:04.93  15.6000A   1.7047   2.0951   2.7189  0.102  1.29  -0.157E+01    650.  0.227E+00  0.232E+00  0.287E-03  0.385E-13   10   0  24    1   8  14   12   8  17    0.000    0.000       0
+    108 00:27:19.33  14.4000A  -1.6865   2.1380   2.6514  0.103  1.29  -0.141E+01    650.  0.224E+00  0.229E+00  0.263E-03  0.352E-13    8  37  33    1   9  14   12   8  17    0.000    0.000       0
+    109 00:27:32.93  13.6000A   1.7101   2.1750   2.6464  0.103  1.29  -0.138E+01    650.  0.221E+00  0.226E+00  0.246E-03  0.328E-13   10  10   4    1   9  14    6  12   0    0.000    0.000       0
+    110 00:27:46.33  13.4000A   1.7037   2.1976   2.6700  0.103  1.32  -0.136E+01    700.  0.218E+00  0.224E+00  0.240E-03  0.319E-13   10  10   4    1   9  14    9  10  11    0.000    0.000       0
+    111 00:27:59.43  13.1000A  -1.7529   2.2135   2.6476  0.103  1.32  -0.139E+01    700.  0.216E+00  0.221E+00  0.233E-03  0.309E-13    5  36  26    1   9  14    9  10  11    0.000    0.000       0
+    112 00:28:12.53  13.1000A  -1.7808   2.2249   2.6316  0.103  1.35  -0.141E+01    750.  0.213E+00  0.219E+00  0.230E-03  0.305E-13    5  36  26    1   9  14    9   9  18    0.000    0.000       0
+    113 00:28:26.13  13.6000A  -1.7741   2.2303   2.6725  0.103  1.35  -0.138E+01    750.  0.211E+00  0.216E+00  0.237E-03  0.313E-13    7  35  25    1   9  14    9   9  18    0.000    0.000       0
+    114 00:28:40.63  14.5000A  -1.8334   2.2235   2.6508  0.103  1.35  -0.137E+01    750.  0.208E+00  0.214E+00  0.250E-03  0.329E-13    7  35  25    1   9  14    9   9  18    0.000    0.000       0
+    115 00:28:55.53  14.9000A  -1.8662   2.2176   2.6175  0.103  1.35  -0.138E+01    750.  0.206E+00  0.212E+00  0.254E-03  0.333E-13    7  35  25    1  10  14   10   9  18    0.000    0.000       0
+    116 00:29:10.03  14.5000A  -1.9194   2.1970   2.5512  0.102  1.35  -0.136E+01    750.  0.203E+00  0.209E+00  0.244E-03  0.320E-13    6  35  25    1  10  14   10  11  11    0.000    0.000       0
+    117 00:29:24.93  14.9000A  -1.9643   2.1498   2.4976  0.102  1.35  -0.137E+01    750.  0.201E+00  0.207E+00  0.248E-03  0.324E-13    6  35  25    1  10  14   13  35   7    0.000    0.000       0
+    118 00:29:40.63  15.7000A  -1.9574   2.1155   2.5043  0.101  1.35  -0.135E+01    750.  0.199E+00  0.205E+00  0.259E-03  0.336E-13    6  35  25    1  10  15   13  35   7    0.000    0.000       0
+    119 00:29:57.63  17.0000A  -1.8768   2.0720   2.4627  0.100  1.35  -0.137E+01    750.  0.197E+00  0.203E+00  0.276E-03  0.356E-13    6  35  25    1  10  16   14  35   7    0.000    0.000       0
+    120 00:30:15.93  18.3000A  -1.8639   2.0444   2.4602  0.099  1.35  -0.139E+01    750.  0.195E+00  0.200E+00  0.293E-03  0.375E-13    5  35  25    1  10  16    9  11   0    0.000    0.000       0
+    121 00:30:33.53  17.6000A  -1.8397   2.0035   2.5563  0.099  1.38  -0.140E+01    800.  0.193E+00  0.199E+00  0.278E-03  0.354E-13    5  35  25    1  10  16    9   0  35    0.000    0.000       0
+    122 00:30:50.83  17.3000A   1.8017   1.9547   2.6063  0.098  1.38  -0.137E+01    800.  0.192E+00  0.197E+00  0.270E-03  0.342E-13    1  22   8    1  10  17    9   0  35    0.000    0.000       0
+    123 00:31:07.63  16.8000A   1.8102   1.9365   2.6768  0.097  1.38  -0.135E+01    800.  0.191E+00  0.196E+00  0.260E-03  0.327E-13    1  22   8    1  10  17   10   0  35    0.000    0.000       0
+    124 00:31:24.03  16.4000A   1.7964   1.9013   2.7394  0.096  1.38  -0.137E+01    800.  0.190E+00  0.195E+00  0.251E-03  0.314E-13    1  22   8    1  10  17   10   0  35    0.000    0.000       0
+    125 00:31:40.33  16.3000A   1.7732   1.9071   2.7604  0.096  1.38  -0.140E+01    800.  0.190E+00  0.194E+00  0.248E-03  0.309E-13    1  22   9    1  10  11   10   0  35    0.000    0.000       0
+    126 00:31:56.63  16.3000A   1.7798   1.9172   2.7606  0.095  1.38  -0.141E+01    800.  0.189E+00  0.193E+00  0.247E-03  0.306E-13    1  22   9    1  10  11   11   0  35    0.000    0.000       0
+    127 00:32:12.93  16.3000A   1.7648   1.9124   2.7688  0.094  1.38  -0.145E+01    800.  0.189E+00  0.193E+00  0.247E-03  0.304E-13    1  22   9    1  11  10   11   0  35    0.000    0.000       0
+    128 00:32:29.23  16.3000A   1.7373   1.9408   2.7543  0.094  1.38  -0.147E+01    800.  0.189E+00  0.193E+00  0.247E-03  0.302E-13    1  22   9    1  11  11   11  39  35    0.000    0.000       0
+    129 00:32:45.33  16.1000A   1.7135   1.9425   2.7997  0.093  1.38  -0.148E+01    800.  0.189E+00  0.193E+00  0.244E-03  0.298E-13    1  21   9    1  11  11   12  39  35    0.000    0.000       0
+    130 00:33:01.23  15.9000A   1.6901   1.9076   2.8378  0.093  1.38  -0.147E+01    800.  0.190E+00  0.193E+00  0.242E-03  0.293E-13    1  21   9    1  12  11   12  39  35    0.000    0.000       0
+    131 00:33:17.03  15.8000A   1.6621   1.9401   2.8413  0.092  1.38  -0.148E+01    800.  0.190E+00  0.194E+00  0.240E-03  0.291E-13    1  22   9    1  12  11   12  39  35    0.000    0.000       0
+    132 00:33:33.03  16.0000A   1.6506   1.9292   2.8142  0.092  1.38  -0.149E+01    800.  0.191E+00  0.194E+00  0.244E-03  0.295E-13    1  22   9    1  12  11   12  39  35    0.000    0.000       0
+    133 00:33:49.13  16.1000A   1.6428   1.8765   2.7978  0.091  1.35  -0.150E+01    750.  0.191E+00  0.195E+00  0.246E-03  0.297E-13    1  22   9    1  12  11   12  39  38    0.000    0.000       0
+    134 00:34:05.13  16.0000A   1.6340   1.8514   2.8169  0.091  1.35  -0.154E+01    750.  0.192E+00  0.196E+00  0.245E-03  0.295E-13    1  22   9    1  13  11   13  39  38    0.000    0.000       0
+    135 00:34:20.53  15.4000A   1.6193   1.8427   2.9195  0.091  1.35  -0.151E+01    750.  0.193E+00  0.196E+00  0.237E-03  0.285E-13    1  22   9    1  13  11   10  25  21    0.000    0.000       0
+    136 00:34:35.73  15.2000A   1.6398   1.8299   2.9692  0.090  1.35  -0.151E+01    750.  0.193E+00  0.197E+00  0.234E-03  0.282E-13    7  28  15    1  13  10   11  25  21    0.000    0.000       0
+    137 00:34:50.83  15.1000A   1.5726   1.8722   2.9776  0.090  1.35  -0.152E+01    750.  0.194E+00  0.197E+00  0.233E-03  0.281E-13    8  28  15    5  28  13   11  25  21    0.000    0.000       0
+    138 00:35:06.23  15.4000A   1.6086   1.8828   2.9160  0.090  1.41  -0.152E+01    850.  0.195E+00  0.198E+00  0.239E-03  0.287E-13    7  29  15    5  28  13   12  25  22    0.000    0.000       0
+    139 00:35:21.33  15.1000A   1.6597   1.8593   2.9792  0.089  1.41  -0.151E+01    850.  0.195E+00  0.198E+00  0.235E-03  0.283E-13    7  29  15    6  29  14   13  24  21    0.000    0.000       0
+    140 00:35:36.13  14.8000A   1.6300   1.8733   3.0332  0.089  1.41  -0.149E+01    850.  0.196E+00  0.199E+00  0.231E-03  0.278E-13    7  29  15    6  29  14   13  24  21    0.000    0.000       0
+    141 00:35:50.63  14.5000A   1.5729   1.8689   3.1043  0.089  1.41  -0.148E+01    850.  0.197E+00  0.200E+00  0.228E-03  0.273E-13    7  30  16    6  29  13   14  24  21    0.000    0.000       0
+    142 00:36:05.23  14.6000A   1.6148   1.9109   3.0728  0.088  1.41  -0.161E+01    850.  0.197E+00  0.200E+00  0.230E-03  0.276E-13    7  30  16    6  29  13   14  24  22    0.000    0.000       0
+    143 00:36:20.23  15.0000A  -1.5863   1.9488   3.0092  0.088  1.41  -0.150E+01    850.  0.198E+00  0.201E+00  0.238E-03  0.285E-13    5  10  18    5  29  12   14  24  22    0.000    0.000       0
+    144 00:36:35.43  15.2000A  -1.6091   1.9169   2.9513  0.088  1.41  -0.154E+01    850.  0.199E+00  0.201E+00  0.242E-03  0.290E-13    5  10  18    5  29  12   15  24  22    0.000    0.000       0
+    145 00:36:50.53  15.1000A  -1.6253   1.8606   2.9742  0.088  1.41  -0.156E+01    850.  0.199E+00  0.202E+00  0.242E-03  0.289E-13    5  10  18    6  29  13   14  24  23    0.000    0.000       0
+    146 00:37:05.33  14.8000A   1.6733   1.8247   3.0435  0.088  1.41  -0.157E+01    850.  0.200E+00  0.202E+00  0.238E-03  0.283E-13   10   2   1    6  30  13   15  24  23    0.000    0.000       0
+    147 00:37:20.43  15.1000A   1.6830   1.8252   2.9896  0.088  1.41  -0.160E+01    850.  0.201E+00  0.203E+00  0.243E-03  0.289E-13   10   2   1    6  29  12   15  24  23    0.000    0.000       0
+    148 00:37:35.73  15.3000A   1.7228   1.8007   2.9340  0.088  1.41  -0.152E+01    850.  0.201E+00  0.203E+00  0.247E-03  0.293E-13   10   2   2    6  30  12   12  26  16    0.000    0.000       0
+    149 00:37:51.03  15.3000A   1.7699   1.7905   2.9399  0.088  1.41  -0.150E+01    850.  0.202E+00  0.204E+00  0.249E-03  0.294E-13   10   2   2    6  30  12   12  26  16    0.000    0.000       0
+    150 00:38:06.43  15.4000A   1.7813   1.6850   2.9239  0.088  1.41  -0.150E+01    850.  0.203E+00  0.205E+00  0.251E-03  0.296E-13   10   2   2   10  31  16   12  26  16    0.000    0.000       0
+    151 00:38:21.93  15.5000A   1.7704   1.7084   2.9002  0.088  1.41  -0.151E+01    850.  0.203E+00  0.205E+00  0.254E-03  0.298E-13   10   2   2   15  31  16   13  26  16    0.000    0.000       0
+    152 00:38:37.63  15.7000A  -1.7634   1.7338   2.8703  0.088  1.38  -0.160E+01    800.  0.204E+00  0.206E+00  0.259E-03  0.303E-13    6  11  17   15  31  16   13  26  16    0.000    0.000       0
+    153 00:38:53.63  16.0000A  -1.7422   1.7364   2.8097  0.089  1.38  -0.156E+01    800.  0.205E+00  0.207E+00  0.265E-03  0.309E-13    6  11  17   15  31  16   13  26  16    0.000    0.000       0
+    154 00:39:09.83  16.2000A  -1.7319   1.7003   2.7847  0.089  1.38  -0.157E+01    800.  0.206E+00  0.207E+00  0.269E-03  0.313E-13    6  10  16   15  31  16   12  19  11    0.000    0.000       0
+    155 00:39:25.83  16.0000A  -1.6881   1.7217   2.8044  0.089  1.38  -0.151E+01    800.  0.207E+00  0.208E+00  0.266E-03  0.309E-13    6  10  16   16  32  16    6  38  37    0.000    0.000       0
+    156 00:39:41.43  15.6000A   1.6873   1.7292   2.8834  0.090  1.38  -0.151E+01    800.  0.208E+00  0.209E+00  0.261E-03  0.302E-13    7  13  11   16  32  16    7  38  37    0.000    0.000       0
+    157 00:39:56.83  15.4000A   1.6931   1.7043   2.9213  0.090  1.38  -0.151E+01    800.  0.209E+00  0.210E+00  0.258E-03  0.299E-13    7  13  11   15  32  17    7  38  37    0.000    0.000       0
+    158 00:40:12.03  15.2000A  -1.7631   1.6749   2.9664  0.090  1.38  -0.152E+01    800.  0.210E+00  0.211E+00  0.257E-03  0.297E-13    7  10  18   15  32  17    8  38  36    0.000    0.000       0
+    159 00:40:26.73  14.7000A  -1.8237   1.7154   3.0609  0.091  1.38  -0.151E+01    800.  0.211E+00  0.212E+00  0.249E-03  0.289E-13    7  10  18   15  33  17    9  38  36    0.000    0.000       0
+    160 00:40:41.03  14.3000A  -1.8252   1.7108   3.1425  0.091  1.38  -0.152E+01    800.  0.212E+00  0.213E+00  0.244E-03  0.283E-13    7  10  18   15  33  17    9  38  36    0.000    0.000       0
+    161 00:40:54.83  13.8000A  -1.8469   1.6486   3.2591  0.091  1.38  -0.148E+01    800.  0.213E+00  0.214E+00  0.237E-03  0.275E-13    7  10  17   15  33  17   10  38  36    0.000    0.000       0
+    162 00:41:08.43  13.6000A  -1.8831   1.6711   3.3158  0.092  1.38  -0.146E+01    800.  0.214E+00  0.216E+00  0.234E-03  0.272E-13    7  10  17   16  34  17   10  38  36    0.000    0.000       0
+    163 00:41:21.93  13.5000A  -1.8756   1.6558   3.3246  0.092  1.38  -0.144E+01    800.  0.216E+00  0.217E+00  0.233E-03  0.271E-13    9   9  17   16  34  17   11  38  36    0.000    0.000       0
+    164 00:41:35.23  13.3000A  -1.9099   1.5999   3.3883  0.092  1.38  -0.145E+01    800.  0.217E+00  0.218E+00  0.230E-03  0.267E-13    9   9  17   16  34  17   11  38  36    0.000    0.000       0
+    165 00:41:48.63  13.4000A  -1.9272   1.6343   3.3563  0.093  1.38  -0.148E+01    800.  0.219E+00  0.220E+00  0.232E-03  0.269E-13    9   9  17   11  17  11   11  37  36    0.000    0.000       0
+    166 00:42:01.73  13.1000A  -1.9304   1.6460   3.4318  0.093  1.38  -0.153E+01    800.  0.220E+00  0.221E+00  0.227E-03  0.263E-13    9   9  17   11  17  11   12  37  36    0.000    0.000       0
+    167 00:42:14.83  13.1000A  -1.9369   1.6727   3.4408  0.093  1.38  -0.153E+01    800.  0.222E+00  0.223E+00  0.227E-03  0.262E-13    9   9  16    6  13  23   12  37  36    0.000    0.000       0
+    168 00:42:27.93  13.1000A  -1.9583   1.7244   3.4359  0.094  1.38  -0.153E+01    800.  0.224E+00  0.225E+00  0.227E-03  0.262E-13    9   9  16    6  13  23   13  38  35    0.000    0.000       0
+    169 00:42:40.83  12.9000A  -1.9542   1.6652   3.4799  0.094  1.38  -0.148E+01    800.  0.226E+00  0.226E+00  0.224E-03  0.257E-13    9   9  16    6  13  23   13  37  35    0.000    0.000       0
+    170 00:42:53.43  12.6000A  -1.9277   1.6484   3.5802  0.094  1.41  -0.149E+01    850.  0.227E+00  0.228E+00  0.218E-03  0.251E-13    9   9  16    7  13  23   13  37  35    0.000    0.000       0
+    171 00:43:05.53  12.1000A  -1.9744   1.7970   3.7067  0.094  1.41  -0.147E+01    850.  0.229E+00  0.230E+00  0.210E-03  0.241E-13   16  38  33    7  13  23   14  37  35    0.000    0.000       0
+    172 00:43:17.53  12.0000A  -1.9712   1.7847   3.7576  0.094  1.41  -0.144E+01    850.  0.231E+00  0.232E+00  0.208E-03  0.239E-13   16  38  32    7  13  23   14  37  35    0.000    0.000       0
+    173 00:43:29.53  12.0000A  -1.9522   1.6696   3.7489  0.095  1.41  -0.143E+01    850.  0.233E+00  0.234E+00  0.207E-03  0.238E-13   16  38  32    6  13  24   13  35  37    0.000    0.000       0
+    174 00:43:41.43  11.9000A  -1.9097   1.6401   3.7698  0.095  1.43  -0.144E+01    900.  0.235E+00  0.236E+00  0.206E-03  0.236E-13   16  37  34    9  37  34   13  35  37    0.000    0.000       0
+    175 00:43:53.43  12.0000A  -1.9635   1.7114   3.7561  0.095  1.43  -0.150E+01    900.  0.237E+00  0.238E+00  0.207E-03  0.238E-13   17  37  33    7  13  24   13  35  37    0.000    0.000       0
+    176 00:44:05.33  11.9000A  -1.9575   1.7504   3.7690  0.095  1.43  -0.148E+01    900.  0.239E+00  0.240E+00  0.206E-03  0.236E-13   17  36  34   13  30  13   14  35  37    0.000    0.000       0
+    177 00:44:17.33  12.0000A  -1.9934  -1.7824   3.7510  0.095  1.43  -0.149E+01    900.  0.241E+00  0.242E+00  0.207E-03  0.237E-13   17  36  34    7   3   0   14  35  37    0.000    0.000       0
+    178 00:44:29.53  12.2000A  -1.9613  -1.7779   3.7032  0.095  1.43  -0.151E+01    900.  0.243E+00  0.244E+00  0.211E-03  0.241E-13   17  35  34    7   3   0   14  35  37    0.000    0.000       0
+    179 00:44:41.93  12.4000A  -2.0902  -1.7806   3.6302  0.095  1.43  -0.150E+01    900.  0.245E+00  0.246E+00  0.215E-03  0.245E-13   18  36  33    7   2   0   14  35  37    0.000    0.000       0
+    180 00:44:54.63  12.7000A  -2.1697   1.7640   3.5391  0.095  1.43  -0.150E+01    900.  0.247E+00  0.248E+00  0.220E-03  0.251E-13   18  36  33   14  30  13   14  35  37    0.000    0.000       0
+    181 00:45:07.73  13.1000A  -2.1866   1.8097   3.4378  0.095  1.43  -0.150E+01    900.  0.250E+00  0.250E+00  0.228E-03  0.258E-13   18  36  32   14  30  13   14  35  37    0.000    0.000       0
+    182 00:45:20.93  13.2000A  -2.1004   1.8603   3.4167  0.095  1.43  -0.148E+01    900.  0.252E+00  0.253E+00  0.230E-03  0.261E-13   18  36  32   14  31  13   12  27  19    0.000    0.000       0
+    183 00:45:33.93  13.0000A  -2.1511   1.8803   3.4715  0.095  1.43  -0.148E+01    900.  0.254E+00  0.255E+00  0.228E-03  0.257E-13   18  35  33   14  31  13   12  27  19    0.000    0.000       0
+    184 00:45:46.73  12.8000A  -2.0333   1.8484   3.5034  0.094  1.41  -0.148E+01    850.  0.257E+00  0.258E+00  0.226E-03  0.254E-13   18  35  33   14  31  13   13  27  19    0.000    0.000       0
+    185 00:45:59.43  12.7000A  -1.9856   1.7671   3.5422  0.094  1.41  -0.147E+01    850.  0.259E+00  0.260E+00  0.225E-03  0.253E-13   19  36  32   14  31  13   13  27  19    0.000    0.000       0
+    186 00:46:12.13  12.7000A  -1.9793   1.8141   3.5472  0.095  1.41  -0.148E+01    850.  0.261E+00  0.262E+00  0.226E-03  0.254E-13   19  36  31   15  31  13   13  27  19    0.000    0.000       0
+    187 00:46:24.83  12.7000A  -2.0450   1.8636   3.5335  0.095  1.41  -0.147E+01    850.  0.264E+00  0.265E+00  0.227E-03  0.255E-13   19  36  31   15  32  12   14  27  19    0.000    0.000       0
+    188 00:46:37.23  12.4000A  -1.9751   1.8908   3.6233  0.095  1.41  -0.146E+01    850.  0.266E+00  0.267E+00  0.223E-03  0.250E-13   19  34  34   15  32  12   12  27  17    0.000    0.000       0
+    189 00:46:49.33  12.1000A  -1.9859   1.8800   3.7117  0.095  1.41  -0.147E+01    850.  0.268E+00  0.270E+00  0.218E-03  0.245E-13   19  34  35   15  32  13   12  27  17    0.000    0.000       0
+    190 00:47:01.23  11.9000A  -1.9649   1.8669   3.7906  0.095  1.43  -0.153E+01    900.  0.270E+00  0.272E+00  0.215E-03  0.243E-13   19  34  35   15  32  13   13  27  17    0.000    0.000       0
+    191 00:47:12.83  11.6000A  -1.9177   1.8155   3.8747  0.095  1.43  -0.152E+01    900.  0.272E+00  0.274E+00  0.210E-03  0.238E-13   19  34  35   15  33  12   13  27  17    0.000    0.000       0
+    192 00:47:24.33  11.5000A  -1.9172   1.8023   3.9215  0.095  1.43  -0.149E+01    900.  0.274E+00  0.276E+00  0.209E-03  0.237E-13   19  33  33   15  33  12   13  27  17    0.000    0.000       0
+    193 00:47:35.63  11.3000A  -1.9288   1.7679   3.9898  0.095  1.43  -0.153E+01    900.  0.275E+00  0.278E+00  0.205E-03  0.234E-13   19  33  34   17  31  16   14  27  17    0.000    0.000       0
+    194 00:47:46.83  11.2000A  -1.9283   1.7499   4.0174  0.095  1.43  -0.158E+01    900.  0.277E+00  0.279E+00  0.203E-03  0.233E-13   19  33  34   17  31  16   14  27  17    0.000    0.000       0
+    195 00:47:57.93  11.1000A  -1.9260   1.8121   4.0520  0.096  1.43  -0.156E+01    900.  0.278E+00  0.280E+00  0.202E-03  0.232E-13   19  33  33   18  31  16   15  27  17    0.000    0.000       0
+    196 00:48:09.03  11.1000A   1.9138   1.8812   4.0589  0.096  1.43  -0.151E+01    900.  0.279E+00  0.282E+00  0.202E-03  0.233E-13   19  26  23   18  31  16   15  27  17    0.000    0.000       0
+    197 00:48:20.23  11.2000A   1.9491   1.8188   4.0098  0.096  1.43  -0.165E+01    900.  0.279E+00  0.283E+00  0.204E-03  0.236E-13   19  26  23   10  18   8   15  27  17    0.000    0.000       0
+    198 00:48:31.53  11.3000A   1.9486   1.8243   3.9990  0.096  1.43  -0.154E+01    900.  0.280E+00  0.283E+00  0.206E-03  0.240E-13   19  26  23   18  32  16   16  27  17    0.000    0.000       0
+    199 00:48:43.03  11.5000A   1.9928   1.8240   3.9242  0.097  1.43  -0.147E+01    900.  0.280E+00  0.284E+00  0.210E-03  0.245E-13   19  25  23   18  32  16   16  27  17    0.000    0.000       0
+    200 00:48:54.93  11.9000A   2.0447  -1.8354   3.7936  0.097  1.43  -0.147E+01    900.  0.281E+00  0.284E+00  0.218E-03  0.255E-13   19  25  23   11  39   1   16  27  17    0.000    0.000       0
+    201 00:49:06.83  11.9000A   2.0545  -1.8822   3.7669  0.097  1.43  -0.154E+01    900.  0.281E+00  0.284E+00  0.219E-03  0.256E-13   19  25  23   11  39   1   16  26  17    0.000    0.000       0
+    202 00:49:18.93  12.1000A   2.0294  -1.8918   3.7184  0.098  1.43  -0.146E+01    900.  0.281E+00  0.284E+00  0.224E-03  0.261E-13   19  25  23   11  39   1   16  26  17    0.000    0.000       0
+    203 00:49:30.93  12.0000A   2.0628  -1.8742   3.7491  0.098  1.43  -0.143E+01    900.  0.281E+00  0.284E+00  0.223E-03  0.260E-13   19  24  23   11  39   1   10  14  22    0.000    0.000       0
+    204 00:49:42.93  12.0000A   2.0971  -1.8691   3.7547  0.098  1.43  -0.143E+01    900.  0.281E+00  0.284E+00  0.225E-03  0.260E-13   19  24  23   12  38   1   10  14  22    0.000    0.000       0
+    205 00:49:55.03  12.1000A  -2.0862  -1.9102   3.7240  0.099  1.43  -0.143E+01    900.  0.281E+00  0.284E+00  0.228E-03  0.263E-13   17  37  32   12  38   1   11  14  22    0.000    0.000       0
+    206 00:50:07.23  12.2000A  -2.0905  -1.9145   3.7023  0.099  1.41  -0.137E+01    850.  0.281E+00  0.284E+00  0.231E-03  0.266E-13   17  37  32   12  38   1   11  14  22    0.000    0.000       0
+    207 00:50:19.73  12.5000A  -2.1722  -1.9829   3.6048  0.099  1.41  -0.135E+01    850.  0.281E+00  0.284E+00  0.238E-03  0.274E-13   17  36  31   16  11  21   12  14  22    0.000    0.000       0
+    208 00:50:32.33  12.6000A  -2.2555  -2.0562   3.5616  0.100  1.41  -0.135E+01    850.  0.281E+00  0.285E+00  0.241E-03  0.277E-13   17  36  31   16  11  21   12  14  22    0.000    0.000       0
+    209 00:50:45.33  13.0000A  -2.2942  -2.0814   3.4709  0.100  1.41  -0.135E+01    850.  0.281E+00  0.285E+00  0.250E-03  0.286E-13   17  36  31   16  11  21   11  14  21    0.000    0.000       0
+    210 00:50:58.23  12.9000A  -2.2962  -2.0729   3.4884  0.101  1.41  -0.134E+01    850.  0.281E+00  0.285E+00  0.249E-03  0.285E-13   17  36  31   16  11  21   13  13  21    0.000    0.000       0
+    211 00:51:11.13  12.9000A  -2.2708  -2.1941   3.4777  0.101  1.41  -0.133E+01    850.  0.282E+00  0.285E+00  0.250E-03  0.285E-13   17  36  31   17  10  21   13  13  21    0.000    0.000       0
+    212 00:51:24.53  13.4000A  -2.2516  -2.2734   3.3579  0.101  1.38  -0.133E+01    800.  0.282E+00  0.286E+00  0.261E-03  0.296E-13   17  36  30   17  10  21   13  13  21    0.000    0.000       0
+    213 00:51:37.63  13.1000A  -2.3359  -2.3003   3.4415  0.101  1.38  -0.132E+01    800.  0.282E+00  0.286E+00  0.255E-03  0.289E-13   16  36  29   17  10  21   10  29  22    0.000    0.000       0
+    214 00:51:50.73  13.1000A  -2.3172  -2.2864   3.4327  0.101  1.38  -0.131E+01    800.  0.283E+00  0.286E+00  0.255E-03  0.289E-13   16  36  29   17  10  20   10  29  22    0.000    0.000       0
+    215 00:52:03.93  13.2000A  -2.3376  -2.2619   3.4039  0.102  1.38  -0.130E+01    800.  0.283E+00  0.287E+00  0.257E-03  0.291E-13   16  36  28   17  10  20   11  15  20    0.000    0.000       0
+    216 00:52:17.33  13.4000A  -2.3173  -2.2088   3.3700  0.102  1.38  -0.129E+01    800.  0.283E+00  0.287E+00  0.261E-03  0.295E-13   16  36  28   18   9  21   11  15  20    0.000    0.000       0
+    217 00:52:30.53  13.2000A  -2.2072  -2.2495   3.4138  0.102  1.38  -0.129E+01    800.  0.284E+00  0.287E+00  0.256E-03  0.291E-13   16  36  28   18   9  21   12  15  20    0.000    0.000       0
+    218 00:52:43.93  13.4000A  -2.1856  -2.2585   3.3599  0.102  1.38  -0.129E+01    800.  0.284E+00  0.288E+00  0.260E-03  0.295E-13   16  35  28   18   9  21   12  15  20    0.000    0.000       0
+    219 00:52:57.13  13.2000A  -2.1054  -2.2597   3.3981  0.101  1.38  -0.130E+01    800.  0.285E+00  0.288E+00  0.256E-03  0.290E-13   16  35  28   18   8  21   13  15  20    0.000    0.000       0
+    220 00:53:10.63  13.5000A  -2.0659  -2.2661   3.3364  0.101  1.38  -0.131E+01    800.  0.285E+00  0.289E+00  0.261E-03  0.296E-13   16  35  27   18   8  21   13  15  20    0.000    0.000       0
+    221 00:53:24.03  13.4000A  -2.0469  -2.2412   3.3682  0.101  1.38  -0.136E+01    800.  0.286E+00  0.289E+00  0.260E-03  0.293E-13   15  36  26   18   8  21   14  15  20    0.000    0.000       0
+    222 00:53:37.63  13.6000A  -2.0226   2.2442   3.3003  0.101  1.38  -0.134E+01    800.  0.287E+00  0.290E+00  0.264E-03  0.298E-13   15  36  26   11   1  33   14  15  20    0.000    0.000       0
+    223 00:53:51.33  13.7000A  -2.0032   2.2579   3.2736  0.100  1.38  -0.134E+01    800.  0.288E+00  0.291E+00  0.267E-03  0.300E-13   15  35  26   11   1  33   15  15  20    0.000    0.000       0
+    224 00:54:05.33  14.0000A  -1.9354   2.2267   3.2037  0.100  1.43  -0.135E+01    900.  0.289E+00  0.292E+00  0.274E-03  0.306E-13   15  35  26   11   1  33   15  15  20    0.000    0.000       0
+    225 00:54:19.33  14.0000A   1.8526   2.1765   3.2042  0.100  1.43  -0.137E+01    900.  0.290E+00  0.292E+00  0.275E-03  0.306E-13   18  37   3   12   2  33    8  24  19    0.000    0.000       0
+    226 00:54:33.33  14.0000A   1.8286   2.1939   3.2118  0.099  1.43  -0.137E+01    900.  0.291E+00  0.293E+00  0.276E-03  0.306E-13   18  37   4   13   2  32    9  24  19    0.000    0.000       0
+    227 00:54:47.03  13.7000A  -1.8187   2.2026   3.2762  0.099  1.43  -0.137E+01    900.  0.291E+00  0.294E+00  0.271E-03  0.300E-13   14  34  34   13   2  32    7  24  19    0.000    0.000       0
+    228 00:55:00.63  13.6000A  -1.9592   2.2377   3.3122  0.099  1.43  -0.141E+01    900.  0.292E+00  0.294E+00  0.269E-03  0.299E-13   14  35  33   12   2  32    7  24  19    0.000    0.000       0
+    229 00:55:13.83  13.2000A  -2.0859   2.2217   3.4151  0.098  1.43  -0.136E+01    900.  0.292E+00  0.295E+00  0.262E-03  0.291E-13   14  35  33   12   2  32    8  24  19    0.000    0.000       0
+    230 00:55:26.83  13.0000A  -2.1155   2.1620   3.4681  0.098  1.43  -0.135E+01    900.  0.293E+00  0.295E+00  0.259E-03  0.287E-13   14  35  33   12   2  32    8  24  19    0.000    0.000       0
+    231 00:55:39.83  13.0000A  -2.1634   2.1815   3.4743  0.097  1.43  -0.138E+01    900.  0.293E+00  0.296E+00  0.259E-03  0.288E-13   14  36  32   12   3  32    8  24  19    0.000    0.000       0
+    232 00:55:52.63  12.8000A  -2.2142   2.2026   3.5212  0.097  1.43  -0.136E+01    900.  0.294E+00  0.296E+00  0.256E-03  0.284E-13   14  36  32   12   3  32    9  24  19    0.000    0.000       0
+    233 00:56:05.13  12.5000A  -2.2226   2.1848   3.5908  0.097  1.43  -0.134E+01    900.  0.294E+00  0.296E+00  0.251E-03  0.278E-13   14  36  32   12   3  32    9  24  19    0.000    0.000       0
+    234 00:56:17.53  12.4000A  -2.1962   2.1370   3.6283  0.096  1.43  -0.135E+01    900.  0.295E+00  0.297E+00  0.250E-03  0.277E-13   14  36  32   12   3  32    9  24  19    0.000    0.000       0
+    235 00:56:29.73  12.2000A  -2.1538   2.0662   3.6908  0.096  1.43  -0.137E+01    900.  0.295E+00  0.297E+00  0.246E-03  0.273E-13   14  37  32   12   3  32   10  24  20    0.000    0.000       0
+    236 00:56:41.63  11.9000A   2.2926   1.9795   3.7819  0.096  1.40  -0.141E+01    850.  0.296E+00  0.298E+00  0.241E-03  0.266E-13   13  23  23   12   4  32   10  24  20    0.000    0.000       0
+    237 00:56:53.53  11.9000A   2.4879  -2.0658   3.7914  0.096  1.40  -0.156E+01    850.  0.296E+00  0.298E+00  0.241E-03  0.266E-13   13  23  23   15  14  18   10  24  20    0.000    0.000       0
+    238 00:57:05.23  11.7000A   2.5308  -2.0787   3.8604  0.095  1.40  -0.139E+01    850.  0.296E+00  0.298E+00  0.237E-03  0.262E-13   13  23  23   15  14  18   11  24  20    0.000    0.000       0
+    239 00:57:16.83  11.6000A   2.4543  -2.0061   3.8780  0.095  1.40  -0.136E+01    850.  0.297E+00  0.299E+00  0.235E-03  0.260E-13   13  23  23   15  14  18   11  24  20    0.000    0.000       0
+    240 00:57:28.53  11.7000A   2.5765  -1.9509   3.8404  0.095  1.43  -0.135E+01    900.  0.297E+00  0.299E+00  0.237E-03  0.261E-13   14  23  24   15  14  17   11  24  20    0.000    0.000       0
+    241 00:57:40.43  11.9000A   2.7843  -1.9425   3.7828  0.094  1.43  -0.135E+01    900.  0.297E+00  0.299E+00  0.241E-03  0.265E-13   14  23  24   15  13  17   11  24  20    0.000    0.000       0
+    242 00:57:52.43  12.0000A   2.7850  -1.9770   3.7573  0.094  1.43  -0.136E+01    900.  0.297E+00  0.299E+00  0.243E-03  0.267E-13   14  23  24   15  13  17   12  24  20    0.000    0.000       0
+    243 00:58:04.43  12.0000A   2.6600  -1.9775   3.7358  0.094  1.43  -0.137E+01    900.  0.296E+00  0.298E+00  0.243E-03  0.267E-13   14  23  25   17  38   6   12  24  20    0.000    0.000       0
+    244 00:58:16.63  12.2000A   2.5855  -1.9419   3.6974  0.093  1.43  -0.138E+01    900.  0.296E+00  0.298E+00  0.247E-03  0.271E-13   15  23  24   17  38   6   12  24  20    0.000    0.000       0
+    245 00:58:28.83  12.2000A   2.6669   1.9855   3.6789  0.093  1.43  -0.138E+01    900.  0.295E+00  0.298E+00  0.247E-03  0.272E-13   15  23  25    9  17  24   11  23  18    0.000    0.000       0
+    246 00:58:40.83  12.0000A   2.5808   2.0433   3.7357  0.092  1.43  -0.140E+01    900.  0.295E+00  0.297E+00  0.243E-03  0.267E-13   15  23  25    9  17  24   11  23  18    0.000    0.000       0
+    247 00:58:52.83  12.0000A   2.4882   2.0625   3.7551  0.092  1.43  -0.142E+01    900.  0.294E+00  0.296E+00  0.243E-03  0.268E-13   15  23  26    9  17  24   11  23  18    0.000    0.000       0
+    248 00:59:04.73  11.9000A   2.3305   2.0541   3.7681  0.092  1.43  -0.142E+01    900.  0.293E+00  0.296E+00  0.241E-03  0.266E-13   15  23  26    9  17  24   12  23  18    0.000    0.000       0
+    249 00:59:16.63  11.9000A  -2.2770   2.0301   3.7872  0.091  1.43  -0.145E+01    900.  0.293E+00  0.295E+00  0.242E-03  0.266E-13   15  23  12    9  17  24   12  23  18    0.000    0.000       0
+    250 00:59:28.33  11.7000A  -2.3207  -2.0141   3.8619  0.091  1.43  -0.143E+01    900.  0.292E+00  0.294E+00  0.238E-03  0.262E-13   17  18  13   16  36   7   11  22  18    0.000    0.000       0
+    251 00:59:39.73  11.4000A  -2.3017   2.0035   3.9583  0.091  1.43  -0.144E+01    900.  0.291E+00  0.293E+00  0.232E-03  0.256E-13   17  18  13    9  18  24   12  22  18    0.000    0.000       0
+    252 00:59:50.93  11.2000A  -2.2327   2.0003   4.0168  0.091  1.43  -0.152E+01    900.  0.290E+00  0.292E+00  0.229E-03  0.252E-13   16  23  11    9  18  24   13  22  18    0.000    0.000       0
+    253 01:00:01.93  11.0000A  -2.2037   1.9839   4.0832  0.090  1.43  -0.147E+01    900.  0.289E+00  0.291E+00  0.225E-03  0.248E-13   16  23  11    9  18  24   13  22  18    0.000    0.000       0

palm/trunk/SOURCE/advec_ws.f90

-                      r1222
+                      r1257
 ! Current revisions:
 ! ------------------
+!
+! accelerator loop directives removed
+!
 ! Former revisions:
 …
 !--    Computation of fluxes and tendency terms
        !$acc kernels present( ddzw, sk, tend, u, v, w, wall_flags_0, wall_flags_00 )
-       !$acc loop
        DO  i = i_left, i_right
           DO  j = j_south, j_north
-             !$acc loop vector( 32 )
              DO  k = nzb+1, nzt
 …
 !--    Computation of fluxes and tendency terms
        !$acc  kernels present( ddzw, tend, u, v, w, wall_flags_0, wall_flags_00 )
-       !$acc  loop
        DO i = i_left, i_right
           DO  j = j_south, j_north
-             !$acc  loop vector( 32 )
              DO  k = nzb+1, nzt
 …
 !--    Computation of fluxes and tendency terms
        !$acc kernels present( ddzw, tend, u, v, w, wall_flags_0, wall_flags_00 )
-       !$acc loop
        DO  i = i_left, i_right
           DO  j = j_south, j_north
-             !$acc loop vector( 32 )
              DO  k = nzb+1, nzt
 …
 !--    Computation of fluxes and tendency terms
        !$acc kernels present( ddzu, tend, u, v, w, wall_flags_0, wall_flags_00 )
-       !$acc loop
        DO i = i_left, i_right
           DO  j = j_south, j_north
-             !$acc loop vector( 32 )
              DO  k = nzb+1, nzt

palm/trunk/SOURCE/boundary_conds.f90

-                      r1242
+                      r1257
 ! Current revisions:
 ! -----------------
+!
+! loop independent clauses added
+!
 ! Former revisions:
 …
     IF ( ibc_pt_b == 0 )  THEN
        !$acc kernels present( nzb_s_inner, pt, pt_p )
+       !$acc loop independent
        DO  i = nxlg, nxrg
+          !$acc loop independent
           DO  j = nysg, nyng
              pt_p(nzb_s_inner(j,i),j,i) = pt(nzb_s_inner(j,i),j,i)
 …
     ELSEIF ( ibc_pt_b == 1 )  THEN
        !$acc kernels present( nzb_s_inner, pt_p )
+       !$acc loop independent
        DO  i = nxlg, nxrg
+          !$acc loop independent
           DO  j = nysg, nyng
              pt_p(nzb_s_inner(j,i),j,i) = pt_p(nzb_s_inner(j,i)+1,j,i)
 …
     IF ( .NOT. constant_diffusion )  THEN
        !$acc kernels present( e_p, nzb_s_inner )
+       !$acc loop independent
        DO  i = nxlg, nxrg
+          !$acc loop independent
           DO  j = nysg, nyng
              e_p(nzb_s_inner(j,i),j,i) = e_p(nzb_s_inner(j,i)+1,j,i)

palm/trunk/SOURCE/buoyancy.f90

-                      r1242
+                      r1257
 ! Currrent revisions:
 ! ------------------
+!
+! vector length (32) removed from openacc clause
+!
 ! Former revisions:
 …
           DO  i = i_left, i_right
              DO  j = j_south, j_north
                 !$acc loop independent vector(32)
+                !$acc loop independent vector
                 DO  k = nzb_s_inner(j,i)+1, nzt-1
                    tend(k,j,i) = tend(k,j,i) + atmos_ocean_sign * g * 0.5 * (  &

palm/trunk/SOURCE/coriolis.f90

-                      r1132
+                      r1257
 ! Current revisions:
 ! -----------------
+!
+! openacc loop and loop vector clauses removed
+!
 ! Former revisions:
 …
           CASE ( 1 )
              !$acc  kernels present( nzb_u_inner, tend, v, vg, w )
-             !$acc  loop
              DO  i = i_left, i_right
                 DO  j = j_south, j_north
-                   !$acc loop vector( 32 )
                    DO  k = 1, nzt
                       IF  ( k > nzb_u_inner(j,i) )  THEN
 …
           CASE ( 2 )
              !$acc  kernels present( nzb_v_inner, tend, u, ug )
-             !$acc  loop
              DO  i = i_left, i_right
                 DO  j = j_south, j_north
-                   !$acc loop vector( 32 )
                    DO  k = 1, nzt
                       IF  ( k > nzb_v_inner(j,i) )  THEN
 …
           CASE ( 3 )
              !$acc  kernels present( nzb_w_inner, tend, u )
-             !$acc  loop
              DO  i = i_left, i_right
                 DO  j = j_south, j_north
-                   !$acc loop vector( 32 )
                    DO  k = 1, nzt
                       IF  ( k > nzb_w_inner(j,i) )  THEN

palm/trunk/SOURCE/diffusion_e.f90

-                      r1182
+                      r1257
 ! Current revisions:
 ! -----------------
+!
+! openacc loop and loop vector clauses removed
+!
 ! Former revisions:
 …
           !$acc kernels present( ddzu, ddzw, dd2zu, diss, e, km, l_grid ) &
           !$acc         present( nzb_s_inner, rif, tend, var, zu, zw )
-          !$acc loop
           DO  i = i_left, i_right
              DO  j = j_south, j_north
-                !$acc loop vector( 32 )
                 DO  k = 1, nzt
 …
           !$acc kernels present( ddzu, ddzw, dd2zu, diss, e, km, l_grid ) &
           !$acc         present( nzb_s_inner, rif, tend, var, zu, zw )
-          !$acc loop
           DO  i = i_left, i_right
              DO  j = j_south, j_north
-                !$acc loop vector( 32 )
                 DO  k = 1, nzt
 …
        IF ( use_sgs_for_particles  .OR.  wang_kernel  .OR.  turbulence )  THEN
           !$acc kernels present( diss, nzb_s_inner )
-          !$acc loop
           DO  i = i_left, i_right
-             !$acc loop vector( 32 )
              DO  j = j_south, j_north
                 diss(nzb_s_inner(j,i),j,i) = diss(nzb_s_inner(j,i)+1,j,i)

palm/trunk/SOURCE/diffusion_s.f90

-                      r1132
+                      r1257
 ! Current revisions:
 ! ------------------
+!
+! openacc loop and loop vector clauses removed
+!
 ! Former revisions:
 …
        !$acc         present( s_flux_b, s_flux_t, tend, wall_s_flux )         &
        !$acc         present( wall_w_x, wall_w_y )
-       !$acc loop
        DO  i = i_left, i_right
           DO  j = j_south, j_north
+!
 !--          Compute horizontal diffusion
-             !$acc loop vector( 32 )
              DO  k = 1, nzt
                 IF ( k > nzb_s_outer(j,i) )  THEN
 …
+!
 !--          Apply prescribed horizontal wall heatflux where necessary
-             !$acc loop vector(32)
              DO  k = 1, nzt
                 IF ( k > nzb_s_inner(j,i)  .AND.  k <= nzb_s_outer(j,i)  .AND. &
 …
 !--          prescribed or computed at bottom and/or top, index k starts/ends at
 !--          nzb+2 or nzt-1, respectively.
-             !$acc loop vector( 32 )
              DO  k = 1, nzt_diff
                 IF ( k >= nzb_diff_s_inner(j,i) )  THEN
 …
 !--          Vertical diffusion at the first computational gridpoint along
 !--          z-direction
-             !$acc loop vector( 32 )
              DO  k = 1, nzt
                 IF ( use_surface_fluxes  .AND.  k == nzb_s_inner(j,i)+1 )  THEN

palm/trunk/SOURCE/diffusion_u.f90

-                      r1132
+                      r1257
 ! Current revisions:
 ! -----------------
+!
+! openacc loop and loop vector clauses removed, declare create moved after
+! the FORTRAN declaration statement
+!
 ! Former revisions:
 …
        REAL    ::  kmym, kmyp, kmzm, kmzp
+       REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) ::  usvs
        !$acc declare create ( usvs )
-       REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) ::  usvs
+!
 …
        !$acc         present ( ddzu, ddzw, fym, fyp, wall_u )           &
        !$acc         present ( nzb_u_inner, nzb_u_outer, nzb_diff_u )
-       !$acc loop
        DO  i = i_left, i_right
           DO  j = j_south, j_north
+!
 !--          Compute horizontal diffusion
-             !$acc loop vector(32)
              DO  k = 1, nzt
                 IF ( k > nzb_u_outer(j,i) )  THEN
 …
+!
 !--          Wall functions at the north and south walls, respectively
-             !$acc loop vector(32)
              DO  k = 1, nzt
                 IF( k > nzb_u_inner(j,i)  .AND.  k <= nzb_u_outer(j,i)  .AND. &
 …
 !--          Compute vertical diffusion. In case of simulating a Prandtl layer,
 !--          index k starts at nzb_u_inner+2.
-             !$acc loop vector(32)
              DO  k = 1, nzt_diff
                 IF ( k >= nzb_diff_u(j,i) )  THEN
 …
        IF ( use_surface_fluxes )  THEN
-          !$acc loop
           DO  i = i_left, i_right
-             !$acc loop vector(32)
              DO  j = j_south, j_north
 …
           k = nzt
-          !$acc loop
           DO  i = i_left, i_right
-             !$acc loop vector(32)
              DO  j = j_south, j_north

palm/trunk/SOURCE/diffusion_v.f90

-                      r1132
+                      r1257
 ! Current revisions:
 ! -----------------
+!
+! openacc loop and loop vector clauses removed, declare create moved after
+! the FORTRAN declaration statement
+!
 ! Former revisions:
 …
        REAL    ::  kmxm, kmxp, kmzm, kmzp
+       REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) ::  vsus
        !$acc declare create ( vsus )
-       REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) ::  vsus
+!
 …
        !$acc         present ( ddzu, ddzw, fxm, fxp, wall_v )           &
        !$acc         present ( nzb_v_inner, nzb_v_outer, nzb_diff_v )
-       !$acc loop
        DO  i = i_left, i_right
           DO  j = j_south, j_north
+!
 !--          Compute horizontal diffusion
-             !$acc loop vector(32)
              DO  k = 1, nzt
                 IF ( k > nzb_v_outer(j,i) )  THEN
 …
+!
 !--          Wall functions at the left and right walls, respectively
-             !$acc loop vector(32)
              DO  k = 1, nzt
                 IF( k > nzb_v_inner(j,i)  .AND.  k <= nzb_v_outer(j,i)  .AND. &
 …
 !--          Compute vertical diffusion. In case of simulating a Prandtl
 !--          layer, index k starts at nzb_v_inner+2.
-             !$acc loop vector(32)
              DO  k = 1, nzt_diff
                 IF ( k >= nzb_diff_v(j,i) )  THEN
 …
        IF ( use_surface_fluxes )  THEN
-          !$acc loop
           DO  i = i_left, i_right
-             !$acc loop vector(32)
              DO  j = j_south, j_north
 …
           k = nzt
-          !$acc loop
           DO  i = i_left, i_right
-             !$acc loop vector(32)
              DO  j = j_south, j_north

palm/trunk/SOURCE/diffusion_w.f90

-                      r1132
+                      r1257
 ! Current revisions:
 ! -----------------
+!
+! openacc loop and loop vector clauses removed, declare create moved after
+! the FORTRAN declaration statement
+!
 ! Former revisions:
 …
        REAL    ::  kmxm, kmxp, kmym, kmyp
+       REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) ::  wsus, wsvs
        !$acc declare create ( wsus, wsvs )
-       REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) ::  wsus, wsvs
+!
 …
        !$acc         present ( ddzu, ddzw, fwxm, fwxp, fwym, fwyp, wall_w_x, wall_w_y )           &
        !$acc         present ( nzb_w_inner, nzb_w_outer )
-       !$acc loop
        DO  i = i_left, i_right
           DO  j = j_south, j_north
-             !$acc loop vector( 32 )
              DO  k = 1, nzt
                 IF ( k > nzb_w_outer(j,i) )  THEN
 …
+!
 !--          Wall functions at all vertical walls, where necessary
-             !$acc loop vector( 32 )
              DO  k = 1,nzt

palm/trunk/SOURCE/exchange_horiz.f90

-                      r1132
+                      r1257
 ! Current revisions:
 ! -----------------
+!
+! openacc loop and loop vector clauses removed
+!
 ! Former revisions:
 …
           DO  i = 0, nbgp_local-1
              DO  j = nys-nbgp_local, nyn+nbgp_local
-                !$acc loop vector( 32 )
                 DO  k = nzb, nzt+1
                    ar(k,j,nxl-nbgp_local+i) = ar(k,j,nxr-nbgp_local+1+i)
 …
        IF ( on_device )  THEN
           !$acc kernels present( ar )
-          !$acc loop
           DO  i = nxl-nbgp_local, nxr+nbgp_local
              !$acc loop independent
              DO  j = 0, nbgp_local-1
                 !$acc loop vector( 32 )
+                !$acc loop independent
                 DO  k = nzb, nzt+1
                    ar(k,nys-nbgp_local+j,i) = ar(k,nyn-nbgp_local+1+j,i)

palm/trunk/SOURCE/fft_xy.f90

-                      r1220
+                      r1257
 ! Current revisions:
 ! -----------------
+!
+! openacc loop and loop vector clauses removed, declare create moved after
+! the FORTRAN declaration statement
+!
 ! Former revisions:
 …
        REAL, DIMENSION(6*(nx+1)) ::  work2
 #elif defined( __cuda_fft )
+       COMPLEX(dpk), DIMENSION(0:(nx+1)/2,nys_x:nyn_x,nzb_x:nzt_x) ::  ar_tmp
        !$acc declare create( ar_tmp )
-       COMPLEX(dpk), DIMENSION(0:(nx+1)/2,nys_x:nyn_x,nzb_x:nzt_x) ::  ar_tmp
 #endif
        REAL, DIMENSION(0:nx,nys_x:nyn_x), OPTIONAL   ::  ar_2d
 …
              !$acc kernels
-             !$acc loop
              DO  k = nzb_x, nzt_x
                 DO  j = nys_x, nyn_x
-                   !$acc loop vector( 32 )
                    DO  i = 0, (nx+1)/2
                       ar(i,j,k)      = REAL( ar_tmp(i,j,k) )  * dnx
                    ENDDO
-                   !$acc loop vector( 32 )
                    DO  i = 1, (nx+1)/2 - 1
                       ar(nx+1-i,j,k) = AIMAG( ar_tmp(i,j,k) ) * dnx
 …
              !$acc data present( ar )
              !$acc kernels
-             !$acc loop
              DO  k = nzb_x, nzt_x
                 DO  j = nys_x, nyn_x
 …
                    ar_tmp(0,j,k) = CMPLX( ar(0,j,k), 0.0 )
-                   !$acc loop vector( 32 )
                    DO  i = 1, (nx+1)/2 - 1
                       ar_tmp(i,j,k) = CMPLX( ar(i,j,k), ar(nx+1-i,j,k) )
 …
        REAL, DIMENSION(6*(ny+1)) ::  work2
 #elif defined( __cuda_fft )
+       COMPLEX(dpk), DIMENSION(0:(ny+1)/2,nxl_y:nxr_y,nzb_y:nzt_y) ::  ar_tmp
        !$acc declare create( ar_tmp )
-       COMPLEX(dpk), DIMENSION(0:(ny+1)/2,nxl_y:nxr_y,nzb_y:nzt_y) ::  ar_tmp
 #endif
        REAL, DIMENSION(0:ny,nxl_y_l:nxr_y_l,nzb_y:nzt_y) ::  ar
 …
              !$acc kernels
-             !$acc loop
              DO  k = nzb_y, nzt_y
                 DO  i = nxl_y, nxr_y
-                   !$acc loop vector( 32 )
                    DO  j = 0, (ny+1)/2
                       ar(j,i,k)      = REAL( ar_tmp(j,i,k) )  * dny
                    ENDDO
-                   !$acc loop vector( 32 )
                    DO  j = 1, (ny+1)/2 - 1
                       ar(ny+1-j,i,k) = AIMAG( ar_tmp(j,i,k) ) * dny
 …
              !$acc data present( ar )
              !$acc kernels
-             !$acc loop
              DO  k = nzb_y, nzt_y
                 DO  i = nxl_y, nxr_y
 …
                    ar_tmp(0,i,k) = CMPLX( ar(0,i,k), 0.0 )
-                   !$acc loop vector( 32 )
                    DO  j = 1, (ny+1)/2 - 1
                       ar_tmp(j,i,k) = CMPLX( ar(j,i,k), ar(ny+1-j,i,k) )

palm/trunk/SOURCE/flow_statistics.f90

-                      r1242
+                      r1257
 ! Current revisions:
 ! -----------------
+!
+! openacc "end parallel" replaced by "end parallel loop"
+!
 ! Former revisions:
 …
           sums_l(k,4,tn) = s3
        ENDDO
        !$acc end parallel
+       !$acc end parallel loop
+!
 …
              sums_l(k,23,tn) = s1
           ENDDO
           !$acc end parallel
+          !$acc end parallel loop
        ENDIF
 …
              sums_l(k,44,tn) = s2
           ENDDO
           !$acc end parallel
+          !$acc end parallel loop
           IF ( cloud_physics )  THEN
 …
                 sums_l(k,43,tn) = s2
              ENDDO
              !$acc end parallel
+             !$acc end parallel loop
           ENDIF
        ENDIF
 …
              sums_l(k,41,tn) = s1
           ENDDO
           !$acc end parallel
+          !$acc end parallel loop
        ENDIF
        !$OMP END PARALLEL
 …
           sums_l(k,38,tn) = s7
        ENDDO
        !$acc end parallel
+       !$acc end parallel loop
        IF ( humidity )  THEN
 …
              sums_l(k,70,tn) = s1
           ENDDO
           !$acc end parallel
+          !$acc end parallel loop
        ENDIF
 …
           ENDDO
        ENDDO
        !$acc end parallel
+       !$acc end parallel loop
        !$acc parallel present( sums_l )
        sums_l(nzb+4,pr_palm,tn) = s1
 …
              sums_l(k,34,tn) = s4
           ENDDO
           !$acc end parallel
+          !$acc end parallel loop
+!
 !--       Total perturbation TKE
 …
           sums_l(k,16,tn) = s3
        ENDDO
        !$acc end parallel
+       !$acc end parallel loop
+!
 …
              sums_l(k,65,tn) = s1
           ENDDO
           !$acc end parallel
+          !$acc end parallel loop
        ENDIF
 …
              sums_l(k,48,tn) = s2
           ENDDO
           !$acc end parallel
+          !$acc end parallel loop
           IF ( cloud_physics ) THEN
 …
                 sums_l(k,51,tn) = s1
              ENDDO
              !$acc end parallel
+             !$acc end parallel loop
           ENDIF
 …
              sums_l(k,48,tn) = s1
           ENDDO
           !$acc end parallel
+          !$acc end parallel loop
        ENDIF
 …
           sums_l(k,37,tn) = s3
        ENDDO
        !$acc end parallel
+       !$acc end parallel loop
+!
 …
                 sums_l(k,66,tn) = s1
              ENDDO
              !$acc end parallel
+             !$acc end parallel loop
           ENDIF
 …
              sums_l(k,71,tn) = s2
           ENDDO
           !$acc end parallel
+          !$acc end parallel loop
        ENDIF
 …
                 sums_l(k,46,tn) = s1
              ENDDO
              !$acc end parallel
+             !$acc end parallel loop
              IF ( .NOT. cloud_droplets )  THEN
 …
                    sums_l(k,52,tn) = s1
                 ENDDO
                 !$acc end parallel
+                !$acc end parallel loop
                 IF ( icloud_scheme == 0  )  THEN
 …
                       sums_l(k,75,tn) = s2
                    ENDDO
                    !$acc end parallel
+                   !$acc end parallel loop
                    IF ( precipitation )  THEN
 …
                          sums_l(k,76,tn) = s3
                       ENDDO
                       !$acc end parallel
+                      !$acc end parallel loop
                    ENDIF
 …
                       sums_l(k,54,tn) = s1
                    ENDDO
                    !$acc end parallel
+                   !$acc end parallel loop
                 ENDIF
 …
                    sums_l(k,54,tn) = s1
                 ENDDO
                 !$acc end parallel
+                !$acc end parallel loop
              ENDIF
 …
                    sums_l(k,46,tn) = s1
                 ENDDO
                 !$acc end parallel
+                !$acc end parallel loop
              ELSEIF ( ws_scheme_sca  .AND.  sr == 0 )  THEN
 …
 .61 * hom(k,1,4,sr) * sums_l(k,49,tn)
                 ENDDO
                 !$acc end parallel
+                !$acc end parallel loop
              ENDIF
 …
              sums_l(k,49,tn) = s1
           ENDDO
           !$acc end parallel
+          !$acc end parallel loop
        ENDIF
 …
              sums_l(k,15,tn) = s1
           ENDDO
           !$acc end parallel
+          !$acc end parallel loop
        ENDIF
 …
              sums_l(k,17,tn) = s1
           ENDDO
           !$acc end parallel
+          !$acc end parallel loop
           IF ( humidity )  THEN
 …
                 sums_l(k,49,tn) = s1
              ENDDO
              !$acc end parallel
+             !$acc end parallel loop
           ENDIF

palm/trunk/SOURCE/modules.f90

-                      r1242
+                      r1257
 ! Current revisions:
 ! ------------------
+!
+! set default values for grid indices of maximum velocity components u|v|w_max_ijk
+!
 ! Former revisions:
 …
     CHARACTER (LEN=40) ::  region(0:9)
     INTEGER ::  pr_palm = 90, statistic_regions = 0
     INTEGER ::  u_max_ijk(3), v_max_ijk(3), w_max_ijk(3)
+    INTEGER ::  u_max_ijk(3) = -1, v_max_ijk(3) = -1, w_max_ijk(3) = -1
     LOGICAL ::  flow_statistics_called = .FALSE.
     REAL ::     u_max, v_max, w_max

palm/trunk/SOURCE/prandtl_fluxes.f90

-                      r1037
+                      r1257
 ! Current revisions:
 ! -----------------
+!
+! openACC "kernels do" replaced by "kernels loop", "loop independent" added
+!
 ! Former revisions:
 …
 !--    for u* use the value from the previous time step
        !$OMP PARALLEL DO
        !$acc kernels do
+       !$acc kernels loop
        DO  i = nxlg, nxrg
           DO  j = nysg, nyng
 …
 !--    (the Richardson number is still the one from the previous time step)
        !$OMP PARALLEL DO PRIVATE( a, b, k, z_p )
        !$acc kernels do
+       !$acc kernels loop
        DO  i = nxlg, nxrg
           DO  j = nysg, nyng
 …
     IF ( .NOT. humidity )  THEN
        !$OMP PARALLEL DO PRIVATE( k, z_p )
        !$acc kernels do
+       !$acc kernels loop
        DO  i = nxlg, nxrg
           DO  j = nysg, nyng
 …
     ELSE
        !$OMP PARALLEL DO PRIVATE( k, z_p )
        !$acc kernels do
+       !$acc kernels loop
        DO  i = nxlg, nxrg
           DO  j = nysg, nyng
 …
 !-- Compute u* at the scalars' grid points
     !$OMP PARALLEL DO PRIVATE( a, b, k, uv_total, z_p )
     !$acc kernels do
+    !$acc kernels loop
     DO  i = nxl, nxr
        DO  j = nys, nyn
 …
 !-- First compute the corresponding component of u* and square it.
     !$OMP PARALLEL DO PRIVATE( a, b, k, rifm, z_p )
     !$acc kernels do
+    !$acc kernels loop
     DO  i = nxl, nxr
        DO  j = nys, nyn
 …
 !-- First compute the corresponding component of u* and square it.
     !$OMP PARALLEL DO PRIVATE( a, b, k, rifm, z_p )
     !$acc kernels do
+    !$acc kernels loop
     DO  i = nxl, nxr
        DO  j = nys, nyn
 …
 !--       For a given water flux in the Prandtl layer:
           !$OMP PARALLEL DO
           !$acc kernels do
+          !$acc kernels loop
           DO  i = nxlg, nxrg
              DO  j = nysg, nyng
 …
           coupled_run = ( coupling_mode == 'atmosphere_to_ocean' .AND. run_coupled )
           !$OMP PARALLEL DO PRIVATE( a, b, k, z_p )
           !$acc kernels do
+          !$acc kernels loop independent
           DO  i = nxlg, nxrg
+             !$acc loop independent
              DO  j = nysg, nyng
 …
     IF ( .NOT. constant_heatflux )  THEN
        !$OMP PARALLEL DO
+       !$acc kernels do
+       DO  i = nxlg, nxrg
+       !$acc kernels loop independent
+       DO  i = nxlg, nxrg
+          !$acc loop independent
           DO  j = nysg, nyng
              shf(j,i) = -ts(j,i) * us(j,i)
 …
     IF ( .NOT. constant_waterflux .AND. ( humidity .OR. passive_scalar ) ) THEN
        !$OMP PARALLEL DO
+       !$acc kernels do
+       DO  i = nxlg, nxrg
+       !$acc kernels loop independent
+       DO  i = nxlg, nxrg
+          !$acc loop independent
           DO  j = nysg, nyng
              qsws(j,i) = -qs(j,i) * us(j,i)
 …
     IF ( ibc_e_b == 2 )  THEN
        !$OMP PARALLEL DO
+       !$acc kernels do
+       DO  i = nxlg, nxrg
+       !$acc kernels loop independent
+       DO  i = nxlg, nxrg
+          !$acc loop independent
           DO  j = nysg, nyng
              e(nzb_s_inner(j,i)+1,j,i) = ( us(j,i) / 0.1 )**2

palm/trunk/SOURCE/pres.f90

-                      r1222
+                      r1257
 ! Current revisions:
 ! ------------------
+!
+! openacc loop and loop vector clauses removed, independent clauses added,
+! end parallel replaced by end parallel loop
+!
 ! Former revisions:
 …
        !$OMP PARALLEL DO SCHEDULE( STATIC )
        !$acc kernels present( d )
-       !$acc loop
        DO  i = nxl, nxr
           DO  j = nys, nyn
-             !$acc loop vector(32)
              DO  k = nzb+1, nzt
                 d(k,j,i) = 0.0
 …
        ENDDO
     ENDDO
     !$acc end parallel
+    !$acc end parallel loop
     localsum = localsum + threadsum * dt_3d * &
                           weight_pres(intermediate_timestep_count)
 …
        !$OMP PARALLEL DO
        !$acc kernels present( d, tend )
-       !$acc loop
        DO  i = nxl, nxr
           DO  j = nys, nyn
-             !$acc loop vector( 32 )
              DO  k = nzb+1, nzt
                 tend(k,j,i) = d(k,j,i)
 …
        !$OMP DO
        !$acc kernels present( p, tend, weight_substep )
        !$acc loop
+       !$acc loop independent
        DO  i = nxl-1, nxr+1
+          !$acc loop independent
           DO  j = nys-1, nyn+1
              !$acc loop vector( 32 )
+             !$acc loop independent
              DO  k = nzb, nzt+1
                 p(k,j,i) = tend(k,j,i) * &
 …
        !$OMP DO
        !$acc kernels present( p, tend, weight_substep )
        !$acc loop
+       !$acc loop independent
        DO  i = nxl-1, nxr+1
+          !$acc loop independent
           DO  j = nys-1, nyn+1
              !$acc loop vector( 32 )
+             !$acc loop independent
              DO  k = nzb, nzt+1
                 p(k,j,i) = p(k,j,i) + tend(k,j,i) * &
 …
     !$OMP DO
     !$acc kernels present( ddzu, nzb_u_inner, nzb_v_inner, nzb_w_inner, tend, u, v, w, weight_pres )
     !$acc loop
+    !$acc loop independent
     DO  i = nxl, nxr
+       DO  j = nys, nyn
+          !$acc loop vector( 32 )
+       !$acc loop independent
+       DO  j = nys, nyn
+          !$acc loop independent
           DO  k = 1, nzt
              IF ( k > nzb_w_inner(j,i) )  THEN
 …
              ENDIF
           ENDDO
           !$acc loop vector( 32 )
+          !$acc loop independent
           DO  k = 1, nzt
              IF ( k > nzb_u_inner(j,i) )  THEN
 …
              ENDIF
           ENDDO
           !$acc loop vector( 32 )
+          !$acc loop independent
           DO  k = 1, nzt
              IF ( k > nzb_v_inner(j,i) )  THEN
 …
        ENDDO
     ENDDO
     !$acc end parallel
+    !$acc end parallel loop
 #endif

palm/trunk/SOURCE/production_e.f90

-                      r1182
+                      r1257
 ! Current revisions:
 ! -----------------
+!
+! openacc loop and loop vector clauses removed, declare create moved after
+! the FORTRAN declaration statement
+!
 ! Former revisions:
 …
                    k1, k2, km_neutral, theta, temp
+       REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) ::  usvs, vsus, wsus, wsvs
        !$acc declare create ( usvs, vsus, wsus, wsvs )
-       REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) ::  usvs, vsus, wsus, wsvs
-!       REAL, DIMENSION(nzb:nzt+1) ::   usvs, vsus, wsus, wsvs
+!
 …
        !$acc         present( shf, tend, tswst, u, v, vpt, w, wall_e_x, wall_e_y )      &
        !$acc         copyin( u_0, v_0 )
-       !$acc loop
        DO  i = i_left, i_right
           DO  j = j_south, j_north
-             !$acc loop vector( 32 )
              DO  k = 1, nzt
 …
 !--       (2) - Will allways be executed.
 !--       'bottom and wall: use u_0,v_0 and wall functions'
-          !$acc loop
           DO  i = i_left, i_right
              DO  j = j_south, j_north
-                !$acc loop vector( 32 )
                 DO  k = 1, nzt
 …
 !--       (1) - will allways be executed.
 !--       'bottom only: use u_0,v_0'
-          !$acc loop
           DO  i = i_left, i_right
              DO  j = j_south, j_north
-                !$acc loop vector( 32 )
                 DO  k = 1, nzt
 …
        ELSEIF ( use_surface_fluxes )  THEN
-          !$acc loop
           DO  i = i_left, i_right
              DO  j = j_south, j_north
+                !$acc loop vector(32)
+                DO  k = 1, nzt
+                 DO  k = 1, nzt
                    IF ( k == nzb_diff_s_outer(j,i)-1 )  THEN
 …
 !--                So far in the ocean no special treatment of density flux
 !--                in the bottom and top surface layer
-                   !$acc loop
                    DO  i = i_left, i_right
                       DO  j = j_south, j_north
-                         !$acc loop vector( 32 )
                          DO  k = 1, nzt
                             IF ( k > nzb_s_inner(j,i) )  THEN
 …
                 ELSE
-                   !$acc loop
                    DO  i = i_left, i_right
                       DO  j = j_south, j_north
-                         !$acc loop vector( 32 )
                          DO  k = 1, nzt_diff
                             IF ( k >= nzb_diff_s_inner(j,i) )  THEN
 …
 !--                So far in the ocean no special treatment of density flux
 !--                in the bottom and top surface layer
-                   !$acc loop
                    DO  i = i_left, i_right
                       DO  j = j_south, j_north
-                         !$acc loop vector( 32 )
                          DO  k = 1, nzt
                             IF ( k > nzb_s_inner(j,i) )  THEN
 …
                 ELSE
-                   !$acc loop
                    DO  i = i_left, i_right
                       DO  j = j_south, j_north
-                         !$acc loop vector( 32 )
                          DO  k = 1, nzt_diff
                             IF( k >= nzb_diff_s_inner(j,i) )  THEN
 …
 !             DO  i = nxl, nxr
 !                DO  j = nys, nyn
 !                   !acc loop vector( 32 )
+!                   !acc loop vector
 !                   DO  k = 1, nzt_diff
+!
 …
 !                DO  i = nxl, nxr
 !                   DO  j = nys, nyn
 !                      !acc loop vector( 32 )
+!                      !acc loop vector
 !                      DO  k = 1, nzt_diff
+!
 …
 !                DO  i = nxl, nxr
 !                   DO  j = nys, nyn
 !                      !acc loop vector( 32 )
+!                      !acc loop vector
 !                      DO  k = 1, nzt
 !                         IF ( k == nzt )  THEN

palm/trunk/SOURCE/prognostic_equations.f90

-                      r1247
+                      r1257
 ! Current revisions:
 ! ------------------
+!
+! openacc loop vector clauses removed, independent clauses added
+!
 ! Former revisions:
 …
 !-- Prognostic equation for u-velocity component
     !$acc kernels present( nzb_u_inner, rdf, tend, tu_m, u, ug, u_p )
     !$acc loop
+    !$acc loop independent
     DO  i = i_left, i_right
+       !$acc loop independent
        DO  j = j_south, j_north
           !$acc loop vector( 32 )
+          !$acc loop independent
           DO  k = 1, nzt
              IF ( k > nzb_u_inner(j,i) )  THEN
 …
 !-- Prognostic equation for v-velocity component
     !$acc kernels present( nzb_v_inner, rdf, tend, tv_m, v, vg, v_p )
     !$acc loop
+    !$acc loop independent
     DO  i = i_left, i_right
+       !$acc loop independent
        DO  j = j_south, j_north
           !$acc loop vector( 32 )
+          !$acc loop independent
           DO  k = 1, nzt
              IF ( k > nzb_v_inner(j,i) )  THEN
 …
 !-- Prognostic equation for w-velocity component
     !$acc kernels present( nzb_w_inner, rdf, tend, tw_m, w, w_p )
     !$acc loop
+    !$acc loop independent
     DO  i = i_left, i_right
+       !$acc loop independent
        DO  j = j_south, j_north
           !$acc loop vector( 32 )
+          !$acc loop independent
           DO  k = 1, nzt-1
              IF ( k > nzb_w_inner(j,i) )  THEN
 …
        !$acc kernels present( nzb_s_inner, rdf_sc, ptdf_x, ptdf_y, pt_init ) &
        !$acc         present( tend, tpt_m, pt, pt_p )
        !$acc loop
+       !$acc loop independent
        DO  i = i_left, i_right
+          !$acc loop independent
           DO  j = j_south, j_north
              !$acc loop vector( 32 )
+             !$acc loop independent
              DO  k = 1, nzt
                 IF ( k > nzb_s_inner(j,i) )  THEN
 …
 !--    value is reduced by 90%.
        !$acc kernels present( e, e_p, nzb_s_inner, tend, te_m )
        !$acc loop
+       !$acc loop independent
        DO  i = i_left, i_right
+          !$acc loop independent
           DO  j = j_south, j_north
              !$acc loop vector( 32 )
+             !$acc loop independent
              DO  k = 1, nzt
                 IF ( k > nzb_s_inner(j,i) )  THEN

palm/trunk/SOURCE/run_control.f90

-                      r1037
+                      r1257
 ! Current revisions:
 ! -----------------
+!
+! output format for theta* changed to avoid output of *****
+!
 ! Former revisions:
 …
               &'------------------'// &
           &'RUN  ITER. HH:MM:SS.SS    DT(E)     UMAX     VMAX     WMAX     U', &
           &'*    W*   THETA*    Z_I     ENERG.   DISTENERG    DIVOLD     DIV', &
           &'NEW     UMAX(KJI)    VMAX(KJI)    WMAX(KJI)   ADVECX   ADVECY   ', &
           &'MGCYC'/                                                            &
+          &'*    W*      THETA*     Z_I     ENERG.   DISTENERG    DIVOLD    ', &
+          &' DIVNEW     UMAX(KJI)    VMAX(KJI)    WMAX(KJI)   ADVECX   ADVEC', &
+          &'Y   MGCYC'/                                                        &
           &'----------------------------------------------------------------', &
           &'----------------------------------------------------------------', &
           &'----------------------------------------------------------------', &
           &'-----')
+          &'---------')
 FORMAT (I3,1X,I6,1X,A8,F3.2,1X,F8.4,A1,1X,F8.4,A1,F8.4,A1,F8.4,2X,F5.3,2X, &
             F4.2, &
 X,F6.3,2X,F6.0,1X,4(E10.3,1X),3(3(I4),1X),F8.3,1X,F8.3,5X,I3)
+X,E10.3,2X,F6.0,1X,4(E10.3,1X),3(3(I4),1X),F8.3,1X,F8.3,5X,I3)
  END SUBROUTINE run_control

palm/trunk/SOURCE/time_integration.f90

-                      r1242
+                      r1257
 ! Current revisions:
 ! ------------------
+!
+! acc-update-host directive for timestep removed
+!
 ! Former revisions:
 …
 !--    Determine size of next time step
        IF ( simulated_time /= 0.0 )  THEN
-          !$acc update host( kh, km, u, v, w )
           CALL timestep
        ENDIF

palm/trunk/SOURCE/timestep.f90

-                      r1093
+                      r1257
 ! Current revisions:
 ! ------------------
+!
+! openacc porting
+! bugfix for calculation of advective timestep in case of vertically stretched
+! grids
+!
 ! Former revisions:
 …
     IMPLICIT NONE
+    INTEGER ::  i, j, k, u_max_cfl_ijk(3), v_max_cfl_ijk(3)
+    REAL ::  div, dt_diff, dt_diff_l, dt_plant_canopy, dt_plant_canopy_l,  &
+             dt_plant_canopy_u, dt_plant_canopy_v, dt_plant_canopy_w,      &
+             dt_u, dt_v, dt_w, u_max_cfl, value, v_max_cfl
+    INTEGER ::  i, j, k
+    REAL ::  div, dt_diff, dt_diff_l, dt_plant_canopy, dt_plant_canopy_l,     &
+             dt_plant_canopy_u, dt_plant_canopy_v, dt_plant_canopy_w,         &
+             dt_u, dt_u_l, dt_v, dt_v_l, dt_w, dt_w_l, u_gtrans_l, u_max_l,   &
+             u_min_l, value, v_gtrans_l, v_max_l, v_min_l, w_max_l, w_min_l
     REAL, DIMENSION(2)         ::  uv_gtrans, uv_gtrans_l
+    REAL, DIMENSION(3)         ::  reduce, reduce_l
     REAL, DIMENSION(nzb+1:nzt) ::  dxyz2_min
 …
+!
 !--       Averaging over the entire model domain.
+          uv_gtrans_l = 0.0
+          u_gtrans_l = 0.0
+          v_gtrans_l = 0.0
+          !$acc parallel present( u, v )
           DO  i = nxl, nxr
              DO  j = nys, nyn
                 DO  k = nzb+1, nzt
                    uv_gtrans_l(1) = uv_gtrans_l(1) + u(k,j,i)
                    uv_gtrans_l(2) = uv_gtrans_l(2) + v(k,j,i)
+                   u_gtrans_l = u_gtrans_l + u(k,j,i)
+                   v_gtrans_l = v_gtrans_l + v(k,j,i)
                 ENDDO
              ENDDO
           ENDDO
+          uv_gtrans_l = uv_gtrans_l / REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb) )
+          !$acc end parallel
+          uv_gtrans_l(1) = u_gtrans_l / REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb) )
+          uv_gtrans_l(2) = v_gtrans_l / REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb) )
 #if defined( __parallel )
           IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
 …
+!
+!-- Determine the maxima of the velocity components.
+!-- Determine the maxima of the velocity components, including their
+!-- grid index positions.
+#if defined( __openacc )
+    IF ( dt_fixed )  THEN  ! otherwise do it further below for better cache usage
+       u_max_l = -999999.9
+       u_min_l =  999999.9
+       v_max_l = -999999.9
+       v_min_l =  999999.9
+       w_max_l = -999999.9
+       w_min_l =  999999.9
+       !$acc parallel present( u, v, w )
+       DO  i = nxl, nxr
+          DO  j = nys, nyn
+             DO  k = nzb+1, nzt
+                u_max_l = MAX( u_max_l, u(k,j,i) )
+                u_min_l = MIN( u_min_l, u(k,j,i) )
+                v_max_l = MAX( v_max_l, v(k,j,i) )
+                v_min_l = MIN( v_min_l, v(k,j,i) )
+                w_max_l = MAX( w_max_l, w(k,j,i) )
+                w_min_l = MIN( w_min_l, w(k,j,i) )
+             ENDDO
+          ENDDO
+       ENDDO
+       !$acc end parallel
+#if defined( __parallel )
+       reduce_l(1) = u_max_l
+       reduce_l(2) = v_max_l
+       reduce_l(3) = w_max_l
+       IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
+       CALL MPI_ALLREDUCE( reduce_l, reduce, 3, MPI_REAL, MPI_MAX, comm2d, ierr )
+       u_max = reduce(1)
+       v_max = reduce(2)
+       w_max = reduce(3)
+       reduce_l(1) = u_min_l
+       reduce_l(2) = v_min_l
+       reduce_l(3) = w_min_l
+       IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
+       CALL MPI_ALLREDUCE( reduce_l, reduce, 3, MPI_REAL, MPI_MIN, comm2d, ierr )
+       IF ( ABS( reduce(1) ) > u_max )  u_max = reduce(1)
+       IF ( ABS( reduce(2) ) > v_max )  v_max = reduce(2)
+       IF ( ABS( reduce(3) ) > w_max )  w_max = reduce(3)
+#else
+       IF ( ABS( u_min_l ) > u_max_l )  THEN
+          u_max = u_min_l
+       ELSE
+          u_max = u_max_l
+       ENDIF
+       IF ( ABS( v_min_l ) > v_max_l )  THEN
+          v_max = v_min_l
+       ELSE
+          v_max = v_max_l
+       ENDIF
+       IF ( ABS( w_min_l ) > w_max_l )  THEN
+          w_max = w_min_l
+       ELSE
+          w_max = w_max_l
+       ENDIF
+#endif
+    ENDIF
+#else
     CALL global_min_max( nzb, nzt+1, nysg, nyng, nxlg, nxrg, u, 'abs', 0.0, &
                          u_max, u_max_ijk )
 …
     CALL global_min_max( nzb, nzt+1, nysg, nyng, nxlg, nxrg, w, 'abs', 0.0, &
                          w_max, w_max_ijk )
+!
+!-- In case of Galilei transformation, the horizontal velocity maxima have
+!-- to be calculated from the transformed horizontal velocities
+    IF ( galilei_transformation )  THEN
+       CALL global_min_max( nzb, nzt+1, nysg, nyng, nxlg, nxrg, u, 'absoff', &
+                            u_gtrans, u_max_cfl, u_max_cfl_ijk )
+       CALL global_min_max( nzb, nzt+1, nysg, nyng, nxlg, nxrg, v, 'absoff', &
+                            v_gtrans, v_max_cfl, v_max_cfl_ijk )
+    ELSE
+       u_max_cfl = u_max
+       v_max_cfl = v_max
+       u_max_cfl_ijk = u_max_ijk
+       v_max_cfl_ijk = v_max_ijk
+    ENDIF
+#endif
     IF ( .NOT. dt_fixed )  THEN
+#if defined( __openacc )
+!
 !--    Variable time step:
+!
+!--    For each component, compute the maximum time step according to the
+!--    CFL-criterion.
+       dt_u = dx / ( ABS( u_max_cfl ) + 1.0E-10 )
+       dt_v = dy / ( ABS( v_max_cfl ) + 1.0E-10 )
+       dt_w = dzu(MAX( 1, w_max_ijk(1) )) / ( ABS( w_max ) + 1.0E-10 )
+!--    Calculate the maximum time step according to the CFL-criterion,
+!--    individually for each velocity component
+       dt_u_l  =  999999.9
+       dt_v_l  =  999999.9
+       dt_w_l  =  999999.9
+       u_max_l = -999999.9
+       u_min_l =  999999.9
+       v_max_l = -999999.9
+       v_min_l =  999999.9
+       w_max_l = -999999.9
+       w_min_l =  999999.9
+       !$acc parallel loop collapse(3) present( u, v, w )
+       DO  i = nxl, nxr
+          DO  j = nys, nyn
+             DO  k = nzb+1, nzt
+                dt_u_l  = MIN( dt_u_l, ( dx     / ( ABS( u(k,j,i) - u_gtrans ) + 1.0E-10 ) ) )
+                dt_v_l  = MIN( dt_v_l, ( dy     / ( ABS( v(k,j,i) - v_gtrans ) + 1.0E-10 ) ) )
+                dt_w_l  = MIN( dt_w_l, ( dzu(k) / ( ABS( w(k,j,i) )            + 1.0E-10 ) ) )
+                u_max_l = MAX( u_max_l, u(k,j,i) )
+                u_min_l = MIN( u_min_l, u(k,j,i) )
+                v_max_l = MAX( v_max_l, v(k,j,i) )
+                v_min_l = MIN( v_min_l, v(k,j,i) )
+                w_max_l = MAX( w_max_l, w(k,j,i) )
+                w_min_l = MIN( w_min_l, w(k,j,i) )
+             ENDDO
+          ENDDO
+       ENDDO
+       !$acc end parallel
+#if defined( __parallel )
+       reduce_l(1) = dt_u_l
+       reduce_l(2) = dt_v_l
+       reduce_l(3) = dt_w_l
+       IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
+       CALL MPI_ALLREDUCE( reduce_l, reduce, 3, MPI_REAL, MPI_MIN, comm2d, ierr )
+       dt_u = reduce(1)
+       dt_v = reduce(2)
+       dt_w = reduce(3)
+       reduce_l(1) = u_max_l
+       reduce_l(2) = v_max_l
+       reduce_l(3) = w_max_l
+       IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
+       CALL MPI_ALLREDUCE( reduce_l, reduce, 3, MPI_REAL, MPI_MAX, comm2d, ierr )
+       u_max = reduce(1)
+       v_max = reduce(2)
+       w_max = reduce(3)
+       reduce_l(1) = u_min_l
+       reduce_l(2) = v_min_l
+       reduce_l(3) = w_min_l
+       IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
+       CALL MPI_ALLREDUCE( reduce_l, reduce, 3, MPI_REAL, MPI_MIN, comm2d, ierr )
+       IF ( ABS( reduce(1) ) > u_max )  u_max = reduce(1)
+       IF ( ABS( reduce(2) ) > v_max )  v_max = reduce(2)
+       IF ( ABS( reduce(3) ) > w_max )  w_max = reduce(3)
+#else
+       dt_u = dt_u_l
+       dt_v = dt_v_l
+       dt_w = dt_w_l
+       IF ( ABS( u_min_l ) > u_max_l )  THEN
+          u_max = u_min_l
+       ELSE
+          u_max = u_max_l
+       ENDIF
+       IF ( ABS( v_min_l ) > v_max_l )  THEN
+          v_max = v_min_l
+       ELSE
+          v_max = v_max_l
+       ENDIF
+       IF ( ABS( w_min_l ) > w_max_l )  THEN
+          w_max = w_min_l
+       ELSE
+          w_max = w_max_l
+       ENDIF
+#endif
+#else
+!
+!--    Variable time step:
+!--    Calculate the maximum time step according to the CFL-criterion,
+!--    individually for each velocity component
+       dt_u_l = 999999.9
+       dt_v_l = 999999.9
+       dt_w_l = 999999.9
+       DO  i = nxl, nxr
+          DO  j = nys, nyn
+             DO  k = nzb+1, nzt
+                dt_u_l = MIN( dt_u_l, ( dx     / ( ABS( u(k,j,i) - u_gtrans ) + 1.0E-10 ) ) )
+                dt_v_l = MIN( dt_v_l, ( dy     / ( ABS( v(k,j,i) - v_gtrans ) + 1.0E-10 ) ) )
+                dt_w_l = MIN( dt_w_l, ( dzu(k) / ( ABS( w(k,j,i) )            + 1.0E-10 ) ) )
+             ENDDO
+          ENDDO
+       ENDDO
+#if defined( __parallel )
+       reduce_l(1) = dt_u_l
+       reduce_l(2) = dt_v_l
+       reduce_l(3) = dt_w_l
+       IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
+       CALL MPI_ALLREDUCE( reduce_l, reduce, 3, MPI_REAL, MPI_MIN, comm2d, ierr )
+       dt_u = reduce(1)
+       dt_v = reduce(2)
+       dt_w = reduce(3)
+#else
+       dt_u = dt_u_l
+       dt_v = dt_v_l
+       dt_w = dt_w_l
+#endif
+#endif
+!
 …
 !$OMP PARALLEL private(i,j,k,value) reduction(MIN: dt_diff_l)
 !$OMP DO
+       !$acc parallel loop collapse(3) present( kh, km )
        DO  i = nxl, nxr
           DO  j = nys, nyn
              DO  k = nzb+1, nzt
+                value = dxyz2_min(k) / ( MAX( kh(k,j,i), km(k,j,i) ) + 1E-20 )
+                dt_diff_l = MIN( value, dt_diff_l )
+                dt_diff_l = MIN( dt_diff_l, dxyz2_min(k) / &
+                                       ( MAX( kh(k,j,i), km(k,j,i) ) + 1E-20 ) )
              ENDDO
           ENDDO
        ENDDO
+       !$acc end parallel
 !$OMP END PARALLEL
 #if defined( __parallel )
 …
                '&dt_diff         = ', dt_diff, ' s',                          &
                '&dt_plant_canopy = ', dt_plant_canopy, ' s',                  &
                '&u_max_cfl   = ', u_max_cfl, ' m/s   k=', u_max_cfl_ijk(1),   &
+               '&u_max           = ', u_max, ' m/s   k=', u_max_ijk(1),       &
                '  j=', u_max_ijk(2), '  i=', u_max_ijk(3),                    &
                '&v_max_cfl   = ', v_max_cfl, ' m/s   k=', v_max_cfl_ijk(1),   &
+               '&v_max           = ', v_max, ' m/s   k=', v_max_ijk(1),       &
                '  j=', v_max_ijk(2), '  i=', v_max_ijk(3),                    &
                '&w_max       = ', w_max, ' m/s   k=', w_max_ijk(1),           &
+               '&w_max           = ', w_max, ' m/s   k=', w_max_ijk(1),       &
                '  j=', w_max_ijk(2), '  i=', w_max_ijk(3)
           CALL message( 'timestep', 'PA0312', 0, 1, 0, 6, 0 )

palm/trunk/SOURCE/transpose.f90

-                      r1217
+                      r1257
 ! Current revisions:
 ! -----------------
+!
+! openacc loop and loop vector clauses removed
+!
 ! Former revisions:
 …
     !$OMP  DO
     !$acc kernels present( f_in, f_inv )
-    !$acc loop
      DO  i = 0, nx
          DO  k = nzb_x, nzt_x
-             !$acc loop vector( 32 )
              DO  j = nys_x, nyn_x
                  f_inv(j,k,i) = f_in(i,j,k)
 …
           ys = 0 + l * ( nyn_x - nys_x + 1 )
           !$acc kernels present( f_out, work )
-          !$acc loop
           DO  i = nxl_y, nxr_y
              DO  k = nzb_y, nzt_y
-                !$acc loop vector( 32 )
                 DO  j = ys, ys + nyn_x - nys_x
                    f_out(j,i,k) = work(j-ys+1,k,i,l)
 …
 !$OMP  DO
        !$acc kernels present( f_inv, f_out )
-       !$acc loop
        DO  k = nzb_y, nzt_y
           DO  i = nxl_y, nxr_y
-             !$acc loop vector( 32 )
              DO  j = 0, ny
                 f_out(j,i,k) = f_inv(j,k,i)
 …
     !$OMP  DO
     !$acc kernels present( f_inv, f_out )
-    !$acc loop
      DO  k = 1, nz
          DO  i = nxl, nxr
-             !$acc loop vector( 32 )
              DO  j = nys, nyn
                  f_out(k,j,i) = f_inv(j,i,k)
 …
           xs = 0 + l * nnx
           !$acc kernels present( f_in, work )
-          !$acc loop
           DO  k = nzb_x, nzt_x
              DO  i = xs, xs + nnx - 1
-                !$acc loop vector( 32 )
                 DO  j = nys_x, nyn_x
                    work(j,i-xs+1,k,l) = f_in(i,j,k)
 …
 !$OMP  DO
        !$acc kernels present( f_in, f_inv )
-       !$acc loop
        DO  i = nxl, nxr
           DO  j = nys, nyn
-             !$acc loop vector( 32 )
              DO  k = 1, nz
                 f_inv(j,i,k) = f_in(i,j,k)
 …
     !$OMP  DO
     !$acc kernels present( f_inv, f_out )
-    !$acc loop
      DO  i = 0, nx
          DO  k = nzb_x, nzt_x
-             !$acc loop vector( 32 )
              DO  j = nys_x, nyn_x
                  f_out(i,j,k) = f_inv(j,k,i)
 …
           ys = 0 + l * ( nyn_x - nys_x + 1 )
           !$acc kernels present( f_in, work )
-          !$acc loop
           DO  i = nxl_y, nxr_y
              DO  k = nzb_y, nzt_y
-                !$acc loop vector( 32 )
                 DO  j = ys, ys + nyn_x - nys_x
                    work(j-ys+1,k,i,l) = f_in(j,i,k)
 …
 !$OMP  DO
        !$acc kernels present( f_in, f_inv )
-       !$acc loop
        DO  i = nxl_y, nxr_y
           DO  k = nzb_y, nzt_y
-             !$acc loop vector( 32 )
              DO  j = 0, ny
                 f_inv(j,k,i) = f_in(j,i,k)
 …
     !$OMP  DO
     !$acc kernels present( f_in, f_inv )
-    !$acc loop
      DO  j = 0, ny
          DO  k = nzb_y, nzt_y
-             !$acc loop vector( 32 )
              DO  i = nxl_y, nxr_y
                  f_inv(i,k,j) = f_in(j,i,k)
 …
 !$OMP  DO
        !$acc kernels present( f_inv, f_out )
-       !$acc loop
        DO  j = 0, ny
           DO  k = nzb_y, nzt_y
-             !$acc loop vector( 32 )
              DO  i = nxl_y, nxr_y
                 f_out(i,j,k) = f_inv(i,k,j)
 …
           zs = 1 + l * ( nzt_y - nzb_y + 1 )
           !$acc kernels present( f_out )
-          !$acc loop
           DO  j = nys_z, nyn_z
              DO  k = zs, zs + nzt_y - nzb_y
-                !$acc loop vector( 32 )
                 DO  i = nxl_z, nxr_z
                    f_out(i,j,k) = work(i,k-zs+1,j,l)
 …
     !$OMP  DO
     !$acc kernels present( f_in, f_inv )
-    !$acc loop
      DO  k = 1,nz
          DO  i = nxl, nxr
-             !$acc loop vector( 32 )
              DO  j = nys, nyn
                  f_inv(j,i,k) = f_in(k,j,i)
 …
 !$OMP  DO
        !$acc kernels present( f_inv, f_out )
-       !$acc loop
        DO  k = 1, nz
           DO  i = nxl, nxr
-             !$acc loop vector( 32 )
              DO  j = nys, nyn
                 f_out(i,j,k) = f_inv(j,i,k)
 …
           xs = 0 + l * nnx
           !$acc kernels present( f_out )
-          !$acc loop
           DO  k = nzb_x, nzt_x
              DO  i = xs, xs + nnx - 1
-                !$acc loop vector( 32 )
                 DO  j = nys_x, nyn_x
                    f_out(i,j,k) = work(j,i-xs+1,k,l)
 …
     !$OMP  DO
     !$acc kernels present( f_inv, f_out )
-    !$acc loop
      DO  k = nzb_y, nzt_y
          DO  j = 0, ny
-             !$acc loop vector( 32 )
              DO  i = nxl_y, nxr_y
                  f_out(j,i,k) = f_inv(i,k,j)
 …
           zs = 1 + l * ( nzt_y - nzb_y + 1 )
           !$acc kernels present( f_in, work )
-          !$acc loop
           DO  j = nys_z, nyn_z
              DO  k = zs, zs + nzt_y - nzb_y
-                !$acc loop vector( 32 )
                 DO  i = nxl_z, nxr_z
                    work(i,k-zs+1,j,l) = f_in(i,j,k)
 …
 !$OMP  DO
        !$acc kernels present( f_in, f_inv )
-       !$acc loop
        DO  k = nzb_y, nzt_y
           DO  j = 0, ny
-             !$acc loop vector( 32 )
              DO  i = nxl_y, nxr_y
                 f_inv(i,k,j) = f_in(i,j,k)

palm/trunk/SOURCE/tridia_solver.f90

-                      r1222
+                      r1257
 ! Current revisions:
 ! ------------------
+!
+! openacc loop and loop vector clauses removed, declare create moved after
+! the FORTRAN declaration statement
+!
 ! Former revisions:
 …
           INTEGER ::  i, j, k, nnxh, nnyh
+          REAL    ::  ll(nxl_z:nxr_z,nys_z:nyn_z)
           !$acc declare create( ll )
-          REAL    ::  ll(nxl_z:nxr_z,nys_z:nyn_z)
 …
           !$acc kernels present( tric )
-          !$acc loop vector( 32 )
           DO  j = nys_z, nyn_z
              DO  i = nxl_z, nxr_z
 …
           ENDDO
-          !$acc loop
           DO  k = 0, nz-1
              DO  j = nys_z, nyn_z
-                !$acc loop vector( 32 )
                 DO  i = nxl_z, nxr_z
                    tric(i,j,k) = ddzuw(k,3) - ll(i,j)
 …
           IF ( ibc_p_b == 1 )  THEN
              !$acc kernels present( tric )
-             !$acc loop
              DO  j = nys_z, nyn_z
                 DO  i = nxl_z, nxr_z
 …
           IF ( ibc_p_t == 1 )  THEN
              !$acc kernels present( tric )
-             !$acc loop
              DO  j = nys_z, nyn_z
                 DO  i = nxl_z, nxr_z
 …
           REAL    ::  ar(nxl_z:nxr_z,nys_z:nyn_z,1:nz)
+          REAL, DIMENSION(nxl_z:nxr_z,nys_z:nyn_z,0:nz-1)   ::  ar1
           !$acc declare create( ar1 )
-          REAL, DIMENSION(nxl_z:nxr_z,nys_z:nyn_z,0:nz-1)   ::  ar1
+!
 …
           DO  k = 0, nz - 1
              !$acc kernels present( ar, tri )
-             !$acc loop
              DO  j = nys_z, nyn_z
                 DO  i = nxl_z, nxr_z
 …
           DO  k = nz-1, 0, -1
              !$acc kernels present( ar, tri )
-             !$acc loop
              DO  j = nys_z, nyn_z
                 DO  i = nxl_z, nxr_z
 …
                    ar(nxl_z,nys_z,k) = 0.0
                 ENDDO
+                !$acc end kernels loop
              ENDIF
           ENDIF

palm/trunk/SOURCE/wall_fluxes.f90

-                      r1154
+                      r1257
 ! Current revisions:
 ! -----------------
+!
+! openacc loop and loop vector clauses removed
+!
 ! Former revisions:
 …
 !--             All subsequent variables are computed for the respective
 !--             location where the respective flux is defined.
                 !$acc loop independent vector( 32 )
+                !$acc loop independent
                 DO  k = nzb_uvw_inner(j,i)+1, nzb_uvw_outer(j,i)
 …
        !$acc kernels present( nzb_diff_s_inner, nzb_diff_s_outer, pt, rif_wall ) &
        !$acc         present( u, v, w, wall, wall_flux, z0 )
-       !$acc loop
        DO  i = i_left, i_right
           DO  j = j_south, j_north
-             !$acc loop vector(32)
              DO  k = min_inner, max_outer
+!

Note: See TracChangeset for help on using the changeset viewer.

Context Navigation

Changeset 1257 for palm/trunk

Legend:

Download in other formats: