!
! $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 dimensions_mod, ONLY: iim, jjm, llm, ndm
USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
          ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
IMPLICIT NONE
  REAL :: lon,lat
  INTEGER :: ilon,jlat
  INTEGER :: i,j



  include "comgeom.h"

  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 SUBROUTINE coordij