1 | ! |
---|
2 | ! $Header$ |
---|
3 | ! |
---|
4 | c----------------------------------------------------------------------- |
---|
5 | c INCLUDE 'fxyprim.h' |
---|
6 | c |
---|
7 | c ................................................................ |
---|
8 | c ................ Fonctions in line ........................... |
---|
9 | c ................................................................ |
---|
10 | c |
---|
11 | REAL fy, fx, fxprim, fyprim |
---|
12 | REAL ri, rj |
---|
13 | c |
---|
14 | c |
---|
15 | fy(rj)=ASIN(1.+2.*((1.-rj)/REAL(jjm))) |
---|
16 | fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj)) |
---|
17 | |
---|
18 | fx ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5* REAL(iim) - 1. ) |
---|
19 | c fx ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5* ( REAL(iim) + 1.) ) |
---|
20 | fxprim( ri ) = 2.*pi/REAL(iim) |
---|
21 | c |
---|
22 | c |
---|
23 | c La valeur de pi est passee par le common/const/ou /const2/ . |
---|
24 | c Sinon, il faut la calculer avant d'appeler ces fonctions . |
---|
25 | c |
---|
26 | c ---------------------------------------------------------------- |
---|
27 | c Fonctions a changer eventuellement, selon x(x) et y(y) choisis . |
---|
28 | c ----------------------------------------------------------------- |
---|
29 | c |
---|
30 | c ..... ici, on a l'application particuliere suivante ........ |
---|
31 | c |
---|
32 | c ************************************** |
---|
33 | c ** x = 2. * pi/iim * X ** |
---|
34 | c ** y = pi/jjm * Y ** |
---|
35 | c ************************************** |
---|
36 | c |
---|
37 | c .................................................................. |
---|
38 | c .................................................................. |
---|
39 | c |
---|
40 | c |
---|
41 | c |
---|
42 | c----------------------------------------------------------------------- |
---|