! $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

  IMPLICIT NONE
  REAL :: lon,lat
  INTEGER :: ilon,jlat
  INTEGER :: i,j

  INCLUDE "dimensions.h"
  INCLUDE "paramet.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 SUBROUTINE coordij