Changeset 5246 for LMDZ6/trunk/libf/grid/fxy_new.h
- Timestamp:
- Oct 21, 2024, 2:58:45 PM (23 hours ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/grid/fxy_new.h
r1907 r5246 2 2 ! $Header$ 3 3 ! 4 c--------------------------------------------------------------5 REALripx6 REALfx,fxprim,fy,fyprim,ri,rj,bigy7 c 8 c....stretching in x...9 c 10 ripx( ri )= (ri-1.0) *2.*pi/REAL(iim)11 fx ( ri )= ripx(ri) + transx +12 *alphax * SIN( ripx(ri)+transx-pxo ) - pi13 fxprim(ri) = 2.*pi/REAL(iim) *14 *( 1.+ alphax * COS( ripx(ri)+transx-pxo ) )4 !-------------------------------------------------------------- 5 REAL :: ripx 6 REAL :: fx,fxprim,fy,fyprim,ri,rj,bigy 7 ! 8 !....stretching in x... 9 ! 10 ripx( ri )= (ri-1.0) *2.*pi/REAL(iim) 11 fx ( ri )= ripx(ri) + transx + & 12 alphax * SIN( ripx(ri)+transx-pxo ) - pi 13 fxprim(ri) = 2.*pi/REAL(iim) * & 14 ( 1.+ alphax * COS( ripx(ri)+transx-pxo ) ) 15 15 16 c....stretching in y...17 c 18 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 ) )16 !....stretching in y... 17 ! 18 bigy(rj) = 2.* (REAL(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 23 24 cfy(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--------------------------------------------------------------24 ! fy(rj)= pyo-pisjjm*(rj-transy)+coefalpha*SIN(depisjm*(rj- 25 ! * transy )) 26 ! fyprim(rj)= pisjjm-pisjjm*coefy2* COS(depisjm*(rj-transy)) 27 !--------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.