dyn3d/Guide/coordij.f

module coordij_m

  IMPLICIT NONE

contains

  SUBROUTINE coordij(lon, lat, ilon, jlat)

    ! From LMDZ4/libf/dyn3d/coordij.F, version 1.1.1.1 2004/05/19 12:53:05

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

    USE dimens_m, only: iim, jjm
    USE comgeom, only: rlonu, rlatv
    USE nr_util, ONLY: pi

    REAL, intent(in):: lon, lat
    INTEGER ilon, jlat
    INTEGER i, j

    DO i = 1, iim + 1
       IF (rlonu(i)>lon) THEN
          ilon = i
          exit
       END IF
    END DO

    j = 0
    DO j = 1, jjm
       IF (rlatv(j)<lat) THEN
          jlat = j
          exit
       END IF
    END DO

    IF (j==0) j = jjm + 1

  END SUBROUTINE coordij

end module coordij_m
1	module coordij_m
2
3	IMPLICIT NONE
4
5	contains
6
7	SUBROUTINE coordij(lon, lat, ilon, jlat)
8
9	! From LMDZ4/libf/dyn3d/coordij.F, version 1.1.1.1 2004/05/19 12:53:05
10
11	! calcul des coordonnees i et j de la maille scalaire dans
12	! laquelle se trouve le point (lon, lat) en radian
13
14	USE dimens_m, only: iim, jjm
15	USE comgeom, only: rlonu, rlatv
16	USE nr_util, ONLY: pi
17
18	REAL, intent(in):: lon, lat
19	INTEGER ilon, jlat
20	INTEGER i, j
21
22	DO i = 1, iim + 1
23	IF (rlonu(i)>lon) THEN
24	ilon = i
25	exit
26	END IF
27	END DO
28
29	j = 0
30	DO j = 1, jjm
31	IF (rlatv(j)<lat) THEN
32	jlat = j
33	exit
34	END IF
35	END DO
36
37	IF (j==0) j = jjm + 1
38
39	END SUBROUTINE coordij
40
41	end module coordij_m