1 |
! |
2 |
! $Header: /home/cvsroot/LMDZ4/libf/grid/fxy_new.h,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $ |
3 |
! |
4 |
c-------------------------------------------------------------- |
5 |
REAL ripx |
6 |
REAL fx,fxprim,fy,fyprim,ri,rj,bigy |
7 |
c |
8 |
c....stretching in x... |
9 |
c |
10 |
ripx( ri )= (ri-1.0) *2.*pi/FLOAT(iim) |
11 |
fx ( ri )= ripx(ri) + transx + |
12 |
* alphax * SIN( ripx(ri)+transx-pxo ) - pi |
13 |
fxprim(ri) = 2.*pi/FLOAT(iim) * |
14 |
* ( 1.+ alphax * COS( ripx(ri)+transx-pxo ) ) |
15 |
|
16 |
c....stretching in y... |
17 |
c |
18 |
bigy(rj) = 2.* (FLOAT(jjp1)-rj ) *pi/jjm |
19 |
fy(rj) = ( bigy(rj) + transy + |
20 |
* alphay * SIN( bigy(rj)+transy-pyo ) ) /2. - pi/2. |
21 |
fyprim(rj) = ( pi/jjm ) * ( 1.+ |
22 |
* alphay * COS( bigy(rj)+transy-pyo ) ) |
23 |
|
24 |
c fy(rj)= pyo-pisjjm*(rj-transy)+coefalpha*SIN(depisjm*(rj- |
25 |
c * transy )) |
26 |
c fyprim(rj)= pisjjm-pisjjm*coefy2* COS(depisjm*(rj-transy)) |
27 |
c-------------------------------------------------------------- |