dyn3d/Guide/coordij.f

!
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/coordij.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
!
      SUBROUTINE coordij(lon,lat,ilon,jlat)

c=======================================================================
c
c   calcul des coordonnees i et j de la maille scalaire dans
c   laquelle se trouve le point (lon,lat) en radian
c
c=======================================================================

      use dimens_m
      use paramet_m
      use comconst
      use comgeom
      use serre
      use nr_util, only: pi
      IMPLICIT NONE
      REAL lon,lat
      INTEGER ilon,jlat
      INTEGER i,j


      real zlon,zlat

      zlon=lon*pi/180.
      zlat=lat*pi/180.

      DO i=1,iim+1
         IF (rlonu(i).GT.zlon) THEN
            ilon=i
            GOTO 10
         ENDIF
      ENDDO
10    CONTINUE

      j=0
      DO j=1,jjm
         IF(rlatv(j).LT.zlat) THEN
            jlat=j
            GOTO 20
         ENDIF
      ENDDO
20    CONTINUE
      IF(j.EQ.0) j=jjm+1

      RETURN
      END
1	!
2	! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/coordij.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
3	!
4	SUBROUTINE coordij(lon,lat,ilon,jlat)
5
6	c=======================================================================
7	c
8	c calcul des coordonnees i et j de la maille scalaire dans
9	c laquelle se trouve le point (lon,lat) en radian
10	c
11	c=======================================================================
12
13	use dimens_m
14	use paramet_m
15	use comconst
16	use comgeom
17	use serre
18	use nr_util, only: pi
19	IMPLICIT NONE
20	REAL lon,lat
21	INTEGER ilon,jlat
22	INTEGER i,j
23
24
25	real zlon,zlat
26
27	zlon=lon*pi/180.
28	zlat=lat*pi/180.
29
30	DO i=1,iim+1
31	IF (rlonu(i).GT.zlon) THEN
32	ilon=i
33	GOTO 10
34	ENDIF
35	ENDDO
36	10 CONTINUE
37
38	j=0
39	DO j=1,jjm
40	IF(rlatv(j).LT.zlat) THEN
41	jlat=j
42	GOTO 20
43	ENDIF
44	ENDDO
45	20 CONTINUE
46	IF(j.EQ.0) j=jjm+1
47
48	RETURN
49	END