trunk/dyn3d/fxysinus.f

!
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/fxysinus.F,v 1.1.1.1 2004/05/19 12:53:06 lmdzadmin Exp $
!
      SUBROUTINE fxysinus (rlatu,yprimu,rlatv,yprimv,rlatu1,yprimu1,
     ,                    rlatu2,yprimu2,
     ,  rlonu,xprimu,rlonv,xprimv,rlonm025,xprimm025,rlonp025,xprimp025)


      use dimens_m
      use paramet_m
      use comconst
      IMPLICIT NONE
c
c     Calcul  des longitudes et des latitudes  pour une fonction f(x,y)
c            avec y = Asin( j )  .
c
c     Auteur  :  P. Le Van
c
c

       INTEGER i,j

       REAL rlatu(jjp1), yprimu(jjp1),rlatv(jjm), yprimv(jjm),
     , rlatu1(jjm), yprimu1(jjm), rlatu2(jjm), yprimu2(jjm)
       REAL rlonu(iip1),xprimu(iip1),rlonv(iip1),xprimv(iip1),
     , rlonm025(iip1),xprimm025(iip1), rlonp025(iip1),xprimp025(iip1)

       include "fxy_sin.h"


c    ......  calcul  des  latitudes  et de y'   .....
c
       DO j = 1, jjm + 1 
          rlatu(j) = fy    ( FLOAT( j )        )
         yprimu(j) = fyprim( FLOAT( j )        )
       ENDDO


       DO j = 1, jjm

         rlatv(j)  = fy    ( FLOAT( j ) + 0.5  )
         rlatu1(j) = fy    ( FLOAT( j ) + 0.25 ) 
         rlatu2(j) = fy    ( FLOAT( j ) + 0.75 ) 

        yprimv(j)  = fyprim( FLOAT( j ) + 0.5  ) 
        yprimu1(j) = fyprim( FLOAT( j ) + 0.25 )
        yprimu2(j) = fyprim( FLOAT( j ) + 0.75 )

       ENDDO

c
c     .....  calcul   des  longitudes et de  x'   .....
c
       DO i = 1, iim + 1
           rlonv(i)     = fx    (   FLOAT( i )          )
           rlonu(i)     = fx    (   FLOAT( i ) + 0.5    )
        rlonm025(i)     = fx    (   FLOAT( i ) - 0.25  )
        rlonp025(i)     = fx    (   FLOAT( i ) + 0.25  )

         xprimv  (i)    = fxprim (  FLOAT( i )          )
         xprimu  (i)    = fxprim (  FLOAT( i ) + 0.5    )
        xprimm025(i)    = fxprim (  FLOAT( i ) - 0.25   )
        xprimp025(i)    = fxprim (  FLOAT( i ) + 0.25   )
       ENDDO

c
       RETURN
       END

1	!
2	! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/fxysinus.F,v 1.1.1.1 2004/05/19 12:53:06 lmdzadmin Exp $
3	!
4	SUBROUTINE fxysinus (rlatu,yprimu,rlatv,yprimv,rlatu1,yprimu1,
5	, rlatu2,yprimu2,
6	, rlonu,xprimu,rlonv,xprimv,rlonm025,xprimm025,rlonp025,xprimp025)
7
8
9	use dimens_m
10	use paramet_m
11	use comconst
12	IMPLICIT NONE
13	c
14	c Calcul des longitudes et des latitudes pour une fonction f(x,y)
15	c avec y = Asin( j ) .
16	c
17	c Auteur : P. Le Van
18	c
19	c
20
21	INTEGER i,j
22
23	REAL rlatu(jjp1), yprimu(jjp1),rlatv(jjm), yprimv(jjm),
24	, rlatu1(jjm), yprimu1(jjm), rlatu2(jjm), yprimu2(jjm)
25	REAL rlonu(iip1),xprimu(iip1),rlonv(iip1),xprimv(iip1),
26	, rlonm025(iip1),xprimm025(iip1), rlonp025(iip1),xprimp025(iip1)
27
28	include "fxy_sin.h"
29
30
31	c ...... calcul des latitudes et de y' .....
32	c
33	DO j = 1, jjm + 1
34	rlatu(j) = fy ( FLOAT( j ) )
35	yprimu(j) = fyprim( FLOAT( j ) )
36	ENDDO
37
38
39	DO j = 1, jjm
40
41	rlatv(j) = fy ( FLOAT( j ) + 0.5 )
42	rlatu1(j) = fy ( FLOAT( j ) + 0.25 )
43	rlatu2(j) = fy ( FLOAT( j ) + 0.75 )
44
45	yprimv(j) = fyprim( FLOAT( j ) + 0.5 )
46	yprimu1(j) = fyprim( FLOAT( j ) + 0.25 )
47	yprimu2(j) = fyprim( FLOAT( j ) + 0.75 )
48
49	ENDDO
50
51	c
52	c ..... calcul des longitudes et de x' .....
53	c
54	DO i = 1, iim + 1
55	rlonv(i) = fx ( FLOAT( i ) )
56	rlonu(i) = fx ( FLOAT( i ) + 0.5 )
57	rlonm025(i) = fx ( FLOAT( i ) - 0.25 )
58	rlonp025(i) = fx ( FLOAT( i ) + 0.25 )
59
60	xprimv (i) = fxprim ( FLOAT( i ) )
61	xprimu (i) = fxprim ( FLOAT( i ) + 0.5 )
62	xprimm025(i) = fxprim ( FLOAT( i ) - 0.25 )
63	xprimp025(i) = fxprim ( FLOAT( i ) + 0.25 )
64	ENDDO
65
66	c
67	RETURN
68	END
69