home *** CD-ROM | disk | FTP | other *** search
/ Simtel MSDOS 1992 September / Simtel20_Sept92.cdr / msdos / spredsht / prof2d.arc / EXAMPLE1.PR2 < prev    next >
Text File  |  1989-02-26  |  3KB  |  107 lines
  1. : Prof2d, The 2D Scientific Spreadsheet 
  2. : by Gertjan L. Ouwerling, Electrical Materials Lab, Delft University
  3. : -------------------------------------------------------------------
  4. : Example 1: Reading, viewing and writing data
  5. :
  6. : Execute this example by typing:  prof2d example1.pr2 <enter>
  7. : 
  8. : NOTE: some steps may be slow without a co-processor!
  9. :       Graphical pictures are terminated by giving <enter>
  10. :
  11. : -------------------------------------------------------------------
  12. :
  13. : Prof2d data types are declared.
  14. : Data is read from a file in matrix form (command mget);
  15. : the alternative is 2D data in list form (x y z1 z2 z3 ... lines)
  16. : with the command get.
  17.  
  18. row xaxis x1 $ column yaxis y1 $ field data z1 $     
  19. mget example1.mat xaxis yaxis data
  20.  
  21. pause
  22.  
  23. :
  24. : A first look is taken at the data using the contour plot preview
  25. : command cview. 
  26. : Give  <enter>  to exit the graphical mode!
  27. :
  28.  
  29. pause
  30. cview xaxis yaxis data
  31.  
  32. :
  33. : This plot does not provide enough information. We will set some 
  34. : contour values by hand and take a second look.
  35. : Again, give  <enter>  to leave the picture.
  36. :
  37.  
  38. setcont 0.0001 0.0005 0.001 0.005 0.01 0.05 0.1 0.5 1 $
  39.  
  40. pause
  41. cview xaxis yaxis data           
  42.  
  43. :
  44. : We decide a plot of the quantity in logaritmic mode would be
  45. : better. This is done by taking the log of the data values into
  46. : a new field data is declared for this purpose.
  47. : The chosen contour values are reset before invoking cview.
  48. :
  49.  
  50. field ldata $
  51. ldata = log(abs(data)+1e-06)
  52. setcont $
  53.  
  54. pause
  55. cview xaxis yaxis ldata
  56.  
  57. :
  58. : To get a different view, the two axes and the field are copied to
  59. : new variables and transposed. Note: the rows and columns are also
  60. : interchanged by the transpose command!
  61. :
  62.  
  63. x1=xaxis
  64. y1=yaxis
  65. z1=ldata
  66.  
  67. transpose x1 y1 z1
  68.  
  69. show rows
  70.  
  71. pause
  72. cview y1 x1 z1
  73.  
  74. :
  75. : We will now write the data and ldata fields to a file in listed
  76. : form using the put command (Note: not mput) that creates list
  77. : format files. You may practise yourself a little by trying to
  78. : read this file again (using get) and making some pictures.
  79. :
  80. : Precision is set to 3 number behind the floating dot.
  81. :
  82. : Note the dollar sign $ to indicate end of field names.
  83. :
  84.  
  85. putprec=3
  86. put example1.lis xaxis yaxis data ldata $
  87.  
  88. :
  89. : Give  quit <enter>  to exit Prof2d
  90. : or continue playing around a little!
  91. :
  92. : Some interesting Prof2d commands not demonstrated in this example
  93. : -----------------------------------------------------------------
  94. :
  95. : help      - read information from the help file
  96. :             give  help help <enter>  for an overview of all commands
  97. :
  98. : view/plot - make 1D cross sections through 2D data
  99. : remx/remy - remove rows/columsn from the data base
  100. : redx/redy - reduce the row/column data density
  101. : mirror    - mirror a row or column with a field
  102. : inte2d    - integrate by 2D trapezium rule (with bilinear interpolation)
  103. :
  104. :
  105.  
  106.  
  107.