source: LMDZ5/trunk/libf/dyn3d_common/fxy.F @ 3866

Last change on this file since 3866 was 2598, checked in by Ehouarn Millour, 8 years ago

Cleanup in the dynamics: turn serre.h into module serre_mod.F90
EM

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.9 KB
RevLine 
[524]1!
[1403]2! $Id: fxy.F 2598 2016-07-22 09:28:39Z fairhead $
[524]3!
4      SUBROUTINE fxy (rlatu,yprimu,rlatv,yprimv,rlatu1,yprimu1,
5     ,                    rlatu2,yprimu2,
6     , rlonu,xprimu,rlonv,xprimv,rlonm025,xprimm025,rlonp025,xprimp025)
7
[2597]8      USE comconst_mod, ONLY: pi
[2598]9      USE serre_mod, ONLY: pxo,pyo,alphax,alphay,transx,transy
10
[524]11      IMPLICIT NONE
12
13c     Auteur  :  P. Le Van
14c
15c     Calcul  des longitudes et des latitudes  pour une fonction f(x,y)
16c           a tangente sinusoidale et eventuellement avec zoom  .
17c
18c
[2598]19      include "dimensions.h"
20      include "paramet.h"
[524]21
22       INTEGER i,j
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)
28
29#include "fxy_new.h"
30
31
32c    ......  calcul  des  latitudes  et de y'   .....
33c
34       DO j = 1, jjm + 1
[1403]35          rlatu(j) = fy    ( REAL( j )        )
36         yprimu(j) = fyprim( REAL( j )        )
[524]37       ENDDO
38
39
40       DO j = 1, jjm
41
[1403]42         rlatv(j)  = fy    ( REAL( j ) + 0.5  )
43         rlatu1(j) = fy    ( REAL( j ) + 0.25 )
44         rlatu2(j) = fy    ( REAL( j ) + 0.75 )
[524]45
[1403]46        yprimv(j)  = fyprim( REAL( j ) + 0.5  )
47        yprimu1(j) = fyprim( REAL( j ) + 0.25 )
48        yprimu2(j) = fyprim( REAL( j ) + 0.75 )
[524]49
50       ENDDO
51
52c
53c     .....  calcul   des  longitudes et de  x'   .....
54c
55       DO i = 1, iim + 1
[1403]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  )
[524]60
[1403]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   )
[524]65       ENDDO
66
67c
68       RETURN
69       END
70
Note: See TracBrowser for help on using the repository browser.