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

animode works

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: 1756
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent 94f5a99b
c this is <gemmodpg.f>
c
c $Id: gemmodpg.f,v 1.3 2005-06-16 11:43:26 tforb Exp $
c $Id: gemmodpg.f,v 1.4 2005-06-16 12:56:14 tforb Exp $
c
c plot GEMINI earth model
c
......@@ -80,6 +80,7 @@ c read file
c----------------------------------------------------------------------
if (debug) then
do i=1,nlayer
print *,rb(i-1)
print *,alpha(i,1),alphah(i,1),beta(i,1),
& betah(i,1),eta(i,1),rho(i,1),qm(i),qk(i)
do j=2,nco(i)
......@@ -88,6 +89,7 @@ c----------------------------------------------------------------------
enddo
print *,' '
enddo
print *,rb(nlayer)
endif
c----------------------------------------------------------------------
c check ranges
......@@ -95,18 +97,18 @@ c check ranges
rmax=min(rmax,sngl(rb(nlayer)))
c----------------------------------------------------------------------
c start plot
if (.false.) then
if (optani) then
call pgp_setdevice(device, 4, 2)
else
call pgp_setdevice(device, 5, 0)
endif
call pgsch(2.3)
c plot velocities
call pgmtxt('T',1.,0.,0.,text)
if (.false.) then
if (optani) then
title='Pv-velocity (km/s)'
call plotcurve(maxlayer, nlayer, rb, alpha, nco, iflso, rmin, rmax,
& title, debug)
call pgmtxt('T',1.,0.,0.,text)
title='Ph-velocity (km/s)'
call plotcurve(maxlayer, nlayer, rb, alphah, nco, iflso, rmin, rmax,
& title, debug)
......@@ -123,6 +125,7 @@ c plot velocities
title='P-velocity (km/s)'
call plotcurve(maxlayer, nlayer, rb, alpha, nco, iflso, rmin, rmax,
& title, debug)
call pgmtxt('T',1.,0.,0.,text)
title='S-velocity (km/s)'
call plotcurve(maxlayer, nlayer, rb, beta, nco, iflso, rmin, rmax,
& title, debug)
......@@ -152,7 +155,7 @@ c
c
integer nsteps, i
parameter(nsteps=400)
real xmax, ymax, xmin, ymin
real xmax, ymax, xmin, ymin,dx
real x(nsteps), y(nsteps)
double precision gemini_par
c
......@@ -171,6 +174,9 @@ c
xmax=max(xmax,x(i))
enddo
if (xmin.eq.xmax) xmin=0.99*xmax
dx=xmax-xmin
xmax=xmax+0.05*dx
xmin=xmin-0.05*dx
if (debug) print *,'DEBUG: go plot xmin/xmax',xmin,xmax
call pgslw(1)
call pgenv(xmin, xmax, ymin, ymax, 0, 2)
......@@ -193,7 +199,7 @@ c
c
integer nsteps, i, gemini_layer, nl
parameter(nsteps=400)
real xmax, ymax, xmin, ymin
real xmax, ymax, xmin, ymin,dx
real x(nsteps), y(nsteps)
c
ymin=rmin
......@@ -212,6 +218,9 @@ c
xmax=max(xmax,x(i))
enddo
if (xmin.eq.xmax) xmin=0.99*xmax
dx=xmax-xmin
xmax=xmax+0.05*dx
xmin=xmin-0.05*dx
if (debug) print *,'DEBUG: go plot xmin/xmax',xmin,xmax
call pgslw(1)
call pgenv(xmin, xmax, ymin, ymax, 0, 2)
......
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