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