|
|
1.1 root 1: integer unit, nx, ny
2: parameter (nx=20,ny=20)
3: real f(nx,ny), fmin, fmax, outsid
4: real pi, z
5: real time
6: unit=30
7: time=0
8: outsid = -99999
9: pi=4*atan(1.)
10: z=2*pi/(nx-1)
11: do 200 k=1,1
12: do 100 i=1,nx
13: do 100 j=1,ny
14: f(i,j) = sqrt(k+5.)*sin((i+2*j+0.2*k-2)*z)
15: $ *cos((2*i-j+0.5*k-1.5)*z)
16: $ + .2
17: 100 continue
18: fmin=0
19: fmax=0
20: call mvefil(unit,time,nx,ny,fmin,fmax,outsid,f)
21: time = time + .1
22: 200 continue
23: end
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.