1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
|
program sample02
use dcl
integer,parameter :: nx=36, ny=36
real,parameter :: xmin=0., xmax=360., ymin=-90., ymax=90.
real,dimension(0:nx,0:ny) :: p
!-- 球面調和関数 ----
do j=0,ny
do i=0,nx
alon = ( xmin + (xmax-xmin)*i/nx ) * DCL_PI/180.
alat = ( ymin + (ymax-ymin)*j/ny ) * DCL_PI/180.
p(i,j) = sqrt(1-sin(alat)**2) * sin(alat) * cos(alon)
end do
end do
!-- グラフ ----
call DclOpenGraphics()
call DclNewFrame
call DclSetWindow( xmin, xmax, ymin, ymax )
call DclSetViewPort( 0.2, 0.8, 0.2, 0.8 )
call DclSetTransFunction
call DclDrawScaledAxis
call DclDrawContour( p )
call DclCloseGraphics
end program
|
