source: LMDZ5/trunk/libf/dyn3dmem/coordij.F @ 1907

Last change on this file since 1907 was 1907, checked in by lguez, 10 years ago

Added a copyright property to every file of the distribution, except
for the fcm files (which have their own copyright). Use svn propget on
a file to see the copyright. For instance:

$ svn propget copyright libf/phylmd/physiq.F90
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

Also added the files defining the CeCILL version 2 license, in French
and English, at the top of the LMDZ tree.
  • 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
File size: 912 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      IMPLICIT NONE
14      REAL lon,lat
15      INTEGER ilon,jlat
16      INTEGER i,j
17
18#include "dimensions.h"
19#include "paramet.h"
20#include "comconst.h"
21#include "comgeom.h"
22#include "serre.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.