libf/dyn3d/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
      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	IMPLICIT NONE
19	REAL lon,lat
20	INTEGER ilon,jlat
21	INTEGER i,j
22
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
35	10 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
44	20 CONTINUE
45	IF(j.EQ.0) j=jjm+1
46
47	RETURN
48	END