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