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 35 36 37 38 39
| m =20
n=20
dx = 1./m
dy = 1./n
do 1 i = 1,m+1
xi(i) = (i-1)*dx
do 3 j = 1,n+1
eta(j) = (j-1)*dy
xg(i,j) = xi(i)
yg(i,j) = y(xi(i),eta(j))
3 continue
1 continue
real function x( xi,eta )
implicit none
real xi,eta
d = 3.5
xo = 0.5
a = exp(d*(xi/xo))-1.
b = exp(d)-1.
x = xo*a/b
end function x(xi,eta)
c
real function y( xi,eta )
implicit none
real xi,eta
d = 3.5
yo = 0.5
h1 = -(2.*xi*xi )
h2 = (2.*xi*xi )
a = exp(d*(eta/yo))-1.
b = exp(d)-1.
y =(h2-h1)*yo*a/b +h1
end function y(xi,eta)
end
return |
Partager