! ! $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 comgeom_mod_h USE comconst_mod, ONLY: pi USE dimensions_mod, ONLY: iim, jjm, llm, ndm USE paramet_mod_h 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 SUBROUTINE coordij