Changeset 5246 for LMDZ6/trunk/libf/dyn3d_common/fxy.f90
- Timestamp:
- Oct 21, 2024, 2:58:45 PM (23 hours ago)
- File:
-
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/dyn3d_common/fxy.f90
r5245 r5246 2 2 ! $Id$ 3 3 ! 4 SUBROUTINE fxy (rlatu,yprimu,rlatv,yprimv,rlatu1,yprimu1, 5 , rlatu2,yprimu2,6 ,rlonu,xprimu,rlonv,xprimv,rlonm025,xprimm025,rlonp025,xprimp025)4 SUBROUTINE fxy (rlatu,yprimu,rlatv,yprimv,rlatu1,yprimu1, & 5 rlatu2,yprimu2, & 6 rlonu,xprimu,rlonv,xprimv,rlonm025,xprimm025,rlonp025,xprimp025) 7 7 8 9 8 USE comconst_mod, ONLY: pi 9 USE serre_mod, ONLY: pxo,pyo,alphax,alphay,transx,transy 10 10 11 11 IMPLICIT NONE 12 12 13 cAuteur : P. Le Van14 c 15 cCalcul des longitudes et des latitudes pour une fonction f(x,y)16 ca tangente sinusoidale et eventuellement avec zoom .17 c 18 c 19 20 13 ! Auteur : P. Le Van 14 ! 15 ! Calcul des longitudes et des latitudes pour une fonction f(x,y) 16 ! a tangente sinusoidale et eventuellement avec zoom . 17 ! 18 ! 19 include "dimensions.h" 20 include "paramet.h" 21 21 22 INTEGERi,j22 INTEGER :: i,j 23 23 24 REAL rlatu(jjp1), yprimu(jjp1),rlatv(jjm), yprimv(jjm),25 ,rlatu1(jjm), yprimu1(jjm), rlatu2(jjm), yprimu2(jjm)26 REAL rlonu(iip1),xprimu(iip1),rlonv(iip1),xprimv(iip1),27 ,rlonm025(iip1),xprimm025(iip1), rlonp025(iip1),xprimp025(iip1)24 REAL :: rlatu(jjp1), yprimu(jjp1),rlatv(jjm), yprimv(jjm), & 25 rlatu1(jjm), yprimu1(jjm), rlatu2(jjm), yprimu2(jjm) 26 REAL :: rlonu(iip1),xprimu(iip1),rlonv(iip1),xprimv(iip1), & 27 rlonm025(iip1),xprimm025(iip1), rlonp025(iip1),xprimp025(iip1) 28 28 29 29 INCLUDE "fxy_new.h" 30 30 31 31 32 c...... calcul des latitudes et de y' .....33 c 34 DO j = 1, jjm + 135 36 37 32 ! ...... calcul des latitudes et de y' ..... 33 ! 34 DO j = 1, jjm + 1 35 rlatu(j) = fy ( REAL( j ) ) 36 yprimu(j) = fyprim( REAL( j ) ) 37 ENDDO 38 38 39 39 40 40 DO j = 1, jjm 41 41 42 43 rlatu1(j) = fy ( REAL( j ) + 0.25 )44 rlatu2(j) = fy ( REAL( j ) + 0.75 )42 rlatv(j) = fy ( REAL( j ) + 0.5 ) 43 rlatu1(j) = fy ( REAL( j ) + 0.25 ) 44 rlatu2(j) = fy ( REAL( j ) + 0.75 ) 45 45 46 yprimv(j) = fyprim( REAL( j ) + 0.5 )47 48 46 yprimv(j) = fyprim( REAL( j ) + 0.5 ) 47 yprimu1(j) = fyprim( REAL( j ) + 0.25 ) 48 yprimu2(j) = fyprim( REAL( j ) + 0.75 ) 49 49 50 50 ENDDO 51 51 52 c 53 c..... calcul des longitudes et de x' .....54 c 55 56 57 58 59 52 ! 53 ! ..... calcul des longitudes et de x' ..... 54 ! 55 DO i = 1, iim + 1 56 rlonv(i) = fx ( REAL( i ) ) 57 rlonu(i) = fx ( REAL( i ) + 0.5 ) 58 rlonm025(i) = fx ( REAL( i ) - 0.25 ) 59 rlonp025(i) = fx ( REAL( i ) + 0.25 ) 60 60 61 62 63 64 65 61 xprimv (i) = fxprim ( REAL( i ) ) 62 xprimu (i) = fxprim ( REAL( i ) + 0.5 ) 63 xprimm025(i) = fxprim ( REAL( i ) - 0.25 ) 64 xprimp025(i) = fxprim ( REAL( i ) + 0.25 ) 65 ENDDO 66 66 67 c 68 69 END 67 ! 68 RETURN 69 END SUBROUTINE fxy 70 70
Note: See TracChangeset
for help on using the changeset viewer.