Changeset 3065


Ignore:
Timestamp:
Jun 12, 2018 7:03:02 AM (3 years ago)
Author:
Giersch
Message:

New vertical stretching procedure has been introduced

Location:
palm/trunk
Files:
22 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/INSTALL/example_cbl_rc

    r2696 r3065  
    11
    22 ******************************    --------------------------------------------
    3  * PALM 4.0  Rev: 2683        *    atmosphere - run without 1D - prerun
     3 * PALM 5.0  Rev: 3058M       *    atmosphere - run without 1D - prerun
    44 ******************************    --------------------------------------------
    55
    6  Date:                 12-12-17    Run:       example_cbl         
    7  Time:                 10:31:35    Run-No.:   00
    8  Run on host:              imuk
    9  Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
     6 Date:                 11-06-18    Run:       example_cbl         
     7 Time:                 15:43:48    Run-No.:   00
     8 Run on host:        imuk_defau
     9 Number of PEs:               1    Processor grid (x,y): (   1,   1) calculated
     10                                   A 1d-decomposition along x is used
    1011 ------------------------------------------------------------------------------
    1112
     
    3536 ----------------------------------
    3637
    37  Grid length:      dx =     50.000 m    dy =     50.000 m    dz =     50.000 m
    38  Domain size:       x =   2000.000 m     y =   2000.000 m  z(u) =   2862.511 m
    39 
    40  dz constant up to   1200.000 m (k=  25), then stretched by factor: 1.080
    41  maximum dz not to be exceeded is dz_max = ********** m
     38 Grid length:      dx =      50.000 m    dy =      50.000 m
     39                dz(1) =      50.000 m
     40                dz(2) =    1000.000 m
     41
     42 Domain size:       x =   2000.000 m     y =   2000.000 m  z(u) =   2997.512 m
     43
     44 Vertical stretching starts at height:    1175.0 m,
     45 Vertical stretching starts at index:           24,
     46 Vertical stretching ends at height:   9999999.9 m,
     47 Vertical stretching ends at index:             41,
     48 Factor used for stretching:                 1.080,
    4249
    4350 Number of gridpoints (x,y,z):  (0:  39, 0:  39, 0:  41)
    44  Subdomain size (x,y,z):        (    20,     20,     42)
     51 Subdomain size (x,y,z):        (    40,     40,     42)
    4552
    4653
     
    109116
    110117    1D-Profiles:
    111        Output every             ******** s
     118       Output every             9999999.90 s
    112119       Time averaged over         600.00 s
    113120       Averaging input every       10.00 s
     
    227234  0      1 00:00:20.00  20.0000D  -0.2125  -0.2273   0.1172  0.009  0.69  -0.132E+02    100.  0.103E-02  0.102E-02  0.570E-05  0.411E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    228235  0      2 00:00:40.00  20.0000D  -0.2117  -0.2252   0.1162  0.009  0.69  -0.133E+02    100.  0.102E-02  0.102E-02  0.935E-05  0.416E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    229   0      3 00:01:00.00  20.0000D  -0.2106  -0.2227   0.1152  0.009  0.69  -0.134E+02    100.  0.101E-02  0.101E-02  0.130E-04  0.421E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     236  0      3 00:01:00.00  20.0000D  -0.2106  -0.2227   0.1152  0.009  0.69  -0.134E+02    100.  0.101E-02  0.101E-02  0.130E-04  0.420E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    230237  0      4 00:01:20.00  20.0000D  -0.2093  -0.2197   0.1141  0.009  0.69  -0.134E+02    100.  0.100E-02  0.997E-03  0.166E-04  0.424E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    231238  0      5 00:01:40.00  20.0000D  -0.2077  -0.2164   0.1129  0.009  0.93  -0.134E+02    250.  0.992E-03  0.988E-03  0.203E-04  0.427E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    232239  0      6 00:02:00.00  20.0000D  -0.2058  -0.2126   0.1117  0.009  0.93  -0.133E+02    250.  0.983E-03  0.979E-03  0.239E-04  0.429E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    233   0      7 00:02:20.00  20.0000D  -0.2037  -0.2086   0.1103  0.009  1.22  -0.131E+02    550.  0.974E-03  0.970E-03  0.275E-04  0.431E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     240  0      7 00:02:20.00  20.0000D  -0.2037  -0.2086   0.1103  0.009  1.22  -0.131E+02    550.  0.974E-03  0.970E-03  0.275E-04  0.430E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    234241  0      8 00:02:40.00  20.0000D   0.3348D  0.3126D  0.1636  0.011  1.22  -0.105E+02    550.  0.202E-02  0.201E-02  0.601E-03  0.856E-13   10   3  39   10   8  27    5  36  28    0.000    0.000       0
    235   0      9 00:03:00.00  20.0000D   0.3330   0.3039   0.1619  0.011  1.22  -0.105E+02    550.  0.199E-02  0.198E-02  0.374E-04  0.830E-15   10   3  39   10   8  27    5  36  28    0.000    0.000       0
     242  0      9 00:03:00.00  20.0000D   0.3330   0.3039   0.1619  0.011  1.22  -0.105E+02    550.  0.199E-02  0.198E-02  0.374E-04  0.829E-15   10   3  39   10   8  27    5  36  28    0.000    0.000       0
    236243  0     10 00:03:20.00  20.0000D   0.3296  -0.2989   0.1601  0.011  1.22  -0.103E+02    550.  0.196E-02  0.195E-02  0.409E-04  0.827E-15   10   3  39    8   4   6    5  36  28    0.000    0.000       0
    237244  0     11 00:03:40.00  20.0000D   0.3249  -0.2986   0.1580  0.011  1.22  -0.101E+02    550.  0.194E-02  0.193E-02  0.444E-04  0.828E-15   10   3  39    8   4   6    5  36  28    0.000    0.000       0
    238245  0     12 00:04:00.00  20.0000D   0.3189  -0.2989   0.1558  0.012  1.22  -0.975E+01    550.  0.191E-02  0.190E-02  0.478E-04  0.831E-15   10   3  39    8   0   7    5  36  28    0.000    0.000       0
    239246  0     13 00:04:20.00  20.0000D   0.3119  -0.3005   0.1535  0.012  1.22  -0.940E+01    550.  0.189E-02  0.188E-02  0.513E-04  0.835E-15   10   3  39    8   0   7    5  36  28    0.000    0.000       0
    240   0     14 00:04:40.00  20.0000D   0.3056  -0.3012   0.1512  0.013  1.22  -0.905E+01    550.  0.187E-02  0.185E-02  0.546E-04  0.842E-15   10   3   0    8   0   7    5  36  28    0.000    0.000       0
     247  0     14 00:04:40.00  20.0000D   0.3056  -0.3012   0.1512  0.013  1.22  -0.905E+01    550.  0.187E-02  0.185E-02  0.546E-04  0.841E-15   10   3   0    8   0   7    5  36  28    0.000    0.000       0
    241248  0     15 00:05:00.00  20.0000D   0.4151D -0.3964D -0.1812  0.013  1.22  -0.851E+01    550.  0.292E-02  0.289E-02  0.625E-03  0.860E-13    8  17   5    8   3   6    8  19   3    0.000    0.000       0
    242249  0     16 00:05:20.00  20.0000D   0.4183  -0.4076  -0.1810  0.014  1.22  -0.811E+01    550.  0.287E-02  0.285E-02  0.640E-04  0.129E-14    8  17   5    8   3   6    8  19   3    0.000    0.000       0
     
    262269  0     36 00:12:00.00  20.0000A  -0.5522  -0.6460   1.4350  0.052  1.18  -0.218E+01    500.  0.920E-02  0.917E-02  0.112E-03  0.522E-14    1   3  23    1  29  17    3   6  12    0.000    0.000       0
    263270  0     37 00:12:20.00  20.0000A   0.6173  -0.7139   1.6131  0.055  1.18  -0.205E+01    500.  0.109E-01  0.108E-01  0.112E-03  0.591E-14    1  35  11    1  29  17    4   6  12    0.000    0.000       0
    264   0     38 00:12:40.00  20.0000A   0.6866  -0.7795   1.7862  0.059  1.18  -0.193E+01    500.  0.130E-01  0.130E-01  0.112E-03  0.673E-14    1  35  11    1  29  17    4   6  12    0.000    0.000       0
     271  0     38 00:12:40.00  20.0000A   0.6866  -0.7795   1.7862  0.059  1.18  -0.193E+01    500.  0.130E-01  0.130E-01  0.112E-03  0.672E-14    1  35  11    1  29  17    4   6  12    0.000    0.000       0
    265272  0     39 00:13:00.00  20.0000A   0.7574  -0.8398   1.9133  0.062  1.18  -0.183E+01    500.  0.158E-01  0.158E-01  0.112E-03  0.767E-14    1  35  11    1  29  17    4   6  12    0.000    0.000       0
    266273  0     40 00:13:20.00  20.0000A   0.8276  -0.8933   2.1010  0.066  1.18  -0.173E+01    500.  0.193E-01  0.193E-01  0.114E-03  0.877E-14    1  35  11    1  29  17    5   6  12    0.000    0.000       0
     
    284291  0     58 00:17:41.39  10.4000A  -1.4206  -1.3889   4.3431  0.094  1.38  -0.122E+01    800.  0.162E+00  0.162E+00  0.144E-03  0.204E-13    5  29  17    1  28  18   11  34  13    0.000    0.000       0
    285292  0     59 00:17:51.79  10.2000A  -1.4342   1.3908   4.3989  0.095  1.38  -0.121E+01    800.  0.170E+00  0.170E+00  0.146E-03  0.208E-13    5  29  17    7  37  19   11  34  13    0.000    0.000       0
    286   0     60 00:18:02.00  10.1000A  -1.4413   1.4902   4.4530  0.095  1.38  -0.120E+01    800.  0.177E+00  0.177E+00  0.148E-03  0.212E-13    4  32  19    7  37  19   12  34  13    0.000    0.000       0
     293  0     60 00:18:02.00  10.1000A  -1.4413   1.4901   4.4530  0.095  1.38  -0.120E+01    800.  0.177E+00  0.177E+00  0.148E-03  0.212E-13    4  32  19    7  37  19   12  34  13    0.000    0.000       0
    287294  0     61 00:18:12.09  10.0000A  -1.4694  -1.5081   4.4813  0.096  1.38  -0.120E+01    800.  0.185E+00  0.185E+00  0.151E-03  0.217E-13    4  32  19    8   7  16   12  34  13    0.000    0.000       0
    288295  0     62 00:18:22.09  10.0000A  -1.4854  -1.6266   4.4842  0.096  1.38  -0.119E+01    800.  0.193E+00  0.192E+00  0.154E-03  0.222E-13    4  32  19    8   7  16   13  34  13    0.000    0.000       0
    289296  0     63 00:18:32.09  10.0000A  -1.4915  -1.7186   4.4785  0.097  1.38  -0.118E+01    800.  0.200E+00  0.200E+00  0.158E-03  0.229E-13    4  32  19    8   7  16   12  34  13    0.000    0.000       0
    290297  0     64 00:18:42.09  10.1000A  -1.5053  -1.7631   4.4346  0.097  1.38  -0.116E+01    800.  0.207E+00  0.207E+00  0.163E-03  0.235E-13    4  32  18    8   7  16   14  34  13    0.000    0.000       0
    291   0     65 00:18:52.19  10.2000A  -1.5591  -1.7454   4.4046  0.098  1.41  -0.116E+01    850.  0.214E+00  0.214E+00  0.168E-03  0.244E-13    4  32  18    8   7  16   13  34  13    0.000    0.000       0
     298  0     65 00:18:52.19  10.2000A  -1.5591  -1.7454   4.4047  0.098  1.41  -0.116E+01    850.  0.214E+00  0.214E+00  0.168E-03  0.244E-13    4  32  18    8   7  16   13  34  13    0.000    0.000       0
    292299  0     66 00:19:02.39  10.4000A  -1.5902  -1.6634   4.3471  0.099  1.41  -0.115E+01    850.  0.221E+00  0.221E+00  0.175E-03  0.253E-13    4  32  18    8   7  16   13  34  13    0.000    0.000       0
    293300  0     67 00:19:12.79  10.5000A  -1.5977  -1.6328   4.2947  0.099  1.41  -0.115E+01    850.  0.228E+00  0.228E+00  0.182E-03  0.264E-13    4  32  18    9   7  16   15  35  12    0.000    0.000       0
     
    297304  0     71 00:19:55.39  11.0000A   1.6043  -1.8306   4.1022  0.102  1.41  -0.112E+01    850.  0.253E+00  0.252E+00  0.204E-03  0.295E-13    7   1  37   10   7  17   12  23  18    0.000    0.000       0
    298305  0     72 00:20:06.39  10.9000A   1.7066  -1.8092   4.1388  0.103  1.41  -0.111E+01    850.  0.258E+00  0.258E+00  0.214E-03  0.308E-13    7  24  30   10   7  17   12  23  18    0.000    0.000       0
    299   0     73 00:20:17.29  10.9000A   1.8396   1.7475   4.1429  0.104  1.41  -0.109E+01    850.  0.262E+00  0.262E+00  0.215E-03  0.309E-13    7  24  30   14  18  28   13  23  18    0.000    0.000       0
     306  0     73 00:20:17.29  10.9000A   1.8395   1.7475   4.1429  0.104  1.41  -0.109E+01    850.  0.262E+00  0.262E+00  0.215E-03  0.309E-13    7  24  30   14  18  28   13  23  18    0.000    0.000       0
    300307  0     74 00:20:28.20  10.7000A   1.9139   1.8281   4.1963  0.105  1.41  -0.107E+01    850.  0.267E+00  0.266E+00  0.219E-03  0.313E-13    7  24  30   14  18  28   13  23  18    0.000    0.000       0
    301308  0     75 00:20:38.90  10.7000A   1.9324  -1.8944   4.1945  0.106  1.41  -0.106E+01    850.  0.270E+00  0.270E+00  0.218E-03  0.311E-13    7  24  30    9   4   0   13  23  17    0.000    0.000       0
     
    305312  0     79 00:21:21.30  10.4000A   2.0495  -2.1041   4.3252  0.110  1.41  -0.102E+01    850.  0.279E+00  0.279E+00  0.223E-03  0.316E-13    7  24  31    9   4   0   14  23  17    0.000    0.000       0
    306313  0     80 00:21:31.70  10.4000A   2.0215  -2.0580   4.3215  0.111  1.41  -0.101E+01    850.  0.281E+00  0.280E+00  0.222E-03  0.315E-13    7  24  31    9   4   0   14  23  17    0.000    0.000       0
    307   0     81 00:21:42.10  10.5000A   1.9801  -1.9716   4.2997  0.112  1.41  -0.100E+01    850.  0.282E+00  0.281E+00  0.224E-03  0.316E-13    7  24  31    9   4   0   15  23  17    0.000    0.000       0
     314  0     81 00:21:42.10  10.5000A   1.9801  -1.9716   4.2996  0.112  1.41  -0.100E+01    850.  0.282E+00  0.281E+00  0.224E-03  0.316E-13    7  24  31    9   4   0   15  23  17    0.000    0.000       0
    308315  0     82 00:21:52.60  10.4000A   1.9327  -1.9463   4.3286  0.113  1.41  -0.982E+00    850.  0.283E+00  0.282E+00  0.227E-03  0.321E-13    7  24  31   10   4   0   15  23  17    0.000    0.000       0
    309316  0     83 00:22:03.00  10.5000A   1.8937  -2.0053   4.3023  0.114  1.41  -0.972E+00    850.  0.283E+00  0.283E+00  0.226E-03  0.320E-13   14   4   0   10   4   0   15  23  17    0.000    0.000       0
     
    312319  0     86 00:22:35.00  11.0000A   1.8894  -1.9584   4.0813  0.117  1.38  -0.945E+00    800.  0.283E+00  0.283E+00  0.237E-03  0.336E-13   14   4   0   10   4   0   16  23  17    0.000    0.000       0
    313320  0     87 00:22:46.00  11.4000A   1.8295  -1.8632   3.9326  0.118  1.38  -0.941E+00    800.  0.283E+00  0.283E+00  0.242E-03  0.344E-13   14   4   0   10   3   0   16  23  17    0.000    0.000       0
    314   0     88 00:22:57.40  12.0000A   1.8992  -1.8015   3.7404  0.119  1.38  -0.929E+00    800.  0.283E+00  0.283E+00  0.252E-03  0.357E-13   15   4   1    7  24   4   16  23  17    0.000    0.000       0
     321  0     88 00:22:57.40  12.0000A   1.8991  -1.8015   3.7404  0.119  1.38  -0.929E+00    800.  0.283E+00  0.283E+00  0.252E-03  0.357E-13   15   4   1    7  24   4   16  23  17    0.000    0.000       0
    315322  0     89 00:23:09.40  12.6000A   1.9542  -1.8235   3.5820  0.119  1.38  -0.923E+00    800.  0.283E+00  0.282E+00  0.265E-03  0.376E-13   15   4   1    8  23   4   17  23  17    0.000    0.000       0
    316323  0     90 00:23:22.00  12.9000A   1.9250  -1.9310   3.4822  0.120  1.38  -0.920E+00    800.  0.282E+00  0.282E+00  0.278E-03  0.395E-13   15   4   1    8  23   4   16  12  29    0.000    0.000       0
    317   0     91 00:23:34.90  13.2000A   1.8498  -2.0264   3.4174  0.121  1.38  -0.928E+00    800.  0.282E+00  0.281E+00  0.284E-03  0.404E-13   15   4   1    8  23   4   11  22  39    0.000    0.000       0
     324  0     91 00:23:34.90  13.2000A   1.8497  -2.0264   3.4174  0.121  1.38  -0.928E+00    800.  0.282E+00  0.281E+00  0.284E-03  0.404E-13   15   4   1    8  23   4   11  22  39    0.000    0.000       0
    318325  0     92 00:23:48.10  13.4000A   1.7617  -2.1071   3.3530  0.121  1.38  -0.909E+00    800.  0.281E+00  0.281E+00  0.290E-03  0.412E-13   15   4   1    8  23   4   12  22  39    0.000    0.000       0
    319   0     93 00:24:01.50  13.7000A   1.7141  -2.1684   3.2965  0.122  1.35  -0.903E+00    750.  0.280E+00  0.280E+00  0.293E-03  0.418E-13    3  22  38    8  23   4   12  22  39    0.000    0.000       0
    320   0     94 00:24:15.20  13.9000A   1.7172  -2.2027   3.2386  0.122  1.35  -0.901E+00    750.  0.280E+00  0.279E+00  0.298E-03  0.426E-13    3  22  38    8  23   4   12  22  39    0.000    0.000       0
    321   0     95 00:24:29.10  14.1000A   1.7197  -2.1962   3.1874  0.123  1.35  -0.896E+00    750.  0.279E+00  0.279E+00  0.300E-03  0.431E-13    3  22  38    8  23   4   12  22  39    0.000    0.000       0
     326  0     93 00:24:01.50  13.7000A   1.7141  -2.1684   3.2966  0.122  1.35  -0.903E+00    750.  0.280E+00  0.280E+00  0.293E-03  0.418E-13    3  22  38    8  23   4   12  22  39    0.000    0.000       0
     327  0     94 00:24:15.20  13.9000A   1.7172  -2.2027   3.2386  0.122  1.35  -0.901E+00    750.  0.280E+00  0.279E+00  0.297E-03  0.426E-13    3  22  38    8  23   4   12  22  39    0.000    0.000       0
     328  0     95 00:24:29.10  14.1000A   1.7197  -2.1962   3.1875  0.123  1.35  -0.896E+00    750.  0.279E+00  0.279E+00  0.300E-03  0.430E-13    3  22  38    8  23   4   12  22  39    0.000    0.000       0
    322329  0     96 00:24:43.20  14.1000A   1.7148  -2.1716   3.1978  0.123  1.35  -0.891E+00    750.  0.279E+00  0.278E+00  0.303E-03  0.435E-13    3  22  38    7  23   5   10  23   1    0.000    0.000       0
    323   0     97 00:24:57.30  13.8000A   1.6954  -2.1302   3.2532  0.123  1.35  -0.890E+00    750.  0.279E+00  0.278E+00  0.302E-03  0.432E-13    3  22  38    7  23   5   10  23   1    0.000    0.000       0
    324   0     98 00:25:11.10  13.8000A  -1.7060  -2.1049   3.2691  0.123  1.35  -0.891E+00    750.  0.278E+00  0.278E+00  0.295E-03  0.420E-13    8  25   1    8  23   5   11  23   1    0.000    0.000       0
    325   0     99 00:25:24.90  13.8000A   1.7235  -2.1802   3.2581  0.123  1.35  -0.892E+00    750.  0.278E+00  0.278E+00  0.295E-03  0.416E-13   10  26  31    8  23   5   11  23   1    0.000    0.000       0
    326   0    100 00:25:38.70  13.8000A   1.7646  -2.1973   3.2551  0.124  1.41  -0.893E+00    850.  0.278E+00  0.277E+00  0.294E-03  0.413E-13   10  26  31    8  23   5   12  23   1    0.000    0.000       0
    327   0    101 00:25:52.50  14.1000A   1.7850  -2.1490   3.2018  0.124  1.38  -0.899E+00    800.  0.278E+00  0.277E+00  0.293E-03  0.409E-13   10  26  31    8  23   5   12  23   1    0.000    0.000       0
    328   0    102 00:26:06.60  14.1000A   1.7815  -2.0922   3.1823  0.124  1.38  -0.897E+00    800.  0.278E+00  0.277E+00  0.299E-03  0.413E-13   10  26  31    8  22   5   13  23   1    0.000    0.000       0
    329   0    103 00:26:20.70  14.4000A  -1.8030  -2.1453   3.1344  0.124  1.38  -0.895E+00    800.  0.277E+00  0.277E+00  0.298E-03  0.409E-13   16  22  38    8  22   5   12  23  19    0.000    0.000       0
    330   0    104 00:26:35.10  14.4000A   1.7946  -2.1085   3.1234  0.124  1.38  -0.894E+00    800.  0.277E+00  0.276E+00  0.303E-03  0.413E-13    1  36  26    8  22   5   12  23  19    0.000    0.000       0
    331   0    105 00:26:49.50  14.4000A   1.8514  -1.9889   3.1209  0.124  1.38  -0.894E+00    800.  0.276E+00  0.276E+00  0.301E-03  0.410E-13    1  37  26    8  22   5    9  38  26    0.000    0.000       0
    332   0    106 00:27:03.90  14.1000A   1.8658   1.8646   3.1965  0.125  1.38  -0.898E+00    800.  0.276E+00  0.275E+00  0.299E-03  0.406E-13    1  37  26    8  35  24   10  38  26    0.000    0.000       0
    333   0    107 00:27:18.00  13.8000A   1.8273   1.8719   3.2544  0.125  1.38  -0.900E+00    800.  0.275E+00  0.274E+00  0.292E-03  0.396E-13    1  37  26    8  35  24   10  38  26    0.000    0.000       0
    334   0    108 00:27:31.80  13.8000A   1.8485   1.8998   3.2650  0.125  1.38  -0.899E+00    800.  0.274E+00  0.274E+00  0.285E-03  0.385E-13    1  16  38    8  35  23   10  38  26    0.000    0.000       0
    335   0    109 00:27:45.60  13.9000A   1.9161   1.9299   3.2455  0.125  1.38  -0.900E+00    800.  0.274E+00  0.273E+00  0.284E-03  0.384E-13    1  16  39    8  35  23   10  38  26    0.000    0.000       0
    336   0    110 00:27:59.50  13.9000A   1.9930   1.9488   3.2423  0.126  1.38  -0.902E+00    800.  0.273E+00  0.272E+00  0.285E-03  0.387E-13    1  16  39    8  36  23   11  38  26    0.000    0.000       0
    337   0    111 00:28:13.40  13.9000A   2.0631   2.0386   3.2310  0.126  1.38  -0.897E+00    800.  0.272E+00  0.271E+00  0.285E-03  0.386E-13    1  16  39    9  36  24   13  37  36    0.000    0.000       0
    338   0    112 00:28:27.30  14.0000A   2.1131   2.0896   3.2197  0.126  1.38  -0.896E+00    800.  0.271E+00  0.270E+00  0.285E-03  0.386E-13    1  16  39    9  36  24   13  37  36    0.000    0.000       0
    339   0    113 00:28:41.30  14.0000A   2.1239   2.1077   3.2062  0.126  1.41  -0.895E+00    850.  0.270E+00  0.269E+00  0.287E-03  0.388E-13    1  16  39    9  36  23   13  37  37    0.000    0.000       0
    340   0    114 00:28:55.30  14.0000A   2.1610   2.1093   3.2184  0.127  1.41  -0.893E+00    850.  0.269E+00  0.268E+00  0.286E-03  0.388E-13    1  16   0    9  36  23   13  37  37    0.000    0.000       0
    341   0    115 00:29:09.30  13.9000A   2.1914   2.1430   3.2457  0.127  1.41  -0.889E+00    850.  0.268E+00  0.267E+00  0.286E-03  0.387E-13    1  16   0    9  37  23   14  37  37    0.000    0.000       0
    342   0    116 00:29:23.20  13.9000A   2.1629   2.1318   3.2336  0.127  1.41  -0.888E+00    850.  0.266E+00  0.266E+00  0.284E-03  0.384E-13    1  16   0    9  37  23   14  37  37    0.000    0.000       0
    343   0    117 00:29:37.10  14.2000A   2.1088   2.0451   3.1754  0.127  1.38  -0.891E+00    800.  0.265E+00  0.264E+00  0.285E-03  0.382E-13    1  16   1    9  37  23   14  37  37    0.000    0.000       0
    344   0    118 00:29:51.30  14.3000A   2.1133   1.9030   3.1387  0.128  1.38  -0.904E+00    800.  0.263E+00  0.262E+00  0.291E-03  0.389E-13    1  16   1    9  37  23   14   9  10    0.000    0.000       0
    345   0    119 00:30:05.60  14.4000A   2.0602  -1.9071   3.1173  0.128  1.38  -0.888E+00    800.  0.261E+00  0.260E+00  0.293E-03  0.390E-13    1  16   1   10   6  21   13   6  16    0.000    0.000       0
    346   0    120 00:30:20.00  14.1000A   2.0168  -2.0168   3.1893  0.128  1.38  -0.885E+00    800.  0.259E+00  0.258E+00  0.294E-03  0.391E-13    1  17   1   10   6  21    7  25  19    0.000    0.000       0
    347   0    121 00:30:34.10  13.8000A  -2.0339  -2.0424   3.2554  0.128  1.38  -0.882E+00    800.  0.256E+00  0.256E+00  0.288E-03  0.381E-13    9  10  18   10   6  21    8  25  19    0.000    0.000       0
    348   0    122 00:30:47.90  13.4000A   1.9253  -2.1019   3.3500  0.129  1.41  -0.884E+00    850.  0.254E+00  0.253E+00  0.281E-03  0.371E-13    1  17   1   10   5  21    8  25  19    0.000    0.000       0
    349   0    123 00:31:01.30  13.2000A  -1.9151  -2.1226   3.3987  0.129  1.41  -0.879E+00    850.  0.252E+00  0.251E+00  0.272E-03  0.358E-13   10   8  19   10   5  21    9  25  19    0.000    0.000       0
    350   0    124 00:31:14.50  13.1000A  -1.9606  -2.0641   3.4475  0.129  1.41  -0.879E+00    850.  0.249E+00  0.249E+00  0.268E-03  0.350E-13   10   9  18   10   5  21    9  25  19    0.000    0.000       0
    351   0    125 00:31:27.60  13.1000A  -1.8503  -1.9343   3.4386  0.128  1.41  -0.885E+00    850.  0.247E+00  0.246E+00  0.265E-03  0.345E-13   10   9  18   10   5  21   10  25  19    0.000    0.000       0
    352   0    126 00:31:40.70  13.1000A  -1.9153  -1.8985   3.4354  0.128  1.41  -0.884E+00    850.  0.244E+00  0.244E+00  0.264E-03  0.343E-13   13  26  20    6  12  21   10  25  19    0.000    0.000       0
    353   0    127 00:31:53.80  13.4000A  -1.9319  -1.8848   3.3614  0.128  1.41  -0.890E+00    850.  0.242E+00  0.241E+00  0.263E-03  0.340E-13   13  26  20    6  12  21   10  25  19    0.000    0.000       0
    354   0    128 00:32:07.20  13.7000A  -1.9610  -1.8736   3.2883  0.128  1.41  -0.894E+00    850.  0.239E+00  0.239E+00  0.268E-03  0.346E-13   13  26  19    1  13  16    9   9  19    0.000    0.000       0
    355   0    129 00:32:20.90  13.9000A  -2.0531  -1.8801   3.2398  0.127  1.41  -0.899E+00    850.  0.237E+00  0.236E+00  0.272E-03  0.351E-13   13  26  19    1  13  16    9   9  19    0.000    0.000       0
    356   0    130 00:32:34.80  14.1000A  -2.0925  -1.8885   3.1975  0.127  1.41  -0.901E+00    850.  0.234E+00  0.233E+00  0.275E-03  0.353E-13   13  26  19    1  13  16   10   9  19    0.000    0.000       0
    357   0    131 00:32:48.90  14.2000A  -2.1181  -1.8974   3.1694  0.127  1.41  -0.903E+00    850.  0.232E+00  0.231E+00  0.277E-03  0.355E-13   13  26  19    1  13  16   10   9  19    0.000    0.000       0
    358   0    132 00:33:03.10  14.6000A  -2.1477  -1.9040   3.0733  0.126  1.41  -0.903E+00    850.  0.229E+00  0.228E+00  0.278E-03  0.354E-13   14  26  19    1  13  16   10   8  19    0.000    0.000       0
    359   0    133 00:33:17.70  14.7000A  -2.2353  -1.9052   3.0618  0.126  1.41  -0.908E+00    850.  0.227E+00  0.226E+00  0.285E-03  0.360E-13   14  26  18    1  13  16   10   8  19    0.000    0.000       0
    360   0    134 00:33:32.40  14.7000A  -2.3231  -1.8976   3.0628  0.125  1.41  -0.911E+00    850.  0.225E+00  0.224E+00  0.286E-03  0.359E-13   14  26  18    1  13  16   11   8  19    0.000    0.000       0
    361   0    135 00:33:47.10  14.6000A  -2.4176  -1.8784   3.0724  0.125  1.41  -0.913E+00    850.  0.223E+00  0.222E+00  0.284E-03  0.355E-13   14  26  18    1  12  16   11   8  19    0.000    0.000       0
    362   0    136 00:34:01.70  14.9000A  -2.5034  -1.8777   3.0298  0.124  1.38  -0.912E+00    800.  0.221E+00  0.220E+00  0.281E-03  0.349E-13   14  26  18    1  12  16   11   8  19    0.000    0.000       0
    363   0    137 00:34:16.60  15.2000A  -2.5428  -1.8641   2.9686  0.124  1.38  -0.911E+00    800.  0.219E+00  0.218E+00  0.286E-03  0.353E-13   14  26  18    1  12  16   11   8  19    0.000    0.000       0
    364   0    138 00:34:31.80  15.3000A  -2.4933  -1.8350   2.9454  0.123  1.38  -0.912E+00    800.  0.218E+00  0.217E+00  0.290E-03  0.357E-13   14  26  18    1  12  16   12   8  19    0.000    0.000       0
    365   0    139 00:34:47.10  15.4000A  -2.3596  -1.7912   2.9241  0.122  1.38  -0.917E+00    800.  0.216E+00  0.215E+00  0.292E-03  0.356E-13   14  27  18    1  12  16   12   8  19    0.000    0.000       0
    366   0    140 00:35:02.50  15.5000A  -2.3919  -1.7717   2.8943  0.122  1.38  -0.917E+00    800.  0.215E+00  0.214E+00  0.293E-03  0.355E-13   15  26  18    1  11  15   12   8  19    0.000    0.000       0
    367   0    141 00:35:18.00  15.8000A  -2.4704  -1.7378   2.8456  0.121  1.38  -0.922E+00    800.  0.214E+00  0.213E+00  0.294E-03  0.354E-13   15  26  17    1  11  15   12   8  19    0.000    0.000       0
    368   0    142 00:35:33.80  15.8000A  -2.4681   1.7976   2.8489  0.121  1.38  -0.929E+00    800.  0.213E+00  0.212E+00  0.299E-03  0.357E-13   15  26  17    8  34  23   14  26  24    0.000    0.000       0
    369   0    143 00:35:49.60  15.7000A  -2.3624  -1.7673   2.8671  0.120  1.38  -0.927E+00    800.  0.213E+00  0.212E+00  0.299E-03  0.354E-13   15  26  17    1   8  35    6  38  32    0.000    0.000       0
    370   0    144 00:36:05.30  15.2000A  -2.2471  -1.7490   2.9525  0.119  1.38  -0.930E+00    800.  0.212E+00  0.211E+00  0.297E-03  0.349E-13   15  26  16    1   7  34    6  38  32    0.000    0.000       0
    371   0    145 00:36:20.50  14.7000A  -2.0989  -1.8137   3.0686  0.119  1.38  -0.939E+00    800.  0.212E+00  0.211E+00  0.288E-03  0.337E-13   15  26  16    1   6  34    7  38  32    0.000    0.000       0
    372   0    146 00:36:35.20  14.2000A  -2.0337  -1.8432   3.1685  0.118  1.38  -0.941E+00    800.  0.212E+00  0.211E+00  0.278E-03  0.325E-13   11  11  13    1   6  34    7  38  32    0.000    0.000       0
    373   0    147 00:36:49.40  13.9000A  -2.0495  -1.8095   3.2329  0.118  1.38  -0.945E+00    800.  0.211E+00  0.211E+00  0.269E-03  0.313E-13   10  12  14    1   6  34    7  38  32    0.000    0.000       0
    374   0    148 00:37:03.30  13.6000A  -2.1261  -1.8874   3.2988  0.117  1.38  -0.949E+00    800.  0.211E+00  0.210E+00  0.264E-03  0.306E-13   10  12  14    1   5  34    8  38  32    0.000    0.000       0
    375   0    149 00:37:16.90  13.4000A  -2.1436  -1.8857   3.3518  0.117  1.38  -0.952E+00    800.  0.211E+00  0.210E+00  0.258E-03  0.299E-13   11  12  13    1   5  34    8  38  32    0.000    0.000       0
    376   0    150 00:37:30.30  13.4000A  -2.1473  -1.8654   3.3558  0.116  1.38  -0.957E+00    800.  0.211E+00  0.210E+00  0.255E-03  0.294E-13   11  12  13    1   6  35    8  38  32    0.000    0.000       0
    377   0    151 00:37:43.70  13.5000A  -2.1150  -1.8966   3.3332  0.116  1.38  -0.961E+00    800.  0.211E+00  0.210E+00  0.255E-03  0.293E-13   12  12  12    1   5  35    9  38  32    0.000    0.000       0
    378   0    152 00:37:57.20  13.6000A  -2.0518  -1.9679   3.3052  0.116  1.38  -0.964E+00    800.  0.211E+00  0.210E+00  0.257E-03  0.293E-13   12  12  12    1   5  35    9  38  32    0.000    0.000       0
    379   0    153 00:38:10.80  13.9000A  -1.9707  -1.9302   3.2273  0.115  1.41  -0.968E+00    850.  0.211E+00  0.210E+00  0.259E-03  0.294E-13   13  12  12    1   5  35    9  38  32    0.000    0.000       0
    380   0    154 00:38:24.70  14.2000A   1.9577  -2.0347   3.1608  0.115  1.41  -0.973E+00    850.  0.211E+00  0.210E+00  0.265E-03  0.298E-13    7  28  21    1   4  35    8  10   7    0.000    0.000       0
    381   0    155 00:38:38.90  13.4000A   1.9885  -2.0495   3.3484  0.115  1.41  -0.974E+00    850.  0.211E+00  0.210E+00  0.271E-03  0.303E-13    7  28  22    1   4  35    8   9   7    0.000    0.000       0
    382   0    156 00:38:52.30  12.8000A   2.0049  -1.9578   3.5085  0.114  1.41  -0.976E+00    850.  0.211E+00  0.210E+00  0.256E-03  0.286E-13    7  28  22    1   4  35    9   9   7    0.000    0.000       0
    383   0    157 00:39:05.10  12.6000A   1.9495  -2.0192   3.5708  0.114  1.41  -0.980E+00    850.  0.211E+00  0.210E+00  0.246E-03  0.273E-13    7  28  22    1   3  35    9   9   7    0.000    0.000       0
    384   0    158 00:39:17.70  12.5000A  -1.8994  -1.9976   3.5936  0.114  1.41  -0.982E+00    850.  0.211E+00  0.210E+00  0.243E-03  0.268E-13   14  13  11    1   3  35   10   9   7    0.000    0.000       0
    385   0    159 00:39:30.20  12.8000A  -1.8343  -1.8951   3.5263  0.114  1.41  -0.982E+00    850.  0.211E+00  0.210E+00  0.241E-03  0.266E-13   14  13  11    1   3  35   10   8   6    0.000    0.000       0
    386   0    160 00:39:43.00  13.0000A  -1.8423   1.8539   3.4660  0.114  1.41  -0.988E+00    850.  0.212E+00  0.211E+00  0.248E-03  0.272E-13    8  13  14   17  15  11   10   8   6    0.000    0.000       0
    387   0    161 00:39:56.00  13.1000A  -1.8734   1.9441   3.4252  0.113  1.41  -0.987E+00    850.  0.212E+00  0.211E+00  0.252E-03  0.276E-13    8  13  14   17  15  11   11   8   6    0.000    0.000       0
    388   0    162 00:40:09.10  13.1000A  -1.8857   1.9860   3.4418  0.113  1.41  -0.989E+00    850.  0.212E+00  0.211E+00  0.255E-03  0.278E-13    8  13  14   17  15  11   10   9   8    0.000    0.000       0
    389   0    163 00:40:22.20  12.9000A  -1.8946   1.9867   3.4953  0.113  1.41  -0.992E+00    850.  0.212E+00  0.211E+00  0.256E-03  0.279E-13    9  13  14   17  16  11   10   9   8    0.000    0.000       0
    390   0    164 00:40:35.10  13.0000A  -1.8988   2.1083   3.4484  0.113  1.41  -0.994E+00    850.  0.213E+00  0.211E+00  0.253E-03  0.275E-13    9  13  14   17  16  11   11   9   8    0.000    0.000       0
    391   0    165 00:40:48.10  13.3000A  -1.9469   2.1590   3.3863  0.113  1.41  -0.994E+00    850.  0.213E+00  0.212E+00  0.256E-03  0.278E-13   18  13   9   17  16  11   11   8   7    0.000    0.000       0
    392   0    166 00:41:01.40  13.5000A  -1.9768   2.1184   3.3307  0.113  1.41  -0.996E+00    850.  0.213E+00  0.212E+00  0.263E-03  0.285E-13   18  13   9   17  16  11   12   8   7    0.000    0.000       0
    393   0    167 00:41:14.90  13.7000A  -1.9837   2.1241   3.2823  0.113  1.41  -0.100E+01    850.  0.213E+00  0.212E+00  0.268E-03  0.290E-13   18  14   9   17  17  11   11   8   8    0.000    0.000       0
    394   0    168 00:41:28.60  13.6000A  -1.9913   2.1833   3.3064  0.113  1.41  -0.100E+01    850.  0.214E+00  0.212E+00  0.274E-03  0.295E-13   18  14   9   17  17  11   12   8   8    0.000    0.000       0
    395   0    169 00:41:42.20  13.4000A  -1.9813   2.1737   3.3472  0.113  1.41  -0.101E+01    850.  0.214E+00  0.213E+00  0.273E-03  0.294E-13   18  14   9   18  18   9   12   8   8    0.000    0.000       0
    396   0    170 00:41:55.60  13.5000A  -1.9944   2.2409   3.3230  0.113  1.38  -0.101E+01    800.  0.214E+00  0.213E+00  0.270E-03  0.291E-13   18  14   8   17  18  10   12   8   8    0.000    0.000       0
    397   0    171 00:42:09.10  13.5000A  -2.0090   2.2335   3.3380  0.113  1.38  -0.101E+01    800.  0.215E+00  0.213E+00  0.272E-03  0.293E-13   18  14   8   17  18  10   13   8   8    0.000    0.000       0
    398   0    172 00:42:22.60  13.6000A  -1.9963   2.1836   3.3089  0.113  1.38  -0.101E+01    800.  0.215E+00  0.214E+00  0.273E-03  0.294E-13   18  14   7   18  19   8   13   8   8    0.000    0.000       0
    399   0    173 00:42:36.20  13.9000A  -1.9847   2.2215   3.2261  0.113  1.38  -0.101E+01    800.  0.216E+00  0.214E+00  0.276E-03  0.297E-13   18  14   7   18  19   8   13   8   8    0.000    0.000       0
    400   0    174 00:42:50.10  14.1000A  -2.0020   2.1920   3.1998  0.113  1.38  -0.102E+01    800.  0.217E+00  0.215E+00  0.282E-03  0.304E-13   18  15   7   17  19  10   14   8   8    0.000    0.000       0
    401   0    175 00:43:04.20  14.3000A  -1.9699   2.2250   3.1385  0.113  1.38  -0.101E+01    800.  0.218E+00  0.216E+00  0.286E-03  0.308E-13   18  15   7   17  20   9   14   8   8    0.000    0.000       0
    402   0    176 00:43:18.50  14.8000A  -1.9215   2.2241   3.0482  0.114  1.38  -0.100E+01    800.  0.219E+00  0.217E+00  0.290E-03  0.313E-13   17  14   4   17  20   9   14   8   8    0.000    0.000       0
    403   0    177 00:43:33.30  15.2000A  -1.9580   2.1654   2.9561  0.114  1.38  -0.100E+01    800.  0.220E+00  0.218E+00  0.300E-03  0.324E-13   17  14   4   17  20   9   14   8   8    0.000    0.000       0
    404   0    178 00:43:48.50  15.5000A  -1.9076   2.1079   2.9114  0.114  1.38  -0.994E+00    800.  0.221E+00  0.219E+00  0.308E-03  0.332E-13   17  14   4   18  18  12   10   9  18    0.000    0.000       0
    405   0    179 00:44:04.00  15.0000A  -1.8769   2.1355   2.9998  0.114  1.38  -0.100E+01    800.  0.222E+00  0.221E+00  0.314E-03  0.339E-13   17  15   4   18  19  11   10   9  18    0.000    0.000       0
    406   0    180 00:44:19.00  14.7000A  -1.8140   2.1157   3.0702  0.115  1.38  -0.992E+00    800.  0.224E+00  0.222E+00  0.305E-03  0.328E-13   16  15   3   18  17  13   11   9  18    0.000    0.000       0
    407   0    181 00:44:33.70  14.6000A  -1.7678   2.1454   3.0915  0.115  1.43  -0.991E+00    900.  0.225E+00  0.223E+00  0.300E-03  0.323E-13   16  15   3   18  17  13   13  11  16    0.000    0.000       0
    408   0    182 00:44:48.30  14.2000A   1.7560   2.1502   3.1645  0.115  1.43  -0.995E+00    900.  0.226E+00  0.225E+00  0.299E-03  0.321E-13    9  27  23   18  17  13   13  11  16    0.000    0.000       0
    409   0    183 00:45:02.50  13.9000A   1.7854   2.1719   3.2326  0.116  1.43  -0.990E+00    900.  0.228E+00  0.226E+00  0.292E-03  0.314E-13    9  27  23   18  18  13   14  11  16    0.000    0.000       0
    410   0    184 00:45:16.40  13.6000A   1.7977   2.1976   3.3144  0.116  1.43  -0.992E+00    900.  0.230E+00  0.228E+00  0.287E-03  0.308E-13    9  27  23   18  18  13   14  11  16    0.000    0.000       0
    411   0    185 00:45:30.00  13.3000A   1.8109   2.2067   3.3797  0.116  1.43  -0.100E+01    900.  0.231E+00  0.230E+00  0.281E-03  0.302E-13    9  27  24   18  18  13   14  11  16    0.000    0.000       0
    412   0    186 00:45:43.30  13.2000A   1.8581   2.2022   3.4134  0.117  1.43  -0.991E+00    900.  0.233E+00  0.231E+00  0.276E-03  0.297E-13    9  27  24   18  18  13   14  11  16    0.000    0.000       0
    413   0    187 00:45:56.50  13.1000A   1.8920   2.1844   3.4439  0.117  1.43  -0.987E+00    900.  0.234E+00  0.233E+00  0.275E-03  0.295E-13    9  27  24   18  18  13   15  11  16    0.000    0.000       0
    414   0    188 00:46:09.60  12.9000A   1.9095   2.1525   3.4749  0.117  1.43  -0.985E+00    900.  0.236E+00  0.235E+00  0.274E-03  0.294E-13    9  27  24   18  18  13   15  11  16    0.000    0.000       0
    415   0    189 00:46:22.50  13.0000A   1.9088   2.1037   3.4713  0.118  1.43  -0.984E+00    900.  0.238E+00  0.237E+00  0.271E-03  0.291E-13    9  27  24   18  18  13   15  11  16    0.000    0.000       0
    416   0    190 00:46:35.50  13.0000A   1.8892   2.0727   3.4501  0.118  1.43  -0.984E+00    900.  0.240E+00  0.238E+00  0.274E-03  0.294E-13    9  27  24   18  19  12   16  11  16    0.000    0.000       0
    417   0    191 00:46:48.50  13.0000A   1.8879   2.0544   3.4586  0.118  1.43  -0.986E+00    900.  0.242E+00  0.240E+00  0.275E-03  0.295E-13    9  26  25   18  19  12   16  11  16    0.000    0.000       0
    418   0    192 00:47:01.50  13.1000A   1.8897   2.0539   3.4431  0.119  1.41  -0.981E+00    850.  0.243E+00  0.242E+00  0.276E-03  0.296E-13    9  27  25    8   4   0   16  11  16    0.000    0.000       0
    419   0    193 00:47:14.60  13.2000A   1.9011   2.0297   3.4101  0.119  1.41  -0.978E+00    850.  0.245E+00  0.244E+00  0.279E-03  0.300E-13    9  27  25    8   4   0   13  10  17    0.000    0.000       0
    420   0    194 00:47:27.80  13.3000A   1.8908   2.0100   3.3868  0.119  1.41  -0.974E+00    850.  0.247E+00  0.246E+00  0.283E-03  0.305E-13    9  27  25   19  18  13   14  10  17    0.000    0.000       0
    421   0    195 00:47:41.10  13.2000A   1.8585   2.0275   3.4008  0.119  1.41  -0.970E+00    850.  0.249E+00  0.248E+00  0.286E-03  0.309E-13    9  27  25   19  18  14   14  10  17    0.000    0.000       0
    422   0    196 00:47:54.30  13.2000A  -1.8883  -1.9919   3.4148  0.119  1.41  -0.968E+00    850.  0.251E+00  0.250E+00  0.286E-03  0.310E-13    6   3   1   13  38  15   11   9  18    0.000    0.000       0
    423   0    197 00:48:07.50  12.9000A  -2.0109  -1.9287   3.4810  0.119  1.41  -0.973E+00    850.  0.253E+00  0.252E+00  0.288E-03  0.313E-13    6   3   1   13  38  15   12   9  18    0.000    0.000       0
    424   0    198 00:48:20.40  12.8000A  -2.0676   1.9262   3.5058  0.119  1.41  -0.971E+00    850.  0.255E+00  0.253E+00  0.283E-03  0.309E-13    6   3   1   19  19  13   12   9  18    0.000    0.000       0
    425   0    199 00:48:33.20  12.8000A  -2.0940   1.8888   3.5269  0.119  1.41  -0.972E+00    850.  0.257E+00  0.255E+00  0.283E-03  0.309E-13    6   4   0   19  19  13   13   9  18    0.000    0.000       0
    426   0    200 00:48:46.00  12.7000A   2.1500   1.8285   3.5411  0.119  1.41  -0.975E+00    850.  0.259E+00  0.257E+00  0.285E-03  0.312E-13    7   2  10   12   4  39   13   9  18    0.000    0.000       0
    427   0    201 00:48:58.70  12.7000A   2.1955   1.8433   3.5368  0.119  1.41  -0.974E+00    850.  0.261E+00  0.259E+00  0.284E-03  0.313E-13    7   2  10   12   4  39   14   9  18    0.000    0.000       0
    428   0    202 00:49:11.40  12.5000A  -2.2761  -1.8582   3.5983  0.119  1.41  -0.976E+00    850.  0.262E+00  0.261E+00  0.286E-03  0.316E-13    7   3   0    3   8   2    9  39   1    0.000    0.000       0
    429   0    203 00:49:23.90  12.9000A  -2.3204  -1.9224   3.4935  0.119  1.41  -0.973E+00    850.  0.264E+00  0.262E+00  0.284E-03  0.315E-13    7   3   0    3   8   2   15   9  18    0.000    0.000       0
    430   0    204 00:49:36.80  12.8000A  -2.2110  -1.9685   3.5230  0.119  1.41  -0.974E+00    850.  0.266E+00  0.264E+00  0.295E-03  0.327E-13    8   2   0    3   8   2    9   0   0    0.000    0.000       0
    431   0    205 00:49:49.60  12.6000A  -2.3818  -2.0373   3.5648  0.119  1.41  -0.975E+00    850.  0.268E+00  0.266E+00  0.294E-03  0.328E-13    8   1   0    3   7   2    9   7  17    0.000    0.000       0
    432   0    206 00:50:02.20  12.4000A  -2.4161  -2.1012   3.6437  0.118  1.41  -0.973E+00    850.  0.269E+00  0.267E+00  0.291E-03  0.325E-13    8   1   0    3   7   2    9   7  17    0.000    0.000       0
    433   0    207 00:50:14.60  12.2000A  -2.3757  -2.1220   3.6851  0.118  1.41  -0.972E+00    850.  0.271E+00  0.269E+00  0.287E-03  0.322E-13    8   1   0    3   7   2   10   7  17    0.000    0.000       0
    434   0    208 00:50:26.80  12.0000A  -2.4802  -2.0799   3.7439  0.118  1.43  -0.972E+00    900.  0.272E+00  0.270E+00  0.284E-03  0.319E-13    9   1  39    3   7   2   10   7  17    0.000    0.000       0
    435   0    209 00:50:38.80  12.0000A  -2.5304  -2.0332   3.7519  0.118  1.43  -0.970E+00    900.  0.273E+00  0.271E+00  0.280E-03  0.315E-13    9   1  39    4   7   3   10   7  17    0.000    0.000       0
    436   0    210 00:50:50.80  12.1000A  -2.4939  -2.0329   3.7144  0.118  1.43  -0.968E+00    900.  0.275E+00  0.272E+00  0.281E-03  0.317E-13    9   1  39    3   6   2   10   7  17    0.000    0.000       0
    437   0    211 00:51:02.90  12.0000A  -2.4209  -1.9740   3.7581  0.118  1.43  -0.968E+00    900.  0.276E+00  0.274E+00  0.284E-03  0.321E-13    9   1  39    3   6   2   14   6  18    0.000    0.000       0
    438   0    212 00:51:14.90  11.7000A  -2.3846  -2.0341   3.8348  0.118  1.43  -0.970E+00    900.  0.277E+00  0.275E+00  0.283E-03  0.320E-13    9   1  38   15  33  38   14   6  18    0.000    0.000       0
    439   0    213 00:51:26.59  11.7000A  -2.3506  -2.1307   3.8576  0.117  1.43  -0.975E+00    900.  0.278E+00  0.276E+00  0.277E-03  0.314E-13   10   1  38   15  33  38   14   6  18    0.000    0.000       0
    440   0    214 00:51:38.29  11.8000A  -2.3205  -2.1303   3.8242  0.117  1.43  -0.977E+00    900.  0.279E+00  0.277E+00  0.277E-03  0.315E-13   10   1  38   15  33  38   14   6  18    0.000    0.000       0
    441   0    215 00:51:50.09  12.0000A  -2.2789  -2.0978   3.7524  0.117  1.43  -0.981E+00    900.  0.280E+00  0.277E+00  0.281E-03  0.319E-13   11   0  38   15  32  38   15   5  18    0.000    0.000       0
    442   0    216 00:52:02.09  12.0000A  -2.2487  -2.0444   3.7616  0.117  1.43  -0.984E+00    900.  0.281E+00  0.278E+00  0.287E-03  0.326E-13   11   0  38   15  32  38   16   5  18    0.000    0.000       0
    443   0    217 00:52:14.09  12.0000A  -2.2217  -1.9955   3.7539  0.117  1.43  -0.980E+00    900.  0.281E+00  0.279E+00  0.288E-03  0.327E-13   11   0  37   15  32  37   16   5  18    0.000    0.000       0
    444   0    218 00:52:26.09  12.2000A  -2.1807  -2.0609   3.6997  0.117  1.43  -0.978E+00    900.  0.282E+00  0.279E+00  0.289E-03  0.329E-13   11   0  37   15  31  37   16   5  18    0.000    0.000       0
    445   0    219 00:52:38.29  12.6000A  -2.0993  -2.1447   3.5845  0.117  1.43  -0.978E+00    900.  0.282E+00  0.280E+00  0.294E-03  0.335E-13   11   0  37   15  31  37   16   5  18    0.000    0.000       0
    446   0    220 00:52:50.89  13.2000A  -1.9941  -2.1586   3.4126  0.116  1.43  -0.979E+00    900.  0.283E+00  0.280E+00  0.304E-03  0.346E-13   13  23   4   15  31  37   16   5  18    0.000    0.000       0
    447   0    221 00:53:04.09  13.6000A  -1.9882  -2.0853   3.3005  0.116  1.43  -0.983E+00    900.  0.283E+00  0.280E+00  0.318E-03  0.362E-13   13  23   4   15  31  37   17   4  18    0.000    0.000       0
    448   0    222 00:53:17.69  13.6000A  -1.9653  -2.0588   3.2981  0.116  1.43  -0.993E+00    900.  0.283E+00  0.280E+00  0.328E-03  0.372E-13   13  23   4   15  30  37    8   4   4    0.000    0.000       0
    449   0    223 00:53:31.29  13.6000A  -1.9332   2.0510   3.3020  0.116  1.43  -0.994E+00    900.  0.283E+00  0.280E+00  0.328E-03  0.370E-13   14  24   3   17  17  10    8  35   5    0.000    0.000       0
    450   0    224 00:53:44.89  13.2000A  -1.9953   2.0070   3.4031  0.116  1.43  -0.987E+00    900.  0.283E+00  0.280E+00  0.327E-03  0.369E-13    9   6   4   17  17  10    9  35   5    0.000    0.000       0
    451   0    225 00:53:58.09  13.1000A  -2.0501   1.9133   3.4417  0.116  1.43  -0.989E+00    900.  0.282E+00  0.280E+00  0.318E-03  0.357E-13    9   6   4   17  17  10    9  35   5    0.000    0.000       0
    452   0    226 00:54:11.19  13.1000A  -2.0938   1.8867   3.4369  0.115  1.43  -0.993E+00    900.  0.282E+00  0.279E+00  0.315E-03  0.352E-13   10   6   4   13   8   4   10  35   5    0.000    0.000       0
    453   0    227 00:54:24.29  13.0000A  -2.1923   1.9127   3.4680  0.115  1.43  -0.100E+01    900.  0.282E+00  0.279E+00  0.315E-03  0.351E-13   10   6   4   13   8   4   10  35   5    0.000    0.000       0
    454   0    228 00:54:37.29  13.1000A  -2.2403   1.9712   3.4314  0.115  1.40  -0.995E+00    850.  0.282E+00  0.279E+00  0.312E-03  0.347E-13   10   6   4   14   9   5   10  35   5    0.000    0.000       0
    455   0    229 00:54:50.39  12.8000A  -2.2360   2.0200   3.5122  0.115  1.40  -0.994E+00    850.  0.282E+00  0.279E+00  0.315E-03  0.349E-13   10   6   4   14   9   5   10  36   6    0.000    0.000       0
    456   0    230 00:55:03.19  12.5000A  -2.2380   2.0059   3.5910  0.115  1.40  -0.996E+00    850.  0.282E+00  0.279E+00  0.308E-03  0.340E-13   11   6   4   14   9   5   10  36   6    0.000    0.000       0
    457   0    231 00:55:15.69  12.4000A  -2.2664   1.9510   3.6299  0.115  1.40  -0.990E+00    850.  0.281E+00  0.279E+00  0.301E-03  0.332E-13   11   6   4   14   9   5   10  36   6    0.000    0.000       0
    458   0    232 00:55:28.09  12.3000A  -2.2644   1.9653   3.6637  0.115  1.40  -0.991E+00    850.  0.281E+00  0.279E+00  0.300E-03  0.329E-13   11   6   4   17  10   6   11  36   6    0.000    0.000       0
    459   0    233 00:55:40.39  12.1000A  -2.2635   1.9614   3.7037  0.115  1.40  -0.992E+00    850.  0.281E+00  0.279E+00  0.298E-03  0.326E-13   11   6   3   17  10   6   11  36   6    0.000    0.000       0
    460   0    234 00:55:52.49  12.1000A  -2.2125   1.9119   3.7274  0.115  1.40  -0.989E+00    850.  0.281E+00  0.278E+00  0.294E-03  0.321E-13   11   6   3   12   3  37   11  36   6    0.000    0.000       0
    461   0    235 00:56:04.59  12.1000A  -2.1711   1.8977   3.7116  0.115  1.40  -0.987E+00    850.  0.281E+00  0.278E+00  0.295E-03  0.321E-13   12   6   3   13   2  38   11  36   6    0.000    0.000       0
    462   0    236 00:56:16.69  12.3000A  -2.1905   1.9584   3.6493  0.115  1.40  -0.985E+00    850.  0.281E+00  0.278E+00  0.295E-03  0.321E-13   12   6   3   13   3  36   11  36   6    0.000    0.000       0
    463   0    237 00:56:28.99  12.3000A  -2.1734   1.9387   3.6577  0.115  1.40  -0.985E+00    850.  0.280E+00  0.278E+00  0.300E-03  0.327E-13   12   6   3   13   3  36   12  36   7    0.000    0.000       0
    464   0    238 00:56:41.29  12.4000A  -2.1181   1.9538   3.6195  0.115  1.40  -0.987E+00    850.  0.280E+00  0.278E+00  0.301E-03  0.327E-13   12   6   3   13   4  36   12  36   7    0.000    0.000       0
    465   0    239 00:56:53.69  12.5000A  -2.0250   1.9152   3.5962  0.115  1.38  -0.988E+00    800.  0.280E+00  0.277E+00  0.304E-03  0.330E-13   12   6   3   13   4  37   13  36   7    0.000    0.000       0
    466   0    240 00:57:06.19  12.7000A  -1.9906   1.9357   3.5442  0.115  1.38  -0.988E+00    800.  0.280E+00  0.277E+00  0.307E-03  0.334E-13   13   6   3   14   4  37   13  36   7    0.000    0.000       0
    467   0    241 00:57:18.89  12.8000A   2.0433   1.9496   3.5178  0.115  1.38  -0.990E+00    800.  0.279E+00  0.277E+00  0.313E-03  0.340E-13   12  37   9   14   5  37   14  35   6    0.000    0.000       0
    468   0    242 00:57:31.69  12.6000A   2.1147   1.9760   3.5706  0.115  1.40  -0.993E+00    850.  0.279E+00  0.277E+00  0.316E-03  0.343E-13   13  37   9   14   5  37   14  35   6    0.000    0.000       0
    469   0    243 00:57:44.29  12.6000A   2.1795   1.9496   3.5723  0.115  1.40  -0.998E+00    850.  0.279E+00  0.276E+00  0.312E-03  0.339E-13   13  37   9   14   5  37   14  35   6    0.000    0.000       0
    470   0    244 00:57:56.89  12.7000A   2.1966   1.9318   3.5502  0.114  1.40  -0.996E+00    850.  0.279E+00  0.276E+00  0.312E-03  0.340E-13   13  37   9   15   5  37   15  35   6    0.000    0.000       0
    471   0    245 00:58:09.59  12.7000A   2.1739   1.9579   3.5513  0.114  1.40  -0.996E+00    850.  0.279E+00  0.276E+00  0.315E-03  0.343E-13   13  37  10   15   5  37   15  35   6    0.000    0.000       0
    472   0    246 00:58:22.29  12.9000A   2.2250   1.9814   3.4946  0.114  1.38  -0.998E+00    800.  0.279E+00  0.276E+00  0.314E-03  0.344E-13   14  37  10   15   6  37   15  35   6    0.000    0.000       0
    473   0    247 00:58:35.19  13.3000A   2.2578   1.9773   3.3900  0.114  1.38  -0.997E+00    800.  0.279E+00  0.276E+00  0.319E-03  0.350E-13   14  37  10   15   6  37   15  35   6    0.000    0.000       0
    474   0    248 00:58:48.49  13.6000A   2.2207   1.9381   3.3184  0.114  1.38  -0.998E+00    800.  0.279E+00  0.276E+00  0.329E-03  0.361E-13   14  37  10   15   6  37   15  34   3    0.000    0.000       0
    475   0    249 00:59:02.09  13.6000A   2.2069   1.9631   3.3056  0.114  1.38  -0.100E+01    800.  0.278E+00  0.276E+00  0.336E-03  0.369E-13   14  37  11   16   6  37   12  37   6    0.000    0.000       0
    476   0    250 00:59:15.69  13.5000A   2.1185   1.9916   3.3423  0.114  1.38  -0.100E+01    800.  0.278E+00  0.276E+00  0.336E-03  0.369E-13   14  37  11   16   7  37    9  37   5    0.000    0.000       0
    477   0    251 00:59:29.19  13.4000A   2.0242   2.0192   3.3527  0.113  1.43  -0.101E+01    900.  0.278E+00  0.276E+00  0.333E-03  0.365E-13   16  36  11   16   7  37   10  37   5    0.000    0.000       0
    478   0    252 00:59:42.59  13.2000A   1.9364   2.0137   3.4014  0.113  1.43  -0.101E+01    900.  0.278E+00  0.276E+00  0.331E-03  0.362E-13   15  37  12   16   7  37   10  37   5    0.000    0.000       0
    479   0    253 00:59:55.79  13.3000A   1.9199   2.0227   3.3908  0.113  1.43  -0.102E+01    900.  0.278E+00  0.276E+00  0.325E-03  0.355E-13   16  36  12   16   7  36   10  37   5    0.000    0.000       0
    480   0    254 01:00:09.09  13.2000A   1.9010   2.0017   3.4011  0.113  1.43  -0.102E+01    900.  0.278E+00  0.276E+00  0.328E-03  0.356E-13   16  36  12   16   7  36   11  37   5    0.000    0.000       0
     330  0     97 00:24:57.30  13.8000A   1.6954  -2.1302   3.2533  0.123  1.35  -0.890E+00    750.  0.279E+00  0.278E+00  0.302E-03  0.432E-13    3  22  38    7  23   5   10  23   1    0.000    0.000       0
     331  0     98 00:25:11.10  13.8000A  -1.7061  -2.1049   3.2692  0.123  1.35  -0.891E+00    750.  0.278E+00  0.278E+00  0.295E-03  0.420E-13    8  25   1    8  23   5   11  23   1    0.000    0.000       0
     332  0     99 00:25:24.90  13.8000A   1.7235  -2.1802   3.2581  0.123  1.35  -0.892E+00    750.  0.278E+00  0.278E+00  0.294E-03  0.416E-13   10  26  31    8  23   5   11  23   1    0.000    0.000       0
     333  0    100 00:25:38.70  13.8000A   1.7646  -2.1973   3.2552  0.124  1.41  -0.893E+00    850.  0.278E+00  0.277E+00  0.294E-03  0.412E-13   10  26  31    8  23   5   12  23   1    0.000    0.000       0
     334  0    101 00:25:52.50  14.1000A   1.7849  -2.1490   3.2018  0.124  1.38  -0.899E+00    800.  0.278E+00  0.277E+00  0.293E-03  0.408E-13   10  26  31    8  23   5   12  23   1    0.000    0.000       0
     335  0    102 00:26:06.60  14.1000A   1.7815  -2.0923   3.1823  0.124  1.38  -0.897E+00    800.  0.278E+00  0.277E+00  0.299E-03  0.413E-13   10  26  31    8  22   5   13  23   1    0.000    0.000       0
     336  0    103 00:26:20.70  14.4000A  -1.8029  -2.1453   3.1344  0.124  1.38  -0.895E+00    800.  0.277E+00  0.277E+00  0.298E-03  0.409E-13   16  22  38    8  22   5   12  23  19    0.000    0.000       0
     337  0    104 00:26:35.10  14.4000A   1.7946  -2.1085   3.1234  0.124  1.38  -0.894E+00    800.  0.277E+00  0.276E+00  0.302E-03  0.413E-13    1  36  26    8  22   5   12  23  19    0.000    0.000       0
     338  0    105 00:26:49.50  14.4000A   1.8514  -1.9889   3.1209  0.124  1.38  -0.894E+00    800.  0.276E+00  0.275E+00  0.300E-03  0.409E-13    1  37  26    8  22   5    9  38  26    0.000    0.000       0
     339  0    106 00:27:03.90  14.1000A   1.8658   1.8647   3.1965  0.125  1.38  -0.898E+00    800.  0.276E+00  0.275E+00  0.299E-03  0.406E-13    1  37  26    8  35  24   10  38  26    0.000    0.000       0
     340  0    107 00:27:18.00  13.8000A   1.8273   1.8719   3.2545  0.125  1.38  -0.900E+00    800.  0.275E+00  0.274E+00  0.291E-03  0.395E-13    1  37  26    8  35  24   10  38  26    0.000    0.000       0
     341  0    108 00:27:31.80  13.8000A   1.8485   1.8999   3.2651  0.125  1.38  -0.899E+00    800.  0.274E+00  0.274E+00  0.284E-03  0.385E-13    1  16  38    8  35  23   10  38  26    0.000    0.000       0
     342  0    109 00:27:45.60  13.9000A   1.9161   1.9300   3.2455  0.125  1.38  -0.900E+00    800.  0.274E+00  0.273E+00  0.283E-03  0.384E-13    1  16  39    8  35  23   10  38  26    0.000    0.000       0
     343  0    110 00:27:59.50  13.9000A   1.9930   1.9489   3.2424  0.126  1.38  -0.902E+00    800.  0.273E+00  0.272E+00  0.285E-03  0.386E-13    1  16  39    8  36  23   11  38  26    0.000    0.000       0
     344  0    111 00:28:13.40  13.9000A   2.0631   2.0386   3.2311  0.126  1.38  -0.897E+00    800.  0.272E+00  0.271E+00  0.284E-03  0.386E-13    1  16  39    9  36  24   13  37  36    0.000    0.000       0
     345  0    112 00:28:27.30  14.0000A   2.1131   2.0896   3.2197  0.126  1.38  -0.896E+00    800.  0.271E+00  0.270E+00  0.284E-03  0.385E-13    1  16  39    9  36  24   13  37  36    0.000    0.000       0
     346  0    113 00:28:41.30  14.0000A   2.1239   2.1078   3.2063  0.126  1.41  -0.895E+00    850.  0.270E+00  0.269E+00  0.286E-03  0.388E-13    1  16  39    9  36  23   13  37  37    0.000    0.000       0
     347  0    114 00:28:55.30  14.0000A   2.1610   2.1093   3.2184  0.127  1.41  -0.893E+00    850.  0.269E+00  0.268E+00  0.286E-03  0.387E-13    1  16   0    9  36  23   13  37  37    0.000    0.000       0
     348  0    115 00:29:09.30  13.9000A   2.1914   2.1431   3.2457  0.127  1.41  -0.889E+00    850.  0.268E+00  0.267E+00  0.286E-03  0.386E-13    1  16   0    9  37  23   14  37  37    0.000    0.000       0
     349  0    116 00:29:23.20  13.9000A   2.1628   2.1319   3.2337  0.127  1.41  -0.888E+00    850.  0.266E+00  0.265E+00  0.284E-03  0.383E-13    1  16   0    9  37  23   14  37  37    0.000    0.000       0
     350  0    117 00:29:37.10  14.2000A   2.1087   2.0452   3.1755  0.127  1.38  -0.891E+00    800.  0.265E+00  0.264E+00  0.284E-03  0.382E-13    1  16   1    9  37  23   14  37  37    0.000    0.000       0
     351  0    118 00:29:51.30  14.3000A   2.1133   1.9031   3.1386  0.128  1.38  -0.904E+00    800.  0.263E+00  0.262E+00  0.290E-03  0.388E-13    1  16   1    9  37  23   14   9  10    0.000    0.000       0
     352  0    119 00:30:05.60  14.4000A   2.0602  -1.9072   3.1172  0.128  1.38  -0.888E+00    800.  0.261E+00  0.260E+00  0.292E-03  0.390E-13    1  16   1   10   6  21   13   6  16    0.000    0.000       0
     353  0    120 00:30:20.00  14.1000A   2.0168  -2.0169   3.1893  0.128  1.38  -0.885E+00    800.  0.259E+00  0.258E+00  0.293E-03  0.390E-13    1  17   1   10   6  21    7  25  19    0.000    0.000       0
     354  0    121 00:30:34.10  13.8000A  -2.0339  -2.0424   3.2555  0.128  1.38  -0.882E+00    800.  0.256E+00  0.256E+00  0.287E-03  0.380E-13    9  10  18   10   6  21    8  25  19    0.000    0.000       0
     355  0    122 00:30:47.90  13.4000A   1.9253  -2.1020   3.3501  0.129  1.41  -0.884E+00    850.  0.254E+00  0.253E+00  0.280E-03  0.370E-13    1  17   1   10   5  21    8  25  19    0.000    0.000       0
     356  0    123 00:31:01.30  13.2000A  -1.9150  -2.1226   3.3988  0.129  1.41  -0.879E+00    850.  0.252E+00  0.251E+00  0.272E-03  0.357E-13   10   8  19   10   5  21    9  25  19    0.000    0.000       0
     357  0    124 00:31:14.50  13.1000A  -1.9605  -2.0641   3.4475  0.129  1.41  -0.879E+00    850.  0.249E+00  0.249E+00  0.267E-03  0.349E-13   10   9  18   10   5  21    9  25  19    0.000    0.000       0
     358  0    125 00:31:27.60  13.1000A  -1.8503  -1.9343   3.4387  0.128  1.41  -0.885E+00    850.  0.247E+00  0.246E+00  0.264E-03  0.344E-13   10   9  18   10   5  21   10  25  19    0.000    0.000       0
     359  0    126 00:31:40.70  13.1000A  -1.9154  -1.8984   3.4355  0.128  1.41  -0.884E+00    850.  0.244E+00  0.244E+00  0.263E-03  0.342E-13   13  26  20    6  12  21   10  25  19    0.000    0.000       0
     360  0    127 00:31:53.80  13.4000A  -1.9320  -1.8848   3.3615  0.128  1.41  -0.890E+00    850.  0.242E+00  0.241E+00  0.262E-03  0.339E-13   13  26  20    6  12  21   10  25  19    0.000    0.000       0
     361  0    128 00:32:07.20  13.7000A  -1.9612  -1.8735   3.2882  0.128  1.41  -0.894E+00    850.  0.239E+00  0.238E+00  0.267E-03  0.344E-13   13  26  19    1  13  16    9   9  19    0.000    0.000       0
     362  0    129 00:32:20.90  13.9000A  -2.0533  -1.8800   3.2398  0.127  1.41  -0.899E+00    850.  0.237E+00  0.236E+00  0.271E-03  0.349E-13   13  26  19    1  13  16    9   9  19    0.000    0.000       0
     363  0    130 00:32:34.80  14.1000A  -2.0927  -1.8885   3.1974  0.127  1.41  -0.901E+00    850.  0.234E+00  0.233E+00  0.274E-03  0.352E-13   13  26  19    1  13  16   10   9  19    0.000    0.000       0
     364  0    131 00:32:48.90  14.2000A  -2.1182  -1.8973   3.1694  0.127  1.41  -0.903E+00    850.  0.231E+00  0.231E+00  0.276E-03  0.354E-13   13  26  19    1  13  16   10   9  19    0.000    0.000       0
     365  0    132 00:33:03.10  14.6000A  -2.1477  -1.9039   3.0732  0.126  1.41  -0.903E+00    850.  0.229E+00  0.228E+00  0.277E-03  0.353E-13   14  26  19    1  13  16   10   8  19    0.000    0.000       0
     366  0    133 00:33:17.70  14.7000A  -2.2354  -1.9052   3.0618  0.126  1.41  -0.908E+00    850.  0.227E+00  0.226E+00  0.283E-03  0.359E-13   14  26  18    1  13  16   10   8  19    0.000    0.000       0
     367  0    134 00:33:32.40  14.7000A  -2.3232  -1.8975   3.0628  0.125  1.41  -0.911E+00    850.  0.224E+00  0.224E+00  0.284E-03  0.358E-13   14  26  18    1  13  16   11   8  19    0.000    0.000       0
     368  0    135 00:33:47.10  14.6000A  -2.4176  -1.8783   3.0724  0.125  1.41  -0.913E+00    850.  0.222E+00  0.222E+00  0.283E-03  0.354E-13   14  26  18    1  12  16   11   8  19    0.000    0.000       0
     369  0    136 00:34:01.70  14.9000A  -2.5034  -1.8776   3.0299  0.124  1.38  -0.912E+00    800.  0.221E+00  0.220E+00  0.280E-03  0.348E-13   14  26  18    1  12  16   11   8  19    0.000    0.000       0
     370  0    137 00:34:16.60  15.2000A  -2.5429  -1.8640   2.9687  0.124  1.38  -0.911E+00    800.  0.219E+00  0.218E+00  0.284E-03  0.351E-13   14  26  18    1  12  16   11   8  19    0.000    0.000       0
     371  0    138 00:34:31.80  15.3000A  -2.4934  -1.8349   2.9455  0.123  1.38  -0.912E+00    800.  0.217E+00  0.216E+00  0.289E-03  0.355E-13   14  26  18    1  12  16   12   8  19    0.000    0.000       0
     372  0    139 00:34:47.10  15.4000A  -2.3593  -1.7912   2.9244  0.122  1.38  -0.917E+00    800.  0.216E+00  0.215E+00  0.290E-03  0.354E-13   14  27  18    1  12  16   12   8  19    0.000    0.000       0
     373  0    140 00:35:02.50  15.5000A  -2.3911  -1.7717   2.8945  0.122  1.38  -0.917E+00    800.  0.215E+00  0.214E+00  0.291E-03  0.353E-13   15  26  18    1  11  15   12   8  19    0.000    0.000       0
     374  0    141 00:35:18.00  15.8000A  -2.4694  -1.7378   2.8459  0.121  1.38  -0.922E+00    800.  0.214E+00  0.213E+00  0.292E-03  0.352E-13   15  26  17    1  11  15   12   8  19    0.000    0.000       0
     375  0    142 00:35:33.80  15.8000A  -2.4673   1.7977   2.8489  0.121  1.38  -0.929E+00    800.  0.213E+00  0.212E+00  0.297E-03  0.355E-13   15  26  17    8  34  23   14  26  24    0.000    0.000       0
     376  0    143 00:35:49.60  15.7000A  -2.3617  -1.7673   2.8672  0.120  1.38  -0.927E+00    800.  0.212E+00  0.212E+00  0.296E-03  0.352E-13   15  26  17    1   8  35    6  38  32    0.000    0.000       0
     377  0    144 00:36:05.30  15.2000A  -2.2459  -1.7490   2.9526  0.119  1.38  -0.930E+00    800.  0.212E+00  0.211E+00  0.295E-03  0.347E-13   15  26  16    1   7  34    6  38  32    0.000    0.000       0
     378  0    145 00:36:20.50  14.7000A  -2.0976  -1.8137   3.0687  0.119  1.38  -0.939E+00    800.  0.212E+00  0.211E+00  0.286E-03  0.335E-13   15  26  16    1   6  34    7  38  32    0.000    0.000       0
     379  0    146 00:36:35.20  14.2000A  -2.0336  -1.8432   3.1686  0.118  1.38  -0.941E+00    800.  0.211E+00  0.210E+00  0.276E-03  0.323E-13   11  11  13    1   6  34    7  38  32    0.000    0.000       0
     380  0    147 00:36:49.40  13.9000A  -2.0496  -1.8095   3.2330  0.118  1.38  -0.945E+00    800.  0.211E+00  0.210E+00  0.267E-03  0.311E-13   10  12  14    1   6  34    7  38  32    0.000    0.000       0
     381  0    148 00:37:03.30  13.6000A  -2.1261  -1.8873   3.2989  0.117  1.38  -0.949E+00    800.  0.211E+00  0.210E+00  0.262E-03  0.304E-13   10  12  14    1   5  34    8  38  32    0.000    0.000       0
     382  0    149 00:37:16.90  13.4000A  -2.1436  -1.8857   3.3518  0.117  1.38  -0.952E+00    800.  0.211E+00  0.210E+00  0.256E-03  0.297E-13   11  12  13    1   5  34    8  38  32    0.000    0.000       0
     383  0    150 00:37:30.30  13.4000A  -2.1473  -1.8654   3.3558  0.116  1.38  -0.957E+00    800.  0.211E+00  0.210E+00  0.253E-03  0.292E-13   11  12  13    1   6  35    8  38  32    0.000    0.000       0
     384  0    151 00:37:43.70  13.5000A  -2.1151  -1.8964   3.3332  0.116  1.38  -0.961E+00    800.  0.211E+00  0.210E+00  0.253E-03  0.291E-13   12  12  12    1   5  35    9  38  32    0.000    0.000       0
     385  0    152 00:37:57.20  13.6000A  -2.0520  -1.9678   3.3051  0.116  1.38  -0.964E+00    800.  0.211E+00  0.210E+00  0.255E-03  0.292E-13   12  12  12    1   5  35    9  38  32    0.000    0.000       0
     386  0    153 00:38:10.80  13.9000A  -1.9710  -1.9302   3.2271  0.115  1.41  -0.968E+00    850.  0.211E+00  0.210E+00  0.257E-03  0.292E-13   13  12  12    1   5  35    9  38  32    0.000    0.000       0
     387  0    154 00:38:24.70  14.2000A   1.9577  -2.0346   3.1608  0.115  1.41  -0.973E+00    850.  0.211E+00  0.210E+00  0.262E-03  0.297E-13    7  28  21    1   4  35    8  10   7    0.000    0.000       0
     388  0    155 00:38:38.90  13.4000A   1.9887  -2.0495   3.3484  0.115  1.41  -0.974E+00    850.  0.211E+00  0.210E+00  0.268E-03  0.301E-13    7  28  22    1   4  35    8   9   7    0.000    0.000       0
     389  0    156 00:38:52.30  12.8000A   2.0050  -1.9579   3.5087  0.114  1.41  -0.976E+00    850.  0.211E+00  0.210E+00  0.254E-03  0.284E-13    7  28  22    1   4  35    9   9   7    0.000    0.000       0
     390  0    157 00:39:05.10  12.6000A   1.9496  -2.0192   3.5707  0.114  1.41  -0.980E+00    850.  0.211E+00  0.210E+00  0.244E-03  0.271E-13    7  28  22    1   3  35    9   9   7    0.000    0.000       0
     391  0    158 00:39:17.70  12.5000A  -1.8998  -1.9977   3.5937  0.114  1.41  -0.982E+00    850.  0.211E+00  0.210E+00  0.240E-03  0.266E-13   14  13  11    1   3  35   10   9   7    0.000    0.000       0
     392  0    159 00:39:30.20  12.8000A  -1.8349  -1.8953   3.5266  0.114  1.41  -0.982E+00    850.  0.211E+00  0.210E+00  0.239E-03  0.264E-13   14  13  11    1   3  35   10   8   6    0.000    0.000       0
     393  0    160 00:39:43.00  13.0000A  -1.8421   1.8536   3.4659  0.114  1.41  -0.988E+00    850.  0.211E+00  0.210E+00  0.245E-03  0.270E-13    8  13  14   17  15  11   10   8   6    0.000    0.000       0
     394  0    161 00:39:56.00  13.1000A  -1.8732   1.9441   3.4251  0.113  1.41  -0.987E+00    850.  0.212E+00  0.210E+00  0.250E-03  0.274E-13    8  13  14   17  15  11   11   8   6    0.000    0.000       0
     395  0    162 00:40:09.10  13.1000A  -1.8856   1.9865   3.4420  0.113  1.41  -0.989E+00    850.  0.212E+00  0.211E+00  0.253E-03  0.276E-13    8  13  14   17  15  11   10   9   8    0.000    0.000       0
     396  0    163 00:40:22.20  12.9000A  -1.8945   1.9871   3.4953  0.113  1.41  -0.992E+00    850.  0.212E+00  0.211E+00  0.254E-03  0.277E-13    9  13  14   17  16  11   10   9   8    0.000    0.000       0
     397  0    164 00:40:35.10  13.0000A  -1.8988   2.1091   3.4483  0.113  1.41  -0.994E+00    850.  0.212E+00  0.211E+00  0.251E-03  0.273E-13    9  13  14   17  16  11   11   9   8    0.000    0.000       0
     398  0    165 00:40:48.10  13.3000A  -1.9469   2.1604   3.3864  0.113  1.41  -0.994E+00    850.  0.212E+00  0.211E+00  0.254E-03  0.276E-13   18  13   9   17  16  11   11   8   7    0.000    0.000       0
     399  0    166 00:41:01.40  13.5000A  -1.9766   2.1203   3.3305  0.113  1.41  -0.996E+00    850.  0.213E+00  0.211E+00  0.261E-03  0.283E-13   18  13   9   17  16  11   12   8   7    0.000    0.000       0
     400  0    167 00:41:14.90  13.7000A  -1.9835   2.1254   3.2826  0.113  1.41  -0.100E+01    850.  0.213E+00  0.212E+00  0.266E-03  0.288E-13   18  14   9   17  17  11   11   8   8    0.000    0.000       0
     401  0    168 00:41:28.60  13.6000A  -1.9910   2.1849   3.3065  0.113  1.41  -0.100E+01    850.  0.213E+00  0.212E+00  0.271E-03  0.293E-13   18  14   9   17  17  11   12   8   8    0.000    0.000       0
     402  0    169 00:41:42.20  13.4000A  -1.9808   2.1753   3.3472  0.113  1.41  -0.101E+01    850.  0.214E+00  0.212E+00  0.270E-03  0.292E-13   18  14   9   18  18   9   12   8   8    0.000    0.000       0
     403  0    170 00:41:55.60  13.5000A  -1.9940   2.2417   3.3229  0.113  1.38  -0.101E+01    800.  0.214E+00  0.213E+00  0.267E-03  0.288E-13   18  14   8   17  18  10   12   8   8    0.000    0.000       0
     404  0    171 00:42:09.10  13.5000A  -2.0081   2.2342   3.3377  0.113  1.38  -0.101E+01    800.  0.214E+00  0.213E+00  0.270E-03  0.291E-13   18  14   8   17  18  10   13   8   8    0.000    0.000       0
     405  0    172 00:42:22.60  13.6000A  -1.9959   2.1839   3.3086  0.113  1.38  -0.101E+01    800.  0.215E+00  0.213E+00  0.270E-03  0.292E-13   18  14   7   18  19   8   13   8   8    0.000    0.000       0
     406  0    173 00:42:36.20  13.9000A  -1.9833   2.2213   3.2258  0.113  1.38  -0.101E+01    800.  0.216E+00  0.214E+00  0.273E-03  0.295E-13   18  14   7   18  19   8   13   8   8    0.000    0.000       0
     407  0    174 00:42:50.10  14.1000A  -2.0015   2.1933   3.1995  0.113  1.38  -0.102E+01    800.  0.216E+00  0.215E+00  0.279E-03  0.301E-13   18  15   7   17  19  10   14   8   8    0.000    0.000       0
     408  0    175 00:43:04.20  14.3000A  -1.9689   2.2253   3.1383  0.113  1.38  -0.101E+01    800.  0.217E+00  0.216E+00  0.283E-03  0.306E-13   18  15   7   17  20   9   14   8   8    0.000    0.000       0
     409  0    176 00:43:18.50  14.8000A  -1.9216   2.2244   3.0482  0.114  1.38  -0.100E+01    800.  0.218E+00  0.217E+00  0.287E-03  0.310E-13   17  14   4   17  20   9   14   8   8    0.000    0.000       0
     410  0    177 00:43:33.30  15.2000A  -1.9577   2.1655   2.9562  0.114  1.38  -0.100E+01    800.  0.219E+00  0.218E+00  0.297E-03  0.321E-13   17  14   4   17  20   9   14   8   8    0.000    0.000       0
     411  0    178 00:43:48.50  15.5000A  -1.9068   2.1078   2.9113  0.114  1.38  -0.994E+00    800.  0.221E+00  0.219E+00  0.305E-03  0.330E-13   17  14   4   18  18  12   10   9  18    0.000    0.000       0
     412  0    179 00:44:04.00  15.0000A  -1.8774   2.1347   2.9999  0.114  1.38  -0.100E+01    800.  0.222E+00  0.220E+00  0.311E-03  0.336E-13   17  15   4   18  19  11   10   9  18    0.000    0.000       0
     413  0    180 00:44:19.00  14.7000A  -1.8143   2.1133   3.0704  0.115  1.38  -0.992E+00    800.  0.223E+00  0.222E+00  0.302E-03  0.326E-13   16  15   3   18  17  13   11   9  18    0.000    0.000       0
     414  0    181 00:44:33.70  14.6000A  -1.7675   2.1428   3.0917  0.115  1.43  -0.991E+00    900.  0.225E+00  0.223E+00  0.297E-03  0.320E-13   16  15   3   18  17  13   13  11  16    0.000    0.000       0
     415  0    182 00:44:48.30  14.2000A   1.7562   2.1477   3.1647  0.115  1.43  -0.994E+00    900.  0.226E+00  0.225E+00  0.296E-03  0.319E-13    9  27  23   18  17  13   13  11  16    0.000    0.000       0
     416  0    183 00:45:02.50  13.9000A   1.7855   2.1698   3.2331  0.116  1.43  -0.990E+00    900.  0.228E+00  0.226E+00  0.288E-03  0.311E-13    9  27  23   18  18  13   14  11  16    0.000    0.000       0
     417  0    184 00:45:16.40  13.6000A   1.7978   2.1955   3.3149  0.116  1.43  -0.992E+00    900.  0.229E+00  0.228E+00  0.283E-03  0.306E-13    9  27  23   18  18  13   14  11  16    0.000    0.000       0
     418  0    185 00:45:30.00  13.3000A   1.8111   2.2046   3.3800  0.116  1.43  -0.100E+01    900.  0.231E+00  0.229E+00  0.278E-03  0.300E-13    9  27  24   18  18  13   14  11  16    0.000    0.000       0
     419  0    186 00:45:43.30  13.2000A   1.8583   2.1999   3.4136  0.117  1.43  -0.991E+00    900.  0.232E+00  0.231E+00  0.273E-03  0.294E-13    9  27  24   18  18  13   14  11  16    0.000    0.000       0
     420  0    187 00:45:56.50  13.1000A   1.8922   2.1820   3.4441  0.117  1.43  -0.987E+00    900.  0.234E+00  0.233E+00  0.272E-03  0.293E-13    9  27  24   18  18  13   15  11  16    0.000    0.000       0
     421  0    188 00:46:09.60  13.0000A   1.9096   2.1503   3.4749  0.117  1.43  -0.985E+00    900.  0.236E+00  0.235E+00  0.271E-03  0.292E-13    9  27  24   18  18  13   15  11  16    0.000    0.000       0
     422  0    189 00:46:22.60  13.0000A   1.9088   2.1015   3.4708  0.118  1.43  -0.984E+00    900.  0.238E+00  0.236E+00  0.270E-03  0.291E-13    9  27  24   18  18  13   15  11  16    0.000    0.000       0
     423  0    190 00:46:35.60  13.0000A   1.8889   2.0704   3.4500  0.118  1.43  -0.984E+00    900.  0.239E+00  0.238E+00  0.271E-03  0.292E-13    9  27  24   18  19  12   16  11  16    0.000    0.000       0
     424  0    191 00:46:48.60  13.0000A   1.8877   2.0526   3.4578  0.118  1.43  -0.985E+00    900.  0.241E+00  0.240E+00  0.272E-03  0.293E-13    9  26  25   18  19  12   16  11  16    0.000    0.000       0
     425  0    192 00:47:01.60  13.1000A   1.8905   2.0533   3.4417  0.119  1.41  -0.981E+00    850.  0.243E+00  0.242E+00  0.273E-03  0.294E-13    9  27  25    8   4   0   16  11  16    0.000    0.000       0
     426  0    193 00:47:14.70  13.2000A   1.9015   2.0285   3.4092  0.119  1.41  -0.979E+00    850.  0.245E+00  0.244E+00  0.276E-03  0.298E-13    9  27  25    8   4   0   13  10  17    0.000    0.000       0
     427  0    194 00:47:27.90  13.3000A   1.8909   2.0073   3.3854  0.119  1.41  -0.974E+00    850.  0.247E+00  0.246E+00  0.280E-03  0.303E-13    9  27  25   19  18  13   14  10  17    0.000    0.000       0
     428  0    195 00:47:41.20  13.2000A   1.8588  -2.0241   3.3991  0.119  1.41  -0.970E+00    850.  0.249E+00  0.247E+00  0.284E-03  0.307E-13    9  27  26   13  38  15   14  10  17    0.000    0.000       0
     429  0    196 00:47:54.40  13.2000A   1.8886  -1.9930   3.4144  0.119  1.41  -0.969E+00    850.  0.251E+00  0.249E+00  0.283E-03  0.308E-13    9  27  26   13  38  15   11   9  18    0.000    0.000       0
     430  0    197 00:48:07.59  12.9000A  -2.0108  -1.9302   3.4806  0.119  1.41  -0.973E+00    850.  0.253E+00  0.251E+00  0.285E-03  0.311E-13    6   3   1   13  38  15   12   9  18    0.000    0.000       0
     431  0    198 00:48:20.50  12.8000A  -2.0666   1.9227   3.5051  0.119  1.41  -0.971E+00    850.  0.254E+00  0.253E+00  0.280E-03  0.307E-13    6   3   1   19  19  13   12   9  18    0.000    0.000       0
     432  0    199 00:48:33.30  12.8000A  -2.0933   1.8878   3.5261  0.119  1.41  -0.972E+00    850.  0.256E+00  0.255E+00  0.280E-03  0.307E-13    6   4   0   19  19  13   13   9  18    0.000    0.000       0
     433  0    200 00:48:46.10  12.7000A   2.1503   1.8281   3.5401  0.119  1.41  -0.975E+00    850.  0.258E+00  0.257E+00  0.282E-03  0.310E-13    7   2  10   12   4  39   13   9  18    0.000    0.000       0
     434  0    201 00:48:58.80  12.7000A   2.1949   1.8433   3.5356  0.119  1.41  -0.974E+00    850.  0.260E+00  0.258E+00  0.281E-03  0.311E-13    7   2  10   12   4  39   14   9  18    0.000    0.000       0
     435  0    202 00:49:11.50  12.5000A  -2.2768  -1.8583   3.5983  0.119  1.41  -0.976E+00    850.  0.262E+00  0.260E+00  0.283E-03  0.314E-13    7   3   0    3   8   2    9  39   1    0.000    0.000       0
     436  0    203 00:49:24.00  12.9000A  -2.3194  -1.9224   3.4924  0.119  1.41  -0.973E+00    850.  0.264E+00  0.262E+00  0.281E-03  0.312E-13    7   3   0    3   8   2   15   9  18    0.000    0.000       0
     437  0    204 00:49:36.90  12.8000A  -2.2131  -1.9682   3.5215  0.119  1.41  -0.975E+00    850.  0.266E+00  0.264E+00  0.291E-03  0.325E-13    8   2   0    3   8   2    9   0   0    0.000    0.000       0
     438  0    205 00:49:49.70  12.6000A  -2.3826  -2.0375   3.5652  0.119  1.41  -0.975E+00    850.  0.267E+00  0.265E+00  0.291E-03  0.325E-13    8   1   0    3   7   2    9   7  17    0.000    0.000       0
     439  0    206 00:50:02.30  12.4000A  -2.4158  -2.1010   3.6437  0.118  1.41  -0.973E+00    850.  0.269E+00  0.267E+00  0.287E-03  0.322E-13    8   1   0    3   7   2    9   7  17    0.000    0.000       0
     440  0    207 00:50:14.70  12.2000A  -2.3751  -2.1213   3.6855  0.118  1.41  -0.972E+00    850.  0.270E+00  0.268E+00  0.284E-03  0.319E-13    8   1   0    3   7   2   10   7  17    0.000    0.000       0
     441  0    208 00:50:26.90  12.0000A  -2.4818  -2.0787   3.7439  0.118  1.43  -0.972E+00    900.  0.272E+00  0.270E+00  0.281E-03  0.316E-13    9   1  39    3   7   2   10   7  17    0.000    0.000       0
     442  0    209 00:50:38.90  12.0000A  -2.5308  -2.0324   3.7515  0.118  1.43  -0.970E+00    900.  0.273E+00  0.271E+00  0.277E-03  0.313E-13    9   1  39    4   7   3   10   7  17    0.000    0.000       0
     443  0    210 00:50:50.90  12.1000A  -2.4937  -2.0321   3.7136  0.118  1.43  -0.968E+00    900.  0.274E+00  0.272E+00  0.278E-03  0.315E-13    9   1  39    3   6   2   10   7  17    0.000    0.000       0
     444  0    211 00:51:03.00  12.0000A  -2.4204  -1.9727   3.7591  0.118  1.43  -0.968E+00    900.  0.275E+00  0.273E+00  0.281E-03  0.319E-13    9   1  39    3   6   2   14   6  18    0.000    0.000       0
     445  0    212 00:51:15.00  11.7000A  -2.3845  -2.0374   3.8360  0.118  1.43  -0.970E+00    900.  0.277E+00  0.274E+00  0.280E-03  0.318E-13    9   1  38   15  33  38   14   6  18    0.000    0.000       0
     446  0    213 00:51:26.69  11.7000A  -2.3517  -2.1313   3.8587  0.117  1.43  -0.975E+00    900.  0.278E+00  0.275E+00  0.274E-03  0.312E-13   10   1  38   15  33  38   14   6  18    0.000    0.000       0
     447  0    214 00:51:38.39  11.8000A  -2.3208  -2.1277   3.8253  0.117  1.43  -0.977E+00    900.  0.279E+00  0.276E+00  0.275E-03  0.313E-13   10   1  38   15  33  38   14   6  18    0.000    0.000       0
     448  0    215 00:51:50.19  12.0000A  -2.2790  -2.0966   3.7533  0.117  1.43  -0.981E+00    900.  0.279E+00  0.277E+00  0.278E-03  0.317E-13   11   0  38   15  32  38   15   5  18    0.000    0.000       0
     449  0    216 00:52:02.19  12.0000A  -2.2486  -2.0403   3.7622  0.117  1.43  -0.984E+00    900.  0.280E+00  0.278E+00  0.284E-03  0.324E-13   11   0  38   15  32  38   16   5  18    0.000    0.000       0
     450  0    217 00:52:14.19  12.0000A  -2.2214  -1.9990   3.7541  0.117  1.43  -0.980E+00    900.  0.281E+00  0.278E+00  0.285E-03  0.326E-13   11   0  37   15  32  37   16   5  18    0.000    0.000       0
     451  0    218 00:52:26.19  12.2000A  -2.1797  -2.0654   3.6993  0.117  1.43  -0.978E+00    900.  0.281E+00  0.279E+00  0.286E-03  0.327E-13   11   0  37   15  31  37   16   5  18    0.000    0.000       0
     452  0    219 00:52:38.39  12.6000A  -2.0979  -2.1467   3.5835  0.117  1.43  -0.978E+00    900.  0.282E+00  0.279E+00  0.291E-03  0.333E-13   11   0  37   15  31  37   16   5  18    0.000    0.000       0
     453  0    220 00:52:50.99  13.2000A  -1.9948  -2.1573   3.4108  0.116  1.43  -0.979E+00    900.  0.282E+00  0.280E+00  0.301E-03  0.344E-13   13  23   4   15  31  37   16   5  18    0.000    0.000       0
     454  0    221 00:53:04.19  13.6000A  -1.9884  -2.0801   3.2998  0.116  1.43  -0.983E+00    900.  0.282E+00  0.280E+00  0.315E-03  0.360E-13   13  23   4   15  31  37   17   4  18    0.000    0.000       0
     455  0    222 00:53:17.79  13.6000A  -1.9649  -2.0529   3.2994  0.116  1.43  -0.992E+00    900.  0.282E+00  0.280E+00  0.324E-03  0.369E-13   13  23   4   15  30  37    8   4   4    0.000    0.000       0
     456  0    223 00:53:31.39  13.6000A  -1.9354   2.0467   3.3034  0.116  1.43  -0.993E+00    900.  0.282E+00  0.280E+00  0.324E-03  0.368E-13   13  23   3   17  17  10    8  35   5    0.000    0.000       0
     457  0    224 00:53:44.99  13.2000A  -1.9955   2.0016   3.4048  0.116  1.43  -0.987E+00    900.  0.282E+00  0.279E+00  0.323E-03  0.366E-13    9   6   4   17  17  10    9  35   5    0.000    0.000       0
     458  0    225 00:53:58.19  13.1000A  -2.0497   1.9075   3.4432  0.116  1.43  -0.988E+00    900.  0.282E+00  0.279E+00  0.314E-03  0.354E-13    9   6   4   17  17  10    9  35   5    0.000    0.000       0
     459  0    226 00:54:11.29  13.1000A  -2.0943   1.8859   3.4378  0.115  1.43  -0.993E+00    900.  0.282E+00  0.279E+00  0.311E-03  0.350E-13   10   6   4   13   8   4   10  35   5    0.000    0.000       0
     460  0    227 00:54:24.39  13.0000A  -2.1924   1.9114   3.4684  0.115  1.43  -0.101E+01    900.  0.282E+00  0.279E+00  0.311E-03  0.348E-13   10   6   4   13   8   4   10  35   5    0.000    0.000       0
     461  0    228 00:54:37.39  13.1000A  -2.2396   1.9699   3.4311  0.115  1.43  -0.995E+00    900.  0.281E+00  0.279E+00  0.309E-03  0.344E-13   10   6   4   14   9   5   10  35   5    0.000    0.000       0
     462  0    229 00:54:50.49  12.8000A  -2.2345   2.0185   3.5139  0.115  1.40  -0.994E+00    850.  0.281E+00  0.279E+00  0.311E-03  0.346E-13   10   6   4   14   9   5   10  36   6    0.000    0.000       0
     463  0    230 00:55:03.29  12.5000A  -2.2386   2.0042   3.5923  0.115  1.40  -0.996E+00    850.  0.281E+00  0.278E+00  0.305E-03  0.338E-13   11   6   4   14   9   5   10  36   6    0.000    0.000       0
     464  0    231 00:55:15.79  12.4000A  -2.2667   1.9491   3.6309  0.115  1.40  -0.990E+00    850.  0.281E+00  0.278E+00  0.298E-03  0.330E-13   11   6   4   14   9   5   10  36   6    0.000    0.000       0
     465  0    232 00:55:28.19  12.3000A  -2.2639   1.9620   3.6661  0.115  1.40  -0.990E+00    850.  0.281E+00  0.278E+00  0.297E-03  0.327E-13   11   6   4   17  10   6   11  36   6    0.000    0.000       0
     466  0    233 00:55:40.49  12.1000A  -2.2609   1.9583   3.7059  0.115  1.40  -0.992E+00    850.  0.281E+00  0.278E+00  0.295E-03  0.324E-13   11   6   3   17  10   6   11  36   6    0.000    0.000       0
     467  0    234 00:55:52.59  12.1000A  -2.2079   1.9095   3.7294  0.115  1.40  -0.989E+00    850.  0.280E+00  0.278E+00  0.291E-03  0.319E-13   11   6   3   12   3  37   11  36   6    0.000    0.000       0
     468  0    235 00:56:04.69  12.1000A  -2.1707   1.8977   3.7131  0.115  1.40  -0.986E+00    850.  0.280E+00  0.278E+00  0.292E-03  0.319E-13   12   6   3   13   2  38   11  36   6    0.000    0.000       0
     469  0    236 00:56:16.79  12.3000A  -2.1897   1.9603   3.6503  0.115  1.40  -0.985E+00    850.  0.280E+00  0.277E+00  0.292E-03  0.320E-13   12   6   3   13   3  36   11  36   6    0.000    0.000       0
     470  0    237 00:56:29.09  12.3000A  -2.1717   1.9372   3.6538  0.115  1.40  -0.985E+00    850.  0.280E+00  0.277E+00  0.297E-03  0.325E-13   12   6   3   13   3  37   12  36   7    0.000    0.000       0
     471  0    238 00:56:41.39  12.4000A  -2.1149   1.9536   3.6151  0.115  1.40  -0.987E+00    850.  0.280E+00  0.277E+00  0.298E-03  0.325E-13   12   6   3   13   4  36   12  36   7    0.000    0.000       0
     472  0    239 00:56:53.79  12.5000A  -2.0201   1.9144   3.5932  0.115  1.40  -0.988E+00    850.  0.279E+00  0.277E+00  0.301E-03  0.328E-13   12   6   3   13   4  37   13  36   7    0.000    0.000       0
     473  0    240 00:57:06.29  12.7000A  -1.9931   1.9369   3.5414  0.115  1.38  -0.988E+00    800.  0.279E+00  0.277E+00  0.304E-03  0.332E-13   13   6   3   14   4  37   13  36   7    0.000    0.000       0
     474  0    241 00:57:18.99  12.8000A   2.0384   1.9510   3.5139  0.115  1.38  -0.991E+00    800.  0.279E+00  0.276E+00  0.310E-03  0.338E-13   12  37   9   14   5  37   14  35   6    0.000    0.000       0
     475  0    242 00:57:31.79  12.6000A   2.1168   1.9781   3.5652  0.115  1.40  -0.993E+00    850.  0.279E+00  0.276E+00  0.313E-03  0.341E-13   13  37   9   14   5  37   14  35   6    0.000    0.000       0
     476  0    243 00:57:44.39  12.6000A   2.1807   1.9523   3.5655  0.114  1.40  -0.998E+00    850.  0.279E+00  0.276E+00  0.308E-03  0.337E-13   13  37   9   14   5  37   14  35   6    0.000    0.000       0
     477  0    244 00:57:56.99  12.7000A   2.1971   1.9347   3.5448  0.114  1.40  -0.997E+00    850.  0.279E+00  0.276E+00  0.308E-03  0.337E-13   13  37   9   15   5  37   15  35   6    0.000    0.000       0
     478  0    245 00:58:09.69  12.7000A   2.1710   1.9616   3.5440  0.114  1.40  -0.997E+00    850.  0.279E+00  0.276E+00  0.311E-03  0.341E-13   13  37  10   15   5  37   15  35   6    0.000    0.000       0
     479  0    246 00:58:22.39  12.9000A   2.2295   1.9868   3.4858  0.114  1.38  -0.998E+00    800.  0.278E+00  0.276E+00  0.311E-03  0.341E-13   14  37  10   15   6  37   15  35   6    0.000    0.000       0
     480  0    247 00:58:35.29  13.3000A   2.2611   1.9834   3.3800  0.114  1.38  -0.997E+00    800.  0.278E+00  0.276E+00  0.315E-03  0.347E-13   14  37  10   15   6  37   15  35   6    0.000    0.000       0
     481  0    248 00:58:48.59  13.6000A   2.2233   1.9450   3.3102  0.114  1.38  -0.999E+00    800.  0.278E+00  0.276E+00  0.325E-03  0.358E-13   14  37  10   15   6  37   15  34   3    0.000    0.000       0
     482  0    249 00:59:02.19  13.6000A   2.2127   1.9673   3.3017  0.114  1.38  -0.100E+01    800.  0.278E+00  0.276E+00  0.333E-03  0.366E-13   14  37  11   16   6  37    9  37   5    0.000    0.000       0
     483  0    250 00:59:15.79  13.5000A   2.1257   1.9989   3.3397  0.113  1.38  -0.100E+01    800.  0.278E+00  0.275E+00  0.333E-03  0.366E-13   14  37  11   16   7  37    9  37   5    0.000    0.000       0
     484  0    251 00:59:29.29  13.4000A   2.0291   2.0298   3.3500  0.113  1.38  -0.101E+01    800.  0.278E+00  0.275E+00  0.330E-03  0.363E-13   16  36  11   16   6  36   10  37   5    0.000    0.000       0
     485  0    252 00:59:42.69  13.2000A   1.9435   2.0214   3.3980  0.113  1.43  -0.101E+01    900.  0.278E+00  0.275E+00  0.328E-03  0.360E-13   15  37  12   16   7  37   10  37   5    0.000    0.000       0
     486  0    253 00:59:55.89  13.3000A   1.9283   2.0368   3.3876  0.113  1.43  -0.102E+01    900.  0.278E+00  0.275E+00  0.322E-03  0.353E-13   16  36  12   16   7  36   10  37   5    0.000    0.000       0
     487  0    254 01:00:09.19  13.2000A   1.9061   2.0144   3.3966  0.113  1.43  -0.102E+01    900.  0.277E+00  0.275E+00  0.325E-03  0.355E-13   16  36  12   16   7  36   11  37   5    0.000    0.000       0
  • palm/trunk/SOURCE/check_parameters.f90

    r3049 r3065  
    2525! -----------------
    2626! $Id$
     27! dz was replaced by dz(1), error message revised
     28!
     29! 3049 2018-05-29 13:52:36Z Giersch
    2730! add variable description
    2831!
     
    41784181    IF ( ( constant_flux_layer .OR.  &
    41794182           INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 )       &
    4180          .AND. roughness_length >= 0.5 * dz )  THEN
     4183         .AND. roughness_length >= 0.5 * dz(1) )  THEN
    41814184       message_string = 'roughness_length must be smaller than dz/2'
    41824185       CALL message( 'check_parameters', 'PA0424', 1, 2, 0, 6, 0 )
     
    41984201!-- Check if vertical grid stretching is switched off in case of complex
    41994202!-- terrain simulations
    4200     IF ( complex_terrain  .AND.  dz_stretch_level < 100000.0_wp )  THEN
     4203    IF ( complex_terrain  .AND.                                                &
     4204         ANY( dz_stretch_level_start /= -9999999.9_wp ) )  THEN
    42014205       message_string = 'Vertical grid stretching is not allowed for ' //      &
    42024206                        'complex_terrain = .T.'
  • palm/trunk/SOURCE/header.f90

    r3045 r3065  
    2525! -----------------
    2626! $Id$
     27! Header output concerning stretching revised
     28!
     29! 3045 2018-05-28 07:55:41Z Giersch
    2730! Error messages revised
    2831!
     
    800803!-- Computational grid
    801804    IF ( .NOT. ocean )  THEN
    802        WRITE ( io, 250 )  dx, dy, dz, (nx+1)*dx, (ny+1)*dy, zu(nzt+1)
    803        IF ( dz_stretch_level_index < nzt+1 )  THEN
    804           WRITE ( io, 252 )  dz_stretch_level, dz_stretch_level_index, &
    805                              dz_stretch_factor, dz_max
    806        ENDIF
     805       WRITE ( io, 250 )  dx, dy
     806       
     807       DO i = 1, number_stretch_level_start+1
     808          WRITE ( io, 253 )  i, dz(i)
     809       ENDDO
     810       
     811       WRITE( io, 251 ) (nx+1)*dx, (ny+1)*dy, zu(nzt+1)
     812       
     813       IF ( ANY( dz_stretch_level_start_index < nzt+1 ) )  THEN
     814          WRITE( io, '(A)', advance='no') ' Vertical stretching starts at height:'
     815          DO i = 1, number_stretch_level_start
     816             WRITE ( io, '(F10.1,A3)', advance='no' )  dz_stretch_level_start(i), ' m,'
     817          ENDDO
     818          WRITE( io, '(/,A)', advance='no') ' Vertical stretching starts at index: '
     819          DO i = 1, number_stretch_level_start
     820             WRITE ( io, '(I12,A1)', advance='no' )  dz_stretch_level_start_index(i), ','
     821          ENDDO
     822          WRITE( io, '(/,A)', advance='no') ' Vertical stretching ends at height:  '
     823          DO i = 1, number_stretch_level_start
     824             WRITE ( io, '(F10.1,A3)', advance='no' )  dz_stretch_level_end(i), ' m,'
     825          ENDDO
     826          WRITE( io, '(/,A)', advance='no') ' Vertical stretching ends at index:   '
     827          DO i = 1, number_stretch_level_start
     828             WRITE ( io, '(I12,A1)', advance='no' )  dz_stretch_level_end_index(i), ','
     829          ENDDO
     830          WRITE( io, '(/,A)', advance='no') ' Factor used for stretching:          '
     831          DO i = 1, number_stretch_level_start
     832             WRITE ( io, '(F12.3,A1)', advance='no' )  dz_stretch_factor_array(i), ','
     833          ENDDO
     834       ENDIF
     835       
    807836    ELSE
    808        WRITE ( io, 250 )  dx, dy, dz, (nx+1)*dx, (ny+1)*dy, zu(0)
    809        IF ( dz_stretch_level_index > 0 )  THEN
    810           WRITE ( io, 252 )  dz_stretch_level, dz_stretch_level_index, &
    811                              dz_stretch_factor, dz_max
    812        ENDIF
    813     ENDIF
    814     WRITE ( io, 254 )  nx, ny, nzt+1, MIN( nnx, nx+1 ), MIN( nny, ny+1 ), &
     837       WRITE ( io, 250 )  dx, dy
     838       DO i = 1, number_stretch_level_start+1
     839          WRITE ( io, 253 )  i, dz(i)
     840       ENDDO
     841       
     842       WRITE ( io, 251 ) (nx+1)*dx, (ny+1)*dy, zu(0)
     843       
     844       IF ( ANY( dz_stretch_level_start_index > 0 ) )  THEN
     845          WRITE( io, '(A)', advance='no') ' Vertical stretching starts at height:'
     846          DO i = 1, number_stretch_level_start
     847             WRITE ( io, '(F10.1,A3)', advance='no' )  dz_stretch_level_start(i), ' m,'
     848          ENDDO
     849          WRITE( io, '(/,A)', advance='no') ' Vertical stretching starts at index: '
     850          DO i = 1, number_stretch_level_start
     851             WRITE ( io, '(I12,A1)', advance='no' )  dz_stretch_level_start_index(i), ','
     852          ENDDO
     853          WRITE( io, '(/,A)', advance='no') ' Vertical stretching ends at height:  '
     854          DO i = 1, number_stretch_level_start
     855             WRITE ( io, '(F10.1,A3)', advance='no' )  dz_stretch_level_end(i), ' m,'
     856          ENDDO
     857          WRITE( io, '(/,A)', advance='no') ' Vertical stretching ends at index:   '
     858          DO i = 1, number_stretch_level_start
     859             WRITE ( io, '(I12,A1)', advance='no' )  dz_stretch_level_end_index(i), ','
     860          ENDDO
     861          WRITE( io, '(/,A)', advance='no') ' Factor used for stretching:          '
     862          DO i = 1, number_stretch_level_start
     863             WRITE ( io, '(F12.3,A1)', advance='no' )  dz_stretch_factor_array(i), ','
     864          ENDDO
     865       ENDIF
     866    ENDIF
     867    WRITE ( io, 254 )  nx, ny, nzt+1, MIN( nnx, nx+1 ), MIN( nny, ny+1 ),      &
    815868                       MIN( nnz+2, nzt+2 )
    816869    IF ( sloping_surface )  WRITE ( io, 260 )  alpha_surface
     
    20832136250 FORMAT (//' Computational grid and domain size:'/ &
    20842137              ' ----------------------------------'// &
    2085               ' Grid length:      dx =    ',F7.3,' m    dy =    ',F7.3, &
    2086               ' m    dz =    ',F7.3,' m'/ &
    2087               ' Domain size:       x = ',F10.3,' m     y = ',F10.3, &
     2138              ' Grid length:      dx =    ',F8.3,' m    dy =    ',F8.3, ' m')
     2139251 FORMAT (  /' Domain size:       x = ',F10.3,' m     y = ',F10.3, &
    20882140              ' m  z(u) = ',F10.3,' m'/)
    2089 252 FORMAT (' dz constant up to ',F10.3,' m (k=',I4,'), then stretched by', &
    2090               ' factor:',F6.3/ &
    2091             ' maximum dz not to be exceeded is dz_max = ',F10.3,' m'/)
    2092 254 FORMAT (' Number of gridpoints (x,y,z):  (0:',I4,', 0:',I4,', 0:',I4,')'/ &
     2141253 FORMAT ( '                dz(',I1,') =    ', F8.3, ' m')
     2142254 FORMAT (//' Number of gridpoints (x,y,z):  (0:',I4,', 0:',I4,', 0:',I4,')'/ &
    20932143            ' Subdomain size (x,y,z):        (  ',I4,',   ',I4,',   ',I4,')'/)
    20942144260 FORMAT (/' The model has a slope in x-direction. Inclination angle: ',F6.2,&
  • palm/trunk/SOURCE/init_grid.f90

    r3051 r3065  
    2525! -----------------
    2626! $Id$
     27! New vertical stretching mechanism introduced
     28!
     29! 3051 2018-05-30 17:43:55Z suehring
    2730! Minor bugfix concerning mapping 3D buildings on top of terrain
    2831!
     
    325328               canyon_height, canyon_wall_left, canyon_wall_south,             &
    326329               canyon_width_x, canyon_width_y, constant_flux_layer,            &
    327                dp_level_ind_b, dz, dz_max, dz_stretch_factor,                  &
    328                dz_stretch_level, dz_stretch_level_index, grid_level,           &
     330               dp_level_ind_b, dz, dz_max, dz_stretch_factor,                  &   
     331               dz_stretch_factor_array, dz_stretch_level, dz_stretch_level_end,&
     332               dz_stretch_level_end_index, dz_stretch_level_start_index,       &
     333               dz_stretch_level_start, grid_level,                             &
    329334               force_bound_l, force_bound_r, force_bound_n, force_bound_s,     &
    330335               ibc_uv_b, inflow_l, inflow_n, inflow_r, inflow_s,               &
    331336               masking_method, maximum_grid_level, message_string,             &
    332                momentum_advec, nest_domain, nest_bound_l, nest_bound_n,        &
    333                nest_bound_r, nest_bound_s, ocean, outflow_l, outflow_n,        &
    334                outflow_r, outflow_s, psolver, scalar_advec, topography,        &
    335                topography_grid_convention, tunnel_height, tunnel_length,       &
    336                tunnel_width_x, tunnel_width_y, tunnel_wall_depth,              &
    337                use_surface_fluxes, use_top_fluxes, wall_adjustment_factor
     337               momentum_advec, nest_domain, nest_bound_l,                      &
     338               nest_bound_n, nest_bound_r, nest_bound_s,                       &
     339               number_stretch_level_end, number_stretch_level_start, ocean,    &
     340               outflow_l, outflow_n, outflow_r, outflow_s, psolver,            &
     341               scalar_advec, topography, topography_grid_convention,           &
     342               tunnel_height, tunnel_length, tunnel_width_x, tunnel_width_y,   &
     343               tunnel_wall_depth, use_surface_fluxes, use_top_fluxes,          &
     344               wall_adjustment_factor
    338345         
    339346    USE grid_variables,                                                        &
     
    361368    IMPLICIT NONE
    362369
    363     INTEGER(iwp) ::  i                !< index variable along x
    364     INTEGER(iwp) ::  j                !< index variable along y
    365     INTEGER(iwp) ::  k                !< index variable along z
    366     INTEGER(iwp) ::  k_top            !< topography top index on local PE
    367     INTEGER(iwp) ::  l                !< loop variable
    368     INTEGER(iwp) ::  nzb_local_max    !< vertical grid index of maximum topography height
    369     INTEGER(iwp) ::  nzb_local_min    !< vertical grid index of minimum topography height
     370    INTEGER(iwp) ::  i                           !< index variable along x
     371    INTEGER(iwp) ::  j                           !< index variable along y
     372    INTEGER(iwp) ::  k                           !< index variable along z
     373    INTEGER(iwp) ::  k_top                       !< topography top index on local PE
     374    INTEGER(iwp) ::  n                           !< loop variable for stretching
     375    INTEGER(iwp) ::  number_dz                   !< number of user-specified dz values       
     376    INTEGER(iwp) ::  nzb_local_max               !< vertical grid index of maximum topography height
     377    INTEGER(iwp) ::  nzb_local_min               !< vertical grid index of minimum topography height
    370378                                     
    371     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nzb_local      !< index for topography top at cell-center
    372     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nzb_tmp        !< dummy to calculate topography indices on u- and v-grid
     379    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nzb_local  !< index for topography top at cell-center
     380    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nzb_tmp    !< dummy to calculate topography indices on u- and v-grid
    373381
    374382    INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE ::  topo !< input array for 3D topography and dummy array for setting "outer"-flags
    375383
     384    REAL(wp) ::  dz_level_end  !< distance between calculated height level for u/v-grid and user-specified end level for stretching
    376385    REAL(wp) ::  dz_stretched  !< stretched vertical grid spacing
     386   
     387    REAL(wp), DIMENSION(:), ALLOCATABLE ::  min_dz_stretch_level_end !< Array that contains all minimum heights where the stretching can end
    377388
    378389
     
    391402!
    392403!-- Compute height of u-levels from constant grid length and dz stretch factors
    393     IF ( dz == -1.0_wp )  THEN
     404    IF ( dz(1) == -1.0_wp )  THEN
    394405       message_string = 'missing dz'
    395406       CALL message( 'init_grid', 'PA0200', 1, 2, 0, 6, 0 )
    396     ELSEIF ( dz <= 0.0_wp )  THEN
    397        WRITE( message_string, * ) 'dz=',dz,' <= 0.0'
     407    ELSEIF ( dz(1) <= 0.0_wp )  THEN
     408       WRITE( message_string, * ) 'dz=',dz(1),' <= 0.0'
    398409       CALL message( 'init_grid', 'PA0201', 1, 2, 0, 6, 0 )
    399410    ENDIF
    400411
    401412!
     413!-- Initialize dz_stretch_level_start with the value of dz_stretch_level
     414!-- if it was set by the user
     415    IF ( dz_stretch_level /= -9999999.9_wp ) THEN
     416       dz_stretch_level_start(1) = dz_stretch_level
     417    ENDIF
     418       
     419!
     420!-- Determine number of dz values and stretching levels specified by the
     421!-- user to allow right controlling of the stretching mechanism and to
     422!-- perform error checks
     423    number_dz = COUNT( dz /= -1.0_wp )
     424    number_stretch_level_start = COUNT( dz_stretch_level_start /=              &
     425                                       -9999999.9_wp )
     426    number_stretch_level_end = COUNT( dz_stretch_level_end /=                  &
     427                                      9999999.9_wp )
     428
     429!
     430!-- The number of specified end levels +1 has to be the same than the number
     431!-- of specified dz values
     432    IF ( number_dz /= number_stretch_level_end + 1 ) THEN
     433       WRITE( message_string, * ) 'The number of values for dz = ',         &
     434                                   number_dz, 'has to be the same than& ',  &
     435                                   'the number of values for ',             &
     436                                   'dz_stretch_level_end + 1 = ',           &
     437                                   number_stretch_level_end+1
     438          CALL message( 'init_grid', 'PA0156', 1, 2, 0, 6, 0 )
     439    ENDIF
     440   
     441!
     442!--    The number of specified start levels has to be the same or one less than
     443!--    the number of specified dz values
     444    IF ( number_dz /= number_stretch_level_start + 1 .AND.                  &
     445         number_dz /= number_stretch_level_start ) THEN
     446       WRITE( message_string, * ) 'The number of values for dz = ',         &
     447                                   number_dz, 'has to be the same or one ', &
     448                                   'more than& the number of values for ',  &
     449                                   'dz_stretch_level_start = ',             &
     450                                   number_stretch_level_start
     451          CALL message( 'init_grid', 'PA0211', 1, 2, 0, 6, 0 )
     452    ENDIF
     453   
     454!--    The number of specified start levels has to be the same or one more than
     455!--    the number of specified end levels
     456    IF ( number_stretch_level_start /= number_stretch_level_end + 1 .AND.   &
     457         number_stretch_level_start /= number_stretch_level_end ) THEN
     458       WRITE( message_string, * ) 'The number of values for ',              &
     459                                  'dz_stretch_level_start = ',              &
     460                                   dz_stretch_level_start, 'has to be the ',&
     461                                   'same or one more than& the number of ', &
     462                                   'values for dz_stretch_level_end = ',    &
     463                                   number_stretch_level_end
     464          CALL message( 'init_grid', 'PA0216', 1, 2, 0, 6, 0 )
     465    ENDIF
     466
     467!
     468!-- Initialize dz for the free atmosphere with the value of dz_max
     469    IF ( dz(number_stretch_level_start+1) == -1.0_wp .AND.                     &
     470         number_stretch_level_start /= 0 ) THEN
     471       dz(number_stretch_level_start+1) = dz_max
     472    ENDIF
     473       
     474!
     475!-- Initialize the stretching factor if (infinitely) stretching in the free
     476!-- atmosphere is desired (dz_stretch_level_end was not specified for the
     477!-- free atmosphere)
     478    IF ( number_stretch_level_start == number_stretch_level_end + 1 ) THEN
     479       dz_stretch_factor_array(number_stretch_level_start) =                   &
     480       dz_stretch_factor
     481    ENDIF
     482   
     483!
     484!-- Allocation of arrays for stretching
     485    ALLOCATE( min_dz_stretch_level_end(number_stretch_level_start) )
     486   
     487!
    402488!-- Define the vertical grid levels
    403489    IF ( .NOT. ocean )  THEN
     490   
     491!
     492!--    The stretching region has to be large enough to allow for a smooth
     493!--    transition between two different grid spacings
     494       DO n = 1, number_stretch_level_start
     495          min_dz_stretch_level_end(n) = dz_stretch_level_start(n) +            &
     496                                        4 * MAX( dz(n),dz(n+1) )
     497       ENDDO
     498
     499       IF ( ANY( min_dz_stretch_level_end > dz_stretch_level_end ) ) THEN
     500             message_string= 'Eeach dz_stretch_level_end has to be larger ' // &
     501                             'than its corresponding value for &' //           &
     502                             'dz_stretch_level_start + 4*MAX(dz(n),dz(n+1)) '//&
     503                             'to allow for smooth grid stretching'
     504             CALL message( 'init_grid', 'PA0224', 1, 2, 0, 6, 0 )
     505       ENDIF
     506       
     507!
     508!--    Stretching must not be applied within the prandtl_layer
     509!--    (first two grid points). For the default case dz_stretch_level_start
     510!--    is negative. Therefore the absolut value is checked here.
     511       IF ( ANY( ABS( dz_stretch_level_start ) < dz(1) * 1.5_wp ) ) THEN
     512          WRITE( message_string, * ) 'Eeach dz_stretch_level_start has to be ',&
     513                                     'larger than ', dz(1) * 1.5
     514             CALL message( 'init_grid', 'PA0226', 1, 2, 0, 6, 0 )
     515       ENDIF
     516
     517!
     518!--    The stretching has to start and end on a grid level. Therefore
     519!--    user-specified values have to ''interpolate'' to the next lowest level
     520       IF ( number_stretch_level_start /= 0 ) THEN
     521          dz_stretch_level_start(1) = INT( (dz_stretch_level_start(1) -        &
     522                                            dz(1)/2.0) / dz(1) )               &
     523                                      * dz(1) + dz(1)/2.0
     524       ENDIF
     525       
     526       IF ( number_stretch_level_start > 1 ) THEN
     527          DO n = 2, number_stretch_level_start
     528             dz_stretch_level_start(n) = INT( dz_stretch_level_start(n) /      &
     529                                              dz(n) ) * dz(n)
     530          ENDDO
     531       ENDIF
     532       
     533       IF ( number_stretch_level_end /= 0 ) THEN
     534          DO n = 1, number_stretch_level_end
     535             dz_stretch_level_end(n) = INT( dz_stretch_level_end(n) /          &
     536                                            dz(n+1) ) * dz(n+1)
     537          ENDDO
     538       ENDIF
     539 
     540!
     541!--    Determine stretching factor if necessary
     542       IF ( number_stretch_level_end >= 1 ) THEN
     543          CALL calculate_stretching_factor( number_stretch_level_end )
     544       ENDIF
     545
    404546!
    405547!--    Grid for atmosphere with surface at z=0 (k=0, w-grid).
     548!--    First compute the u- and v-levels. In case of dirichlet bc for u and v
     549!--    the first u/v- and w-level (k=0) are defined at same height (z=0).
    406550!--    The second u-level (k=1) corresponds to the top of the
    407551!--    Prandtl-layer.
    408 
    409552       IF ( ibc_uv_b == 0 .OR. ibc_uv_b == 2 ) THEN
    410553          zu(0) = 0.0_wp
    411       !    zu(0) = - dz * 0.5_wp
    412554       ELSE
    413           zu(0) = - dz * 0.5_wp
    414        ENDIF
    415        zu(1) =   dz * 0.5_wp
    416 
    417        dz_stretch_level_index = nzt+1
    418        dz_stretched = dz
     555          zu(0) = - dz(1) * 0.5_wp
     556       ENDIF
     557         
     558       zu(1) =   dz(1) * 0.5_wp
     559       
     560!
     561!--    Determine u and v height levels considering the possibility of grid
     562!--    stretching in several heights.
     563       n = 1
     564       dz_stretch_level_start_index = nzt+1
     565       dz_stretch_level_end_index = nzt+1
     566       dz_stretched = dz(1)
     567
     568!--    The default value of dz_stretch_level_start is negative, thus the first
     569!--    condition is always true. Hence, the second condition is necessary.
    419570       DO  k = 2, nzt+1
    420           IF ( dz_stretch_level <= zu(k-1)  .AND.  dz_stretched < dz_max )  THEN
    421              dz_stretched = dz_stretched * dz_stretch_factor
    422              dz_stretched = MIN( dz_stretched, dz_max )
    423              IF ( dz_stretch_level_index == nzt+1 ) dz_stretch_level_index = k-1
    424           ENDIF
     571          IF ( dz_stretch_level_start(n) <= zu(k-1) .AND.                      &
     572               dz_stretch_level_start(n) /= -9999999.9_wp ) THEN
     573             dz_stretched = dz_stretched * dz_stretch_factor_array(n)
     574             
     575             IF ( dz(n) > dz(n+1) ) THEN
     576                dz_stretched = MAX( dz_stretched, dz(n+1) ) !Restrict dz_stretched to the user-specified (higher) dz
     577             ELSE
     578                dz_stretched = MIN( dz_stretched, dz(n+1) ) !Restrict dz_stretched to the user-specified (lower) dz
     579             ENDIF
     580             
     581             IF ( dz_stretch_level_start_index(n) == nzt+1 )                         &
     582             dz_stretch_level_start_index(n) = k-1
     583             
     584          ENDIF
     585         
    425586          zu(k) = zu(k-1) + dz_stretched
     587         
     588!
     589!--       Make sure that the stretching ends exactly at dz_stretch_level_end
     590          dz_level_end = ABS( zu(k) - dz_stretch_level_end(n) )
     591         
     592          IF ( dz_level_end  < dz(n+1)/3.0 ) THEN
     593             zu(k) = dz_stretch_level_end(n)
     594             dz_stretched = dz(n+1)
     595             dz_stretch_level_end_index(n) = k
     596             n = n + 1             
     597          ENDIF
    426598       ENDDO
    427599
     
    438610
    439611    ELSE
     612
     613!
     614!--    The stretching region has to be large enough to allow for a smooth
     615!--    transition between two different grid spacings
     616       DO n = 1, number_stretch_level_start
     617          min_dz_stretch_level_end(n) = dz_stretch_level_start(n) -            &
     618                                        4 * MAX( dz(n),dz(n+1) )
     619       ENDDO
     620       
     621       IF ( ANY( min_dz_stretch_level_end < dz_stretch_level_end ) ) THEN
     622             message_string= 'Eeach dz_stretch_level_end has to be less ' //   &
     623                             'than its corresponding value for &' //           &
     624                             'dz_stretch_level_start - 4*MAX(dz(n),dz(n+1)) '//&
     625                             'to allow for smooth grid stretching'
     626             CALL message( 'init_grid', 'PA0224', 1, 2, 0, 6, 0 )
     627       ENDIF
     628       
     629!
     630!--    Stretching must not be applied within the prandtl_layer
     631!--    (first two grid points). For the default case dz_stretch_level_start
     632!--    is negative. Therefore the absolut value is checked here.
     633       IF ( ANY( dz_stretch_level_start > dz(1) * 1.5_wp ) ) THEN
     634          WRITE( message_string, * ) 'Eeach dz_stretch_level_start has to be ',&
     635                                     'less than ', dz(1) * 1.5
     636             CALL message( 'init_grid', 'PA0226', 1, 2, 0, 6, 0 )
     637       ENDIF
     638
     639!
     640!--    The stretching has to start and end on a grid level. Therefore
     641!--    user-specified values have to ''interpolate'' to the next highest level
     642       IF ( number_stretch_level_start /= 0 ) THEN
     643          dz_stretch_level_start(1) = INT( (dz_stretch_level_start(1) +        &
     644                                            dz(1)/2.0) / dz(1) )               &
     645                                      * dz(1) - dz(1)/2.0
     646       ENDIF
     647       
     648       IF ( number_stretch_level_start > 1 ) THEN
     649          DO n = 2, number_stretch_level_start
     650             dz_stretch_level_start(n) = INT( dz_stretch_level_start(n) /      &
     651                                              dz(n) ) * dz(n)
     652          ENDDO
     653       ENDIF
     654       
     655       IF ( number_stretch_level_end /= 0 ) THEN
     656          DO n = 1, number_stretch_level_end
     657             dz_stretch_level_end(n) = INT( dz_stretch_level_end(n) /          &
     658                                            dz(n+1) ) * dz(n+1)
     659          ENDDO
     660       ENDIF
     661       
     662!
     663!--    Determine stretching factor if necessary
     664       IF ( number_stretch_level_end >= 1 ) THEN
     665          CALL calculate_stretching_factor( number_stretch_level_end )
     666       ENDIF
     667
    440668!
    441669!--    Grid for ocean with free water surface is at k=nzt (w-grid).
     
    444672!--    w-level are defined at same height, but staggered from the second level.
    445673!--    The second u-level (k=1) corresponds to the top of the Prandtl-layer.
    446        zu(nzt+1) =   dz * 0.5_wp
    447        zu(nzt)   = - dz * 0.5_wp
    448 
    449        dz_stretch_level_index = 0
    450        dz_stretched = dz
     674!--    z values are negative starting from z=0 (surface)
     675       zu(nzt+1) =   dz(1) * 0.5_wp
     676       zu(nzt)   = - dz(1) * 0.5_wp
     677
     678!
     679!--    Determine u and v height levels considering the possibility of grid
     680!--    stretching in several heights.
     681       n = 1
     682       dz_stretch_level_start_index = 0
     683       dz_stretch_level_end_index = 0
     684       dz_stretched = dz(1)
     685
    451686       DO  k = nzt-1, 0, -1
    452 !
    453 !--       The default value of dz_stretch_level is positive, thus the first
    454 !--       condition is always true. Hence, the second condition is necessary.
    455           IF ( dz_stretch_level >= zu(k+1)  .AND.  dz_stretch_level <= 0.0  &
    456                .AND.  dz_stretched < dz_max )  THEN
    457              dz_stretched = dz_stretched * dz_stretch_factor
    458              dz_stretched = MIN( dz_stretched, dz_max )
    459              IF ( dz_stretch_level_index == 0 ) dz_stretch_level_index = k+1
    460           ENDIF
     687         
     688          IF ( dz_stretch_level_start(n) >= zu(k+1) ) THEN
     689             dz_stretched = dz_stretched * dz_stretch_factor_array(n)
     690
     691             IF ( dz(n) > dz(n+1) ) THEN
     692                dz_stretched = MAX( dz_stretched, dz(n+1) ) !Restrict dz_stretched to the user-specified (higher) dz
     693             ELSE
     694                dz_stretched = MIN( dz_stretched, dz(n+1) ) !Restrict dz_stretched to the user-specified (lower) dz
     695             ENDIF
     696             
     697             IF ( dz_stretch_level_start_index(n) == 0 )                             &
     698             dz_stretch_level_start_index(n) = k+1
     699             
     700          ENDIF
     701         
    461702          zu(k) = zu(k+1) - dz_stretched
     703         
     704!
     705!--       Make sure that the stretching ends exactly at dz_stretch_level_end
     706          dz_level_end = ABS( zu(k) - dz_stretch_level_end(n) )
     707         
     708          IF ( dz_level_end  < dz(n+1)/3.0 ) THEN
     709             zu(k) = dz_stretch_level_end(n)
     710             dz_stretched = dz(n+1)
     711             dz_stretch_level_end_index(n) = k
     712             n = n + 1             
     713          ENDIF
    462714       ENDDO
    463 
     715       
    464716!
    465717!--    Compute the w-levels. They are always staggered half-way between the
     
    468720!--    same height. The top w-level (nzt+1) is not used but set for
    469721!--    consistency, since w and all scalar variables are defined up tp nzt+1.
    470        zw(nzt+1) = dz
     722       zw(nzt+1) = dz(1)
    471723       zw(nzt)   = 0.0_wp
    472724       DO  k = 0, nzt
     
    8551107 END SUBROUTINE init_grid
    8561108
     1109
     1110! Description:
     1111! -----------------------------------------------------------------------------!
     1112!> Calculation of the stretching factor through an iterative method. Ideas were
     1113!> taken from the paper "Regional stretched grid generation and its application
     1114!> to the NCAR RegCM (1999)". Normally, no analytic solution exists because the
     1115!> system of equations has two variables (r,l) but four requirements
     1116!> (l=integer, r=[0,88;1,2], Eq(6), Eq(5) starting from index j=1) which
     1117!> results into an overdetermined system.
     1118!------------------------------------------------------------------------------!
     1119 SUBROUTINE calculate_stretching_factor( number_end )
     1120 
     1121    USE control_parameters,                                                    &
     1122        ONLY:  dz, dz_stretch_factor, dz_stretch_factor_array,                 &   
     1123               dz_stretch_level_end, dz_stretch_level_start, message_string
     1124 
     1125    USE kinds
     1126   
     1127    IMPLICIT NONE
     1128   
     1129    INTEGER(iwp) ::  iterations  !< number of iterations until stretch_factor_lower/upper_limit is reached 
     1130    INTEGER(iwp) ::  l_rounded   !< after l_rounded grid levels dz(n) is strechted to dz(n+1) with stretch_factor_2
     1131    INTEGER(iwp) ::  n           !< loop variable for stretching
     1132   
     1133    INTEGER(iwp), INTENT(IN) ::  number_end !< number of user-specified end levels for stretching
     1134       
     1135    REAL(wp) ::  delta_l               !< absolute difference between l and l_rounded
     1136    REAL(wp) ::  delta_stretch_factor  !< absolute difference between stretch_factor_1 and stretch_factor_2
     1137    REAL(wp) ::  delta_total_new       !< sum of delta_l and delta_stretch_factor for the next iteration (should be as small as possible)
     1138    REAL(wp) ::  delta_total_old       !< sum of delta_l and delta_stretch_factor for the last iteration
     1139    REAL(wp) ::  distance              !< distance between dz_stretch_level_start and dz_stretch_level_end (stretching region)
     1140    REAL(wp) ::  l                     !< value that fulfil Eq. (5) in the paper mentioned above together with stretch_factor_1 exactly
     1141    REAL(wp) ::  numerator             !< numerator of the quotient
     1142    REAL(wp) ::  stretch_factor_1      !< stretching factor that fulfil Eq. (5) togehter with l exactly
     1143    REAL(wp) ::  stretch_factor_2      !< stretching factor that fulfil Eq. (6) togehter with l_rounded exactly
     1144   
     1145    REAL(wp), PARAMETER ::  stretch_factor_interval = 1.0E-06  !< interval for sampling possible stretching factors
     1146    REAL(wp), PARAMETER ::  stretch_factor_lower_limit = 0.88  !< lowest possible stretching factor
     1147    REAL(wp), PARAMETER ::  stretch_factor_upper_limit = 1.12  !< highest possible stretching factor
     1148 
     1149 
     1150       l = 0
     1151       DO  n = 1, number_end
     1152       
     1153          iterations = 1
     1154          stretch_factor_1 = 1.0
     1155          stretch_factor_2 = 1.0
     1156          delta_total_old = 1.0
     1157         
     1158          IF ( dz(n) > dz(n+1) ) THEN
     1159             DO WHILE ( stretch_factor_1 >= stretch_factor_lower_limit )
     1160               
     1161                stretch_factor_1 = 1.0 - iterations * stretch_factor_interval
     1162                distance = ABS( dz_stretch_level_end(n) -                   &
     1163                           dz_stretch_level_start(n) )
     1164                numerator = distance*stretch_factor_1/dz(n) +               &
     1165                            stretch_factor_1 - distance/dz(n)
     1166               
     1167                IF ( numerator > 0.0 ) THEN
     1168                   l = LOG( numerator ) / LOG( stretch_factor_1 ) - 1.0
     1169                   l_rounded = NINT( l )
     1170                   delta_l = ABS( l_rounded - l ) / l
     1171                ENDIF
     1172               
     1173                stretch_factor_2 = EXP( LOG( dz(n+1)/dz(n) ) / (l_rounded) )
     1174               
     1175                delta_stretch_factor = ABS( stretch_factor_1 -              &
     1176                                            stretch_factor_2 ) /            &
     1177                                       stretch_factor_2
     1178               
     1179                delta_total_new = delta_l + delta_stretch_factor
     1180
     1181!
     1182!--                stretch_factor_1 is taken to guarantee that the stretching
     1183!--                procedure ends as close as possible to dz_stretch_level_end.
     1184!--                stretch_factor_2 would guarantee that the stretched dz(n) is
     1185!--                equal to dz(n+1) after l_rounded grid levels.
     1186                IF (delta_total_new < delta_total_old) THEN
     1187                   dz_stretch_factor_array(n) = stretch_factor_1
     1188                   delta_total_old = delta_total_new
     1189                ENDIF
     1190               
     1191                iterations = iterations + 1
     1192               
     1193             ENDDO
     1194               
     1195          ELSEIF ( dz(n) < dz(n+1) ) THEN
     1196             DO WHILE ( stretch_factor_1 <= stretch_factor_upper_limit )
     1197                       
     1198                stretch_factor_1 = 1.0 + iterations * stretch_factor_interval
     1199                distance = ABS( dz_stretch_level_end(n) -                   &
     1200                           dz_stretch_level_start(n) )
     1201                numerator = distance*stretch_factor_1/dz(n) +               &
     1202                            stretch_factor_1 - distance/dz(n)
     1203               
     1204                l = LOG( numerator ) / LOG( stretch_factor_1 ) - 1.0
     1205                l_rounded = NINT( l )
     1206                delta_l = ABS( l_rounded - l ) / l
     1207               
     1208                stretch_factor_2 = EXP( LOG( dz(n+1)/dz(n) ) / (l_rounded) )
     1209
     1210                delta_stretch_factor = ABS( stretch_factor_1 -              &
     1211                                           stretch_factor_2 ) /             &
     1212                                           stretch_factor_2
     1213               
     1214                delta_total_new = delta_l + delta_stretch_factor
     1215               
     1216!
     1217!--                stretch_factor_1 is taken to guarantee that the stretching
     1218!--                procedure ends as close as possible to dz_stretch_level_end.
     1219!--                stretch_factor_2 would guarantee that the stretched dz(n) is
     1220!--                equal to dz(n+1) after l_rounded grid levels.
     1221                IF (delta_total_new < delta_total_old) THEN
     1222                   dz_stretch_factor_array(n) = stretch_factor_1
     1223                   delta_total_old = delta_total_new
     1224                ENDIF
     1225               
     1226                iterations = iterations + 1
     1227             ENDDO
     1228             
     1229          ELSE
     1230             message_string= 'Two adjacent values of dz must be different'
     1231             CALL message( 'init_grid', 'PA0228', 1, 2, 0, 6, 0 )
     1232             
     1233          ENDIF
     1234       ENDDO
     1235       
     1236 END SUBROUTINE calculate_stretching_factor
     1237 
     1238 
    8571239! Description:
    8581240! -----------------------------------------------------------------------------!
     
    17032085!--       Tunnel-wall depth
    17042086          IF ( tunnel_wall_depth == 9999999.9_wp )  THEN 
    1705              td = MAX ( dx, dy, dz )
     2087             td = MAX ( dx, dy, dz(1) )
    17062088          ELSE
    17072089             td = tunnel_wall_depth
  • palm/trunk/SOURCE/init_masks.f90

    r3049 r3065  
    2525! -----------------
    2626! $Id$
     27! dz_stretch_level was replaced by dz_stretch_level_start
     28!
     29! 3049 2018-05-29 13:52:36Z Giersch
    2730! Error messages revised
    2831!
     
    133136        ONLY:  constant_diffusion, cloud_droplets, cloud_physics,              &
    134137               data_output_masks, data_output_masks_user,                      &
    135                doav, doav_n, domask, domask_no, dz, dz_stretch_level, humidity,&
    136                mask, masks, mask_scale, mask_i,                                &
     138               doav, doav_n, domask, domask_no, dz, dz_stretch_level_start,    &
     139               humidity, mask, masks, mask_scale, mask_i,                      &
    137140               mask_i_global, mask_j, mask_j_global, mask_k, mask_k_global,    &
    138141               mask_loop, mask_size, mask_size_l, mask_start_l, mask_x,        &
     
    141144               microphysics_morrison, microphysics_seifert, passive_scalar,    &
    142145               ocean, varnamelength
     146               
    143147
    144148    USE grid_variables,                                                        &
     
    490494       CALL set_mask_locations( 1, dx, 'dx', nx, 'nx', nxl, nxr )
    491495       CALL set_mask_locations( 2, dy, 'dy', ny, 'ny', nys, nyn )
    492        CALL set_mask_locations( 3, dz, 'dz', nz, 'nz', nzb, nzt )
     496       CALL set_mask_locations( 3, dz(1), 'dz', nz, 'nz', nzb, nzt )
    493497!
    494498!--    Set global masks along all three dimensions (required by
     
    727731!--          The following line assumes a constant vertical grid spacing within
    728732!--          the vertical mask range; it fails for vertical grid stretching.
    729 !--          Maybe revise later. Issue warning but continue execution.
     733!--          Maybe revise later. Issue warning but continue execution. ABS(...)
     734!--          within the IF statement is necessary because the default value of
     735!--          dz_stretch_level_start is -9999999.9_wp.
    730736             loop_stride = NINT( mask_loop(mid,dim,3) * mask_scale(dim) * ddxyz )
    731737
    732              IF ( mask_loop(mid,dim,2) * mask_scale(dim) > dz_stretch_level )  &
    733                   THEN
     738             IF ( mask_loop(mid,dim,2) * mask_scale(dim) >                     &
     739                  ABS( dz_stretch_level_start(1) ) )  THEN
    734740                WRITE ( message_string, '(A,I3,A,I1,A,F9.3,A,F8.2,3A)' )       &
    735741                     'mask_loop(',mid,',',dim,',2)=', mask_loop(mid,dim,2),    &
    736                      ' exceeds dz_stretch_level=',dz_stretch_level,            &
     742                     ' exceeds dz_stretch_level=',dz_stretch_level_start(1),   &
    737743                     '.&Vertical mask locations will not ',                    &
    738744                     'match the desired heights within the stretching ',       &
  • palm/trunk/SOURCE/lpm_advec.f90

    r2969 r3065  
    2525! -----------------
    2626! $Id$
     27! dz values were replaced by dzw or dz(1) to allow for right vertical stretching
     28!
     29! 2969 2018-04-13 11:55:09Z thiele
    2730! Bugfix in Interpolation indices.
    2831!
     
    158161
    159162    USE arrays_3d,                                                             &
    160         ONLY:  de_dx, de_dy, de_dz, diss, e, km, u, v, w, zu, zw
     163        ONLY:  de_dx, de_dy, de_dz, diss, dzw, e, km, u, v, w, zu, zw
    161164
    162165    USE cpulog
     
    409412                            + ( gg-cc ) * u(k+1,j+1,i) + ( gg-dd ) *           &
    410413                            u(k+1,j+1,i+1) ) / ( 3.0_wp * gg ) - u_gtrans
    411                 u_int(n) = u_int_l + ( zv(n) - zu(k) ) / dz *                  &
     414                u_int(n) = u_int_l + ( zv(n) - zu(k) ) / dzw(k) *               &
    412415                           ( u_int_u - u_int_l )
    413416             ENDIF
     
    503506                          + ( gg-cc ) * v(k+1,j+1,i) + ( gg-dd ) * v(k+1,j+1,i+1) &
    504507                          ) / ( 3.0_wp * gg ) - v_gtrans
    505                 v_int(n) = v_int_l + ( zv(n) - zu(k) ) / dz * &
     508                v_int(n) = v_int_l + ( zv(n) - zu(k) ) / dzw(k) *              &
    506509                                  ( v_int_u - v_int_l )
    507510             ENDIF
     
    540543                            ( gg-dd ) * w(k+1,j+1,i+1) &
    541544                          ) / ( 3.0_wp * gg )
    542                 w_int(n) = w_int_l + ( zv(n) - zw(k) ) / dz * &
     545                w_int(n) = w_int_l + ( zv(n) - zw(k) ) / dzw(k) *              &
    543546                           ( w_int_u - w_int_l )
    544547             ENDIF
     
    611614                               ( gg - dd ) * e(k+1,j+1,i+1) &
    612615                            ) / ( 3.0_wp * gg )
    613                    e_int(n) = e_int_l + ( zv(n) - zu(k) ) / dz * &
     616                   e_int(n) = e_int_l + ( zv(n) - zu(k) ) / dzw(k) *            &
    614617                                     ( e_int_u - e_int_l )
    615618                ENDIF
     
    637640                                   ( gg - dd ) * de_dx(k+1,j+1,i+1) &
    638641                                  ) / ( 3.0_wp * gg )
    639                    de_dx_int(n) = de_dx_int_l + ( zv(n) - zu(k) ) / dz * &
     642                   de_dx_int(n) = de_dx_int_l + ( zv(n) - zu(k) ) / dzw(k) *    &
    640643                                              ( de_dx_int_u - de_dx_int_l )
    641644                ENDIF
     
    655658                                   ( gg - dd ) * de_dy(k+1,j+1,i+1) &
    656659                                  ) / ( 3.0_wp * gg )
    657                       de_dy_int(n) = de_dy_int_l + ( zv(n) - zu(k) ) / dz * &
     660                      de_dy_int(n) = de_dy_int_l + ( zv(n) - zu(k) ) / dzw(k) * &
    658661                                                 ( de_dy_int_u - de_dy_int_l )
    659662                ENDIF
     
    661664!
    662665!--             Interpolate the TKE gradient along z
    663                 IF ( zv(n) < 0.5_wp * dz )  THEN
     666                IF ( zv(n) < 0.5_wp * dz(1) )  THEN
    664667                   de_dz_int(n) = 0.0_wp
    665668                ELSE
     
    678681                                      ( gg - dd ) * de_dz(k+1,j+1,i+1) &
    679682                                     ) / ( 3.0_wp * gg )
    680                       de_dz_int(n) = de_dz_int_l + ( zv(n) - zu(k) ) / dz * &
     683                      de_dz_int(n) = de_dz_int_l + ( zv(n) - zu(k) ) / dzw(k) * &
    681684                                                 ( de_dz_int_u - de_dz_int_l )
    682685                   ENDIF
     
    699702                                  ( gg - dd ) * diss(k+1,j+1,i+1) &
    700703                                 ) / ( 3.0_wp * gg )
    701                    diss_int(n) = diss_int_l + ( zv(n) - zu(k) ) / dz * &
     704                   diss_int(n) = diss_int_l + ( zv(n) - zu(k) ) / dzw(k) *      &
    702705                                            ( diss_int_u - diss_int_l )
    703706                ENDIF
  • palm/trunk/SOURCE/lpm_exchange_horiz.f90

    r3049 r3065  
    2525! -----------------
    2626! $Id$
     27! dz was replaced by dz(1) to allow for right vertical stretching
     28!
     29! 3049 2018-05-29 13:52:36Z Giersch
    2730! Error messages revised
    2831!
     
    910913       ip = particle_array(n)%x * ddx
    911914       jp = particle_array(n)%y * ddy
    912        kp = particle_array(n)%z / dz + 1 + offset_ocean_nzt
     915       kp = particle_array(n)%z / dz(1) + 1 + offset_ocean_nzt
    913916!
    914917!--    In case of grid stretching a particle might be above or below the 
  • palm/trunk/SOURCE/lpm_init.f90

    r3049 r3065  
    2525! -----------------
    2626! $Id$
     27! dz was replaced by dzw or dz(1) to allow for right vertical stretching
     28!
     29! 3049 2018-05-29 13:52:36Z Giersch
    2730! Error messages revised
    2831!
     
    794797                               ip = tmp_particle%x * ddx
    795798                               jp = tmp_particle%y * ddy
    796                                kp = tmp_particle%z / dz + 1 + offset_ocean_nzt                               
     799                               kp = tmp_particle%z / dz(1) + 1 + offset_ocean_nzt                               
    797800                               DO WHILE( zw(kp) < tmp_particle%z )
    798801                                  kp = kp + 1
     
    969972                                     pdz(particles(n)%group)
    970973                      particles(n)%z = particles(n)%z +                        &
    971                               MERGE( rand_contr, SIGN( dz, rand_contr ),       &
    972                                      ABS( rand_contr ) < dz                    &
     974                              MERGE( rand_contr, SIGN( dzw(kp), rand_contr ),  &
     975                                     ABS( rand_contr ) < dzw(kp)               &
    973976                                   )
    974977                   ENDIF
     
    987990                   i = particles(n)%x * ddx
    988991                   j = particles(n)%y * ddy
    989                    k = particles(n)%z / dz + 1 + offset_ocean_nzt
     992                   k = particles(n)%z / dz(1) + 1 + offset_ocean_nzt
    990993                   DO WHILE( zw(k) < particles(n)%z )
    991994                      k = k + 1
  • palm/trunk/SOURCE/lpm_set_attributes.f90

    r2718 r3065  
    2525! -----------------
    2626! $Id$
     27! dz was replaced by dzw to allow for right vertical stretching
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    8083
    8184    USE arrays_3d,                                                             &
    82         ONLY:  pt, u, v, w, zu, zw
     85        ONLY:  dzw, pt, u, v, w, zu, zw
    8386
    8487    USE control_parameters,                                                    &
    85         ONLY:  u_gtrans, v_gtrans, dz
     88        ONLY:  u_gtrans, v_gtrans
    8689
    8790    USE cpulog,                                                                &
     
    207210                                     ( gg - dd ) * u(k+1,j+1,i+1)                &
    208211                                   ) / ( 3.0_wp * gg ) - u_gtrans
    209                          u_int(n) = u_int_l + ( zv(n) - zu(k) ) / dz *         &
     212                         u_int(n) = u_int_l + ( zv(n) - zu(k) ) / dzw(k) *      &
    210213                                           ( u_int_u - u_int_l )
    211214                      ENDIF
     
    240243                                      ( gg - dd ) * v(k+1,j+1,i+1)                &
    241244                                    ) / ( 3.0_wp * gg ) - v_gtrans
    242                          v_int(n) = v_int_l + ( zv(n) - zu(k) ) / dz *         &
     245                         v_int(n) = v_int_l + ( zv(n) - zu(k) ) / dzw(k) *      &
    243246                                           ( v_int_u - v_int_l )
    244247                      ENDIF
     
    344347                                 ) / ( 3.0_wp * gg ) - sums(k,4)
    345348
    346                       pt_int = pt_int_l + ( zv(n) - zu(k) ) / dz *    &
     349                      pt_int = pt_int_l + ( zv(n) - zu(k) ) / dzw(k) *          &
    347350                                          ( pt_int_u - pt_int_l )
    348351
  • palm/trunk/SOURCE/modules.f90

    r3045 r3065  
    2525! -----------------
    2626! $Id$
     27! Variables concerning stretching introduced or revised
     28!
     29! 3045 2018-05-28 07:55:41Z Giersch
    2730! z_max_do2d removed
    2831!
     
    11481151    INTEGER(iwp) ::  dp_level_ind_b = 0                !< lowest grid index for external pressure gradient forcing
    11491152    INTEGER(iwp) ::  dvrp_filecount = 0                !< parameter for dvr visualization software
    1150     INTEGER(iwp) ::  dz_stretch_level_index            !< vertical grid level index above which the vertical grid spacing is stretched
    11511153    INTEGER(iwp) ::  ensemble_member_nr = 0            !< namelist parameter
    11521154    INTEGER(iwp) ::  gamma_mg                          !< switch for steering the multigrid cycle: 1: v-cycle, 2: w-cycle
     
    11901192    INTEGER(iwp) ::  num_var_fl                        !< number of sampling/output variables in virtual flight measurements
    11911193    INTEGER(iwp) ::  num_var_fl_user=0                 !< number of user-defined sampling/output variables in virtual flight measurements
     1194    INTEGER(iwp) ::  number_stretch_level_start        !< number of user-specified start levels for stretching
     1195    INTEGER(iwp) ::  number_stretch_level_end          !< number of user-specified end levels for stretching
    11921196    INTEGER(iwp) ::  nz_do3d = -9999                   !< namelist parameter
    11931197    INTEGER(iwp) ::  prt_time_count = 0                !< number of output intervals for particle data output
     
    11991203    INTEGER(iwp) ::  timestep_count = 0                !< number of timesteps carried out since the beginning of the initial run
    12001204    INTEGER(iwp) ::  y_shift = 0                       !< namelist parameter
     1205   
    12011206    INTEGER(iwp) ::  dist_nxl(0:1)                               !< left boundary of disturbance region
    12021207    INTEGER(iwp) ::  dist_nxr(0:1)                               !< right boundary of disturbance region
     
    12111216    INTEGER(iwp) ::  domask_no(max_masks,0:1) = 0                !< number of masked output quantities
    12121217    INTEGER(iwp) ::  domask_time_count(max_masks,0:1)            !< number of output intervals for masked data
     1218    INTEGER(iwp) ::  dz_stretch_level_end_index(9)               !< vertical grid level index until which the vertical grid spacing is stretched
     1219    INTEGER(iwp) ::  dz_stretch_level_start_index(9)             !< vertical grid level index above which the vertical grid spacing is stretched
    12131220    INTEGER(iwp) ::  mask_size(max_masks,3) = -1                 !< size of mask array per mask and dimension (for netcdf output)
    12141221    INTEGER(iwp) ::  mask_size_l(max_masks,3) = -1               !< subdomain size of mask array per mask and dimension (for netcdf output)
     
    14161423    REAL(wp) ::  dt_spinup = 60.0_wp                           !< namelist parameter
    14171424    REAL(wp) ::  dt_3d = 1.0_wp                                !< time step
    1418     REAL(wp) ::  dz = -1.0_wp                                  !< namelist parameter
    1419     REAL(wp) ::  dz_max = 9999999.9_wp                         !< namelist parameter
     1425    REAL(wp) ::  dz_max = 1000.0_wp                            !< namelist parameter
    14201426    REAL(wp) ::  dz_stretch_factor = 1.08_wp                   !< namelist parameter
    1421     REAL(wp) ::  dz_stretch_level = 100000.0_wp                !< namelist parameter
     1427    REAL(wp) ::  dz_stretch_level = -9999999.9_wp              !< namelist parameter
    14221428    REAL(wp) ::  e_init = 0.0_wp                               !< namelist parameter
    14231429    REAL(wp) ::  e_min = 0.0_wp                                !< namelist parameter
     
    15291535    REAL(wp) ::  dpdxy(1:2) = 0.0_wp                               !< namelist parameter
    15301536    REAL(wp) ::  dt_domask(max_masks) = 9999999.9_wp               !< namelist parameter
     1537    REAL(wp) ::  dz(10) = -1.0_wp                                  !< namelist parameter
     1538    REAL(wp) ::  dz_stretch_level_start(9) = -9999999.9_wp         !< namelist parameter
     1539    REAL(wp) ::  dz_stretch_level_end(9) = 9999999.9_wp            !< namelist parameter
     1540    REAL(wp) ::  dz_stretch_factor_array(9) = 1.08_wp              !< namelist parameter
    15311541    REAL(wp) ::  mask_scale(3)                                     !< collective array for mask_scale_x/y/z
    15321542    REAL(wp) ::  pt_vertical_gradient(10) = 0.0_wp                 !< namelist parameter
  • palm/trunk/SOURCE/parin.f90

    r3049 r3065  
    2525! -----------------
    2626! $Id$
     27! New initialization parameters added
     28!
     29! 3049 2018-05-29 13:52:36Z Giersch
    2730! Error messages revised
    2831!
     
    509512             dp_external, dp_level_b, dp_smooth, dpdxy, dry_aerosol_radius,    &
    510513             dt, dt_pr_1d, dt_run_control_1d, dt_spinup, dx, dy, dz, dz_max,   &
    511              dz_stretch_factor, dz_stretch_level, end_time_1d,                 &
    512              ensemble_member_nr, e_init, e_min, fft_method,                    &
    513              flux_input_mode, flux_output_mode, forcing,                       &
     514             dz_stretch_factor, dz_stretch_level, dz_stretch_level_start,      &
     515             dz_stretch_level_end, end_time_1d, ensemble_member_nr, e_init,    &
     516             e_min, fft_method, flux_input_mode, flux_output_mode, forcing,    &
    514517             galilei_transformation, humidity,                                 &
    515518             inflow_damping_height, inflow_damping_width,                      &
     
    579582             dp_external, dp_level_b, dp_smooth, dpdxy, dry_aerosol_radius,    &
    580583             dt, dt_pr_1d, dt_run_control_1d, dt_spinup, dx, dy, dz, dz_max,   &
    581              dz_stretch_factor, dz_stretch_level, end_time_1d,                 &
    582              ensemble_member_nr, e_init, e_min, fft_method,                    &
    583              flux_input_mode, flux_output_mode, forcing,                       &
     584             dz_stretch_factor, dz_stretch_level, dz_stretch_level_start,      &
     585             dz_stretch_level_end, end_time_1d, ensemble_member_nr, e_init,    &
     586             e_min, fft_method, flux_input_mode, flux_output_mode, forcing,    &
    584587             galilei_transformation, humidity,                                 &
    585588             inflow_damping_height, inflow_damping_width,                      &
  • palm/trunk/SOURCE/plant_canopy_model_mod.f90

    r3049 r3065  
    2525! -----------------
    2626! $Id$
     27! dz was replaced by the help of zw to allow for vertical stretching
     28!
     29! 3049 2018-05-29 13:52:36Z Giersch
    2730! Error messages revised
    2831!
     
    550553
    551554       USE control_parameters,                                                 &
    552            ONLY: dz, passive_scalar
     555           ONLY: passive_scalar
    553556
    554557
     
    568571       REAL(wp) ::  canopy_height       !< canopy height (in m)
    569572       
    570        canopy_height = pch_index * dz
     573       canopy_height = zw(pch_index)
    571574
    572575       WRITE ( io, 1 )  canopy_mode, canopy_height, pch_index,                 &
     
    668671
    669672       USE control_parameters,                                                 &
    670            ONLY: dz, humidity, io_blocks, io_group, message_string, ocean,     &
     673           ONLY: humidity, io_blocks, io_group, message_string, ocean,         &
    671674                 passive_scalar, urban_surface
    672675
     
    761764!--       Use beta function for lad-profile construction
    762765          int_bpdf = 0.0_wp
    763           canopy_height = pch_index * dz
     766          canopy_height = zw(pch_index)
    764767
    765768          DO k = 0, pch_index
  • palm/trunk/SOURCE/pmc_interface_mod.f90

    r3049 r3065  
    2525! -----------------
    2626! $Id$
     27! dz was replaced by dz(1)
     28!
     29! 3049 2018-05-29 13:52:36Z Giersch
    2730! Error messages revised
    2831!
     
    913916          parent_grid_info_real(5) = lower_left_coord_x + ( nx + 1 ) * dx
    914917          parent_grid_info_real(6) = lower_left_coord_y + ( ny + 1 ) * dy
    915           parent_grid_info_real(7) = dz
     918          parent_grid_info_real(7) = dz(1)
    916919
    917920          parent_grid_info_int(1)  = nx
     
    13181321       fval(3) = dx
    13191322       fval(4) = dy
    1320        fval(5) = dz
     1323       fval(5) = dz(1)
    13211324
    13221325       IF ( myid == 0 )  THEN
  • palm/trunk/SOURCE/radiation_model_mod.f90

    r3049 r3065  
    2828! -----------------
    2929! $Id$
     30! dz was replaced by dz(1), error message concerning vertical stretching was
     31! added 
     32!
     33! 3049 2018-05-29 13:52:36Z Giersch
    3034! Error messages revised
    3135!
     
    330334! ------------
    331335!> Radiation models and interfaces
     336!> @todo Replace dz(1) appropriatly to account for grid stretching
    332337!> @todo move variable definitions used in radiation_init only to the subroutine
    333338!>       as they are no longer required after initialization.
     
    43714376     IF ( plant_canopy )  THEN
    43724377         pchf_prep(:) = r_d * (hyp(nzub:nzut) / 100000.0_wp)**0.286_wp &
    4373                      / (cp * hyp(nzub:nzut) * dx*dy*dz) !< equals to 1 / (rho * c_p * Vbox * T)
     4378                     / (cp * hyp(nzub:nzut) * dx*dy*dz(1)) !< equals to 1 / (rho * c_p * Vbox * T)
    43744379         pctf_prep(:) = r_d * (hyp(nzub:nzut) / 100000.0_wp)**0.286_wp &
    43754380                     / (l_v * hyp(nzub:nzut) * dx*dy*dz)
     
    43964401!--         each dimension, so that this new direction vector will allow us
    43974402!--         to traverse the ray path within grid coordinates directly
    4398          sunorig_grid = (/ sunorig(1)/dz, sunorig(2)/dy, sunorig(3)/dx /)
     4403         sunorig_grid = (/ sunorig(1)/dz(1), sunorig(2)/dy, sunorig(3)/dx /)
    43994404!--         sunorig_grid = sunorig_grid / norm2(sunorig_grid)
    44004405         sunorig_grid = sunorig_grid / SQRT(SUM(sunorig_grid**2))
     
    44034408!--            precompute effective box depth with prototype Leaf Area Density
    44044409            pc_box_dimshift = MAXLOC(ABS(sunorig), 1) - 1
    4405             CALL box_absorb(CSHIFT((/dz,dy,dx/), pc_box_dimshift),      &
     4410            CALL box_absorb(CSHIFT((/dz(1),dy,dx/), pc_box_dimshift),      &
    44064411                                60, prototype_lad,                          &
    44074412                                CSHIFT(ABS(sunorig), pc_box_dimshift),      &
     
    48014806        IF ( idir(d) == 0 ) facearea(d) = facearea(d) * dx
    48024807        IF ( jdir(d) == 0 ) facearea(d) = facearea(d) * dy
    4803         IF ( kdir(d) == 0 ) facearea(d) = facearea(d) * dz
     4808        IF ( kdir(d) == 0 ) facearea(d) = facearea(d) * dz(1)
    48044809     ENDDO
    48054810!
     
    50145019    SUBROUTINE radiation_interaction_init
    50155020
     5021       USE control_parameters,                                                 &
     5022           ONLY:  dz_stretch_level_start
     5023           
    50165024       USE netcdf_data_input_mod,                                              &
    50175025           ONLY:  leaf_area_density_f
     
    50955103       nzpt = nzptl
    50965104#endif
     5105!
     5106!--    Stretching (non-uniform grid spacing) is not considered in the radiation
     5107!--    model. Therefore, vertical stretching has to be applied above the area
     5108!--    where the parts of the radiation model which assume constant grid spacing
     5109!--    are active. ABS (...) is required because the default value of
     5110!--    dz_stretch_level_start is -9999999.9_wp (negative).
     5111       IF ( ABS( dz_stretch_level_start(1) ) <= zw(nzut) ) THEN
     5112          WRITE( message_string, * ) 'The lowest level where vertical ',       &
     5113                                     'stretching is applied have to be ',      &
     5114                                     'greater than ', zw(nzut)
     5115          CALL message( 'init_grid', 'PA0496', 1, 2, 0, 6, 0 )
     5116       ENDIF
    50975117!
    50985118!--    global number of urban and plant layers
     
    53905410            IF ( idir(d) == 0 ) facearea(d) = facearea(d) * dx
    53915411            IF ( jdir(d) == 0 ) facearea(d) = facearea(d) * dy
    5392             IF ( kdir(d) == 0 ) facearea(d) = facearea(d) * dz
     5412            IF ( kdir(d) == 0 ) facearea(d) = facearea(d) * dz(1)
    53935413        ENDDO
    53945414
     
    55125532
    55135533!--                 unit vector source -> target
    5514                     uv = (/ (ta(1)-sa(1))*dz, (ta(2)-sa(2))*dy, (ta(3)-sa(3))*dx /)
     5534                    uv = (/ (ta(1)-sa(1))*dz(1), (ta(2)-sa(2))*dy, (ta(3)-sa(3))*dx /)
    55155535                    sqdist = SUM(uv(:)**2)
    55165536                    uv = uv / SQRT(sqdist)
     
    56775697
    56785698!--               unit vector source -> target
    5679                   uv = (/ (ta(1)-sa(1))*dz, (ta(2)-sa(2))*dy, (ta(3)-sa(3))*dx /)
     5699                  uv = (/ (ta(1)-sa(1))*dz(1), (ta(2)-sa(2))*dy, (ta(3)-sa(3))*dx /)
    56805700                  sqdist = SUM(uv(:)**2)
    56815701                  uv = uv / SQRT(sqdist)
     
    61276147        ENDIF
    61286148        uvect(:) = delta(:) / distance
    6129         realdist = SQRT(SUM( (uvect(:)*(/dz,dy,dx/))**2 ))
     6149        realdist = SQRT(SUM( (uvect(:)*(/dz(1),dy,dx/))**2 ))
    61306150
    61316151        lastdist = 0._wp
     
    65096529                     zb1 = CEILING(zexit * zsgn - .5_wp) - 1  ! because it must be smaller than exit
    65106530                     nz = MAX(zb1 - zb0 + 3, 2)
    6511                      rt2_dist(nz) = SQRT(((zexit-zorig)*dz)**2 + dxxyy)
     6531                     rt2_dist(nz) = SQRT(((zexit-zorig)*dz(1))**2 + dxxyy)
    65126532                     qdist = rt2_dist(nz) / (zexit-zorig)
    65136533                     rt2_dist(2:nz-1) = (/( ((REAL(l, wp) + .5_wp) * zsgn - zorig) * qdist , l = zb0, zb1 )/)
     
    65676587                  zb1 = CEILING(zexit * zsgn - .5_wp) - 1  ! because it must be smaller than exit
    65686588                  nz = MAX(zb1 - zb0 + 3, 2)
    6569                   rt2_dist(nz) = SQRT(((zexit-zorig)*dz)**2 + dxxyy)
     6589                  rt2_dist(nz) = SQRT(((zexit-zorig)*dz(1))**2 + dxxyy)
    65706590                  qdist = rt2_dist(nz) / (zexit-zorig)
    65716591                  rt2_dist(2:nz-1) = (/( ((REAL(l, wp) + .5_wp) * zsgn - zorig) * qdist , l = zb0, zb1 )/)
  • palm/trunk/SOURCE/read_restart_data_mod.f90

    r3056 r3065  
    2525! -----------------
    2626! $Id$
     27! New parameters concerning vertical grid stretching have been added
     28!
     29! 3056 2018-06-04 07:49:35Z Giersch
    2730! found variable has to be set to false inside overlap loop
    2831!
     
    407410             CASE ( 'dz_stretch_factor' )
    408411                READ ( 13 )  dz_stretch_factor
     412             CASE ( 'dz_stretch_factor_array' )
     413                READ ( 13 )  dz_stretch_factor_array
    409414             CASE ( 'dz_stretch_level' )
    410415                READ ( 13 )  dz_stretch_level
     416             CASE ( 'dz_stretch_level_end' )
     417                READ ( 13 )  dz_stretch_level_end
     418             CASE ( 'dz_stretch_level_start' )
     419                READ ( 13 )  dz_stretch_level_start
    411420             CASE ( 'e_min' )
    412421                READ ( 13 )  e_min
     
    12511260    THEN
    12521261       WRITE( message_string, * ) 'number of PEs or virtual PE-grid changed ', &
    1253                         'in restart run PE', myid, ' will read from files ', &
     1262                        'in restart run & PE', myid, ' will read from files ', &
    12541263                         file_list(1:files_to_be_opened)
    12551264       CALL message( 'rrd_local', 'PA0285', 0, 0, 0, 6, 0 )
  • palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90

    r3051 r3065  
    2525! -----------------
    2626! $Id$
     27! Error message related to vertical stretching has been added, dz was replaced
     28! by dz(1)
     29!
     30! 3051 2018-05-30 17:43:55Z suehring
    2731! Bugfix in calculation of initial Reynolds-stress tensor.
    2832!
     
    313317
    314318    USE control_parameters,                                                    &
    315         ONLY:  bc_lr, bc_ns, forcing, nest_domain, rans_mode, turbulent_inflow
     319        ONLY:  bc_lr, bc_ns, forcing, nest_domain, number_stretch_level_start, &
     320               rans_mode, turbulent_inflow
    316321
    317322    USE pmc_interface,                                                         &
     
    365370                              'is not allowed'
    366371          CALL message( 'stg_check_parameters', 'PA0039', 1, 2, 0, 6, 0 )
     372       ENDIF
     373       
     374       IF ( number_stretch_level_start > 0 )  THEN
     375          message_string = 'Using synthetic turbulence generator ' //          &
     376                           'in combination with stretching is not allowed'
     377          CALL message( 'stg_check_parameters', 'PA0420', 1, 2, 0, 6, 0 )
    367378       ENDIF
    368379
     
    607618
    608619!
    609 !--       Convert length scales from meter to number of grid points
     620!--       Convert length scales from meter to number of grid points. Attention:
     621!--       Does not work if grid stretching is used
    610622          nuy(k) = INT( luy * ddy )
    611           nuz(k) = INT( luz / dz  )
     623          nuz(k) = INT( luz / dz(1)  )
    612624          nvy(k) = INT( lvy * ddy )
    613           nvz(k) = INT( lvz / dz  )
     625          nvz(k) = INT( lvz / dz(1)  )
    614626          nwy(k) = INT( lwy * ddy )
    615           nwz(k) = INT( lwz / dz  )
     627          nwz(k) = INT( lwz / dz(1)  )
    616628!
    617629!--       Workaround, assume isotropic turbulence
  • palm/trunk/SOURCE/urban_surface_mod.f90

    r3049 r3065  
    2828! -----------------
    2929! $Id$
     30! Unused array dxdir was removed, dz was replaced by dzu to consider vertical
     31! grid stretching
     32!
     33! 3049 2018-05-29 13:52:36Z Giersch
    3034! Error messages revised
    3135!
     
    282286#if ! defined( __nopointer )
    283287    USE arrays_3d,                                                             &
    284         ONLY:  hyp, zu, pt, pt_1, pt_2, p, u, v, w, hyp, tend
     288        ONLY:  dzu, hyp, zu, pt, pt_1, pt_2, p, u, v, w, hyp, tend
    285289#endif
    286290
     
    292296   
    293297    USE control_parameters,                                                    &
    294         ONLY:  coupling_start_time, dz, topography, dt_3d,                     &
     298        ONLY:  coupling_start_time, topography, dt_3d,                         &
    295299               intermediate_timestep_count, initializing_actions,              &
    296300               intermediate_timestep_count_max, simulated_time, end_time,      &
     
    72077211        REAL(wp), DIMENSION(nzb:nzt)          :: exn                !< value of the Exner function in layers
    72087212       
    7209         REAL(wp), DIMENSION(0:4)              :: dxdir              !< surface normal direction gridbox length
    72107213        REAL(wp)                              :: dtime              !< simulated time of day (in UTC)
    72117214        INTEGER(iwp)                          :: dhour              !< simulated hour of day (in UTC)
     
    72137216
    72147217
    7215         dxdir = (/dz,dy,dy,dx,dx/)
    72167218#if ! defined( __nopointer )
    72177219        exn(nzb:nzt) = (hyp(nzb:nzt) / 100000.0_wp )**0.286_wp          !< Exner function
     
    76347636                         (dtime/3600.0_wp-REAL(dhour,wp))*aheatprof(k,dhour+1)
    76357637                 IF ( aheat(k,j,i) > 0.0_wp )  THEN
    7636                     pt(k,j,i) = pt(k,j,i) + aheat(k,j,i)*acoef*dt_3d/(exn(k)*rho_cp*dz)
     7638                    pt(k,j,i) = pt(k,j,i) + aheat(k,j,i)*acoef*dt_3d/(exn(k)*rho_cp*dzu(k))
    76377639                 ENDIF
    76387640              ENDIF
  • palm/trunk/SOURCE/user_init_grid.f90

    r2718 r3065  
    2525! -----------------
    2626! $Id$
     27! dz was replaced by dz(1)
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    108111!--       topography, bit is 1 for atmospheric grid point.
    109112!--       The following example shows how to prescribe sine-like topography
    110 !--       along x-direction with amplitude of 10 * dz and wavelength 10 * dy.
     113!--       along x-direction with amplitude of 10 * dz(1) and wavelength 10 * dy.
    111114!           DO  i = nxlg, nxrg
    112 !              h_topo = 10.0_wp * dz * (SIN(3.14_wp*0.5_wp)*i*dx / ( 5.0_wp * dy ) )**2
     115!              h_topo = 10.0_wp * dz(1) * (SIN(3.14_wp*0.5_wp)*i*dx / ( 5.0_wp * dy ) )**2
    113116!
    114117!              k_topo = MINLOC( ABS( zw - h_topo ), 1 ) - 1
  • palm/trunk/SOURCE/vertical_nesting_mod.f90

    r3049 r3065  
    2626! -----------------
    2727! $Id$
     28! dz was replaced by dz(1), error messages related to vertical grid stretching
     29! have been added
     30!
     31! 3049 2018-05-29 13:52:36Z Giersch
    2832! Error messages revised
    2933!
     
    7377!> after spin-up of the CG
    7478!>
     79!> @todo Replace dz(1) appropriatly to account for grid stretching
    7580!> @todo Ensure that code can be compiled for serial and parallel mode. Please
    7681!>       check the placement of the directive "__parallel".
     
    36393644
    36403645       USE control_parameters,                                                    &
    3641            ONLY:  coupling_mode, coupling_mode_remote, coupling_topology, dz
     3646           ONLY:  coupling_mode, coupling_mode_remote, coupling_topology, dz,     &
     3647                  dz_stretch_level_start, message_string
    36423648
    36433649       USE grid_variables,                                                        &
     
    36783684           dxc = dx
    36793685           dyc = dy
    3680            dzc = dz
     3686           dzc = dz(1)
    36813687           cg_nprocs = numprocs
    36823688
    36833689           IF ( myid == 0 )  THEN
    36843690
    3685                CALL MPI_SEND( nxc, 1, MPI_INTEGER  , numprocs, 1,  comm_inter,  &
     3691               CALL MPI_SEND( nxc, 1, MPI_INTEGER  , numprocs, 1,  comm_inter, &
    36863692                   ierr )
    3687                CALL MPI_SEND( nyc, 1, MPI_INTEGER  , numprocs, 2,  comm_inter,  &
     3693               CALL MPI_SEND( nyc, 1, MPI_INTEGER  , numprocs, 2,  comm_inter, &
    36883694                   ierr )
    3689                CALL MPI_SEND( nzc, 1, MPI_INTEGER  , numprocs, 3,  comm_inter,  &
     3695               CALL MPI_SEND( nzc, 1, MPI_INTEGER  , numprocs, 3,  comm_inter, &
    36903696                   ierr )
    3691                CALL MPI_SEND( dxc, 1, MPI_REAL     , numprocs, 4,  comm_inter,  &
     3697               CALL MPI_SEND( dxc, 1, MPI_REAL     , numprocs, 4,  comm_inter, &
    36923698                   ierr )
    3693                CALL MPI_SEND( dyc, 1, MPI_REAL     , numprocs, 5,  comm_inter,  &
     3699               CALL MPI_SEND( dyc, 1, MPI_REAL     , numprocs, 5,  comm_inter, &
    36943700                   ierr )
    3695                CALL MPI_SEND( dzc, 1, MPI_REAL     , numprocs, 6,  comm_inter,  &
     3701               CALL MPI_SEND( dzc, 1, MPI_REAL     , numprocs, 6,  comm_inter, &
    36963702                   ierr )
    3697                CALL MPI_SEND( pdims, 2, MPI_INTEGER, numprocs, 7,  comm_inter,  &
     3703               CALL MPI_SEND( pdims, 2, MPI_INTEGER, numprocs, 7,  comm_inter, &
    36983704                   ierr )
    36993705               CALL MPI_SEND( cg_nprocs, 1, MPI_INTEGER, numprocs, 8,  comm_inter,  &
    37003706                   ierr )
    3701                CALL MPI_RECV( nxf, 1, MPI_INTEGER,   numprocs, 21, comm_inter,  &
     3707               CALL MPI_RECV( nxf, 1, MPI_INTEGER,   numprocs, 21, comm_inter, &
    37023708                   status, ierr )
    3703                CALL MPI_RECV( nyf, 1, MPI_INTEGER,   numprocs, 22, comm_inter,  &
     3709               CALL MPI_RECV( nyf, 1, MPI_INTEGER,   numprocs, 22, comm_inter, &
    37043710                   status, ierr )
    3705                CALL MPI_RECV( nzf, 1, MPI_INTEGER,   numprocs, 23, comm_inter,  &
     3711               CALL MPI_RECV( nzf, 1, MPI_INTEGER,   numprocs, 23, comm_inter, &
    37063712                   status, ierr )
    3707                CALL MPI_RECV( dxf, 1, MPI_REAL,      numprocs, 24, comm_inter,  &
     3713               CALL MPI_RECV( dxf, 1, MPI_REAL,      numprocs, 24, comm_inter, &
    37083714                   status, ierr )
    3709                CALL MPI_RECV( dyf, 1, MPI_REAL,      numprocs, 25, comm_inter,  &
     3715               CALL MPI_RECV( dyf, 1, MPI_REAL,      numprocs, 25, comm_inter, &
    37103716                   status, ierr )
    3711                CALL MPI_RECV( dzf, 1, MPI_REAL,      numprocs, 26, comm_inter,  &
     3717               CALL MPI_RECV( dzf, 1, MPI_REAL,      numprocs, 26, comm_inter, &
    37123718                   status, ierr )
    3713                CALL MPI_RECV( pdims_partner, 2, MPI_INTEGER,                     &
     3719               CALL MPI_RECV( pdims_partner, 2, MPI_INTEGER,                   &
    37143720                   numprocs, 27, comm_inter, status, ierr )
    3715                CALL MPI_RECV( fg_nprocs, 1, MPI_INTEGER,                     &
     3721               CALL MPI_RECV( fg_nprocs, 1, MPI_INTEGER,                       &
    37163722                   numprocs, 28, comm_inter, status, ierr )
    37173723           ENDIF
     
    37253731           CALL MPI_BCAST( pdims_partner, 2, MPI_INTEGER, 0, comm2d, ierr )
    37263732           CALL MPI_BCAST( fg_nprocs,  1, MPI_INTEGER, 0, comm2d, ierr )
     3733           
     3734!
     3735!--        Check if stretching is used within the nested domain. ABS(...) is
     3736!--        necessary because of the default value of -9999999.9_wp (negative)
     3737           IF ( ABS( dz_stretch_level_start(1) ) <= (nzf+1)*dzf )  THEN       
     3738               message_string = 'Stretching in the parent domain is '//        &
     3739                                'only allowed above the nested domain'
     3740               CALL message( 'vertical_nesting_mod', 'PA0497', 1, 2, 0, 6, 0 )
     3741           ENDIF
    37273742
    37283743       ELSEIF ( coupling_mode ==  'vnested_fine' )  THEN
     
    37333748           dxf = dx
    37343749           dyf = dy
    3735            dzf = dz
     3750           dzf = dz(1)
    37363751           fg_nprocs = numprocs
    37373752
     
    37743789
    37753790       ENDIF
    3776  
     3791       
    37773792       ngp_c = ( nxc+1 + 2 * nbgp ) * ( nyc+1 + 2 * nbgp )
    37783793       ngp_f = ( nxf+1 + 2 * nbgp ) * ( nyf+1 + 2 * nbgp )
     
    39103925 
    39113926#if defined( __parallel )
    3912           USE arrays_3d,                                                             &
     3927          USE arrays_3d,                                                       &
    39133928              ONLY:  zu, zw
    39143929             
    3915           USE control_parameters,                                                    &
    3916               ONLY:  coupling_mode
     3930          USE control_parameters,                                              &
     3931              ONLY:  coupling_mode, message_string, number_stretch_level_start
    39173932             
    3918           USE indices,                                                               &
     3933          USE indices,                                                         &
    39193934              ONLY:  nzt
    39203935         
     
    39253940          IMPLICIT NONE
    39263941     
    3927           !-- Allocate and Exchange zuc and zuf, zwc and zwf
     3942!
     3943!--       Allocate and Exchange zuc and zuf, zwc and zwf
    39283944          IF ( coupling_mode(1:8)  == 'vnested_' )  THEN
    39293945         
     
    39323948         
    39333949             IF ( coupling_mode ==  'vnested_crse' )  THEN
    3934          
    3935                    zuc = zu
    3936                    zwc = zw
     3950               
     3951                zuc = zu
     3952                zwc = zw
     3953               
    39373954                IF ( myid == 0 )  THEN
    39383955         
     
    39543971             ELSEIF ( coupling_mode ==  'vnested_fine' )  THEN
    39553972         
    3956                    zuf = zu
    3957                    zwf = zw
     3973!
     3974!--             Check if stretching is used within the nested domain
     3975                IF ( number_stretch_level_start > 0 )  THEN
     3976                   message_string = 'Stretching in the nested domain is not '//&
     3977                           'allowed'
     3978                   CALL message( 'vertical_nesting_mod', 'PA0498', 1, 2, 0, 6, 0 )
     3979                ENDIF
     3980               
     3981                zuf = zu
     3982                zwf = zw
     3983               
    39583984                IF ( myid == 0 )  THEN
    39593985         
  • palm/trunk/SOURCE/virtual_flight_mod.f90

    r3049 r3065  
    2525! -----------------
    2626! $Id$
     27! IF statement revised to consider new vertical grid stretching
     28!
     29! 3049 2018-05-29 13:52:36Z Giersch
    2730! Error messages revised
    2831!
     
    734737
    735738       USE control_parameters,                                                 &
    736            ONLY:  dz, dz_stretch_level   
     739           ONLY:  dz, dz_stretch_level_start   
    737740 
    738741      USE grid_variables,                                                     &
     
    777780!--    if flight position is above the vertical grid-stretching level.
    778781!--    fac=1 if variable is on scalar grid level, fac=0 for w-component.
    779        IF ( z_pos(l) < dz_stretch_level )  THEN
    780           k = ( z_pos(l) + fac * 0.5_wp * dz ) / dz
     782       IF ( z_pos(l) < dz_stretch_level_start(1) )  THEN
     783          k = ( z_pos(l) + fac * 0.5_wp * dz(1) ) / dz(1)
    781784       ELSE
    782785!
  • palm/trunk/SOURCE/wind_turbine_model_mod.f90

    r3049 r3065  
    2626! -----------------
    2727! $Id$
     28! dz was replaced by dz(1), error message concerning grid stretching has been
     29! introduced
     30!
     31! 3049 2018-05-29 13:52:36Z Giersch
    2832! Error messages revised
    2933!
     
    132136!> automatically).
    133137!>
     138!> @todo Replace dz(1) appropriatly to account for grid stretching
    134139!> @todo Revise code according to PALM Coding Standard
    135140!> @todo Implement ADM and ALM turbine models
     
    740745       delta_r_factor = segment_width
    741746       delta_t_factor = segment_length
    742        delta_r_init   = delta_r_factor * MIN( dx, dy, dz)
     747       delta_r_init   = delta_r_factor * MIN( dx, dy, dz(1))
    743748
    744749       DO inot = 1, nturbines
     
    937942
    938943   
     944       USE control_parameters,                                                 &
     945           ONLY:  dz_stretch_level_start
     946   
    939947       IMPLICIT NONE
    940948
     
    9921000                 pi**( 1.0_wp / 6.0_wp ) * eps_kernel
    9931001!
     1002!--    Stretching (non-uniform grid spacing) is not considered in the wind
     1003!--    turbine model. Therefore, vertical stretching has to be applied above
     1004!--    the area where the wtm is active. ABS (...) is required because the
     1005!--    default value of dz_stretch_level_start is -9999999.9_wp (negative).
     1006       IF ( ABS( dz_stretch_level_start(1) ) <= MAXVAL(rcz) + MAXVAL(rr) +     &
     1007       eps_min)  THEN
     1008          WRITE( message_string, * ) 'The lowest level where vertical ',       &
     1009                                     'stretching is applied &have to be ',     &
     1010                                     'greater than ', MAXVAL(rcz) +            &
     1011                                      MAXVAL(rr) + eps_min
     1012          CALL message( 'init_grid', 'PA0484', 1, 2, 0, 6, 0 )
     1013       ENDIF
     1014!
    9941015!--    Square of eps_min:
    9951016       eps_min2 = eps_min**2
     
    10221043          i_hub(inot) = INT(   rcx(inot)                 / dx )
    10231044          j_hub(inot) = INT( ( rcy(inot) + 0.5_wp * dy ) / dy )
    1024           k_hub(inot) = INT( ( rcz(inot) + 0.5_wp * dz ) / dz )
     1045          k_hub(inot) = INT( ( rcz(inot) + 0.5_wp * dz(1) ) / dz(1) )
    10251046
    10261047!
     
    10311052          i_smear(inot) = CEILING( ( rr(inot) + eps_min ) / dx )
    10321053          j_smear(inot) = CEILING( ( rr(inot) + eps_min ) / dy )
    1033           k_smear(inot) = CEILING( ( rr(inot) + eps_min ) / dz )
     1054          k_smear(inot) = CEILING( ( rr(inot) + eps_min ) / dz(1) )
    10341055       
    10351056       ENDDO
     
    11211142!--       surface. All points between z=0 and z=dz/s would already be contained
    11221143!--       in grid box 1.
    1123           index_nacb(inot) = INT( ( rcz(inot) - rnac(inot) ) / dz ) + 1
    1124           index_nact(inot) = INT( ( rcz(inot) + rnac(inot) ) / dz ) + 1
     1144          index_nacb(inot) = INT( ( rcz(inot) - rnac(inot) ) / dz(1) ) + 1
     1145          index_nact(inot) = INT( ( rcz(inot) + rnac(inot) ) / dz(1) ) + 1
    11251146
    11261147!
     
    11481169                            IF ( j == tower_s )  THEN
    11491170                               tow_cd_surf(k,j,i) = ( rcz(inot) -              &
    1150                                     ( k_hub(inot) * dz - 0.5_wp * dz ) )  *    & ! extension in z-direction
     1171                                    ( k_hub(inot) * dz(1) - 0.5_wp * dz(1) ) )*& ! extension in z-direction
    11511172                                  ( ( tower_s + 1.0_wp + 0.5_wp ) * dy    -    &
    11521173                                    ( rcy(inot) - 0.5_wp * dtow(inot) ) ) *    & ! extension in y-direction
     
    11541175                            ELSEIF ( j == tower_n )  THEN
    11551176                               tow_cd_surf(k,j,i) = ( rcz(inot)            -   &
    1156                                     ( k_hub(inot) * dz - 0.5_wp * dz ) )  *    & ! extension in z-direction
     1177                                    ( k_hub(inot) * dz(1) - 0.5_wp * dz(1) ) )*& ! extension in z-direction
    11571178                                  ( ( rcy(inot) + 0.5_wp * dtow(inot) )   -    &
    11581179                                    ( tower_n + 0.5_wp ) * dy )           *    & ! extension in y-direction
     
    11631184                            ELSE
    11641185                               tow_cd_surf(k,j,i) = ( rcz(inot) -              &
    1165                                     ( k_hub(inot) * dz - 0.5_wp * dz ) )  *    &
     1186                                    ( k_hub(inot) * dz(1) - 0.5_wp * dz(1) ) )*&
    11661187                                    dy * turb_cd_tower(inot)
    11671188                            ENDIF
     
    11691190!--                      tower lies completely within one grid box:
    11701191                         ELSE
    1171                             tow_cd_surf(k,j,i) = ( rcz(inot)                 - &
    1172                                        ( k_hub(inot) * dz - 0.5_wp * dz ) ) *  &
     1192                            tow_cd_surf(k,j,i) = ( rcz(inot) - ( k_hub(inot) * &
     1193                                       dz(1) - 0.5_wp * dz(1) ) ) *            &
    11731194                                       dtow(inot) * turb_cd_tower(inot)
    11741195                         ENDIF
     
    11821203!--                         leftmost and rightmost grid box:
    11831204                            IF ( j == tower_s )  THEN                         
    1184                                tow_cd_surf(k,j,i) = dz * (                     &
     1205                               tow_cd_surf(k,j,i) = dz(1) * (                  &
    11851206                                      ( tower_s + 1 + 0.5_wp ) * dy         -  &
    11861207                                      ( rcy(inot) - 0.5_wp * dtow(inot) )      &
    11871208                                                        ) * turb_cd_tower(inot)
    11881209                            ELSEIF ( j == tower_n )  THEN
    1189                                tow_cd_surf(k,j,i) = dz * (                     &
     1210                               tow_cd_surf(k,j,i) = dz(1) * (                  &
    11901211                                      ( rcy(inot) + 0.5_wp * dtow(inot) )   -  &
    11911212                                      ( tower_n + 0.5_wp ) * dy                &
     
    11951216!--                         (where tow_cd_surf = grid box area):
    11961217                            ELSE
    1197                                tow_cd_surf(k,j,i) = dz * dy * turb_cd_tower(inot)
     1218                               tow_cd_surf(k,j,i) = dz(1) * dy *               &
     1219                                                    turb_cd_tower(inot)
    11981220                            ENDIF
    11991221!
    12001222!--                         tower lies completely within one grid box:
    12011223                         ELSE
    1202                             tow_cd_surf(k,j,i) = dz * dtow(inot) *             &
     1224                            tow_cd_surf(k,j,i) = dz(1) * dtow(inot) *          &
    12031225                                                turb_cd_tower(inot)
    12041226                         ENDIF ! end if larger than grid box
     
    13021324       ENDDO   ! end of loop over turbines
    13031325
    1304        tow_cd_surf   = tow_cd_surf   / ( dx * dy * dz )      ! Normalize tower drag
    1305        nac_cd_surf = nac_cd_surf / ( dx * dy * dz )      ! Normalize nacelle drag
     1326       tow_cd_surf   = tow_cd_surf   / ( dx * dy * dz(1) )  ! Normalize tower drag
     1327       nac_cd_surf = nac_cd_surf / ( dx * dy * dz(1) )      ! Normalize nacelle drag
    13061328
    13071329       CALL wtm_read_blade_tables
     
    17461768                   ii =   rbx(ring,rseg) * ddx
    17471769                   jj = ( rby(ring,rseg) - 0.5_wp * dy ) * ddy
    1748                    kk = ( rbz(ring,rseg) + 0.5_wp * dz ) / dz
     1770                   kk = ( rbz(ring,rseg) + 0.5_wp * dz(1) ) / dz(1)
    17491771!
    17501772!--                Interpolate only if all required information is available on
     
    17761798
    17771799                         u_int_1_l(inot,ring,rseg) = u_int_l          +        &
    1778                                      ( rbz(ring,rseg) - zu(kk) ) / dz *        &
     1800                                     ( rbz(ring,rseg) - zu(kk) ) / dz(1) *     &
    17791801                                     ( u_int_u - u_int_l )
    17801802
     
    17911813                   ii = ( rbx(ring,rseg) - 0.5_wp * dx ) * ddx
    17921814                   jj =   rby(ring,rseg)                 * ddy
    1793                    kk = ( rbz(ring,rseg) + 0.5_wp * dz ) / dz
     1815                   kk = ( rbz(ring,rseg) + 0.5_wp * dz(1) ) / dz(1)
    17941816!
    17951817!--                Interpolate only if all required information is available on
     
    18211843
    18221844                         v_int_1_l(inot,ring,rseg) = v_int_l +                 &
    1823                                      ( rbz(ring,rseg) - zu(kk) ) / dz *        &
     1845                                     ( rbz(ring,rseg) - zu(kk) ) / dz(1) *     &
    18241846                                     ( v_int_u - v_int_l )
    18251847
     
    18361858                   ii = ( rbx(ring,rseg) - 0.5_wp * dx ) * ddx
    18371859                   jj = ( rby(ring,rseg) - 0.5_wp * dy ) * ddy
    1838                    kk =   rbz(ring,rseg)                 / dz
     1860                   kk =   rbz(ring,rseg)                 / dz(1)
    18391861!
    18401862!--                Interpolate only if all required information is available on
     
    18661888
    18671889                         w_int_1_l(inot,ring,rseg) = w_int_l +                 &
    1868                                      ( rbz(ring,rseg) - zw(kk) ) / dz *        &
     1890                                     ( rbz(ring,rseg) - zw(kk) ) / dz(1) *     &
    18691891                                     ( w_int_u - w_int_l )
    18701892                      ELSE
     
    22732295                            dist_u_3d = ( i * dx               - rbx(ring,rseg) )**2 + &
    22742296                                        ( j * dy + 0.5_wp * dy - rby(ring,rseg) )**2 + &
    2275                                         ( k * dz - 0.5_wp * dz - rbz(ring,rseg) )**2
     2297                                        ( k * dz(1) - 0.5_wp * dz(1) - rbz(ring,rseg) )**2
    22762298                            dist_v_3d = ( i * dx + 0.5_wp * dx - rbx(ring,rseg) )**2 + &
    22772299                                        ( j * dy               - rby(ring,rseg) )**2 + &
    2278                                         ( k * dz - 0.5_wp * dz - rbz(ring,rseg) )**2
     2300                                        ( k * dz(1) - 0.5_wp * dz(1) - rbz(ring,rseg) )**2
    22792301                            dist_w_3d = ( i * dx + 0.5_wp * dx - rbx(ring,rseg) )**2 + &
    22802302                                        ( j * dy + 0.5_wp * dy - rby(ring,rseg) )**2 + &
    2281                                         ( k * dz               - rbz(ring,rseg) )**2
     2303                                        ( k * dz(1)               - rbz(ring,rseg) )**2
    22822304
    22832305!
  • palm/trunk/SOURCE/write_restart_data_mod.f90

    r3004 r3065  
    2525! -----------------
    2626! $Id$
     27! New parameters concerning vertical grid stretching have been added
     28!
     29! 3004 2018-04-27 12:33:25Z Giersch
    2730! precipitation_rate_av removed
    2831!
     
    380383       CALL wrd_write_string( 'dz' )
    381384       WRITE ( 14 )  dz
    382 
     385       
    383386       CALL wrd_write_string( 'dz_max' )
    384387       WRITE ( 14 )  dz_max
     
    386389       CALL wrd_write_string( 'dz_stretch_factor' )
    387390       WRITE ( 14 )  dz_stretch_factor
    388 
     391       
     392       CALL wrd_write_string( 'dz_stretch_factor_array' )
     393       WRITE ( 14 )  dz_stretch_factor_array
     394       
    389395       CALL wrd_write_string( 'dz_stretch_level' )
    390396       WRITE ( 14 )  dz_stretch_level
    391397
     398       CALL wrd_write_string( 'dz_stretch_level_end' )
     399       WRITE ( 14 )  dz_stretch_level_end
     400       
     401       CALL wrd_write_string( 'dz_stretch_level_start' )
     402       WRITE ( 14 )  dz_stretch_level_start
     403       
    392404       CALL wrd_write_string( 'e_min' )
    393405       WRITE ( 14 )  e_min
Note: See TracChangeset for help on using the changeset viewer.