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