source: LMDZ.3.3/trunk/libf/dyn3d/fxy.F @ 1611

Last change on this file since 1611 was 2, checked in by lmdz, 25 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 1.8 KB
Line 
1      SUBROUTINE fxy (rlatu,yprimu,rlatv,yprimv,rlatu1,yprimu1,
2     ,                    rlatu2,yprimu2,
3     , rlonu,xprimu,rlonv,xprimv,rlonm025,xprimm025,rlonp025,xprimp025)
4
5      IMPLICIT NONE
6
7c     Auteur  :  P. Le Van
8c
9c     Calcul  des longitudes et des latitudes  pour une fonction f(x,y)
10c           a tangente sinusoidale et eventuellement avec zoom  .
11c
12c
13#include "dimensions.h"
14#include "paramet.h"
15#include "serre.h"
16#include "comconst.h"
17
18       INTEGER i,j
19
20       REAL rlatu(jjp1), yprimu(jjp1),rlatv(jjm), yprimv(jjm),
21     , rlatu1(jjm), yprimu1(jjm), rlatu2(jjm), yprimu2(jjm)
22       REAL rlonu(iip1),xprimu(iip1),rlonv(iip1),xprimv(iip1),
23     , rlonm025(iip1),xprimm025(iip1), rlonp025(iip1),xprimp025(iip1)
24
25#include "fxy_new.h"
26
27
28c    ......  calcul  des  latitudes  et de y'   .....
29c
30       DO j = 1, jjm + 1
31          rlatu(j) = fy    ( FLOAT( j )        )
32         yprimu(j) = fyprim( FLOAT( j )        )
33       ENDDO
34
35
36       DO j = 1, jjm
37
38         rlatv(j)  = fy    ( FLOAT( j ) + 0.5  )
39         rlatu1(j) = fy    ( FLOAT( j ) + 0.25 )
40         rlatu2(j) = fy    ( FLOAT( j ) + 0.75 )
41
42        yprimv(j)  = fyprim( FLOAT( j ) + 0.5  )
43        yprimu1(j) = fyprim( FLOAT( j ) + 0.25 )
44        yprimu2(j) = fyprim( FLOAT( j ) + 0.75 )
45
46       ENDDO
47
48c
49c     .....  calcul   des  longitudes et de  x'   .....
50c
51       DO i = 1, iim + 1
52           rlonv(i)     = fx    (   FLOAT( i )          )
53           rlonu(i)     = fx    (   FLOAT( i ) + 0.5    )
54        rlonm025(i)     = fx    (   FLOAT( i ) - 0.25  )
55        rlonp025(i)     = fx    (   FLOAT( i ) + 0.25  )
56
57         xprimv  (i)    = fxprim (  FLOAT( i )          )
58         xprimu  (i)    = fxprim (  FLOAT( i ) + 0.5    )
59        xprimm025(i)    = fxprim (  FLOAT( i ) - 0.25   )
60        xprimp025(i)    = fxprim (  FLOAT( i ) + 0.25   )
61       ENDDO
62
63c
64       RETURN
65       END
66
Note: See TracBrowser for help on using the repository browser.