1 |
! |
2 |
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/coordij.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $ |
3 |
! |
4 |
SUBROUTINE coordij(lon,lat,ilon,jlat) |
5 |
|
6 |
c======================================================================= |
7 |
c |
8 |
c calcul des coordonnees i et j de la maille scalaire dans |
9 |
c laquelle se trouve le point (lon,lat) en radian |
10 |
c |
11 |
c======================================================================= |
12 |
|
13 |
use dimens_m |
14 |
use paramet_m |
15 |
use comconst |
16 |
use comgeom |
17 |
use serre |
18 |
IMPLICIT NONE |
19 |
REAL lon,lat |
20 |
INTEGER ilon,jlat |
21 |
INTEGER i,j |
22 |
|
23 |
|
24 |
real zlon,zlat |
25 |
|
26 |
zlon=lon*pi/180. |
27 |
zlat=lat*pi/180. |
28 |
|
29 |
DO i=1,iim+1 |
30 |
IF (rlonu(i).GT.zlon) THEN |
31 |
ilon=i |
32 |
GOTO 10 |
33 |
ENDIF |
34 |
ENDDO |
35 |
10 CONTINUE |
36 |
|
37 |
j=0 |
38 |
DO j=1,jjm |
39 |
IF(rlatv(j).LT.zlat) THEN |
40 |
jlat=j |
41 |
GOTO 20 |
42 |
ENDIF |
43 |
ENDDO |
44 |
20 CONTINUE |
45 |
IF(j.EQ.0) j=jjm+1 |
46 |
|
47 |
RETURN |
48 |
END |