Commit 25c68e9d authored by thomas.forbriger's avatar thomas.forbriger Committed by thomas.forbriger
Browse files

correction

This is a legacy commit from before 2015-03-01.
It may be incomplete as well as inconsistent.
See COPYING.legacy and README.history for details.


SVN Path:     http://gpitrsvn.gpi.uni-karlsruhe.de/repos/TFSoftware/trunk
SVN Revision: 2945
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent efa3a9d9
c this is file <grepg.f>
c======================================================================
c $Id: grepg.f,v 1.39 2009-11-27 11:05:53 tforb Exp $
c $Id: grepg.f,v 1.40 2009-12-04 11:57:24 tforb Exp $
c
c GREPG.F
c
......@@ -114,6 +114,7 @@ c so we may change the current dimension without
c affecting the CVS
c V2.38 provide adaptive scaling along time axis
c (option -ts, flag adaptivetimescaling)
c V2.39 04/12/09 use correct DIN notation for units
c
c======================================================================
program grepg
......@@ -127,11 +128,11 @@ c
character titleformat*130
c program version:
character*79 version
parameter(version='GREPG V2.38 plot green amplitudes')
parameter(version='GREPG V2.39 plot green amplitudes')
character*120 CVS_VERSION
parameter(CVS_VERSION=
& '$Id: grepg.f,v 1.39 2009-11-27 11:05:53 tforb Exp $')
& '$Id: grepg.f,v 1.40 2009-12-04 11:57:24 tforb Exp $')
c declare variables for io
character*80 filename
......@@ -1134,20 +1135,20 @@ c ----------------------------
c
if (isgrereso) then
titleformat='("phase-slowness/test-slowness - plot of ",a,a)'
xlabel='wavefield phase-slowness (s/km)'
ylabel='resulting phase-slowness (s/km)'
xlabel='wavefield phase-slowness / s km\u-1'
ylabel='resulting phase-slowness / s km\u-1'
elseif (isgabor) then
titleformat='("gabor-matrix (spectrogram) - plot of ",a,a)'
xlabel='frequency (Hz)'
ylabel='time (s)'
xlabel='frequency / Hz'
ylabel='time / s'
elseif (iskspec) then
titleformat='("wavenumber/frequency - plot of ",a,a)'
xlabel='frequency (Hz)'
ylabel='wavenumber (1/km)'
xlabel='frequency / Hz'
ylabel='wavenumber / km\u-1'
else
titleformat='("phase-slowness/frequency - plot of ",a,a)'
xlabel='frequency (Hz)'
ylabel='phase-slowness (s/km)'
xlabel='frequency / Hz'
ylabel='phase-slowness / s km\u-1'
endif
write(title, fmt=titleformat)
& filename(1:index(filename,' ')-1), ' amplitude'
......@@ -1458,7 +1459,7 @@ c plot cross section
endif
call pgenv(xcross(1),xcross(nslo),mincross,maxcross,0,2)
write(title, 50) cux
call pglab('p [s/km]','Re[A(p)]',title)
call pglab('p / s km\u-1','Re[A(p)]',title)
call pgline(nslo, xcross, ycross)
do islo=1,nslo
xcross(islo)=(islo-1)*(maxs-mins)/(nslo-1)+mins
......@@ -1476,7 +1477,7 @@ c plot cross section
mincross=-1.e-6
endif
call pgenv(xcross(1),xcross(nslo),mincross,maxcross,0,2)
call pglab('p [s/km]','Im[A(p)]',title)
call pglab('p / s km\u-1','Im[A(p)]',title)
call pgline(nslo, xcross, ycross)
call pgclos
call pgslct(maindev)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment