PaulLaidler Site Admin
Joined: 21 Feb 2005 Posts: 8011 Location: Salford, UK
|
Posted: Wed Jan 04, 2023 9:51 pm Post subject: |
|
|
This information is from Bill Bardsley who can be contacted via bill.bardsley@simfit.org.uk. (I assume that there is no automatic response to mouse events.)
Three dimensional surface plots and two dimensional contour diagrams using simdem.
The routine surd2s is called as shown below by option 29 of simdem70 with user-supplied coordinates and gives the option to plot either a 3D surface, a 2D contour, or a 3D surface projected onto a 2D contour, all of which can be rotated in multiples of quarter pi and then edited interactively before saving a hardcopy.
Code: |
!
! 3D surface/bar chart/contours
!
i1 = 10
i2 = nsmax
if (nsur.gt.i2) nsur = i2
call getjm1 (i1, nsur, i2, 'No. of divisions required')
xmin = - one
xmax = one
ymin = - one
ymax = one
delta = (xmax - xmin)/(dble(nsur - 1))
atemp = xmin - delta
do j = 1, nsur
atemp = atemp + delta
btemp = ymin - delta
do i = 1, nsur
btemp = btemp + delta
z(i,j) = btemp**2 - atemp**2
enddo
enddo
isend = 4
call surd2s (isend, nsmax, nsur, nsur, &
vector, xmax, xmin, ymax, ymin, z, &
unused)
|
|
|