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