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 coordonnées ilon et jlat de la maille scalaire dans
    ! laquelle se trouve le point (lon, lat).

    USE comgeom, only: rlonu, rlatv
    USE dimens_m, only: iim, jjm

    REAL, intent(in):: lon, lat ! in rad
    INTEGER, intent(out):: ilon, jlat

    !----------------------------------------------------------

    ilon = 1
    do while (ilon <= iim .and. rlonu(ilon) <= lon)
       ilon = ilon + 1
    end do

    jlat = 1
    do while (jlat <= jjm - 1 .and. rlatv(jlat) >= lat)
       jlat = jlat + 1
    end do

  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	132	! Calcul des coordonnées ilon et jlat de la maille scalaire dans
12			! laquelle se trouve le point (lon, lat).
13	guez	3
14	guez	132	USE comgeom, only: rlonu, rlatv
15	guez	127	USE dimens_m, only: iim, jjm
16	guez	3
17	guez	132	REAL, intent(in):: lon, lat ! in rad
18			INTEGER, intent(out):: ilon, jlat
19	guez	3
20	guez	132	!----------------------------------------------------------
21	guez	81
22	guez	132	ilon = 1
23			do while (ilon <= iim .and. rlonu(ilon) <= lon)
24			ilon = ilon + 1
25			end do
26	guez	81
27	guez	132	jlat = 1
28			do while (jlat <= jjm - 1 .and. rlatv(jlat) >= lat)
29			jlat = jlat + 1
30			end do
31	guez	81
32	guez	127	END SUBROUTINE coordij
33	guez	81
34	guez	127	end module coordij_m