forums.silverfrost.com

[JohnCampbell]

After watching for a long time, I have had a need to look at many (400) time history plots that I have generated in my FE analysis.
I thought I should put a %PL routine in my results analysis to see how it was going.
I found an example of Ken's posts from Oct 2017 and give it a go.

To do and review 400 plots, I converted Ken's program to a subroutine which would be called multiple times.

It has worked very well, but I wanted to improve, by giving each plot a title and dumping to a .png file.

I included %ws for the title (first at the top, then at the bottom) and
%mn for a screen_dump_func to dump the graphics region to a .png file.

The problem I have is export_image@ only appears to grab the %PL region, but not the added text above/below. Puzzling !

How can get the %ws (or text) to be included in the export_image@ region?
Basically I want to annotate the plot.

[PaulLaidler]

John

Try using DRAW_CHARACTERS@ and its associated routines that set the font etc. Also GET_PLOT_POINT@ can be used to get the pixel coordinates of a particular point relative to the axes.

[JohnCampbell]

Thanks Paul,

I did do some more reading of the documentation and found that using %^PL is demonstrated as a way to include a call to draw_characters@ ( trim(plot_desc), lx, ly, 0 ) in the region, via a call-back "Legend".
I created a module for key data, the title and integer function draw_legend ().

So I think I have a solution.

I also call plot test about 80 times, so this repeated definition of a clearwin %PL region does not cause any noticeable problems.

I have workable code for applying the title and a screen dump, if needed.

[Kenneth_Smith]

John,

If a single title above the plot is sufficient, the quickest solution would be something like:

[JohnCampbell]

Hi Ken,

Thanks for your advice and your earlier posted examples, which were very helpful.

I was unsure where the latest documentation is located. I found recent via Plato Help > ClearWin+.

%PL is not included in "Alphabetical index of format codes" ?
My understanding is that :
FTN95 /64 %PL is the "native mode" with the recently developed options, while,
FTN95 (/32) %PL uses the Simpleplot .dll version.

Having now found "Native graph plotting" and code example of %^PL, I can now include annotation in the %PL plotting region and dump my enhanced plots to a .png file for inclusion in reports.

It was helpful to easily include this plot capability in a Fortran program that generates 400 force functions so they can be systematically reviewed, rather than attacking the logistics of loading them all into Excel.

[PaulLaidler]

John

The %pl "native" options are also available for 32 bit programming.

[JohnCampbell]

Paul, thanks for the update , I will have to understand this better.
I thought %pl "native" was a Silverfrost replacement for changes to SimplePlot, which was no longer being developed.
Can you direct me to where I might be able to read more details ?
Although I am only using FTN95 /64 for this graphics.

[PaulLaidler]

John

The native %pl was introduced because SIMPLEPLOT is a third party 32 bit library. So the development was motivated by the move to 64 bits. But the native %pl soon had more options that work equally well for 32 bits.

Documentation can be found the the FTN95 help file and in Silverfrost\FTN95\doc\cwplus.enh both of which can be accessed from the Plato Help menu.

[JohnCampbell]

Thanks for the advice, however I am now wanting to plot multiple (4) lines on a chart.
I am assumed this could be achieved by providing a 2d Y array, but can't find it described ? ie N_GRAPH=4:

integer npts, ng
real*8 x(npts), y(npts,ng)
...
call winop@ ('%pl[native,x_array,N_GRAPHS=4]')
call winop@ ('%pl[y_min=-12,y_max=2,x_min=0,x_max=7.2]')
...
call winop@ ('%pl[link=lines, colour=blue, symbol=0, pen_style=0]')

i = winio@ ('%`bg[white]&')
i = winio@ ('%^pl',x_size,y_size, npts,x,y, draw_legend)

[JohnCampbell]

I got the following to work.

Is it correct, but is there a neater way ?
Is there any advice on the order required for winop@ and winio@ ?

[Kenneth_Smith]

For N_GRAPHS=4, use the [independent] option, now npts is an array npts=[n1,n2,n3,n4] where n# is the number of points per curve.

You must now supply a pair of arrays (x(n1) and y(n1) etc.) for each line.

You can hide any of the individual plots by setting the corresponding element of npts to zero.

So to draw only the first curve npts=[n1,0,0,0]. i.e you don't need two different calls to %pl to set up the two different plots.

[Kenneth_Smith]

This demonstrates the basic idea.

[Kenneth_Smith]

And it then becomes very easy to make the selection of which graphs to plot interactive:

[JohnCampbell]

Hi Ken,

Thanks for these examples.
I am struggling with a large set of results, so I am looking to give Excel away and write a Fortran program that can better manage all the different results.

I am modelling a train crossing a set of girders, so have an analysis that generates a time varying force function for each rail fastener.
I then take each fastener time force function and introduce this into a time history analysis for all these force functions.
Conceptually it is fairly simple, but logistically the data transfer files are huge, so I want a more systematic approach than what Excel can provide.

At the end, I have to pick out the important results from many, for which %PL may be challenging, if I want to compare "n" graphs !

Thanks for your assistance so far and I will keep asking the questions.

[Kenneth_Smith]

John,

If n_graphs is going to become a very large number, then you should look at the %pl option [link=user] and set_plot_mode@ which allow the graphs to be drawn via the call back.