C  Instructions for running the program
C  ====================================
C
C  1) Go to line 220, change NPTS to the number of data points to simulate. For each
C     minute of operation, 60 points are generated.
C
C  2) Go to line 226, change SSPTS to the number of data points to simulate in steady
C     state operation before implementing the disturbance.
C
C  3) Go to line 367, implement any of the 21 programmed disturbances. For example, to
C     implement disturbance 2, type IDV(2)=1 .
C
C  4) The program will generate 15 output files and all data are recorded every
C     180 seconds, see Table 1 for details.  The default path is the home directory.
C     To change the file name and path, modify lines 346-360 accordingly.  
C     To overwrite the files that already existed, change STATUS='new' to
C     STATUS='old' from lines 346-360.
C              
C
C                     Table 1: Content of the output files
C
C      File Name                                Content
C      ---------                                -------
C    TE_data_inc.dat                             Time (in seconds)
C    TE_data_mv1.dat           Measurements for manipulated variables 1 to 4
C    TE_data_mv2.dat           Measurements for manipulated variables 5 to 8
C    TE_data_mv3.dat           Measurements for manipulated variables 9 to 12
C    TE_data_me01.dat           Measurements for measurement variables 1 to 4
C    TE_data_me02.dat           Measurements for measurement variables 5 to 8
C    TE_data_me03.dat           Measurements for measurement variables 9 to 12
C    TE_data_me04.dat           Measurements for measurement variables 13 to 16
C    TE_data_me05.dat           Measurements for measurement variables 17 to 20
C    TE_data_me06.dat           Measurements for measurement variables 21 to 24
C    TE_data_me07.dat           Measurements for measurement variables 25 to 28
C    TE_data_me08.dat           Measurements for measurement variables 29 to 32
C    TE_data_me09.dat           Measurements for measurement variables 33 to 36
C    TE_data_me10.dat           Measurements for measurement variables 37 to 40
C    TE_data_me11.dat           Measurements for measurement variable 41
C
C  5) To ensure the randomness of the measurement noises, the random number
C     G in the sub program (teprob.f, line 1187) has to be changed each time before
C     running 'temain_mod.f'.
C
C  6) Save the changes in 'temain_mod.f' and 'teprob.f' and compile the program in
C      unix by typing
C      f77 temain_mod.f teprob.f
C
C  7) Run the program by typing
C      a.out