TweetFollow Us on Twitter

Plotting
Volume Number:9
Issue Number:6
Column Tag:Jörg's Folder

Plotting Packages for Fortran

Inspecting TSI Graphics and SuperPlot

By Jörg Langowski, MacTech Magazine Regular Contributing Author

Note: Source code files accompanying article are located on MacTech CD-ROM or source code disks.

Sorry to say for you C++ and Forth lovers, my desk recently gets swamped with goodies from both main players in the Macintosh Fortran game. Absoft has just sent a big parcel with version 3.2 of MacFortran II. Now, the deadline is again approaching in giant leaps, and I have no time to review it for this column; but I promise you a detailed overview in the next one. Here I’ll only mention that Absoft Fortran’s Macintosh run time environment - one of the main shortcomings of previous versions - has finally been improved. Apple Events are supported, as is publish/subscribe, you can use multiple windows, change fonts, change menus, etc. Just from looking at the manual it seems that they offer as many possibilities as LS Fortran, if not more. First tests with the Linpack benchmark indicate that the speed it still about the same (on my old MacII), which means the usual 20-30% faster than LS Fortran (I guess the compiler that makes my MacII turn as fast as a Cray still has to be invented). So stay tuned for a review.

TSI Graphics

Unless you are content with interfacing your Fortran output to a spreadsheet (you saw in the last column how that goes), you might have struggled with writing your own code for graphing the results of a Fortran program. If you have done that, you know how much time you can waste reinventing the wheel. Because there are good plotting libraries around, and one of them is TSI Graphics (Techno-Sciences Inc., 10001 Derekwood Lane, Suite 204, Lanham, MD 20706; phone (301) 577-6000). Language Systems is promoting this package, which works with their Fortran and contains some very powerful tools.

TSI Graphics does several things that are quite tedious to implement if you do it on your own. First of all, it creates one or several - automatically updated - graphics windows. Each window can contain several plot regions in which to display your data. The beautiful thing is that plots created using the TSI plotting routines are automatically adjusted to fit into the active plot area, so you don’t have to worry about setting the length of an axis, the rectangle of a frame, etc. If you have written code that plots correctly in a certain region, the plot will still have the same aspect when you resize the plot area.

The types of plots that are possible are line, bar, pie, X-Y, contour, and 3D mesh surface plots. You can add error bars to plots, create cumulative plots from groups of data, of course use linear or logarithmic scales, set the ranges of the plot axes automatically, label the plot with text in the font, size and style you want, print the plot or save it to a file, and plot in any color (using your own color tables). In fact, you have most of the functions available that a good data graphing program will give you, including the possibility to offset slices of a pie chart to emphasize them. I found one important feature missing - you cannot draw dashed or dotted lines, or change the line thickness. For the rest of us who does not have color printers available, this would be a very important addition - the only way you can distinguish line plots on a black-and-white printer so far is by adding symbols to data points. Adding different line styles to version 2.0 is strongly suggested.

Let us now see how TSI graphics does its work. For the example program, I have chosen the contour plot and the 3D mesh plot, both things you wouldn’t wish you had to write your own code for. We create a 2-dimensional data array and display it both ways in one window, using two different plot areas. The function we plot is

z = 10.0 * exp (-r/5.) * sin(r/2.) *  cos(2.*atan2(x,y))

where r is the distance from the origin x=0, y=0. The output of the program (see listing) is shown in Figure 1.

Fig. 1: Output of the plot program

How is TSI Graphics used? The example shows that you first call a setup routine, InitGraphicsPackage(), set the size of the output window, and then, after creating the data, call the actual plotting routine DrawInNewWindow. This last routine takes one parameter, a pointer to a subroutine that does the drawing. All this happens in the main Fortran program; as you know if you are familiar with the LS Fortran environment, the main program is often used only for setup purposes. After the main program has exited, the run time environment lets you call subroutines to do the actual work by connecting them to menu items. TSI Graphics works closely connected to this LS Fortran run time environment and hooks itself into the main event loop (at least I think so, I haven’t analyzed the code), so that update events to the graphics window are intercepted and the drawing subroutine is called, even after the main program has exited.

This gives you the advantage that you can easily write very modular code. You define your graphics window(s) with the associated drawing routines and your menus with the routines that do the calculation in the main program, and then exit to the run time environment. Then you can initiate your calculations from the menus, change the data that is to be plotted, and have the graphics window automatically update itself (assuming you have written the plotting subroutine so that it takes into account the changed data correctly). Even if you resize the graphics window, all plot areas will be recalculated to fit into the new window.

Look at the example again: all drawing is done through the routine CreatePlot, which receives its parameters through a common block (you can’t pass subroutines with parameters as arguments to other subroutines in Fortran). The first part of the drawing routine creates the 3D mesh plot. Two plotting areas are defined in the window: the two calls to SetPercentPlotRegion define the left and right halves of the window as plotting areas 1 and 2. The coordinates of the areas are given in percent of the total window height and width. We then set the plotting area to 1 (the left half), and inset the plot by 10% from each side of the area (the call to SetPlotRatio). The rectangle to be displayed in the mesh plot goes from -10 to +10 in both x and y directions (the call to Set3DViewingRect), and we look at this rectangle from a certain point in space (the call to Set3DEyePosition). After these setup calls, we are ready to call MeshPlot, which gives us the complicated plot that you see in the left half of the figure, and label it with some text.

The contour plot is then created by the second half of the program. We set the new plotting area, define the x and y scaling (both from -15 to 15), and draw axes on all four sides of the plot. Then we call ContourPlot and are done! You see the result in the right half of the figure. Quite impressive, and I challenge you to write such a routine yourself, you’ll spend a lot of time on it.

SuperPlot

Another philosophy for displaying data is followed by a plotting package that is compatible with both Language Systems and Absoft Fortran: SuperPlot. It is available through SuperSoft, 498 E Robin Rd., Orem, UT 84057, (801) 225-4356. SuperPlot is at an intermediate position between choosing a separate graphing program like Excel for displaying the Fortran data, and an integrated plotting library like TSI Graphics. When you link the SuperPlot routines into your Fortran application, the program behaves as if you had a complete 2D plotting application embedded within your program. When you invoke SuperPlot it replaces the application’s menu bar with one of its own and, once the plot has been drawn, it allows you to edit, save or print the graph, change the axis scaling, line size and styles, add arrows or text, etc. You have to select “Quit to application” from SuperPlot’s menu bar to return control to the FORTRAN program. Multiple calls to SuperPlot can be made from a single program and up to 10 sets of data may be overlaid on the same graph.

In my opinion, SuperPlot is not so much for writing Fortran programs for general distribution (except if you write a general-purpose graphing program, in which case SuperPlot has already done most of the job), but for the research/engineering type user who has to display data from a program and wants to do so without having to save the data first, open a graphics program, read the data in, create the graph etc. If you write a program in which a graphics window is an integral part of the application and the graphics output needs to be updated without user intervention, you need TSI Graphics. Also, the contour and 3D mesh plot options are very powerful and not available in SuperPlot. On the other hand, for creating nicely-looking black and white graphs, SuperPlot is much superior (of course it can do color, too), allowing you to do dashed lines of varying thickness and easily change the plot style without having to recompile the Fortran program.

Next month - a review of Absoft Fortran (as promised), and more.

Example: Window with double plot area, containing a contour and a 3D mesh 
plot

      Program JLDemo
      External CreatePlot
      INTEGER maxContourX, maxContourY
      REAL XYZdata(31,31)
      integer contourColors(128)
      Common /theData/ XYZdata, contourcolors 

!Initialize TSI Graphics
      Call InitGraphicsPackage()

!Create data and load into XYZdata array
      CALL createData(XYZdata)

!Change Output window 
   Call SetOutputGraphWindowSize (20,20,600,300)
   
!Create a graphics window and hook up our subroutine
      Call DrawInNewWindow(CreatePlot)

!Move the FORTRAN text window off-screen
      Call MoveOutWindow(10,800,200,900)

!Exit and let the FORTRAN menus take over
      End
   
   
      Subroutine CreatePlot
   include 'tsiconstants.fi'
      INTEGER i
      REAL XYZdata(31,31)
      integer contourColors(128)
      Common /theData/ XYZdata, contourcolors 

!Set the contour colors
      DO i = 1, 128
        contourColors(i) = mod(i, 7)
      END DO


!Define two plot windows for the contour plot 
!and the 3D plot
      CALL SetPercentPlotRegion(1, 0.02,0.05, 0.48,0.95)
      CALL SetPercentPlotRegion(2, 0.52,0.05, 0.97,0.95)

!3D mesh plot
      CALL ChangePlotRegion(1)
      CALL SetPlotRatio(1, 0.1, 0.1, 0.1, 0.1)

!Set the size of the viewing rectangle
      Call Set3DViewingRect(-10.0,-10.0,10.0,10.0)

!Set the view point location
      Call Set3DEyePosition(15.0, 7.5, -10.0)

!Create the plot using a blue outline and yellow fill
      Call SetColor(blue)
      Call MeshPlot(XYZdata, 31, 31, -15., -15., 
 1.0, 1.0, yellow)

!  Switch to black and label the mesh using 
!18pt Palatino Bold-Italic
      Call SetColor(black)
      Call SetTextStyle('Palatino', 1+2, 0, 12)
      Call LabelGraph(500.0, 900.0, 
 'exp(r/5)sin(r/2)cos(2*phi)', 1, 0)

!Select region for the contour plot
      CALL ChangePlotRegion(2)
      CALL SetPlotRatio(2, 0.1, 0.1, 0.1, 0.1)

!Set plotting range scale 
      CALL SetXYPlottingRange( -15.0, -15.0, 15.0, 15.0)

!Change to black and draw axes
      CALL SetColor(black)

!Manually draw and label the lower pair of axes
      CALL SetXYIntercepts(-15.0,-15.0)
      CALL DrawYAxis(0.5,1)
      CALL LabelYAxis(10, 1)
      CALL DrawXAxis(0.5,1)
      CALL LabelXAxis(10,1)

!Manually draw and label the upper pair of axes
      CALL SetXYIntercepts(15.0,15.0)
      CALL DrawYAxis(0.5,0)
      CALL LabelYAxis(10,0)
      CALL DrawXAxis(0.5,0)
      CALL LabelXAxis(10,0)

!Border the plot area
      CALL BorderPlotRegion(black)
   
!Plot the data
      CALL ContourPlot(XYZdata,31,31,-15.0,-15.0, 
 1.0,1.0, 1.0,contourColors)

      Return
      END


      SUBROUTINE CreateData(XYZdata)
      INTEGER i,j
      REAL x,y,z,r
      REAL XYZdata(31,31)

!The subroutine for creating the data
      DO i = 1, 31
         x = REAL(i-16)
         DO j = 1, 31
             y = REAL(j-16)
  r = sqrt(x*x+y*y)
             z = 10.0 * exp (-r/5.) * sin(r/2.) *  
 cos(2.*atan2(x,y)) 
             XYZdata(i,j) = z
         END DO
      END DO
      END

 

Community Search:
MacTech Search:

Software Updates via MacUpdate

BusyCal 3.7.2 - Powerful calendar app wi...
BusyCal is the powerful, flexible, reliable calendar app for macOS. It's packed with time-saving features and compatible with all leading cloud services including iCloud, Google, Exchange and more.... Read more
BusyContacts 1.4.2 - Fast, efficient con...
BusyContacts is a contact manager for OS X that makes creating, finding, and managing contacts faster and more efficient. It brings to contact management the same power, flexibility, and sharing... Read more
iShowU Instant 1.3.0 - Full-featured scr...
iShowU Instant gives you real-time screen recording like you've never seen before! It is the fastest, most feature-filled real-time screen capture tool from shinywhitebox yet. All of the features you... Read more
Capo 3.8 - Slow down and learn to play y...
Capo lets you slow down your favorite songs so you can hear the notes and learn how they are played. With Capo, you can quickly tab out your songs atop a highly-detailed OpenCL-powered spectrogram... Read more
Viber 11.7.0 - Send messages and make fr...
Viber lets you send free messages and make free calls to other Viber users, on any device and network, in any country! Viber syncs your contacts, messages and call history with your mobile device, so... Read more
Drive Genius 5.3.1 - $79.00
Drive Genius features a comprehensive Malware Scan. Automate your malware protection. Protect your investment from any threat. The Malware Scan is part of the automated DrivePulse utility. DrivePulse... Read more
CorelDRAW 21.2.0.708 - Graphic design so...
CorelDRAW - professional graphic design software for vector illustration, layout, and so much more. Get started quickly and easily with a wealth of intuitive tools, built-in learning materials,... Read more
Navicat Premium Essentials 12.1.27 - Pro...
Navicat Premium Essentials is a compact version of Navicat which provides basic and necessary features you will need to perform simple administration on a database. It supports the latest features... Read more
OmniFocus 3.4.3 - GTD task manager with...
OmniFocus is an organizer app. It uses projects to organize tasks naturally, and then add tags to organize across projects. Easily enter tasks when you’re on the go, and process them when you have... Read more
Microsoft OneNote 16.30 - Free digital n...
OneNote is your very own digital notebook. With OneNote, you can capture that flash of genius, that moment of inspiration, or that list of errands that's too important to forget. Whether you're at... Read more

Latest Forum Discussions

See All

Upcoming mobile MMO TeeTINY Online’s sec...
TeeTINY Online is an upcoming anime-themed MMO which boasts striking visuals and an ambitious open world to explore. Following a successful CBT (closed beta test) on Android, developer TEENY.studio has now released a follow-up test which adds an... | Read more »
Breakout: Dark Prison is a fast-paced ac...
Breakout: Dark Prison is an action RPG from LaterSoft. Set in the aftermath of a deadly virus outbreak your daughter has been taken from you because she has an immunity to the illness in her DNA. Not being a fan of experimentation on children –... | Read more »
Apple Arcade in review
This weekend, Apple Arcade will officially be one month old. That means anyone who signed up for the free trial on day one has a decision to make: Stick with the service and shell out $5 a month, or cancel and go about your merry way. | Read more »
Alluris is a choose-your-own adventure g...
Alluris is an RPG that the developer's are calling a swipe-your-own adventure game. This is because the game incorporates a Reigns-style - swiping left or right - selection mechanic to make all the decisions you'd usually expect to make across... | Read more »
Hello Hero All Stars receives update wit...
The first Hello Hero game hit global platforms in 2013 and proved a huge success, with developer Fincon adding two more entries to this popular series of casual RPG games since. Released in June this year, Hello Hero All Stars brings many of the... | Read more »
Zombieland: Double Tapper, a cartoon idl...
Zombieland: Double Tapper is the idle RPG tie-in to the upcoming Zombieland: Double Tap. Oddly, it's one of two different Zombieland games launching today, with the other being the Switch title Zombieland: Double Tap - Road Trip. [Read more] | Read more »
Rusty Lake's The White Door launche...
Rusty Lake and Second Maze's intriguing point-and-click adventure game, The White Door, is now up for pre-order on the App Store. This one sees you playing as Robert Hill, a mental health patient who is suffering from severe memory loss. The game... | Read more »
Hellrule is an auto-runner inspired by G...
Hellrule is an upcoming auto-runner game from independent developer Pedrocorp where players will take control of a dapperly dressed gentlemen who comes equipped with a razor-sharp umbrella for slicing up his foes. The game will be available for... | Read more »
Grobo is a gravity bending puzzle platfo...
Grobo is a 2D puzzle platformer that marks the first release from developers Hot Chocolate Games. You'll find yourself manipulating gravity as you make your through this title that's available now for iOS and Android. [Read more] | Read more »
Adrenaline, Compulsive Entertainment’s h...
Compulsive Entertainment’s high-octane arcade racer, Adrenaline, has now made its way to the App Store following a successful launch on Google Play. It’s a ton of challenging, fast-paced fun, boasting easy-to-learn controls and a varied selection... | Read more »

Price Scanner via MacPrices.net

11″ WiFi iPad Pros on sale today for up to $2...
Amazon has new 2018 Apple 11″ WiFi iPad Pros in stock today and on sale for up to $250 off Apple’s MSRP. These are the same iPad Pros sold by Apple in its retail and online stores. Be sure to select... Read more
Apple offers clearance 2018 15″ 6-Core MacBoo...
Apple has clearance 2018 15″ 6-Core Touch Bar MacBook Pros, Certified Refurbished, available starting at only $1999. Each model features a new outer case, shipping is free, and an Apple 1-year... Read more
Apple continues to offer refurbished 2019 21″...
Apple has Certified Refurbished 2019 21″ & 27″ iMacs available starting at $929 and up to $350 off the cost of new models. Apple’s one-year warranty is standard, shipping is free, and each iMac... Read more
Verizon offers free iPhone 7 for customers wh...
Verizon is offering a free 32GB iPhone 7 for new or existing customers who open a new line of service, no trade-in required. Cost of the phone is credited to your account monthly for 24 months. The... Read more
13″ 2.4GHz 4-Core MacBook Pros remain availab...
Apple has a full line of Certified Refurbished 2019 13″ 2.4GHz 4-Core Touch Bar MacBook Pros available starting at $1529 and up to $300 off MSRP. Apple’s one-year warranty is included, shipping is... Read more
Lease one iPhone Xs, Xr, or X at Sprint and g...
Purchase an Apple iPhone X, Xr, Xs, or Xs Max at Sprint, and get a 64GB iPhone Xr for free. Requires 2 new lines or 1 upgrade-eligible line and 1 new line. The fine print: “iPhone Xs (64 GB) $37.50/... Read more
How to use your Apple Education discount to s...
Purchase a new MacBook Pro or MacBook Air using Apple’s Education discount, and take up to $200 off MSRP. All teachers, students, and staff of any educational institution with a .edu email address... Read more
Base 2019 13″ 1.4GHz 4-Core MacBook Pro on sa...
Amazon has the new 2019 13″ 1.4GHz/128GB 4-Core Space Gray Touch Bar MacBook Pro on sale for $100 off Apple’s MSRP. This is the same MacBook Pro sold by Apple in its retail and online stores: – 2019... Read more
Save up to $30 on Apple’s AirPods at these re...
Amazon is offering discounts on new 2019 Apple AirPods ranging up to $30 off MSRP. Shipping is free: – AirPods with Charging Case: $144 $15 off MSRP – AirPods with Wireless Charging Case: $169 $30... Read more
Save $15 on Apple Watch Series 5 models on Wa...
Walmart has new Apple Watch Series 5 models on sale for $15 off Apple’s MSRP on their online store. Choose free shipping or free local store pickup (if available). Sale prices for online orders only... Read more

Jobs Board

Best Buy *Apple* Computing Master - Best Bu...
**734688BR** **Job Title:** Best Buy Apple Computing Master **Job Category:** Store Associates **Location Number:** 000297-Reston-Store **Job Description:** The Core Read more
Best Buy *Apple* Computing Master - Best Bu...
**738529BR** **Job Title:** Best Buy Apple Computing Master **Job Category:** Sales **Location Number:** 000258-Lewisville-Store **Job Description:** **What does a Read more
*Apple* Mobility Pro - Best Buy (United Stat...
**741888BR** **Job Title:** Apple Mobility Pro **Job Category:** Store Associates **Location Number:** 000870-N Scottsdale-Store **Job Description:** At Best Buy, Read more
*Apple* Mobile App Developer - eiWorkflow So...
…eiWorkflow Solutions, LLC is currently looking for a consultant for the following role. Apple Mobile App Developer Tasks the role will be performing: ? Mobile App Read more
Essbase Developer - *Apple* - Theorem, LLC...
Job Summary Apple is seeking an experienced, detail-minded Essbase developer to join our worldwide business development and strategy team. If you are someone who Read more
All contents are Copyright 1984-2011 by Xplain Corporation. All rights reserved. Theme designed by Icreon.