! $Header$

SUBROUTINE coordij(lon,lat,ilon,jlat)

  !=======================================================================

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

  !=======================================================================

  USE comconst_mod, ONLY: pi
  USE lmdz_comgeom

  USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
  USE lmdz_paramet
  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)>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 SUBROUTINE coordij