source: trunk/MESOSCALE/LMDZ.MARS/libf_gcm/grid/fxy_new.h @ 1242

Last change on this file since 1242 was 57, checked in by aslmd, 14 years ago

mineur LMD_MM_MARS: ajout du GCM ancienne physique, systeme maintenant complet sur SVN (ne manque que la base de donnees d'etats initiaux)

File size: 749 bytes
Line 
1c--------------------------------------------------------------
2         REAL ripx
3         REAL fx,fxprim,fy,fyprim,ri,rj,bigy
4c
5c....stretching in x...
6c
7        ripx(  ri )= (ri-1.0) *2.*pi/FLOAT(iim) 
8        fx  (  ri )= ripx(ri) + transx  +
9     *         alphax * SIN( ripx(ri)+transx-pxo ) - pi
10        fxprim(ri) = 2.*pi/FLOAT(iim)  *
11     *        ( 1.+ alphax * COS( ripx(ri)+transx-pxo ) )
12
13c....stretching in y...
14c
15        bigy(rj)   = 2.* (FLOAT(jjp1)-rj) *pi/jjm
16        fy(rj)     =  ( bigy(rj) + transy  +
17     *        alphay * SIN( bigy(rj)+transy-pyo ) ) /2.  - pi/2.
18        fyprim(rj) = ( pi/jjm ) * ( 1.+
19     *           alphay * COS( bigy(rj)+transy-pyo ) )
20
21c--------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.