[524] | 1 | ! |
---|
| 2 | ! $Header$ |
---|
| 3 | ! |
---|
[5246] | 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)) |
---|
[524] | 17 | |
---|
[5246] | 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 | !----------------------------------------------------------------------- |
---|