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

 
AAPL
$97.67
Apple Inc.
+0.64
MSFT
$44.50
Microsoft Corpora
+0.10
GOOG
$589.02
Google Inc.
-4.33

MacTech Search:
Community Search:

Software Updates via MacUpdate

TotalFinder 1.6.2 - Adds tabs, hotkeys,...
TotalFinder is a universally acclaimed navigational companion for your Mac. Enhance your Mac's Finder with features so smart and convenient, you won't believe you ever lived without them. Tab-based... Read more
Vienna 3.0.0 RC 2 :be5265e: - RSS and At...
Vienna is a freeware and Open-Source RSS/Atom newsreader with article storage and management via a SQLite database, written in Objective-C and Cocoa, for the OS X operating system. It provides... Read more
VLC Media Player 2.1.5 - Popular multime...
VLC Media Player is a highly portable multimedia player for various audio and video formats (MPEG-1, MPEG-2, MPEG-4, DivX, MP3, OGG, ...) as well as DVDs, VCDs, and various streaming protocols. It... Read more
Default Folder X 4.6.7 - Enhances Open a...
Default Folder X attaches a toolbar to the right side of the Open and Save dialogs in any OS X-native application. The toolbar gives you fast access to various folders and commands. You just click... Read more
TinkerTool 5.3 - Expanded preference set...
TinkerTool is an application that gives you access to additional preference settings Apple has built into Mac OS X. This allows to activate hidden features in the operating system and in some of the... Read more
Audio Hijack Pro 2.11.0 - Record and enh...
Audio Hijack Pro drastically changes the way you use audio on your computer, giving you the freedom to listen to audio when you want and how you want. Record and enhance any audio with Audio Hijack... Read more
Intermission 1.1.1 - Pause and rewind li...
Intermission allows you to pause and rewind live audio from any application on your Mac. Intermission will buffer up to 3 hours of audio, allowing users to skip through any assortment of audio... Read more
Autopano Giga 3.6 - Stitch multiple imag...
Autopano Giga allows you to stitch 2, 20, or 2,000 images. Version 3.0 integrates impressive new features that will definitely make you adopt Autopano Pro or Autopano Giga: Choose between 9... Read more
Airfoil 4.8.7 - Send audio from any app...
Airfoil allows you to send any audio to AirPort Express units, Apple TVs, and even other Macs and PCs, all in sync! It's your audio - everywhere. With Airfoil you can take audio from any... Read more
Microsoft Remote Desktop 8.0.8 - Connect...
With Microsoft Remote Desktop, you can connect to a remote PC and your work resources from almost anywhere. Experience the power of Windows with RemoteFX in a Remote Desktop client designed to help... Read more

Latest Forum Discussions

See All

This Week at 148Apps: July 21-25, 2014
Another Week of Expert App Reviews   At 148Apps, we help you sort through the great ocean of apps to find the ones we think you’ll like and the ones you’ll need. Our top picks become Editor’s Choice, our stamp of approval for apps with that little... | Read more »
Reddme for iPhone - The Reddit Client (...
Reddme for iPhone - The Reddit Client 1.0 Device: iOS iPhone Category: News Price: $.99, Version: 1.0 (iTunes) Description: Reddme for iPhone is an iOS 7-optimized Reddit client that offers a refreshing new way to experience Reddit... | Read more »
Jacob Jones and the Bigfoot Mystery : Ep...
Jacob Jones and the Bigfoot Mystery : Episode 2 1.0 Device: iOS Universal Category: Games Price: $1.99, Version: 1.0 (iTunes) Description: Jacob Jones is back in Episode 2 of one of Apples 'Best of 2013' games and an App Store... | Read more »
New Trailer For Outcast Odyssey, A New K...
New Trailer For Outcast Odyssey, A New Kind of Card Battler Posted by Jennifer Allen on July 25th, 2014 [ permalink ] Out this Fall is a new kind of card battle game: Outcast Odyssey. | Read more »
Hay Day – Tip, Tricks, Strategies, and C...
Recently got into Supercell’s other huge hit, Hay Day and could do with some advice on what to do? We’ve got you covered with some helpful trips and tricks to bear in mind! Ticking Along One of the key things to keep in mind while building up that... | Read more »
Monster Head Review
Monster Head Review By Nadia Oxford on July 25th, 2014 Our Rating: :: FEEDING TIMEUniversal App - Designed for iPhone and iPad Racking up a high score with Monster Head is trickier than it first appears. The appeal wears out fairly... | Read more »
Garfield: Survival of the Fattest Coming...
Garfield: Survival of the Fattest Coming to iOS this Fall Posted by Jennifer Allen on July 25th, 2014 [ permalink ] Who loves lasagna? Me. Also everyone’s favorite grumpy fat cat, Garfield. | Read more »
Happy Flock Review
Happy Flock Review By Andrew Fisher on July 25th, 2014 Our Rating: :: HERD IT ALL BEFOREUniversal App - Designed for iPhone and iPad Underneath the gloss of Happy Flock’s visuals is a game of very little substance. It’s cute, but... | Read more »
Square Register Updates Adds Offline Pay...
Square Register Updates Adds Offline Payments Posted by Ellis Spice on July 25th, 2014 [ permalink ] Universal App - Designed for iPhone and iPad | Read more »
Looking For Group – Hearthstone’s Curse...
For the first time since its release (which has thankfully been a much shorter window for iPad players than their PC counterparts), Blizzard’s wildly successful Hearthstone: Heroes of Warcraft CCG is sporting some brand new content: the single... | Read more »

Price Scanner via MacPrices.net

13-inch 2.5GHz MacBook Pro on sale for $1099,...
Best Buy has the 13″ 2.5GHz MacBook Pro available for $1099.99 on their online store. Choose free shipping or free instant local store pickup (if available). Their price is $100 off MSRP. Price is... Read more
Roundup of Apple refurbished MacBook Pros, th...
The Apple Store has Apple Certified Refurbished 13″ and 15″ MacBook Pros available for up to $400 off the cost of new models. Apple’s one-year warranty is standard, and shipping is free. Their prices... Read more
Record Mac Shipments In Q2/14 Confound Analys...
A Seeking Alpha Trefis commentary notes that Apple’s fiscal Q3 2014 results released July 22, beat market predictions on earnings, although revenues were slightly lower than anticipated. Apple’s Mac’... Read more
Intel To Launch Core M Silicon For Use In Not...
Digitimes’ Monica Chen and Joseph Tsai, report that Intel will launch 14nm-based Core M series processors specifically for use in fanless notebook/tablet 2-in-1 models in Q4 2014, with many models to... Read more
Apple’s 2014 Back to School promotion: $100 g...
 Apple’s 2014 Back to School promotion includes a free $100 App Store Gift Card with the purchase of any new Mac (Mac mini excluded), or a $50 Gift Card with the purchase of an iPad or iPhone,... Read more
iMacs on sale for $150 off MSRP, $250 off for...
Best Buy has iMacs on sale for up to $160 off MSRP for a limited time. Choose free home shipping or free instant local store pickup (if available). Prices are valid for online orders only, in-store... Read more
Mac minis on sale for $100 off MSRP, starting...
Best Buy has Mac minis on sale for $100 off MSRP. Choose free shipping or free instant local store pickup. Prices are for online orders only, in-store prices may vary: 2.5GHz Mac mini: $499.99 2.3GHz... Read more
Global Tablet Market Grows 11% in Q2/14 Notwi...
Worldwide tablet sales grew 11.0 percent year over year in the second quarter of 2014, with shipments reaching 49.3 million units according to preliminary data from the International Data Corporation... Read more
New iPhone 6 Models to Have Staggered Release...
Digitimes’ Cage Chao and Steve Shen report that according to unnamed sources in Apple’s upstream iPhone supply chain, the new 5.5-inch iPhone will be released several months later than the new 4.7-... Read more
New iOS App Helps People Feel Good About thei...
Mobile shoppers looking for big savings at their favorite stores can turn to the Goodshop app, a new iOS app with the latest coupons and deals at more than 5,000 online stores. In addition to being a... Read more

Jobs Board

*Apple* Solutions Consultant (ASC) - Apple (...
**Job Summary** The ASC is an Apple employee who serves as an Apple brand ambassador and influencer in a Reseller's store. The ASC's role is to grow Apple Read more
WW Sales Program Manager, *Apple* Online St...
**Job Summary** Imagine what you could do here. At Apple , great ideas have a way of becoming great products, services, and customer experiences very quickly. Bring Read more
Lead Software Engineer, *Apple* Online Stor...
**Job Summary** Imagine what you could do here. At Apple , great ideas have a way of becoming great products, services, and customer experiences very quickly. Bring Read more
Manager, *Apple* Fullfillment Operation (AF...
…cross-functional teams to drive the highest level of program management quality for Apple . You will help plan launch strategy and demand generation programs with Read more
Project Manager / Business Analyst, WW *Appl...
…a senior project manager / business analyst to work within our Worldwide Apple Fulfillment Operations and the Business Process Re-engineering team. This role will work Read more
All contents are Copyright 1984-2011 by Xplain Corporation. All rights reserved. Theme designed by Icreon.