CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Using gnuplot to plot probe data

Register Blogs Community New Posts Updated Threads Search

Like Tree5Likes
  • 2 Post By tomf
  • 3 Post By SHUBHAM9595

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 13, 2019, 03:36
Default Using gnuplot to plot probe data
  #1
New Member
 
Fab
Join Date: Apr 2019
Posts: 27
Rep Power: 7
Rotidpor is on a distinguished road
Hello everyone.


I started using gnuplot because it is a convenient tool to plot data. Basically, you just have to type a line or two in a terminal to have your data plotted.
Yet I don't understand the syntax, I found a way to plot my residuals using an openfoam command which extracts the data in my log file and puts it in a single column text file. This way I can tell gnuplot to plot it correctly using simple commands.
This is what I do :


foamLog -n "name of my log file"
gnuplot
plot './Uz_0'


It works well for single column data files created by foamLog. Now my problem is that I want to do the same with probes. Often, my probe files are two columns files, time/iteration versus quantity of interest. But for the velocity, the file is much more complicated with "(" and ")" and four columns.


Is there a simple way to tell gnuplot how to plot a two column file ? Is there an openFoam utility that can convert the velocity probe file into something more manageable like several one column files, as foamLog does for the logfile ?


Thank you all for any input.
Rotidpor is offline   Reply With Quote

Old   June 14, 2019, 04:34
Default
  #2
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Breda, Netherlands
Posts: 647
Rep Power: 32
tomf will become famous soon enoughtomf will become famous soon enough
Send a message via MSN to tomf Send a message via Skype™ to tomf
Hi,

Not an OpenFOAM utility, but sed would work for this:

Code:
sed 's/[()]/ /g' postProcessing/probes/0/U > postProcessing/probes/0/U_noBrackets
Or if you just want to keep the file, you can do replacement inside the file itself as:

Code:
sed -i 's/[()]/ /g' postProcessing/probes/0/U
And than plot on the U_noBrackets or new U file.

Regards,
Tom
raj kumar saini and Rotidpor like this.
tomf is offline   Reply With Quote

Old   June 17, 2019, 03:29
Default
  #3
New Member
 
Fab
Join Date: Apr 2019
Posts: 27
Rep Power: 7
Rotidpor is on a distinguished road
Thank you for your answer tomf, it was helpful.


For those who will read this, here is the gnuplot syntax to plot a file with multiple columns : https://stackoverflow.com/questions/...-all-on-y-axis


Another thing : I realised that gnuplot can actually handle parenthesis. You don't need to remove them before trying to plot. Looks like it just ignores them, which is practical.

Last edited by Rotidpor; June 18, 2019 at 05:02.
Rotidpor is offline   Reply With Quote

Old   March 9, 2022, 06:44
Default
  #4
Member
 
MNM
Join Date: Aug 2017
Posts: 69
Rep Power: 9
SHUBHAM9595 is on a distinguished road
To all the FOAMers...

This is a small attempt to save ur time.....The following script can be be used to monitor as well as export the probes of pressure, velocity and temperature. The script is written for 3 probe points....you might need to edit it accordingly as per the number of probes and fields that u r observing

Code:
#! /bin/sh

set terminal pngcairo enhanced font "Century Schoolbook,20.0" size 1500,1100
set output "U1.png"
set grid
set autoscale
set title "Velocity Probe-1"
set ylabel 'magnitude'
set xlabel 'Time'
set style line 1 
plot "< cat postProcessing/probes/0/U | tr -d '()' " using 1:2 lw 3 title '1probeUx' with lines,\
	 "< cat postProcessing/probes/0/U | tr -d '()'" using 1:3 lw 3 title '1probeUy' with lines,\
     "< cat postProcessing/probes/0/U | tr -d '()'" using 1:4 lw 3 title '1probeUz' with lines,
    set terminal x11 0
	set output
	replot
    pause 1

set terminal pngcairo enhanced font "Century Schoolbook,20.0" size 1500,1100
set output "U2.png"
set grid
set autoscale
set title "Velocity Probe-2"
set ylabel 'magnitude'
set xlabel 'Time'
set style line 1 
plot "< cat postProcessing/probes/0/U | tr -d '()' " using 1:5 lw 3 title '2probeUx' with lines,\
	  "< cat postProcessing/probes/0/U | tr -d '()'" using 1:6 lw 3 title '2probeUy' with lines,\
     "< cat postProcessing/probes/0/U | tr -d '()'" using 1:7 lw 3 title '2probeUz' with lines,
    set terminal x11 1
	set output
	replot
    pause 1
    
set terminal pngcairo enhanced font "Century Schoolbook,20.0" size 1500,1100
set output "U3.png"
set grid
set autoscale
set title "Velocity Probe-3"
set ylabel 'magnitude'
set xlabel 'Time'
set style line 1 
plot "< cat postProcessing/probes/0/U | tr -d '()' " using 1:8 lw 3 title '3probeUx' with lines,\
	  "< cat postProcessing/probes/0/U | tr -d '()'" using 1:9 lw 3 title '3probeUy' with lines,\
     "< cat postProcessing/probes/0/U | tr -d '()'" using 1:10 lw 3 title '3probeUz' with lines,
    set terminal x11 2
	set output
	replot
    pause 1
 
set terminal pngcairo enhanced font "Century Schoolbook,20.0" size 1500,1100
set output "pressure_probes.png"
set grid
set autoscale
set title "Pressure Probes"
set ylabel 'magnitude'
set xlabel 'Time'
set style line 1 
plot "< cat postProcessing/probes/0/p | tr -d '()' " using 1:2 lw 3 title 'Pressure probe 1' with lines,\
	 "< cat postProcessing/probes/0/p | tr -d '()'" using 1:3 lw 3 title 'Pressure probe 2' with lines,\
     "< cat postProcessing/probes/0/p | tr -d '()'" using 1:4 lw 3 title 'Pressure probe 3' with lines,
    set terminal x11 3
	set output
	replot
    pause 1
 
set terminal pngcairo enhanced font "Century Schoolbook,20.0" size 1500,1100
set output "Thermal_probes.png"
set grid
set autoscale
set title "Temperature Probes"
set ylabel 'magnitude'
set xlabel 'Time'
set style line 1 
plot "< cat postProcessing/probes/0/T | tr -d '()' " using 1:2 lw 3 title 'Thermal probe 1' with lines,\
	 "< cat postProcessing/probes/0/T | tr -d '()'" using 1:3 lw 3 title 'Thermal probe 2' with lines,\
     "< cat postProcessing/probes/0/T | tr -d '()'" using 1:4 lw 3 title 'Thermal probe 3' with lines,
    set terminal x11 4
	set output
	replot
    pause 60
SHUBHAM9595 is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[General] Extracting ParaView Data into Python Arrays Jeffzda ParaView 30 November 6, 2023 22:00
CFD by anderson, chp 10.... supersonic flow over flat plate varunjain89 Main CFD Forum 18 May 11, 2018 08:31
plot probe data ehsanshams OpenFOAM Post-Processing 15 August 10, 2017 07:00
paraview - cell data or point data on plot over line bye bye my blue OpenFOAM 0 December 13, 2016 07:07
probe and sample point in same place producing different data odellar OpenFOAM Pre-Processing 3 May 2, 2016 03:35


All times are GMT -4. The time now is 09:44.