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 |
265 |
USE dimensions, only: iim, jjm |
15 |
guez |
139 |
USE dynetat0_m, only: rlonu, rlatv |
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 |