treatment dummy variables with  toxlevel (0,1,2) for control,low,high  basins 88
 (1X,I5,(T10,17F4.0))                                                            
         4
     1     0.  0.  1.  0.
     2     0.  0.  1.  0.
     3     0.  0.  1.  0.
     4     0.  0.  1.  0.
     5     0.  0.  1.  0.
     6     0.  0.  1.  0.
     7     0.  0.  1.  0.
     8     0.  0.  1.  0.
     9     0.  0.  1.  0.
    10     0.  0.  1.  0.
    11     0.  0.  1.  0.
    12     0.  0.  1.  0.
    13     0.  0.  1.  0.
    14     0.  0.  1.  0.
    15     0.  0.  1.  0.
    16     0.  0.  1.  0.
    17     0.  0.  1.  0.
    18     0.  0.  1.  0.
    19     0.  0.  1.  0.
    20     1.  0.  0.  1.
    21     0.  1.  0.  2.
    22     0.  0.  1.  0.
    23     1.  0.  0.  1.
    24     0.  1.  0.  2.
    25     0.  0.  1.  0.
    26     1.  0.  0.  1.
    27     0.  1.  0.  2.
    28     0.  0.  1.  0.
    29     1.  0.  0.  1.
    30     0.  1.  0.  2.
     0     0.  0.  0.  0.
low     high    control toxlevel
c_1     l_1     h_1     c_2     l_2     h_2     c_3     l_3     h_3     c_4     
l_4     h_4     c_5     l_5     h_5     c_6     l_6     h_6     c_7     l_7     
h_7     c_8     l_8     h_8     c_9     l_9     h_9     c_10    l_10    h_10    


Coding sample namen: c_7: control pont timepoint 7          

Before treatment : timepoint  1-6
After treatment  : timepoint  7-10