! $Header$

      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 comconst_mod, ONLY: pi
      
      IMPLICIT NONE
      REAL lon,lat
      INTEGER ilon,jlat
      INTEGER i,j

      include "dimensions.h"
      include "paramet.h"
      include "comgeom.h"

      real zlon,zlat

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

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

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

      RETURN
      END