source: LMDZ6/branches/SETHET_DECOUPLE/libf/dyn3d_common/coordij.F @ 5440

Last change on this file since 5440 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: 926 bytes
Line 
1!
2! $Header$
3!
4      SUBROUTINE coordij(lon,lat,ilon,jlat)
5
6c=======================================================================
7c
8c   calcul des coordonnees i et j de la maille scalaire dans
9c   laquelle se trouve le point (lon,lat) en radian
10c
11c=======================================================================
12
13      USE comconst_mod, ONLY: pi
14     
15      IMPLICIT NONE
16      REAL lon,lat
17      INTEGER ilon,jlat
18      INTEGER i,j
19
20      include "dimensions.h"
21      include "paramet.h"
22      include "comgeom.h"
23
24      real zlon,zlat
25
26      zlon=lon*pi/180.
27      zlat=lat*pi/180.
28
29      DO i=1,iim+1
30         IF (rlonu(i).GT.zlon) THEN
31            ilon=i
32            GOTO 10
33         ENDIF
34      ENDDO
3510    CONTINUE
36
37      j=0
38      DO j=1,jjm
39         IF(rlatv(j).LT.zlat) THEN
40            jlat=j
41            GOTO 20
42         ENDIF
43      ENDDO
4420    CONTINUE
45      IF(j.EQ.0) j=jjm+1
46
47      RETURN
48      END
Note: See TracBrowser for help on using the repository browser.